Abstract

In this work, we construct an interacting model of the Rényi holographic dark energy in the Brans-Dicke theory of gravity using Rényi entropy in a spatially flat Friedmann-Lemaître-Robertson-Walker Universe considering the infrared cut-off as the Hubble horizon. In this setup, we then study the evolutionary history of some important cosmological parameters, in particular, deceleration parameter, Hubble parameter, equation of state parameter, and Rényi holographic dark energy density parameter in both nonflat Universe and flat Universe scenarios and also observe satisfactory behaviors of these parameters in the model. We find that during the evolution, the present model can give rise to a late-time accelerated expansion phase for the Universe preceded by a decelerated expansion phase for both flat and nonflat cases. Moreover, we obtain as , which indicates that this model behaves like the cosmological constant at the future. The stability analysis for the distinct estimations of the Rényi parameter and coupling coefficient has been analyzed. The results indicate that the model is stable at the late time.

1. Introduction

The Virial theorem (1930s) which provided the Coma galaxy cluster mass [1, 2], accompanied by the galaxy rotation curve study (1970) [3] and the two different research groups’ observational results in the 1990s [4, 5], has uncovered one of the most interesting issues of cosmology at present: the dark sector. It is suggested by the researchers that the five percent of the present energy content of the cosmos is composed of the radiation and the ordinary matter (baryons); the remaining ninety-five percent is dominated by this dark component to clarify the late accelerated expansion of the Universe. It is believed that this dark sector of the Universe mainly includes two constituents: dark energy (DE) and dark matter (DM). Both are important and significant to understand the phenomena of scales and nature. The significance of DM lies primarily in the structure formation, for instance, to permit baryonic structures to become nonlinear in the wake of decoupling from the photons. Interestingly, dark energy is the subject of study to answer the late-time accelerated expansion for the observable Universe [6]. Also, the DM is narrated as cold dark matter (CDM), and dark energy is portrayed by the cosmological constant () in the standard cosmological scenario. The dark component of the Universe with radiation and baryons combined the CDM model. Also, despite the fact that the CDM model appreciates an impressive observational achievement [79], there are still a number of hypothetical and observational focuses that have the right to be completely researched [10]. The greatest test lies in understanding the crucial idea of these dark sectors from the theoretical perspective [6]. In 2004, Li [11] proposed the idea of holographic dark energy (HDE) which is also used to explain the DE scenario to explain the late-time accelerated expansion of the Universe inspired by the holographic principle [1218]. Right after a paper by Li, the most complete generalization which includes all known HDE models were suggested [19]. Furthermore, it is shown that the Nojiri-Odintsov HDE describes also covariant theories unlike Li’s HDE [20].

Recently, inspired by the holographic principle and using the Rényi entropy [21], a new dark energy model has been proposed by Moradpour et al. [22] named the Rényi holographic dark energy (RHDE) model for the cosmological and gravitational investigations. Generalizing one of the entropy or gravity, as entropy-area connection relies on the gravity hypothesis, will change the corresponding one. It is proposed that by using the Rényi entropy, the modified Friedmann equations can be obtained [2325]. Ghaffari et al. [26] proposed that inflation may be found in the Rényi formalism. The RHDE models have been explored with IR cut-off as the particle and future event horizons [27]. The spatially homogeneous and anisotropic Bianchi Universe filled with RHDE with Granda-Oliveros and Hubble horizons as the IR cut-off has been investigated in general relativity [28]. Recently, Sharma et al. [29, 30] discriminated the RHDE model from the CDM model by using different diagnostic tools such as statefinder diagnostic and statefinder hierarchy in ample details. Also, the RHDE model has been compared with the holographic and Tsallis holographic dark energy through the statefinder diagnostic tool [31].

Indeed, all the above attempts claim that, at least mathematically, the DE density profile introduced under the shadow of the RHDE hypothesis has considerable potential for modelling the DE behavior, and thus, more studies on this density profile are motivated. Further, at large scales, the models presenting interaction fare well when confronted with observational outcomes from the CMB [32] and matter distribution [33]. Therefore, the interaction between DE and DM must be handled seriously. Then again, there exist limits for the quality of this association for different setups [3448]. This newly proposed Rényi HDE has also been examined by many researchers by considering the interaction between DE and DM to explain the accelerated expansion of the Universe with different IR cut-offs in general relativity, braneworld, loop quantum cosmology, and modified gravity [4952]. Sharma and Dubey [53] investigated the Rényi HDE model in the Friedmann-Lemaître-Robertson-Walker (FLRW) Universe considering different parametrizations of the interaction between the DM and DE.

On the other side, the modified gravity theories have been broadly applied in cosmology [5456]. The modified theories of gravity are not new and have a long history. A well-known modified gravity theory, the Brans-Dicke gravity [57], is also a choice to general relativity to explain the accelerated expansion of the cosmos [58] and also can pass the experimental tests from the solar system [59]. Theoretically, the value of the Brans-Dicke coupling parameter has a smaller value than observed by the observational data, which encouraged physicists to suggest various DE scenarios to describe the present cosmos in the Brans-Dicke formalism [5860]. Using the holographic principle, Gong [61] investigated the holographic bound in the Jordan and Einstein frames to the Brans-Dicke gravity, and for larger , the similar results were proposed as those in general relativity. The similar problem was studied in [62], by considering Bianchi identity as a consistency condition. For the IR cut-off as a future event horizon, the importance of Brans-Dicke gravity for the dust matter and the HDE has been explored in [63]. It is proposed that with the Hubble radius as infrared cut-off, the standard HDE may not produce the accelerated expansion Universe in the Brans-Dicke gravity, but the suitable description of the Universe can be obtained by taking the IR cut-off as a future event horizon [64]. Therefore, many other works proposed that the Brans-Dicke gravity is suitable for the examination in the holographic dark energy scenario [6571]. Observational constraints also have been proposed for a sign-changeable interaction among the Universe sectors [7274]. Considering different IR cut-offs, the noninteracting and interacting Tsallis HDE and their cosmological consequences are explored in the Brans-Dicke theory [7577].

Very recently, the authors constructed the noninteracting RHDE model in the Brans-Dicke theory taking the Hubble horizon as the IR cut-off [78]. While in this work, we propose the interacting RHDE model in the framework of the Brans-Dicke theory in both flat and nonflat Universes. The paper is organized as follows: we explored the interacting RHDE model and physical parameters in the Brans-Dicke theory in Section 2. We study the stability of the RHDE model in Section 3. The conclusion is given in the last section.

Throughout the text, an “overdot” represents a derivative with respect to cosmic time.

2. Interacting Rényi Holographic Dark Energy in the Brans-Dicke Theory

We consider a homogeneous and isotropic FLRW Universe which is described by the line element where is the scale factor of the Universe, is the cosmic time, and the curvature constant corresponds to closed, flat, and open Universes, respectively. The coordinates , , and are known as comoving coordinates.

In BD theory, the action is given by [57, 79] where is the BD scalar field, is the Ricci scalar, is the BD parameter, and is the Lagrangian matter. Here, the gravitational constant () takes the place of the time-dependent scalar field , which is inversely proportional to , i.e., . If we assume the matter field to consist of a perfect fluid, then the BD field equations from the variation of action (2) and for the FLRW space-time are obtained as [79] where is the Hubble parameter, is the matter energy density, is the RHDE density, and is the RHDE pressure. The BD scalar field evolution equation is

2.1. Rényi Entropy and HDE

It is important to mention here that it seems there is a deep connection between quantum gravity and generalized entropy scenarios, and indeed, quantum aspects of gravity may also be considered as another motivation for considering generalized entropies [22, 80]. Tsallis entropy is one of the generalized entropy measures which lead to acceptable results in the gravitational and different cosmological setups [23, 25, 43, 8189]. Usually, Tsallis entropy is defined as [86] for a system consisting of discrete states. In the above equation, is the ordinary probability of accessing state , and is a real parameter which may be originated from the nonextensive features of the system such as the long-range nature of gravity [22, 83, 86]. In fact, the concept of nonextensivity is more complex than that of nonadditivity [87]. For example, the well-known Bekenstein entropy is nonadditive and nonextensive simultaneously (for details, see Refs. [84, 85]). It is proposed recently that the Bekenstein entropy (, where and is the IR cut-off) is actually a Tsallis entropy leading to for the Rényi entropy content of the system [21, 22]. Here, is a free parameter and known in the current literature as the real nonextensive parameter that quantifies the degree of nonextensibility [22, 86, 87]. It is proposed in [90] that the parameter affects the energy balance of the Universe. When , the gravitational field is strong enough in such a way that we need only a small quantity of DE and DM to construct the observable Universe. On the other hand, when , the gravitational field is weak in such a way that we need, contrarily to the case, a larger quantity of DE and DM. To sum up, implies less DE and implies more DE than we would have if we consider the standard Boltzmann-Gibbs scenario [90, 91]. In [22], the authors used the value of from to . There are wide ranges for which can produce desired results, while we have taken the values of from to . Authors investigated late-time acceleration for a spatially flat dust filled Universe in the Brans-Dicke theory in the presence of a positive cosmological constant , where the value for the Brans-Dicke-coupling constant is taken as [92]. Authors have studied the Tsallis holographic dark energy in the Brans-Dicke framework using and [77]. The primary focus is in [93] on the FLRW Universe specified by WMAP data. The role of dark energy is played by the vacuum energy density in this model, that is, one had . With the assumption [22] and (i.e., Hubble horizon) and using Equation (7), we obtained energy density for the RHDE as where is a numeric constant. We used and ; relations to get this equation corroborated in a flat FLRW Universe [94]. One can get without , which is in complete agreement with the standard HDE [1418]. It deserves mention here that the apparent horizon is a proper causal boundary for the cosmos in agreement with the thermodynamics laws. Besides, in a flat FLRW Universe, Friedmann equations indicate that whenever DE is dominant in the cosmos, its energy density will scale with (for details, see [22, 94]). Therefore, from a thermodynamic point of view, a HDE model in the flat FLRW Universe, for which the radii of the apparent horizon and the Hubble horizon are the same, will be more compatible with the thermodynamics laws, if it can provide a proper description for the Universe by using the Hubble horizon as its IR cut-off. Following [68], we assume that , i.e., the power law of scale factor in this case to the BD scalar field . One can now easily obtain and hence,

The Rényi HDE density with the Hubble horizon as the IR cut-off is given as

Here, the holographic principle [17] is used, and the effective gravitational constant is given by . The gravitational constant may be found from as a limit. The RHDE energy density can be recovered in the fundamental cosmology [22]. The Holographic DE can also be found in the Brans-Dicke gravity for the case [64]. The dimensionless density parameters are defined as

Our main goal of this work is to build a cosmological model of late acceleration based on the BD theory of gravity and on the assumption that the RHDE and the pressureless dark matter do not conserve separately. Therefore, we assume that both components—the RHDE and the pressureless matter—interact with each other, i.e., one component may grow at the expense of the other. Hence, the energy conservation equations for them are given as follows: and where represents the Rényi HDE equation of state (EoS) parameter and denotes the interaction term. Clearly, for (), there is an energy flow from pressureless matter (RHDE) to RHDE (pressureless matter). We assume the form of interaction as [41, 42, 74], in which is the coupling constant and denotes the deceleration parameter. Here, the main ingredient is the deceleration parameter in the interaction term , and hence, can change its sign when the expansion of our Universe changes from the early decelerated () phase to the late-time accelerated () phase. So, the above interacting term deserves further investigation in the present context. Now, taking derivatives with respect to time of Equation (11), we get combined with relation to obtain where the prime denotes derivative with . Now, taking the derivative with respect to time of Equation (3) and substituting the value of , , , and from Equations (9), (10), (14), and (15), respectively, we get

Defining, as usual, the deceleration parameter as and using Equation (17), we obtain

The evolutionary behavior of the deceleration parameter is plotted for the interacting Rényi HDE model versus redshift by finding its numerical solution using the initial values and , for both flat Universe and nonflat Universe . It is proposed by different observations that the Universe is in an accelerated expansion phase, and the value of the deceleration parameter lies in the range . Also, we have used [77] for all plots. All the physical parameters are examined through and coupling coefficient because they play a crucial role in the evolution of dynamical parameters of the RHDE. From Figure 1, we see the evolutionary behavior of for interacting RHDE in BD gravity, for distinct estimations of and in the nonflat Universe (lower two panels) and flat Universe (upper two panels). We can observe from Figure 1 that the RHDE model shows the transition from an early decelerated stage to a current accelerated stage for both cases for distinct estimations of and parameter . In this context, it is worthwhile mentioning that the standard HDE in the framework of BD theory can explain the accelerated expansion if the event horizon is taken as the role of the IR cut-off [64]. Such a scenario also predicts no acceleration if the Hubble horizon is considered as the IR cut-off. Therefore, the novelty of the present work is that it can explain the current accelerated phase of the Universe if we choose the IR cut-off to be the Hubble horizon.

Combining Equations (13), (15), and (17) with each other, the EoS parameter is obtained as

We have graphed the behavior of EoS parameter of our derived interacting RHDE model for both (two upper panels) and (two lower panels) cases, in Figure 2 for distinct values of parameter and coupling coefficient . According to this figure, it can be seen that of the RHDE model varies from quintessence to the phantom region (). Moreover, we can observe that the EoS parameter approaches CDM model () for all values of and in the future, which is in agreement with cosmological observations. We also noted that the evolution of the EoS parameter at an early time in both flat and nonflat Universes is more distinct for different values of in comparison to coupling coefficient .

By putting Equation (17) in Equation (16), we also obtain the evolution of dimensionless RHDE density parameter as

We have shown the behavior of interacting RHDE density parameter in Figure 3 for both (two upper panels) and , (two lower panels) cases, for distinct values of coupling coefficient and . The thermal history of the Universe, in particular, the successive sequence of matter and DE era, can be observed from these figures for different values of and in both nonflat and flat Universes. We also observed that the RHDE density parameter is consistent with cosmological observations [7], and our results are consistent.

We have plotted the behavior of the Hubble parameter of our derived interacting RHDE model for both the (two upper panels) and , (two lower panels) cases, in Figure 4 for distinct values of parameter and coupling coefficient . It depicts that the variation of affects the behavior of , while different values of coupling coefficient do not affect it. The value of decreases and approaches to a positive value near in the far future.

3. Stability

In this section, we shall discuss the stability of the interacting RHDE model through the squared sound speed in both flat and nonflat Universes. The (the real value of speed), shows a regular propagating mode for a density perturbation. For , the perturbation becomes an irregular wave equation. Hence the negative squared speed (imaginary value of speed) shows an exponentially growing mode for a density perturbation. That is, an increasing density perturbation induces a lowering pressure, supporting the emergence of instability [95].

The squared sound speed is given as [96, 97]

Now, inserting the result of Equation (15) in Equation (22), we get

We solve the Equation (23) numerically by using Mathematica package NDSolve and plotted the squared sound speed versus redshift in Figure 5, for both (two upper panels) and , (two lower panels) cases for distinct values of the parameter and coupling coefficient .From Figure 5, we observe that the RHDE model is not stable initially by taking different values of and in both flat and nonflat Universes, while for in both flat and nonflat Universes, the value of the squared sound speed diverges. By taking different values of in both flat and nonflat Universes, the RHDE model becomes stable at the late time. By analyzing all these plots, we can say that values of and have qualitative effects on the nature of the squared sound speed in both the nonflat and the flat cosmos. The inset plot of Figure 5 shows a close-up of the outer plot around in which the difference can be seen. They are not exactly identical but the difference is very small.

4. Concluding Remarks

In this work, we explored the role of the interacting FLRW cosmos to model dark energy in the Brans-Dicke theory framework using RHDE by taking an infrared cut-off as the Hubble radius in both nonflat and flat Universes. The pressureless matter is assumed to interact with the RHDE through a sign-changeable interaction. In this analysis, we have used the initial values , , , and [77] for both flat Universe () and nonflat Universe (). It has been found that for different values of the Rényi parameter and the coupling coefficient , the interacting RHDE model produces the suitable behavior for the deceleration parameter (), the EoS parameter (), the RHDE density parameter (), and the Hubble parameter, in both the cases (see Figures 14). The effect of different values of and is only quantitative on these parameters.

As discussed earlier, the Brans-Dicke theory in the framework of HDE can explain the accelerated expansion if we choose the IR cut-off to be the event horizon. The theory also predicts no acceleration if we choose the IR cut-off to be the Hubble horizon. However, in our case, the deceleration parameter shows a smooth transition from the decelerated phase () early to the accelerated phase () at a later time. Hence, a remarkable feature of this model is that RHDE in the framework of the Brans-Dicke theory explains the accelerated expansion if we choose the IR cut-off to be the Hubble horizon. It has also been found that the EoS parameter varies from a quintessence () to the phantom region (), and the RHDE model transits decelerating to an accelerating stage of the Universe and approach to as , which implies that the RHDE model imitates the cosmological constant at a far future. It is observed that the RHDE density parameter becomes as . Moreover, it is observed that the variation of affects the behavior of the Hubble parameter , while different values of do not affect it. Also, the value of decreases and approaches to a value near in the far future. Furthermore, we have investigated the classical stability of our model by analyzing the squared sound speed . It has been found that the stability of our model crucially depends on the choices of the parameter in both flat and nonflat Universes (see Figure 5).

As we showed, the present model exhibits more interesting phenomenology comparing to the standard scenario, and hence, it can be a candidate for the description of nature. In a follow-up study, we would like to perform an observational analysis to constrain the parameter .

Data Availability

This manuscript has no associated data or the data will not be deposited. (Authors’ comment: data sharing is not applicable to this article as no new data were created or analyzed in this study).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

VCD and UKS are thankful to GLA University, India, for providing the support and help to carry out this research work.