The model of holographic dark energy in which dark energy interacts with dark matter is investigated in this paper. In particular, we consider the interacting holographic dark energy model in the context of a perturbed universe, which was never investigated in the literature. To avoid the large-scale instability problem in the interacting dark energy cosmology, we employ the generalized version of the parameterized post-Friedmann approach to treat the dark energy perturbations in the model. We use the current observational data to constrain the model. Since the cosmological perturbations are considered in the model, we can then employ the redshift-space distortions (RSD) measurements to constrain the model, in addition to the use of the measurements of expansion history, which was either never done in the literature. We find that, for both the cases with $Q=\beta H\rho_{\rm c}$ and $Q=\beta H_0\rho_{\rm c}$, the interacting holographic dark energy model is more favored by the current data, compared to the holographic dark energy model without interaction. It is also found that, with the help of the RSD data, a positive coupling $\beta$ can be detected at the $2.95\sigma$ statistical significance for the case of $Q=\beta H_0\rho_{\rm c}$.

We apply the statefinder hierarchy plus the fractional growth parameter to explore the extended Ricci dark energy (ERDE) model, in which there are two independent coefficients $\alpha$ and $\beta$. By adjusting them, we plot evolution trajectories of some typical parameters, including Hubble expansion rate $E$, deceleration parameter $q$, the third and fourth order hierarchy $S_3^{(1)}$ and $S_4^{(1)}$ and fractional growth parameter $\epsilon$, respectively, as well as several combinations of them. For the case of variable $\alpha$ and constant $\beta$, in the low-redshift region the evolution trajectories of $E$ are in high degeneracy and that of $q$ separate somewhat. However, the $\Lambda$CDM model is confounded with ERDE in both of these two cases. $S_3^{(1)}$ and $S_4^{(1)}$, especially the former, perform much better. They can differentiate well only varieties of cases within ERDE except $\Lambda$CDM in the low-redshift region. For high-redshift region, combinations $\{S_n^{(1)},\epsilon\}$ can break the degeneracy. Both of $\{S_3^{(1)},\epsilon\}$ and $\{S_4^{(1)},\epsilon\}$ have the ability to discriminate ERDE with $\alpha=1$ from $\Lambda$CDM, of which the degeneracy cannot be broken by all the before-mentioned parameters. For the case of variable $\beta$ and constant $\alpha$, $S_3^{(1)}(z)$ and $S_4^{(1)}(z)$ can only discriminate ERDE from $\Lambda$CDM. Nothing but pairs $\{S_3^{(1)},\epsilon\}$ and $\{S_4^{(1)},\epsilon\}$ can discriminate not only within ERDE but also ERDE from $\Lambda$CDM. Finally we find that $S_3^{(1)}$ is surprisingly a better choice to discriminate within ERDE itself, and ERDE from $\Lambda$CDM as well, rather than $S_4^{(1)}$.

We apply the statefinder diagnostic to the extended holographic Ricci dark energy (ERDE) model without and with interaction to study their behaviors. We plot the trajectories of various parameters for different cases. It is shown that the non-interacting model does not reach the LCDM point $\{1,0\}$ and the interacting one is favored, because the interaction makes the evolution of the statefinder pair $\{r,s\}$ quite different.

So far, there have been no theories or observational data that deny the presence of interaction between dark energy and dark matter. We extend naturally the holographic dark energy (HDE) model, proposed by Granda and Oliveros, in which the dark energy density includes not only the square of the Hubble scale, but also the time derivative of the Hubble scale to the case with interaction and the analytic forms for the cosmic parameters are obtained under the specific boundary conditions. The various behaviors concerning the cosmic expansion depend on the introduced numerical parameters which are also constrained. The more general interacting model inherits the features of the previous ones of HDE, keeping the consistency of the theory.

We employ the metric of Schwarzschild space surrounded by quintessential matter to study the trajectories of test masses on the motion of a binary system. The results, which are obtained through the gradually approximate approach, can be used to search for dark energy via the difference of the azimuth angle of the pericenter. The classification of the motion is discussed.