Reexploration of interacting holographic dark energy model: cases of interaction term excluding the Hubble parameter

In this paper, we make a deep analysis for the five typical interacting holographic dark energy models with the interaction terms Q=3βH0ρde\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=3\beta H_{0}\rho _\mathrm{{de}}$$\end{document}, Q=3βH0ρc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=3\beta H_{0}\rho _\mathrm{{c}}$$\end{document}, Q=3βH0(ρde+ρc)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=3\beta H_{0}(\rho _\mathrm{{de}}+\rho _\mathrm{c})$$\end{document}, Q=3βH0ρdeρc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=3\beta H_{0}\sqrt{\rho _\mathrm{{de}}\rho _\mathrm{c}}$$\end{document}, and Q=3βH0ρdeρcρde+ρc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=3\beta H_{0}\frac{\rho _\mathrm{{de}}\rho _{c}}{\rho _\mathrm{{de}}+\rho _\mathrm{c}}$$\end{document}, respectively. We obtain observational constraints on these models by using the type Ia supernova data (the Joint Light-Curve Analysis sample), the cosmic microwave background data (Planck 2015 distance priors), the baryon acoustic oscillations data, and the direct measurement of the Hubble constant. We find that the values of χmin2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _\mathrm{min}^2$$\end{document} for all the five models are almost equal (around 699), indicating that the current observational data equally favor these IHDE models. In addition, a comparison with the cases of an interaction term involving the Hubble parameter H is also made.


Introduction
The accelerated expansion of the universe has been discovered by the type Ia supernova observations [1,2] and further confirmed by various other cosmological observations [3][4][5][6]. Dark energy that has negative pressure has been proposed to explain the phenomenon of cosmic acceleration [7][8][9][10][11][12][13][14][15]. In the present universe, dark energy contributes about 70% of the cosmic energy density, and thus it is now dominating the evolution of the universe. The study of dark energy has become one of the most important issues in theoretical physics and modern cosmology. Although enormous efforts have been made to investigate dark energy, its nature is still in the dark.
The primary candidate of dark energy is the so-called "cosmological constant" (denoted as ), which is equivalent to the density of vacuum energy and thus is a constant in space and time. The cosmological constant has a constant equation-of-state parameter (EoS) w ≡ p /ρ = −1. The a e-mail: zhangxin@mail.neu.edu.cn cosmological model with and cold dark matter (CDM) is called the CDM model, which can explain the current various cosmological observations quite well [16]. However, the cosmological constant always suffers from serious theoretical challenges, i.e., the so-called "fine-tuning" and "cosmic coincidence" puzzles [17,18]. The value of the vacuum energy density calculated by quantum field theory is higher than the fit value of the cosmological constant by cosmological observations by about the 120 orders of magnitude, and so a bare cosmological constant needs to be introduced to make an offset, leading to the fine-tuning problem. The coincidence problem asks why the densities of vacuum energy and matter are in the same order today, although their evolutionary histories differ enormously. These two puzzles have been frustrating the CDM cosmology in theoretical aspect.
Actually, there are many other candidates for dark energy, for which the vast majority believes that dark energy has dynamics, often realized by some scalar field [19][20][21][22][23][24][25][26][27][28][29]. However, more theoretically, it is believed that the dark energy problems are closely related to the theory of quantum gravity in nature. It is actually obvious that has a quantum origin and at the same time it yields repulsive gravity leading the current universe to accelerate. Therefore, it is of great interest to explore the nature of dark energy from the perspective of quantum gravity. In the current circumstance that we have no a complete theory of quantum gravity, we have to appeal to the holographic principle of quantum gravity for an effective theory of dark energy.
Together with the effective quantum field theory, the holographic principle leads to a model of dark energy, named "holographic dark energy" (HDE) model, which can solve the two theoretical puzzles of cosmological constant at the same time to some extent [30] and can explain the current cosmological observations well (see Ref. [31] for a recent review). By considering the holographic principle, it is required that in a spatial region there is an upper limit for the number of degrees of freedom involved in it, due to the gravita-tional effects of them (the condition of black hole formation sets such a bound) [32]. That is to say, in this theory, the infrared (IR) cutoff (with length scale L) is related to the ultraviolet (UV) cutoff (with energy scale k max ). Recall that for the vacuum energy density we have the evaluation ρ vac k 4 max /(16π 2 ). Thus, one finds that by such a theoretical consideration the density of dark energy can be decided by some IR cutoff length scale of the universe. The holographic reasoning gives the density of dark energy of the form [30] where is the reduced Planck mass and c is a dimensionless constant characterizing ambiguous factors in the effective theory. In the HDE model, Li [30] argues that L should be chosen as the future event horizon of the universe, which can lead to a cosmic acceleration. Namely, in the HDE model, we have The HDE model has been widely studied in depth [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52].
There are also some variants of this kind [53][54][55][56][57][58][59][60][61][62]. In the HDE model, it is found that the parameter c solely determines the evolution of dark energy, by solving a differential equation (see, e.g., Ref. [30]). When c > 1, the dark energy has w > −1 (in the case of c = 1, w will eventually evolve to get −1); when c < 1, the EoS of dark energy w will cross the phantom divide −1 from w → −1/3 to w → −1/3−2/(3c) [63]. The cosmological constraints show that c is around 0.7 (see, e.g., Refs. [64][65][66]). A recent work [67] on comparing popular dark energy models shows that the HDE model is still a competitive model in the aspect of fitting the current cosmological observations. On the other hand, there might be some direct interaction between dark energy and dark matter, which is capable of helping resolve (or alleviate) the coincidence problem of dark energy [68][69][70]. Therefore, in addition to probing the dynamics of dark energy, another important mission for the investigation of dark energy is to detect such a "fifth force" between dark energy and dark matter by accurate cosmological observations. The interacting dark energy models have been widely studied . The interacting models in the framework of holographic dark energy have also been deeply explored (see, e.g., Refs. [110][111][112]).
In interacting models of dark energy, one considers that there is an energy transfer between dark energy and dark matter in the background universe (and in a perturbed universe there is also a momentum transfer between them). In a concrete model, a form of the energy (density) transfer rate (denoted as Q) should be assumed. Usually, consulting from the theories of nuclear decay and inflationary reheating, the form of Q is assumed to be proportional to the density of dark energy or dark matter, i.e., Q = 3β Hρ de or Q = 3β Hρ c , where ρ de and ρ c are the densities of dark energy and cold dark matter, respectively, H is the Hubble parameter, and β denotes the dimensionless coupling between dark energy and dark matter. Note that here 3H appears only for mathematical convenience. In our recent work [112], Feng and Zhang explored the interacting models in the framework of holographic dark energy and made a comparison for five interacting cases (Q = 3β Hρ de , Q = 3β Hρ c , Q = 3β H (ρ de +ρ c ), Q = 3β H √ ρ de ρ c , and Q = 3β H ρ de ρ c ρ de +ρ c ), according to the constraint results of current observations. However, in the research area of interacting dark energy, there is another perspective that Q should not involve the Hubble parameter H because the local interaction should not depend on the global expansion of the universe [102]. According to this perspective, one should write down the form of Q as, e.g., Q = 3β H 0 ρ de or Q = 3β H 0 ρ c , where the appearance of H 0 is only for a dimensional consideration. In this paper, we will revisit the exploration of interacting holographic dark energy models by adopting this perspective. We will consider the five cases with We constrain the models by using the current cosmological observations, and we report the results and make an analysis for them. This paper is organized as follows. In Sect. 2, we briefly describe the interacting holographic dark energy model. In Sect. 3, we present the analysis method and the observational data used in this paper. In Sect. 4, we report the constraint results and make a deep discussion for them. Conclusion is given in Sect. 5.

The interacting model of holographic dark energy
In this section, we briefly derive the equations describing the interacting holographic dark energy (IHDE) model for a flat cosmology.
In the context of the flat Friedmann-Roberston-Walker universe, the Friedmann equation can be written as where 3M 2 pl H 2 is the critical density of the universe, ρ c , ρ b , ρ r , and ρ de represent the energy densities of cold dark matter, baryon, radiation, and dark energy, respectively. For convenience, we define the fractional energy densities of various components as By definition, we have In the IHDE model, there is some direct, non-gravitational interaction between holographic dark energy and dark matter, and thus we have the following continuity equations for the various components: where Q is the phenomenological interaction term [113][114][115][116][117][118][119][120][121][122], denoting the energy transfer rate between dark energy and dark matter. In this paper, we consider the following five cases in the IHDE model: As has been mentioned above, β is a dimensionless coupling parameter describing the strength of interaction between dark energy and dark matter. Combining Eqs. (3) and (5)- (9), we obtain Substituting Eq. (15) into Eq. (7), we obtain Here, for convenience, following Ref. [110] we define From Eq. (1) (i.e., the definition of density of holographic dark energy), we can get the relation We now write the IR cutoff length scale L as the form Combining Eqs. (18) and (19), we get Combining Eqs. (2) and (20), we have the relation Taking the derivative of Eq. (21) with respect to t, we can get the equatioṅ Combining Eqs. (16) and (22), we get the following two differential equations governing the dynamical evolution of dark energy in the IHDE model for a flat cosmology: where E(z) = H (z)/H 0 is the dimensionless Hubble expansion rate, de (z) is the fractional density of dark energy, and In this paper, for convenience, we occasionally call the cases with Q 1 -Q 5 [described by Eqs. (10)- (14)] the IHDE1-IHDE5 models, respectively.

Method and data
In a flat universe, the IHDE models have four free parameters, c, h, m0 , and β. We will use the current observational data to constrain the models.
We use the χ 2 statistic to estimate the model parameters. The form of χ 2 function is as follows: where ξ th is the theoretically predicted value for the observable ξ , ξ obs is the corresponding experimentally measured value, and σ ξ is the standard deviation. The total χ 2 is the sum of all χ 2 ξ , The observational data we use in this paper include the type Ia supernova (SN) data, the cosmic microwave background (CMB) anisotropy data, the baryon acoustic oscilla-tion (BAO) data, and the direct measurement of the Hubble constant H 0 . Thus the total χ 2 function is Since these IHDE models have the same parameter number, we can direct compare them with their χ 2 values. But when we compare them with the CDM model and the HDE model (without interaction), the χ 2 comparison becomes unfair because their parameter numbers are different. We thus employ the Akaike information criterion (AIC) [123] and the Bayesian information criterion (BIC) [124] to do the model comparison in this situation. By definition, we have AIC = −2 ln L max + 2k and BIC = −2 ln L max + k ln N , where k is the number of parameters, and N is the number of data points. Since we wish to measure the difference between models, we are more interested in the relative values of them. In this work, we choose the CDM model as a reference model, and then calculate AIC = χ 2 min + 2 k and BIC = χ 2 min + k ln N . A model with a lower value of AIC or BIC is believed to be more favored by data.

Type Ia supernovae
We use the JLA compilation of type Ia supernovae [125] in this work. It is from a joint analysis of type Ia supernova observations. The JLA compilation consists of 740 Ia supernovae data points, obtained by the SDSS-II and SNLS collaborations. The distance modulus of a SN Ia iŝ where m * B is the observed peak magnitude, M B is the absolute magnitude, X 1 is the time stretching of the light curve, and C is the supernova color at maximum brightness. The luminosity distance d L of a supernova in a spatially flat FRW universe is defined as where z and z hel are the CMB frame and heliocentric redshifts, respectively. The χ 2 function for the JLA SN observation is where C SN is the covariance matrix of the JLA SN observation and μ th denotes the theoretical distance modulus, defined as

Cosmic microwave background
For the CMB data, we use the "Planck distance priors" from the Planck 2015 data [126]. The distance priors contain the shift parameter R, the "acoustic scale" A , and the baryon density ω b ≡ b0 h 2 . R and A are defined as where m0 is the present-day fractional energy density of matter, D A (z * ) is the proper angular diameter distance at the redshift of the decoupling epoch of photons z * , and r s (z * ) is the comoving size of the sound horizon at z * . In a flat universe, D A can be expressed as and r s (z) can be expressed as where b0 and γ 0 are the present-day fractional energy densities of baryons and photons, respectively. From the measurement of CMB, we have 3 b0 /4 γ 0 = 31500 b0 h 2 (T cmb /2.7 K) −4 , where T cmb = 2.7255 K. The fitting formula of z * is given by [127] z where Using the Planck TT + LowP data, the values of the three quantities are obtained: R = 1.7488 ± 0.0074, A = 301.76 ± 0.14, and b h 2 = 0.02228 ± 0.00023. The inverse covariance matrix for them, Cov −1 CMB , can be found in Ref. [126], The χ 2 function for CMB is thus given by where p 1 = A , p 2 = R, and p 3 = ω b .

Baryon acoustic oscillations
The BAO data can be used to measure the angular diameter distance D A (z) and the expansion rate of the universe H (z). The BAO measurements can provide the ratio of the effective distance measure D V (z) and the comoving sound horizon size r s (z d ) for us [i.e., ξ(z) = D V (z)/r s (z d )]. The expression of D V (z) from the spherical average is where where We use four BAO points from the six-degree-field galaxy survey (6dFGS) at z eff = 0.106 [129], the SDSS main galaxy sample (MGS) at z eff = 0.15 [130], the baryon oscillation spectroscopic survey (BOSS) "LOWZ" at z eff = 0.32 [131], and the BOSS CMASS at z eff = 0.57 [131]. The χ 2 function for BAO is given by

The Hubble constant
The direct measurement of the Hubble constant we use in this work is given by Efstathiou [132], H 0 = 70.6 ± 3.3 km s −1 Mpc −1 . It is a re-analysis of the Cepheid data of Riess et al. [133]. The χ 2 function of the Hubble constant measurement is

Results and discussion
In this section, we report the fitting results of the IHDE models and discuss the implications of these results. We use the observational data combination SN + CMB + BAO + H 0 to constrain the models (the CDM model, the HDE model, and the IHDE1-5 models). The fitting results are summarized in Tables 1 and 2.
In Table 1, we give the values of χ 2 min , AIC, and BIC for these models. We find that, among these models, the CDM model is still the best one in fitting the current observational data. The CDM model has the least number of parameters, but it gets the smallest χ 2 min value in this fit. The HDE model has one more parameter than the CDM model, but it yields a greater χ 2 min value, by χ 2 ∼ 5. The IHDE models have two more parameters than the CDM model, but they only yield similar χ 2 min values (all around 699) to that of CDM. So, although the CDM model has been facing the severe theoretical problems, it is the simplest dark energy theoretical model and can explain the observations best. The HDE model indeed can provide an attractive theoretical scheme for avoiding the cosmological constant problems, but it performs worse than the CDM model in fitting  Fig. 1 Graphical representation of the results of AIC and BIC for the HDE model and the IHDE models the observational data. It seems that one should explore more possible factors in the framework of holographic dark energy for seeking for a better model of dark energy in the sense of quantum gravity.

ΔAIC ΔBIC
Using the information criteria to assess the models (the CDM model is selected as a reference model in this work), we can see that the HDE model has AIC = 7.2 and BIC = 11.8, and the IHDE1-5 models have AIC ∼ 4.3 and BIC ∼ 13.6. The comparison of HDE and IHDE shows that, only considering the factor of number of parameters (i.e., AIC), the IHDE performs better, but when further considering the factor of number of data points (i.e., BIC), the HDE performs better. We show the graphical representation of the results of AIC and BIC for the HDE model and the IHDE models in Fig. 1. We also find that the values of χ 2 min (also, AIC and BIC) for all the five IHDE models are almost equal, indicating that the current observational data equally favor these IHDE models (see also Fig. 1).
In the previous study [112], Feng and Zhang investigated the IHDE models with interaction terms involving the Hubble parameter H (also Q1-Q5, but with the form like Q = 3β Hρ). We would like to make a comparison of our results in the present paper with those in Ref. [112]. We will occasionally use the names like Q = 3β H 0 ρ models and Q = 3β Hρ models to distinguish the models in this paper and those in Ref. [112]. It was shown in Ref. [112] that, for the Q = 3β Hρ models, according to the same data sets as this work, the IHDE5 model is the best one, the IHDE1 model is the next best one, and the IHDE2 model is the worst one. Namely, in the framework of holographic dark energy with Q = 3β Hρ, the Q = 3β H ρ de ρ c ρ de +ρ c model is most favored by the current data, the Q = 3β Hρ de model is also a good model in the sense of fitting data, and the Q = 3β Hρ c model is relatively not favored by the current data. However, in the present work, we find that all the IHDE models with Q = 3β H 0 ρ are equally favored by the current data.
From Table 2, we find that the fitting values of c in both HDE and IHDE cases are all around 0.7. For the HDE model, we obtain c = 0.73 (the best-fit value); and for the IHDE models, we obtain c = 0.69-0.71 (the best-fit values). The c value in the IHDE models is slightly smaller than that in the HDE model. For the coupling parameter β in the IHDE models, we find that in all the cases β > 0 is favored at the more than 1σ level, indicating that the decay of dark energy into cold dark matter is favored at the more than 1σ statistical significance by the current data. For the IHDE1, IHDE2, and IHDE4 models, we have β ∼ 0.03 and σ β ∼ 0.02; for the IHDE3 model, we have β ∼ 0.01 and σ β ∼ 0.01; and for the IHDE5 model, we have β ∼ 0.06 and σ β ∼ 0.04-0.06.
In Figs. 2, 3 and 4, we show the 1σ and 2σ confidence level contours in the m0 -c, m0 -β, and c-β planes, respectively. The blue contours and the pink contours correspond to the Q = 3β H 0 ρ models (in the present work) and the Q = 3β Hρ models (in Ref. [112]), respectively. In Fig. 2, we find that the Q = 3β H 0 ρ models systematically move towards lower left, relative to the Q = 3β Hρ models, in the m0 -c plane, indicating that for the Q = 3β H 0 ρ models both m0 and c are smaller. In Figs. 3 and 4, we find that, for the IHDE1 and IHDE5 cases, the β value in the Q = 3β H 0 ρ models, relative to the Q = 3β Hρ models, is evidently smaller. In this work, we find that β > 0 is favored at more than 1σ level but less than 2σ level, for all the cases. But for the Q = 3β Hρ models investigated in Ref. [112], it was shown that the interaction between dark energy and dark matter can be detected at more than 2σ significance; for example, for the IHDE1 model, β > 0 is favored at the 2.3σ level, and for the IHDE5 model, β > 0 is favored at the 2.1σ level.
From Fig. 4, we find that for all the IHDE models β and c are in positive correlation. A positive β means that dark energy decays into dark matter, and thus the result of β > 0 will affect the parameter estimation of c. The positive correlation between β and c implies that the decay of dark energy into dark matter will decrease the possibility of a big rip in a finite future.
In Fig. 5, we show the one-dimensional marginalized posterior distributions of β (left panel) and c (right panel) for the HDE model and the IHDE models from the current observations. The blue dashed lines in the figure denote the cases of β = 0 (left panel) and c = 1 (right panel). We can see from the right panel that the c values for the IHDE models are almost the same and slightly smaller than the c value of the HDE model. But, although the c value of the IHDE models is smaller, thanks to the positive β, the risk of a big rip in the IHDE models, compared to the HDE model, is still decreased. (Note that, in the model of holographic dark energy, c < 1 will lead to a late-time phantom and thus a big rip in the finite future; see, e.g., Ref. [134].) In Fig. 6  we show the reconstructed evolution of w (with 1-3σ errors) for the HDE model and the IHDE models. The red dashed line in the figure denotes the cosmological constant boundary w = −1. We can clearly see that, compared to the HDE model, the risk in the IHDE models is indeed decreased.

Conclusion
In this paper, we have investigated the interacting holographic dark energy models in which the interaction term Q does not involve the Hubble parameter H . We consider five typical IHDE models with the interaction terms Q = 3β H 0 ρ de , Q = 3β H 0 ρ c , Q = 3β H 0 (ρ de +ρ c ), Q = 3β H 0 √ ρ de ρ c , and Q = 3β H 0 ρ de ρ c ρ de +ρ c , respectively. We use the current observational data, including SN (JLA) data, CMB (Planck 2015 distance priors) data, BAO data, and H 0 measurement, to constrain these models.
We find that the current observational data equally favor these IHDE models. We also find that in all the cases the coupling parameter β > 0 is favored at more than 1σ level (but less than 2σ level), indicating that the current observations slightly favor the decay of dark energy into dark matter in the current framework of IHDE (with Q excluding H ).
We have made a comparison of our results in the present work with those in the previous work [112] in which the IHDE models with Q involving H are investigated. In Ref.  [112], it was shown that the IHDE5 (Q = 3β H ρ de ρ c ρ de +ρ c ) model and the IHDE1 (Q = 3β Hρ de ) model are most favored by the current data, and the IHDE2 (Q = 3β Hρ c ) model is relatively not favored by the current data; and, in some cases, the coupling of β > 0 can be detected at more than 2σ level (e.g., 2.3σ in the IHDE1 model and 2.1σ in the IHDE5 model). The comparison is shown in Figs. 2, 3 and 4. We also show that for all the cases β and c are in positive correlation, and as a result in the IHDE models the risk of a big rip is decreased, compared to the HDE model.