Holographic dark energy from the laws of thermodynamics with Rényi entropy

This article investigates the relationship between the holographic principle and the laws of thermodynamics in explaining the late-time acceleration of the universe. First, we explore the possibilities of generating the standard holographic dark energy (SHDE) from the laws of horizon thermodynamics. Except for entropies that follow an exponent stretched area law, unless we redefine the horizon temperature, we found it challenging to construct a one-to-one correspondence between the dark energies defined by the holographic principle and the laws of thermodynamics. Secondly, in SHDE models, unless we invoke some phenomenological interactions, it is impossible to explain the late-time cosmic acceleration with the Hubble horizon as the IR cutoff. On the other hand, it is possible to induce dark energy as an integration constant using the laws of thermodynamics on the Hubble horizon. These motivated us to explore a feasible way to invoke the holographic principle from the laws of horizon thermodynamics. We show that the additional terms that appear in the modified Friedmann equations on using entropies other than the Bekenstein–Hawking entropy in the first law of thermodynamics can behave like a dynamic holographic dark energy (HDE). We study the features of such an HDE with Rényi entropy as the choice without considering any non-standard interactions. Interestingly, the resulting form of dark energy reduces to the standard cosmological constant when Rényi entropy reduces to the Bekenstein–Hawking entropy. By examining different parameters, we affirm the validity of our approach to dark energy, which respects both holographic principle and thermodynamics.


Introduction
Explaining the recent accelerated expansion of the universe is one of the significant theoretical endeavours in modern cosmology [1,2].A possible explanation requires a cosmic component called the "dark energy" with negative pressure, whose true nature is still obscure.The simplest candidate for dark energy is the cosmological constant, advocated by the ΛCDM model, which should be same as the energy density of the vacuum [3,4].However, the vacuum energy density predicted by the standard quantum field theory (QFT) and the observed dark energy density have a discrepancy of ∼ 120 orders in magnitude [5].This mismatch leads to the finetuning problem.Further, there is no conclusive explanation for why the present value of dark energy density has the same order of magnitude as that of the matter component, other than being a coincidence [6].A practical solution to these problems is to replace the cosmological constant with a dynamical dark energy [7,8].There are many efforts to explain the observed data with several dynamical dark energy models [9,10], and more conclusions will follow as we improve our observations [11].
Among numerous dynamical dark energy models, a promising approach comes from the holographic principle.In its most general sense, the holographic principle refers to the duality between theories of the bulk and its boundary [12][13][14][15][16].This notion of dimensional reduction proposed by 't Hooft [12], based on the non-extensive scaling of black hole entropy, prevents one from over-counting the degrees of freedom.The disparity mentioned above between the measured and the predicted values of the dark energy density is because local quantum field theory over-counts the degrees of freedom.Accounting the holographic principle, Cohen, Kaplan and Nelson [17] conjectured that an effective field theory which connects the energy density and length scale through a saturation entropy could adequately describe the observed dark energy density ( Λ ).More precisely, their conjecture reads, Here,  is the entropy of the cosmological horizon, proportional to its surface area as motivated by black hole physics [18][19][20][21][22][23][24][25], and  is the length scale called the infra-red (IR) cutoff.The IR cutoff is generally a function of the Hubble parameter and its time derivatives [26].Thus, the nonextensive entropy associated with the horizon connects to the entropy of the "dark energy" by the notion of holography.Many authors refer to this relation as the Cohen-Kaplan-Nelson (CKN) bound or simply the "holographic principle" [17,27,28].Although, the term "holographic principle" is widely used in various areas, here we refer to the CKN relation given above.
In standard holographic dark energy (HDE) models, one identifies the dark energy density as the vacuum energy density of the underlying effective field theory and proposes a dynamical expression from the above CKN bound as, The advantage of the HDE model is that it favours a very small value for the current vacuum energy density [29] and is also stable against divergent quantum corrections [30].Lately, Banks and Draper [27] suggested modifications to the CKN relation in cosmology based on the Nariai limit.Originally, the CKN bound was an attempt to establish a correlation between the IR and UV cutoffs in quantum field theories.However, there are new proposals which give different interpretation to the same.For example, Blinov and Draper gives an interpretation in which the QFT degrees of freedom are depleted as a function of scale [28].Their interpretation is consistent with the apparent success of standard QFT in particle physics and its failure at the cosmic scale.It is important to note that Cohen et al. [17] proposed the CKN bound for a 3 + 1 dimensional spacetime with a horizon entropy following the Bekenstein-Hawking area law [18,21,22].Reference [31] proposes a generalization of this bound to discrete higher dimensions.Many authors investigated the implications of this holographic principle in cosmology [32][33][34][35][36][37][38][39][40][41] and other areas of physics [42,43].The above references give solid arguments to adapt the holographic principle and the corresponding dynamical dark energy  Λ ∝ / 4 in cosmology.For instance, in [44], Saridakis demonstrated the validity of the principle with Barrow entropy as the choice.
The connection between horizon thermodynamics and holography [45] led Moradpour et al. [46] to reason that a proper connection between  Λ ,  and  should also come from the laws of thermodynamics.Recently, reference [47] pointed out a thermodynamic origin for the CKN bound based on the principle of free energy minimization.To propose the thermodynamic connection, Moradpour et al. assumed that the total energy (  ) during the late phase of the universe mostly consists of dark energy ( Λ ) alone.Then, using the relation  Λ ∼  Λ    ∝  , where  is the horizon temperature and  is the areal volume bounded by the horizon, they proposed, as the expression for dark energy density.The immediate question is whether there is a one-to-one correspondence between this expression and the HDE density from the CKN bound.If yes, then equation ( 3) and ( 2) must be at least proportional to each other or there should exist some function relation between them.In this article, we show that one cannot confirm such a proportionality or some functional correspondence for all choices of entropy with  Λ coming from the first law of thermodynamics.
Prior to Moradpour et al.'s proposal, Luongo proposed a thermodynamic approach to holographic dark energy [48].They related the holographic minimal information density to the de Broglie wavelength at a given temperature.They demonstrate the existence of an energy density that appears in the energy-momentum tensor, which they interpret differently from the vacuum energy.The essence of their approach lies in finding a proper thermal cutoff associated with massive and massless particles.Although Luongo's approach is based on the holographic principle and laws of thermodynamics, it is different from what we follow in this manuscript, which is based on the conjecture proposed by Cohen et al. [17].
In proposing the relation (3), the authors in [46] assumed a flat 3 + 1 dimensional FLRW universe and the validity of the standard Friedmann equations.However, Golanbari et al. [49] refuted the cosmological model proposed in [46] and pointed out that one must use the modified Friedmann equations when using any generalised entropies.Although reference [49] considered a modified version of Friedmann equations to study cosmic evolution, they used the same definition for dark energy proposed in [46], which according to our calculations, already assumes the validity of the standard Friedmann equations.Together, references [49] and [46] illustrates the inconsistency between the conventional holographic and the thermodynamic approaches in cosmology.This disparity led us to reconsider the following questions.
"Is it possible to model dark energy density from the laws of thermodynamics, which is consistent with the holographic principle?If yes, how? and if no, why not?"Our answer is "yes", and the rest of this manuscript answers "how".
Further, there is no argument that the ΛCDM model is the simplest model that can very well explain the observational data, and so far, no dynamical dark energy models have shown strong evidence against it.Then, it is only reasonable to assume that, if at all dark energy is dynamic, any new approach must closely resemble the ΛCDM model.We show that such a model is possible by combining the notion of holography with the laws of thermodynamics.Additionally, one can see that Moradpour et al.'s proposal is very similar to the entropic force formalism by Verlinde [50].It is thus essential to investigate the possible correlations between entropic dark energy, ΛCDM and HDE models from the perspective of thermodynamics [51].Considering all these statements, we found that the additional entropic functions that appear in the modified Friedmann equations can take up the role of holographic dark energy (HDE) to explain cosmic evolution.We study the features of such an HDE with Rényi entropy as the choice without considering any non-standard interactions.Interestingly, our form of dark energy reduces to the standard cosmological constant when Rényi entropy reduces to the Bekenstein-Hawking entropy.By examining different parameters, we affirm the validity of our approach to dark energy, which respects holography and thermodynamics.
The structure of the article is as follows.First we review the notion of standard HDE and its thermodynamic analogue proposed by Moradpour et al..We then compare and contrast between these two notions using different choices of entropies.Since we found difference between the SHDE and Moradpour et al.'s proposal, we investigated the underlying reasons for the differences.Further, we explored the possible expressions for dark energy from the Clausius relation and the unified first law.We found that the new expression is also different from Moradpour et al.'s proposal and the CKN bound.The calculation confirms that the thermodynamic picture of dark energy density is not the same as the conventional HDE density.We then illustrate specific issues when we blindly assume the standard HDE model with Hubble horizon as IR cutoff.To resolve these inconsistencies, we then propose a new approach to extract the dark energy component from the laws of thermodynamics, which is consistent with the holographic bound.Further, we demonstrated the validity of our approach using Rényi entropy, and investigated the late time cosmological behaviour.We also studied the entropy evolution and found that the new model respects the second law of thermodynamics.Finally we summarise our findings.

Standard holographic dark energy and the thermodynamic analogue
There is no conclusive evidence that quantum fields fluctuate independently over infinitely large cosmological scales [28].Hence, instead of a local QFT, Cohen, Kaplan and Nelson [17] conjectured that an effective QFT with a UV-IR connection bounded by the horizon entropy is a feasible description.In holographic dark energy models, one identifies the dynamical dark energy from the CKN bound as, Here C is a constant,  is the horizon entropy, and  is the IR cutoff (all in natural units)."The sole purpose of this approach is to replace the vacuum energy density predicted by the standard QFT with an appropriate one".Here, the dynamical dark energy given in eq. ( 4) simultaneously resolves the fine-tuning and coincidence problems.In the standard HDE approach, although  Λ depends on the choice of entropy, the Friedmann equations remain the same as in the 3 + 1 Einstein's gravity.In other words, HDE modifies the stress tensor (  ) rather than the Einstein tensor (  ).
At this juncture, we must remember that the CKN relation is a bound that predicts a consistent value of the cosmological constant, rather than a dynamical function.The holographic principle thus only demands the dark energy density to be / 4 , not necessarily equal or proportional to it.Further, the notion of dark energy density to vary as a function of ∼  −2 also comes from the renormalization group flow expression given as [52], Here,  0 is a constant of integration, and the first order term is similar to the HDE with Hubble horizon as IR cutoff.
If this running vacuum approach respects the holographic principle, the motivation for HDE must be to identify a term that scales like ∼  −2 up to the first order, not necessarily equal or proportional to  −2 .Now, why should one approach the HDE models from the laws of thermodynamics?Since the discovery of black hole thermodynamics, many seminal works illustrated the deeper connection between the equations of gravity and thermodynamics.Works by Jacobson [53], Cai and Kim [54], Padmanabhan [55], Akbar and Cai [56], Hayward et al. [57,58] are some notable ones.Along this line, one can study the cosmic evolution using the laws of thermodynamics and extend it to Gauss-Bonnet and more general Lovelock theories of gravity [54,56].Since horizon entropy connects holography and thermodynamics, it is natural to anticipate that any quantity described in one framework should also have a description in the other.
Considering the thermodynamic properties of the cosmological horizon [45,53,56], Moradpour et al. [46] proposed a thermodynamic analogue of the standard HDE density.Assuming dark energy as the dominant component in the late epoch of cosmic evolution, they argued that the relation    Λ ∼  Λ  is a good approximation, with   as the total energy,  Λ as the dark energy,  Λ as the dark energy density and  as the areal volume.Then, they proposed a dark energy density which takes the form, Here C  is the proportionality constant.Similar to the standard HDE approach, the authors assumed the validity of standard Friedmann equations and only accounted for modifications to the vacuum energy density.
Before going into the details, or verifying the validity of Moradpour et al.'s proposal, let us contrast it with the standard HDE approach.To compare, we will use different nonextensive entropies for  and use the apparent horizon radius, r , as the IR cutoff ().Using  = r allows one to use the standard definition of horizon temperature,  = 1/(2r  ), in eq. ( 6).Although there are several choices for entropies, here we will restrict ourselves to Bekenstein-Hawking, Barrow [59], Tsallis [60] and Rényi [61] entropies.Other entropies such as Sharma -Mittal entropy [62], Kaniadakis entropy [63] or Nojiri -Odintsov -Faraoni entropy [64] represents different generalization schemes, which we will address in future works.

Standard HDE versus Moradpour et al.'s proposal
In this section we compare the Standard HDE with Moradpour et al.'s proposal.To compare, we consider two definitions for dark energy density,  Λ1 from the CKN bound and  Λ2 from Moradpour et al.'s proposal.Each of which respectively reads, In the previous section we discussed the motivation for the above expressions.The entropy choices to contrast the above definitions are (in natural units), Here, Δ,  and  are the non-extensive parameters for Barrow, Tsallis and Rényi entropies respectively.Although some of the above entropies might look similar, the motivations behind each are different.For more details on different nonextensive entropies, one can consult references [64,65] and the references therein.Here we have,

Barrow entropy:
In 2020, John Barrow considered fractal structure to the black hole horizon and introduced a toy model to account for the effects of quantum gravity on spacetime [59].Due to the fractal structure, the horizon that bounds a finite volume can have a finite or an infinite surface area.
Here, 0 ≤ Δ ≤ 1 is the permissible range of the non-extensive parameter Δ.When Δ = 0, we recover a smooth horizon that satisfies the standard area law.Although Δ can be a varying function [66], here we consider it as an unknown constant.

Tsallis-Cirto entropy:
One of the biggest difference between the conventional Boltzmann-Gibbs entropy and the Bekenstein-Hawking entropy is how they scale.While the former scales as ∼  3 , the latter scales like ∼  2 .In order to describe them in a unified framework, Tsallis and Cirto [67] proposed a generalised version of Tsallis statistics [60] to replace the Boltzmann-Gibbs statistics.The new generalised entropy reads, where  is the Tsallis parameter,  is the Tsallis-Cirto parameter,  is the total number of internal configurations and {  } is the probability distribution.Here log  () ≡  1− − 1 /(1 − ).For  =  = 1, the above expression is additive, and non-additive otherwise.The new composition rule reproduces the ordinary and the black hole entropy for different values of non-extensive parameters  and .In the context of black holes and cosmology, it takes the form given in equation (8) as mentioned in [68].For a system that obeys the area law, Total number of internal configuration  is a function of  given as, ln  () ∝  −1 , ( > 1).
Here  is the length of the system such that   is the volume in -dimension.Then Tsallis and Cirto showed that,  1, =/(−1) ∝  /(−1)

BH
, where  BH is the Bekenstein-Hawking area law.In other words,  TC =  1,  ∝   BH .Thus, in the context of cosmology we consider only a single parameter Tsallis-Cirto entropy by fixing the Tsallis parameter to unity.Here we assume  as an unknown constant and  is a numerical factor introduced to make the expression dimensionally consistent in appropriate units (natural units in this case).Additionally, the value of  also depends on .One may also assume a dynamic , which we do not consider in this manuscript.

Rényi entropy:
Another generalized version of entropy is the Rényi entropy.It also appears as the formal logarithm of the Tsallis entropy, which allows one to confirm the compatibility of non-extensive thermodynamics and the zeroth law [69].Here,  is the non-extensive parameter introduced by Rényi in [61].Unlike Tsallis-Cirto entropy discussed above, Rényi entropy is a single parameter generalization of the standard Boltzmann-Gibbs entropy.In the context of cosmology and black hole horizons, one can write the Rényi entropy as given in equation (8).Czinner and Iguchi [70] introduced this expression to study the thermodynamic stability of black holes.This follows form the fact that the Rényi entropy is the formal logarithm of Tsallis entropy and Bekenstein-Hawking entropy is a special case of Tsallis-Cirto entropy.In the context of quantum information theory, Rényi entropy is a suitable measure of entanglement between quantum states other than the standard von Neumann entropy.In the classical information theory, it is the generalization of the renowned Shannon entropy [71].The definition of Tsallis-Cirto and Rényi entropy relies on the assumption that, ln  ∝ Area is correct.And it is this proportionality with the area, instead of volume, that makes these entropies non-extensive.Now, to compare equation ( 7a) and (7b), let us consider the simplest 3 + 1 dimensional flat FLRW universe, where r = 1/.Now using different choices of entropies given in ( 8), we compute the corresponding expressions for dark energy densities, and we summarize the results in table (1).In the calculations, we assumed  = 1/(2r  ) as the horizon temperature and  = (4/3)r 3 as the areal volume.For a constant non-extensive parameter, one can easily see that both definitions appear to be proportional to each other for Bekenstein-Hawking, Barrow and Tsallis entropies.However, the case differs for Rényi entropy when  ≠ 0. This result indicates that, although both definitions give similar scaling, they are not consistent with each other for all entropies.
This inconsistency could arise due to many reasons, one of which might be the simple proportionality argument used in the relation (7b).In the standard thermodynamic approach, the choice of entropy modifies the Friedmann equations, and this modification affects the relation between the deceleration parameter () and the total equation of state parameter ().In the following sections, we will derive an expression for dark energy density in the late epoch from the laws of thermodynamics and show that, even in the simplest case of a flat 3 + 1 dimensional FLRW universe, it is impossible to ensure the proportionality for all entropy choices.
Another reason for the above inconsistency might be the standard horizon temperature and volume in eq.(7b).In black hole thermodynamics, the Hawking temperature is,  H = (/ ) −1  , , where ,  and  are the black hole's mass, angular momentum and electric charge respectively [23].If we modify the functional form of  without modifying , this relation will yield different expressions for  H . Nojiri et al. [72] noted this critical aspect and illustrated how the use of different entropies leads to inconsistencies in black hole thermodynamics.Thus, when we replace Bekenstein-Hawking entropy with other non-extensive entropies, even if  has the dimension of [r −1  ], it is not obvious that  ∝ 1/r  is always valid.Further, if we follow the arguments by Cai and Kim [54] and assume that the temperature corresponding to the apparent horizon is 1/(2r  ), we cannot claim that equations (7a) and (7b) represent the same quantity.
Furthermore, in the standard thermodynamic approach, one can only obtain the standard Friedmann equations starting from the Bekenstein-Hawking entropy.In the process, we will get an integration constant, which many authors consider as the cosmological constant [54].Thus, in the thermodynamic approach, adding additional cosmic components to explain the accelerated expansion is generally not a trend.Restricted covariant theories also give the same explanation, where the cosmological constant depends on the initial conditions [73].
Finally, the inconsistency might be due to all of the above reasons.Then, it is not trivial to bring a simple connection between holography and thermodynamics.The argument that any modified version of HDE density, either from equation (7a) or from relations similar to equation (7b), reducing to the standard HDE model with Hubble scale is not pressing enough, as such an HDE model cannot explain the late time acceleration.Instead, if an expression reduces to the ΛCDM model, that would be more convincing.
To derive the expression for dark energy density from the laws of thermodynamics we need certain tools, which we will introduce here.For this purpose, we will assume a general  + 1 dimensional FLRW spacetime with the metric given as, Here, ,  ∈ (0, 1), , with scale factor '()', comoving distance '',  − 1 spatial surface metric 'Ω −1 ' and ℎ 10 = ℎ 01 = 0. Further,  can take values 1, 0 and -1 for closed, flat and open universes, respectively.It is convenient to discuss further in terms of distance r = (), as we are interested in a spherically symmetric apparent horizon which satisfies the condition ℎ    r  r = 0.One can satisfy this condition for r = r given as, Here, r is the radius of the apparent horizon.For a flat universe, with  = 0, the above relation reduces to r = 1/ with  as the Hubble parameter.Additionally, is the area of the apparent horizon.The temperature of apparent horizon, as motivated by black hole thermodynamics is,  = /(2), where,  = (2 is the surface gravity.A straightforward calculation using the FLRW metric yields, Using the above definitions, we can derive the energy density bound from the Clausius relation and the unified first law.Before that, let us see the problem with the proportionality argument used in [46] when we start from the first law of thermodynamics. 4 Is   ∝  (/) consistent with the first law?
In this section, we will reconsider the simplest case of a flat 3 + 1 dimensional FLRW universe and assume that the total energy density is  ≈  Λ .We then check whether the relation  Λ ∝  (/) is consistent with the first law.Now from the Clausius relation we have, − =   and the explicit expression reads [54], With  = (4/3)r 3  , we have  = 4r 2  r  and from r = 1/, we have r  / = − / 2 =⇒  = ( 2 / )r  .Thus the above expression becomes, Since pressure  := , we finally get, The above expression is an exact equation for the "total" energy density.Since we assumed a dark energy-dominated universe, we can follow Moradpour et al.'s argument and write, Comparing the above expression with Dividing the equation (18b) by equation (18a) we get, Thus, for Bekenstein-Hawking entropy, Moradpour et al.'s proposal appear to be valid.We again emphasise that, in the thermodynamic approach, one can only derive the standard Friedmann equations using the Bekenstein-Hawking entropy.
Let us now consider other non-extensive entropies.Since Barrow and Tsallis-Cirto entropies are mathematically similar, we will consider Barrow-modified Friedmann equations.With the same background assumptions, the Barrow corrected Friedmann equations read [74] (in natural units), With  =  +  2 , one can easily show that, Since we started with a constant Δ, this again is a constant.Hence, for Barrow and Tsallis-Cirto entropies, Moradpour et al.'s proposal appear to be valid.This is why we get similar expressions for dark energy density from equation (7b) and eq.(7a).(See table (1)).
In the same manner, considering the Rényi entropy, we get the corresponding Rényi modified Friedmann equations as, (See Appendix A.) In this case, we get, which is a function of  for non-zero . the first law of thermodynamics, the proposal in ref. [46] is only consistent with the standard HDE when we consider the standard Friedmann equations or entropies that follow exponent stretched area laws.One can also prove this by equating equation (7a) and (7b), and then solve for .In doing so, we get, Here, K accounts for the integration constant.When 8 2 C/C  = 2, we recover the standard area law, and any deviation might point towards Barrow or Tsallis-Cirto type area stretched entropies.This is why we say that the simple proportionality argument used in [46] may not be generally valid.The other possibility is to redefine the horizon temperature or the areal volume as in [72].Redefining the temperature essentially redefines the surface gravity and thus the underlying metric.Since that would lead to a universe different from what we expect from an FLRW spacetime, redefining the temperature appears premature.
Till now, we considered the first law of thermodynamics and the Friedmann equations (standard/modified).However, one needs to consider the continuity relation to get the Friedmann equations from the first law.In the following subsections, we use the continuity equation and derive an expression for energy density without explicitly using the standard (modified) Friedmann equations and study whether it can resolve the inconsistency we encountered earlier.For simplicity, we do not consider any interaction between dark energy and matter.It will be interesting to investigate such situations, which we will consider elsewhere.

Dark energy density from the Clausius relation
In the previous section we checked whether  Λ ∝  (/) is consistent with the first law or not.It is possible to reverse the question and ask "what form of  Λ one can obtain from the first law?".In this section we derive an expression for the dark energy density from the Clausius relation.One of the seminal works that first illustrated the deep connection between gravity and thermodynamics belongs to Jacobson [53], where he obtained Einstein's field equation from the Clausius relation,  =  .In cosmology, we can take  = − as the heat flow across the apparent horizon, and the Clausius relation reads [54], ( Here,  is the total energy density,  is the pressure,  is the Hubble parameter,  is the number of spatial dimensions, and the over-dot represents the derivative with respect to cosmic time.One may also write this relation using the equation of state parameter defined as  := /.Now, using the expressions for ,  and the continuity equation, we arrive at, Plugging this expression back to the Clausius relation given in equation ( 25), we get, Now we approximate the expression for  in equation ( 12) as, Although this is an approximation, it is valid, as the first law holds the same for an infinitesimal change in temperature.Thus we get the total energy density as We will show that the above result remains unchanged even when considering the unified first law, where we do not take the above approximation in temperature.Now, the total energy density () consist of different cosmic components such as matter (including dark matter) (  ), radiation (  ), dark energy ( Λ ) etc.In this article, we restrict ourselves to a period around transition redshift to late acceleration.Since   dilutes much faster than any other components, we may ignore its effect in our analysis and assume the relation   Λ +   for an almost flat universe.Then the dark energy density is, In the asymptotic de Sitter limit, when redshift  → −1, the matter density, which scales as   ∝ (1 + ) 3 , dilutes away, and we get a "limit" where, In [46], the authors considered this bound as the source for HDE.Then the problem is that this expression is valid only in the limit where no matter component exists.One may also assume such a limit for a matter-dominated era, and we will end up in a cyclic loop of arguments [75].However, if we demand to define a thermodynamic analogue of HDE density from this limit as in [46], it can be of the form, Here,  2 is the proportionality constant.We will now derive the same from the unified first law for completeness and then discuss its features.

Dark energy density from the unified first law
Similar to the previous section, here we derive a possible expression for dark energy density from the unified first law.The notable work by Akbar and Cai [56] used the unified first law to discuss Einstein, Gauss-Bonnet and Lovelock gravities in a unified framework.In the previous section, we used an approximation in the temperature for our calculations, while in the unified first law, there is no need to take such assumptions.Here, the unified first law of thermodynamics reads, where,  = Ω  is the volume,  = 1 2 ( − ) is the work density and  =  is the total gravitating energy bounded by the horizon [76].Substituting for  in equation (34) we get, Now, using the continuity equation given in equation ( 26) and  = Ω  r−1  r  to equation (35), we arrive at, Now the proper definition of temperature in terms of surface gravity,  = /(2) reads, Then immediately, we get, Once again, in the asymptotic de Sitter limit, when redshift  → −1, the matter density, which scales as   ∝ (1 + ) 3 , dilutes away, and we get the limit, Interestingly, for Rényi entropy, equation ( 33) is also different from Moradpour et al.'s definition given in equation (7b).Why should the laws of thermodynamics give distinct expressions for the same quantity?This could be because, we considered only the validity of generalized continuity equation, while equation (7b) demanded the validity of equation (17).However, one cannot guarantee the validity of equation Reference [77] gives a detailed account of the Misner-Sharp energy ( ), and references therein point out various other energy candidates.

Entropy
Λ using equation ( 33) Different entropies and corresponding dark energy densities using the equation ( 33).
(17) from the Clausius relation or the unified first law of thermodynamics.Further, it is not immediately clear from the arguments in [46], how one can assume  Λ ∼  Λ  (the Misner-Sharp energy) and consider  Λ =  Λ .
As already mentioned, standard Friedmann equations need Bekenstein-Hawking entropy as the entropy choice.Then, it is unclear how one can choose one entropy to define the equations of motion and another to define dark energy density from the same governing laws.For example, it is not logical to construct modified Friedmann equations using Rényi entropy and use Barrow entropy to define the dark energy density from the same laws of thermodynamics.A possible justification might be that the equations of motion come from Einstein's field equation and the expression for dark energy density has its roots in holography or thermodynamics.However, in the thermodynamic approach, Einstein's gravity, which yields the standard Friedmann equations, is only a particular case when we use Bekenstein-Hawking entropy.This leads us to investigate the possibility of proposing a dark energy component using the laws of thermodynamics and the notion of holography.

Standard Rényi HDE and its inconsistencies
Before we jump into a new definition, let us see some inconsistencies in the standard approach to HDE, especially with Rényi entropy.Here, we consider the Rényi entropy and propose Rényi HDE density using the equation (7a).It is interesting to investigate this, as there is no such work in the literature to our knowledge that uses equation (7a) for Rényi HDE density.This might be because of the lack of simple analytical solutions for the Hubble parameter.Here the Rényi HDE density from the CKN bound reads, We will consider the Hubble horizon as the IR cutoff for our discussion.Now, the standard Friedmann equations are, For simplicity, in our analysis, we do not consider any interaction between HDE and matter.Thus each component obeys independent continuity equations, and for the matter sector, we get, Here  is the redshift.Since it is convenient to deal with density parameters rather than actual densities we use, Thus, using the expression for dark energy density, matter density and  () =  0 ℎ(), we get, where Ω 0 and  0 are the present value (where  = 0) of the matter density parameter and Hubble parameter.Now to obtain the expression for the free parameter C, we consider the situation when  = 0 for which we have ℎ( = 0) = 1.
Then we get, Now the previous expression in terms of  =  () becomes, This equation makes it difficult to get an exact analytical form for  as a function of .Thus, we resort to numerical methods to further investigate the model.It is noteworthy to recognize Moradpour et al.'s results, as it is possible to get an exact expression for  in terms of  from their proposal.In figure .(2), we plot the evolution of the Hubble parameter as a function of redshift for different values of .The model does not exhibit a late time acceleration for zero or large positive values of .However, similar to Moradpour  [78]).OHD corresponds to the Observational Hubble Data set with 57 entries [79] et al.'s results, we get a late time acceleration for negative values of , which could explain the observation.The limit in which  → 0, the expression reduces to ordinary HDE model and cannot explain the cosmic evolution, as we did not consider any interactions.
To study the features of the cosmic evolution, we now investigate the deceleration parameter  for  that closely follows the data points.The definition of  is, Evaluating the previous expression, we plot the evolution for  as a function of redshift in figure (3) for different values of .The model shows late time acceleration with transition redshift close to  = 0.5.Similar to ΛCDM, this model also exhibits a final de Sitter limit since  → −1 implies  → 0, which then implies a constant Hubble parameter.Moreover, we can evaluate the behaviour of the total equation of state parameter ().For this, we first use, The validity of the above expression depends on whether we consider the standard Friedmann equations or not.The behaviour of  is similar to , as the above equation indicates a linear relationship between them.we have, One can see from the behaviour of  Λ in figure (4) that the dark energy equation of state crosses the phantom divide and approach the phantom divide as  → −1.It is interesting to note that  never crosses the phantom divide ( = −1) while  Λ does.This is due to the linear relation between  and , as  never crosses the phantom divide, nor does .This feature is the consequence of assuming the validity of standard Friedmann equations.On the other hand, for calculating  Λ , we used the continuity equation instead of the second Friedmann equation.This is why we could see the Phantom behaviour of the dark energy density.This is the first logical inconsistency we encounter in all standard HDE models.To clarify, let us calculate the total equation of state parameter using the continuity relation.When we consider the exact relation  =  Λ /(  +  Λ ), without assuming the validity of the standard Friedmann equations, we cannot arrive at the relation (48).To see this, consider the expression,  Λ = (  Λ /(−3)) −  Λ from the continuity relation of dark energy.Now the new expression for total equation of state parameter ( w) becomes, This expression crosses the phantom divide which one can see from figure (5).Further, the total equation of state parameter gives a higher value in the matter-dominated era, where we expect it to be zero.Thus we end up in a situation where we get two different expressions for the same total equation of state parameter.Thus, it is clear that there are inconsistencies in assuming the validity of standard Friedmann equations and the continuity equation simultaneously.Now, why does the total equation of state parameter goes above zero in figure ( 5)?This must be because of the nature of dark energy itself.If the matter is not the dominant component in the early (not as early as the radiation dominant period) epoch, then the behaviour of dark energy affects the total equation of state.To see this, we study the behaviour of the dark energy density parameter given as Unlike the ΛCDM model, here, the value of Ω Λ do not tend to zero for high redshift.Instead, it settles at a nonzero value (See figure (6)).From the above expression, for  → ∞, we have,  Thus for a fixed value of , Ω Λ remains a constant for large .This effectively gives the wrong total equation of state, when we consider the relation  =  Λ /(  +  Λ ).Because in the matter-dominated era, total pressure must be close to zero, and matter density must dominate the total energy density.Thus, even if w goes above zero, it should at least come down to zero for large values of .Here, since Ω Λ never approaches zero, it contributes to the total equation of state parameter.Thus, simultaneously considering the validity of Friedmann equations and the validity of continuity equations for HDE models may lead to counter-intuitive results.All of this demands us to reconsider the standard HDE approach.So far, we did not bother about the possible values of .In the equation ( 40), the dependency on  comes from the factor log 1 . Now, by the conventional definition of Rényi entropy, the value of  must lie between 1 and −∞ to be Schur concave.A stronger concavity condition demands  to lie between 1 and zero.However, positive values, such that  > 1 is in principle allowed because of the skew symmetry of Rényi entropy about  = 1/2.An interesting situation will arise in our case when / 2 < −1.In such cases, the dark energy density defined using equation ( 40) will be complex.Thus, the value of  cannot be arbitrarily negative.Thus, mathematically, the allowed values for  must be between +∞ to − 2 /.However, in Moradpour et al.'s approach, this problem never arises, as there is no log function present in their expression (See equation (14) in [46]).Article [80] also demands a very tiny positive value for  in the context of holographic equipartition due to Padmanabhan [81].

Thermodynamics near Phantom divide
In the model discussed above, for certain values of , the dark energy equation of state appears to cross the Phantom divide and then reaches the final de Sitter state in the asymptotic future.Since Phantom models are not ruled out by current observations, these results needs some extra attentions.Since our primary motive is to have a consistent thermodynamic picture, let us explore the thermodynamics at the Phantom divide by asking ourself certain questions.

What happens to the temperature at 𝑤 = −1?
In this article we adapt two definitions for horizon temperature, depending whether we use the Clausius relation or unified first law.While using the Clausius relation (− =  ) we have,  = 1/ 2r .Thus, by definition, the temperature does not become zero or negative throughout the entire course of evolution.This temperature also goes by the name "Cai and Kim temperature", as it was introduced in [54].
When the model is at the Phantom divide, i.e.  = −1, we have  :=  = − and we get − = 0.This result implies that either  or  has to be zero.Since by definition  is always positive-non-zero (except for the limit r → ∞), the only choice is to assume that  = 0. Thus, the third law of thermodynamics is perfectly respected and the second law approaches its limit.When the model crosses the Phantom divide (i.e. < −1), we will have a situation in which the entropy decreases.This situation is because  is positive for the outward flux and negative for the inward flux by convention [82].With  as a positive quantity,  becomes negative and the entropy decreases.Hence, the second law of thermodynamics is not respected, when  < −1.
If we were to start from the unified first law, which states that  =   +  , with  = ( − )/2, the temperature from the induced flat FLRW metric is, In the standard HDE adapted above, we assumed the validity of the standard Friedmann equations.Thus, one can rewrite the above temperature in terms of  as, Then, at the phantom divide ( = −1) we have  = −1/(2r  ) and  = .Now the unified first law reduces to,  =   + .Since,  =  (the Misner-Sharp energy) and  =   + , we have,   =  .In the dark energy dominated epoch (i.e. ≈  Λ ) the density  is a constant when  = −1 which implies  = 0. Thus we have   = 0. Based on the previous arguments, temperature, by definition, remains non-zero while entropy is at its extremum.Further, if  becomes less than unity, entropy decreases and the second law no longer holds.However, unlike some Phantom models with  < −1 for  → −1, the above model shows a final de Sitter space.Thus all the laws of thermodynamics are respected in the far future.This might indicate a deviation from the standard equilibrium thermodynamics [83] on the course of cosmic evolution.
Another interesting observation would be for the radiation dominant era.Here, the equation of state will be  = 1/3, and the Cai and Kim temperature posses no problem, as it is always well defined.However, the temperature defined using the induced flat FLRW metric, equation ( 54), becomes zero.Thus, there is a sign flip for the temperature at  = 1/3.At this point (also known as the Hayward zero temperature divide), r is a null surface that divides the spacelike and timeline regions, and it does not respect the third law of thermodynamics [82].The systematic study of radiation dominated epoch in the concerned HDE is kept for elsewhere.

Effects of interacting dark energy
An important caveat in the above discussion was the absence of matter -dark energy interaction.So far there are no observational evidence which disprove the existence of non-standard interaction at cosmic scale.In fact this gives us the freedom to have scalar field, vector fields etc. as our dark energy candidates.Further, as the nature of interaction is unknown, nothing stops us from considering linear as well as non-linear interaction between dark matter and dark energy, as long as we can explain the late time acceleration and structure formations [84,85].
Here, instead of the standard conservation equations, one has,   + 3 (1 +   )   = , where  is the interaction term, which can be either linear or non-linear functions of densities along with unknown constants.Once there is interactions, it is possible to set  = −1, by which the model reduces to ΛCDM when  = 0 (see [86] for example).Thus the model never crosses the Phantom divide by definition.One can also rule out certain interactions, as it cannot explain the power spectrum [87].One way to test Phantom models is to investigate the density perturbations and check for the stability.This depends on the form of interaction and is subjected to its compatibility with energy conditions [88].Reviews by Wand et.al [89] and Bolotin et.al [90] gives a comprehensive treatment of the subject.A detailed analysis including various interaction will not be in the rest of the article.

A holographic dark energy from entropic function
So far, the holographic principle and the laws of thermodynamics appears far from each other, and there are many inconsistencies in the standard HDE approach.In this section we propose a novel way to model a dark energy density from the entropy correction factors appearing in the modified Friedmann equations.Our approach will respect both holographic bound and the laws of thermodynamics.
Clearly, we get additional terms in the modified Friedmann equations when we consider entropies other than the Bekenstein-Hawking entropy.One can interpret these additional factors as the source of holographic dark energy.Very recently, authors in references [64,65,91] adopted a similar approach to explain the cosmic evolution.They equated the entropic correction factor with the standard HDE density and redefined the IR cutoff.They illustrate the possibility of explaining the cosmic evolution with and without the integration constant, which is essentially the cosmological constant.Here, instead of equating the extra entropic function with the standard HDE density, we assume its proportionality with the CKN bound to define a new HDE density.In the process, we drop the integration constant, which one can recover as a limiting case.
From the laws of thermodynamics (both the Clausius relation and the unified first law), we originally had, This is essentially the first Friedmann equation.Upon integrating the right-hand side, we get the standard ∼  2 term, additional factors, and an integration constant.In the ΛCDM model, this integration constant is responsible for the late time acceleration.Thus, to recover the factor responsible for the cosmic acceleration, we must remove the standard Bekenstein-Hawking factor from the integral.In a more general setting, we argue that one must remove the contribution from the standard area law to extract the explicit contribution from the entropic correction factor.Using this notion, we propose that the dark energy must be of the form, Here, ( − 1)/(16r 2  ) removes the terms coming from the Bekenstein-Hawking entropy.This proposal is dimensionally consistent with the holographic dark energy and gives the cosmological constant in the ΛCDM model when  = r 2  .Besides, the holographic principle is a bound rather than a functional form.Since we remove the contribution from the Bekenstein-Hawking term, the effective functional form in equation ( 56) is well under the allowed CKN bound.Thus respecting the holographic principle and thermodynamics simultaneously.Other HDE models also show ΛCDM-like behaviour; however, such models assumed  = −1 from the beginning, which is a trademark feature of ΛCDM [86,92].Now, using the notion of holography, we define the dark energy density as, For entropies other than Bekenstein-Hawking area law, we get additional factors, and it is then possible to set the integration constant to zero by defining  2 using the conservation relation, for a flat universe ( = 0).Here, Ω  = 8  /(3 2 ) and Ω Λ = 8 Λ /(3 2 ) are the dimensionless density parameters.It is useless if our HDE model cannot explain cosmic evolution.Hence, in the next section, we discuss the cosmic evolution by using the Rényi entropy as the choice.

New holographic dark energy with Rényi entropy
Out of all HDE densities discussed earlier, Rényi entropy was exceptional, and the standard approach to Rényi HDE is not satisfactory.Here, we use our definition of dark energy density, given in equation (57), introduced in the previous section to define an HDE density.Here, we consider the validity of the continuity equation and the first law of thermodynamics.From equation ( 57), we have the definition of new holographic dark energy as Considering a flat 3 + 1 dimensional FLRW universe with Rényi entropy as the choice, we get, With the integration constant Λ, we will get acceleration, where the extra factor acts like a small correction.Here, we set Λ = 0 and explain the cosmic evolution from the correction factor alone.To do so, we use the relation Ω Λ + Ω  = 1, and with r = 1/, the expression for  becomes, Here, Ω 0 and  0 are the present value of the matter density parameter and the Hubble parameter, respectively.Thus the expression connecting Hubble parameter, energy densities, and redshift becomes, Interestingly, in the limit  → 0, we get, ℎ() 2 = Ω 0 ( + 1) 3 + (1 − Ω 0 ), which is precisely the ΛCDM model.Once again, it is not easy to obtain an analytical solution for  ≠ 0; hence we rely on numerical techniques.Figure (7) shows the behaviour of the Hubble pa- Here we took  = −19.366as the nuisance parameter [93] to estimate the apparent magnitude (with Ω 0 = 0.315 and  0 = 67.4km/s/Mpc [78]).Observational Hubble Data from [79] and Pantheon data from [94].
rameter and the apparent magnitude as a function of redshift.Here, the expression for apparent magnitude is, () = 5 log 10 (1 + ) , where  is the nuisance parameter.Clearly, for a very small value of , the model behaves very close to the ΛCDM model.As  increases, the model deviates from the ΛCDM, and for a very high value of , it will approach a final de Sitter universe.

Value of '𝜆':-
A complete data analysis using the available cosmological observation is a rigorous work that needs more attention beyond the scope and length of the current manuscript.However, we expect  to be very small based on the following reasons.
By the definition of Rényi entropy,  must lie between 1 and −∞ to be Schur concave.An even stronger concavity condition from quantum information theory demands  to be between 1 and 0 [71].Positive values greater than one are mathematically allowed due to the skew symmetry of Rényi entropy.Although there is no physical motivation to assume so.Even if we relax the stronger concavity condition and assume negative values, for / 2 < −1, the Rényi entropy becomes complex-valued.Such entropies are of interest in information theory [95,96], and it will be interesting to investigate them in the context of cosmology and black hole physics.In summary, mathematically the uniform prior range of  can be between +∞ to − 2 / if we restrict to real valued entropy, and physically, based on quantum information theory, the uniform prior range with stronger concavity is between 0 and 1.A freehanded range would be − 2 / to 1, thus making  → 0 a special case rather than a limit.Since such a lower bound is time dependent, we did not used this as we already assume  as an unknown constant.A time dependent non-extensive parameter is interesting to investigate and will be done elsewhere.N. Komatsu [80] illustrated a behaviour similar to ours with Rényi entropy in the context of Padmanabhan's law of emergence [81].Although different from our approach, they also preferred a very small positive value for .evolution of dark energy equation of state parameter as a function of  in figure (8).Clearly, for  ≤ 1, the model never crosses the phantom divide.Unlike the previous case of the standard Rényi HDE model, the dark energy density parameter Ω Λ = 8 Λ /(3 2 ) reduces to zero for large , indicating a proper matter dominated early phase.Here we have, For  → 0, the above expression reduces to  2 0 (1 − Ω 0 )/ 2 , which is the ΛCDM behaviour.Further, it is evident from the above expression and the figure (9) that Ω Λ goes to zero, rather than settling in a non-zero value for large .
Now we evaluate the stability of the model by calculating the sound speed square ( 2  ) given as,  For a stable model, the value of  2  lie between the region zero and 1.Values beyond this region indicate instabilities of various kinds.We can see from the figure (10) that the model goes from an unstable state to a final stable de Sitter limit.
To distinguish our model from the ΛCDM model, we use the statefinder diagnostic.This geometric diagnostic tool has two parameters {, } defined as, Here,  is the jerk parameter.For the ΛCDM model, these parameters take a fixed value {, } = {1, 0}, while for a varying dark energy model, it traces a parametric path depending on the evolution of the density parameters, the equation of state parameter and its time derivatives.In terms of scaled Hubble parameter ℎ() =  ()/ 0 , the statefinder parameter {, } takes the form Where,  = ln() and  = 1/(1 + ) is the scale factor.From the analysis, our model is distinguishable from the ΛCDM model.Figure (11) shows the parametric {, } plot for different values of .For  ≠ 0, the model lies in the quintessence region, where  < 1 and  > 0 and goes to ΛCDM in the future.
The jerk parameter is also denoted by ' ' instead of '' in the literature [111].However, we used '' to quickly compare it with the statefinder diagnostics illustrated by V. Sahni et al. [112].In our analysis, for the given parameter space of  ∈ [0, 1], our model can explain the late time acceleration with a quintessence type dark energy, where  is positive by construction.Previously, Luongo [113] showed that a positive jerk parameter can explain the late time acceleration and departure from the ΛCDM model.Their calculations suggest that the cosmological constant can be interpreted as a limiting case of a more general dark energy model.Here in our model, the ΛCDM precisely appears as the limiting case when  → 0, and also in the asymptotic future, when  ≠ 0.

Entropy evolution and the second law of thermodynamics
Finally, we explore the thermodynamic evolution of our new model.The entropy of the universe is the sum of the entropy of the horizon and the entropy of everything inside the horizon.Since the total entropy inside the horizon is very less when compared to the horizon entropy, one may approximate the total entropy as the entropy of the horizon [103].Also, according to the second law of thermodynamics, the entropy of an isolated system must be an increasing function of time.Here, we have Rényi entropy as the horizon entropy, and it reads, For clarity, we plot entropy as a function of the scale factor ().From figure (12), we can see that the entropy tends to saturate at a maximum value as we go forward in time.To make this evident, we calculate the first and second derivatives of entropy with respect to the scale factor.The first derivative reads, Since / is a negative function, as  decreases with respect to time in our model,  is a positive valued function.
Further, as  tends to attain a maximum, the first derivative must tend to zero.These features are evident from the figure (13).To confirm the maximization of horizon entropy, we evaluate the second derivative with respect to the scale factor, Here,  → 0 implies ΛCDM behaviour (in natural units).
and it reads, Here, the last two terms are positive, while the first term is negative for a model where  continuously decreases in the future direction.In the long run, the last two terms dilute away faster, and the first term becomes dominant.Thus, in the long run, the second derivative is negative, and we can confirm that the horizon entropy saturates at a maximum value (See figure (13)).Thus our model is thermodynamically consistent and respects the second law of thermodynamics.

Conclusions
The holographic principle and the laws of thermodynamics are well-established frameworks to explain the evolution of the universe.Since the thermodynamic properties of the cosmic horizon are the fundamental basis of the holographic principle, we explored the possible connection between holographic dark energy and the laws of thermodynamics.First, we explored Moradpour et al.'s idea in [46], where they proposed a holographic dark energy from the principles of thermodynamics.If we were to start from the first law of horizon thermodynamics, we found that one cannot satisfy the proposed proportionality by Moradpour et al. for all choices of entropies.Their proposal appears consistent with the standard holographic approach only when the entropy follows an exponent stretched area law.One can artificially work around this problem by redefining the horizon temperature.Although such a strategy appears inconsistent with the standard laws of thermodynamics and is premature in the current stage of research, there are promising attempts via modified gravity theories [104].Similar to the standard holographic dark energy approach, Moradpour et al.'s proposal assume the validity of the standard Friedmann equations.However, one can only derive the standard Friedmann equations from the first law of thermodynamics if we consider Bekenstein-Hawking entropy as the horizon entropy.Thus, as Golanbari et al. point out in [49], it is crucial to consider modified Friedmann equations when using entropies other than the standard area law.In other words, from the thermodynamic perspective, it is inconsistent with using the standard Friedmann equations and the continuity relation simultaneously for entropies other than the Bekenstein-Hawking entropy.Now, to see if this disparity between the definition of dark energy densities by holography and thermodynamics is due to the assumption of standard Friedmann equations, we derived the expression for dark energy density from the Clausius relation and the unified first law for a dark energy dominated FLRW universe.Instead of assuming the validity of standard Friedmann equations, we took the continuity relation and the first law as the principle equations.Once again, we found it challenging to construct a one-to-one correspondence between the standard HDE and the dark energy stipulated by the first law of thermodynamics.Additionally, although the Hubble horizon is a well-established thermodynamic surface using which one can explain the cosmic evolution from the laws of thermodynamics, in the standard HDE model, the Hubble horizon cannot explain the effects of dark energy unless we invoke some phenomenological interactions.In a recent paper, R. G. Landim [87] points out that the Hubble scale cutoff with certain phenomenological interactions fails to explain the matter and the CMB power spectrums.Thus, the standard holographic and thermodynamics approaches may not go hand in hand.
Further, the standard HDE approach with the Hubble horizon has additional issues.We first discuss the Rényi HDE using the conventional holographic approach to see these problems explicitly.Although the model could explain latetime accelerated expansion, it suffers from inconsistencies due to the assumption of the standard Friedmann equations and the continuity equation simultaneously with the Hubble horizon as the IR cutoff.We saw that the total equation of state parameter has different expressions when considering the validity of standard Friedmann equations and that of the continuity equation.Further, evaluating the dark energy equation of state parameter, we found that it crosses the phantom divide, which may not be a favourable behaviour [105], as we do not consider any interaction between the dark energy and matter.
Additionally, the model cannot explain the matterdominated epoch when we use the continuity equation to define the equation of state parameter.A possible reason might be the non-vanishing behaviour of the dark energy density parameter.We found that for the above Rényi HDE model, the dark energy density parameter settles at a constant value instead of diluting to zero for large values of .This behaviour effectively generates a non-zero pressure, resulting in the wrong equation of state parameter.
On the other hand, in the simplest scenario, it is possible to induce dark energy as an integration constant using the laws of thermodynamics on the Hubble horizon.Thus, the thermodynamic approach does not demand any additional dark energy component.These observations led us to reconsider the standard HDE approach from the laws of thermodynamics.
In the thermodynamic approach, when considering entropies other than the Bekenstein-Hawking entropy, we can set this integration constant to zero by the proper definition of a dark energy component.We remove the contribution from the standard area law from the total energy density to achieve this dark energy density.This definition is con-sistent with the original CKN "bound", as we remove the contribution by the Bekenstein-Hawking term.Interestingly this definition reduces to the cosmological constant in the ΛCDM model when the generalized entropy reduces to the Bekenstein-Hawking entropy.Thus, if dark energy is dynamical, our model gives a simple extension to the ΛCDM model from the laws of thermodynamics, which is consistent with the holographic principle.
We then demonstrate the validity of our model with Rényi entropy as our choice by investigating different cosmological parameters.Our model can explain the late-time acceleration and gives a final de Sitter epoch.In our model, the dark energy equation of state parameter never crosses the Phantom divide, and the dark energy density parameter reduces to zero for large .This behaviour thus guarantees the proper matter-dominated and dark energy-dominated epochs.By analysing the sound speed squared, we found that the model goes to a final stable state, and the statefinder diagnostic tool can distinguish the model from the ΛCDM model.The numerical estimates of various parameters are well within the observational bounds for  ∈ [0, 1].Finally, the model is consistent with the second law of thermodynamics, by which the entropies tend to a maximum towards the future.
In future, we plan to investigate the effects of phenomenological interaction terms in the continuity relation and the ability of the model to explain the observed power spectrums.It is also promising to look for solutions to problems such as the Hubble tension with our model or some extended versions of the same, as several varying dark energy models could give a reasonable explanation to it [106].Whether this demands a phantom or Quintessence behaviour is still an open area of research [107,108].Finally, it will be interesting to look for a possible connection between our approach and other holographic approaches, such as Padmanabhan's law of emergence of cosmic space [81], as the footings of his idea are the laws of thermodynamics.

Figure ( 4
) illustrates the behaviour of  for different values of .Additional to the total equation of state parameter, one can also evaluate the equation of state parameter for the dark energy component using the relation  Λ =  Λ / Λ .In a 3 + 1 dimensional spacetime, from the continuity equation for dark energy density,

10
Late time cosmological behaviourIn this section we study various cosmological parameters to understand the late time behaviour of the proposed dark energy model.We will focus on the behaviour of the deceleration parameter (), dark energy equation of state parameter ( Λ ), dark energy density parameter (Ω Λ ), sound speed squared ( 2  ), and the statefinder parameters {, }.Now, using the expression  = −1+ 1+  ()  ()  , we can evaluate the deceleration parameter.From the figure (8), the model shows late time acceleration around redshift  ∼ 0.5 and tends to a final de Sitter phase for a very small positive value of .For very large values of , the model could reenter the decelerated phase in the future.Thus, the model prefers a very small value of , indicating evolution close to the ΛCDM model.Using the standard definition  Λ =  Λ / Λ , where  Λ = −(  Λ /(3)) −  Λ , we numerically evaluate and plot the

Fig. 11
Fig. 11 Statefinder diagnostic for different values of .The sectors divided by dotted lines about the markers represents the past and the future (with Ω 0 = 0.315 and  0 = 67.4km/s/Mpc [78]).The markers represents the present time for different .

Table 1
Different entropies and corresponding dark energy densities using the standard HDE density in equation (7a) and Moradpour et al.'s proposal in equation (7b).
(33) is the same as the previous result from the Clausius relation.Either way, we do get the same bound for the dark energy density as proposed by Moradpour et al.If we consider equation(33)as our definition for HDE density from thermodynamics, we get expressions similar to  Λ2 given in table(1), except for Rényi entropy.Table(2)summarise the expressions for dark energy density using equation(33).It is again evident that equation(33)generates HDE proportional to the HDE from the CKN bound, for Bekenstein-Hawking, Barrow and Tsallis entropies.However, similar to Moradpour et al.'s results, Rényi entropy gives a different result.

Table 3
[78] support our estimates[99- 102].The following table tabulates various estimates for different values of  based on the model discussed above.Here we used Ω 0 = 0.315 and  0 = 67.4km/s/Mpcfrom the Planck 2018 result[78].Estimates of various parameters for different values of  ∈ [0, 1] with Ω 0 = 0.315 and  0 = 67.4km/s/Mpc[78].Here  is the model parameter,  0 is the present value of deceleration parameter,   is the transition redshift,  Λ0 is the present value of dark energy equation of state and  is the jerk parameter (also the 1  statefinder parameter).