The main goal of galaxy surveys is to map the distribution of the galaxies, for the purpose of understanding the properties of this distribution and its implications for the content and evolution of the universe. However, in order to realise the potential of these surveys, we need to ensure that we are using the correct analysis: the relativistic analysis, which has been widely studied recently. In this work, the known relativistic overdensity of galaxy surveys is re-examined. Subtle, yet crucial parameters which appear to have been missed by previous works are uncovered. The possible implication of these parameters on the observed galaxy power spectrum is demonstrated for a generic survey, in the cosmological concordance model. The results show that these parameters can alter the predictions of galaxy clustering on all scales; hence, the relativistic effects (on ultra-large scales). In particular, the results show that ignoring these parameters - as in the analysis in previous works - will lead to a false overestimation of the galaxy clustering strength in structure formation.

The beyond-Horndeski gravity has recently been reformulated in the dark energy paradigm - which has been dubbed, Unified Dark Energy (UDE). The evolution equations for the given UDE appear to correspond to a non-conservative dark energy scenario, in which the total energy-momentum tensor is not conserved. We investigate both the background cosmology and, the large-scale imprint of the UDE by probing the angular power spectrum of galaxy number counts, on ultra-large scales; taking care to include the full relativistic corrections in the observed overdensity. The background evolution shows that only an effective mass smaller than the Planck mass is needed in the early universe in order for predictions in the given theory to match current observational constraints. We found that the effective mass-evolution-rate parameter, which drives the evolution of the UDE, acts to enhance the observed power spectrum and, hence, relativistic effects (on ultra-large scales) by enlarging the UDE sound horizon. Conversely, both the (beyond) Horndeski parameter and the kineticity act to diminish the observed power spectrum, by decreasing the UDE sound horizon. Our results show that, in a universe with UDE, a multi-tracer analysis will be needed to detect the relativistic effects in the large-scale structure. In the light of a multi-tracer analysis, the various relativistic effects hold the potential to distinguish different gravity models. Moreover, while the Doppler effect will remain significant at all epochs and, thus can not be ignored, the integrated Sachs-Wolfe, the time-delay and the potential (difference) effects, respectively, will only become significant at epochs near z=3 and beyond, and may be neglected at late epochs. In the same vein, the Doppler effect alone can serve as an effective cosmological probe for the large-scale structure or gravity models, in the angular power spectrum - at all z.

Apart from the known weak gravitational lensing effect, the cosmic magnification acquires relativistic corrections owing to Doppler, integrated Sachs-Wolfe, time-delay and other (local) gravitational potential effects, respectively. These corrections grow on very large scales and high redshifts z, which will be the reach of forthcoming surveys. In this work, these relativistic corrections are investigated in the magnification angular power spectrum, using both (standard) noninteracting dark energy (DE), and interacting DE (IDE). It is found that for noninteracting DE, the relativistic corrections can boost the magnification large-scale power by ~ 40% at z = 3, and increases at lower z. It is also found that the IDE effect is sensitive to the relativistic corrections in the magnification power spectrum, particularly at low z---which will be crucial for constraints on IDE. Moreover, the results show that if relativistic corrections are not taken into account, this may lead to an incorrect estimate of the large-scale imprint of IDE in the cosmic magnification; including the relativistic corrections can enhance the true potential of the cosmic magnification as a cosmological probe.

Thus far, there does not appear to be an agreed (or adequate) definition of homogeneous dark energy (DE). This paper seeks to define a valid, adequate homogeneity condition for DE. Firstly, it is shown that as long as $w_x \neq -1$, DE must have perturbations. It is then argued, independent of $w_x$, that a correct definition of homogeneous DE is one whose density perturbation vanishes in comoving gauge: and hence, in the DE rest frame. Using phenomenological DE, the consequence of this approach is then investigated in the observed galaxy power spectrum - with the power spectrum being normalized on small scales, at the present epoch z = 0. It is found that for high magnification bias, relativistic corrections in the galaxy power spectrum are able to distinguish the concordance model from both a homogeneous DE and a clustering DE - on super-horizon scales.

The observed galaxy power spectrum acquires relativistic corrections from lightcone effects, and these corrections grow on very large scales. Future galaxy surveys in optical, infrared and radio bands will probe increasingly large wavelength modes and reach higher redshifts. In order to exploit the new data on large scales, an accurate analysis requires inclusion of the relativistic effects. This is especially the case for primordial non-Gaussianity and for extending tests of dark energy models to horizon scales. Here we investigate the latter, focusing on models where the dark energy interacts non-gravitationally with dark matter. Interaction in the dark sector can also lead to large-scale deviations in the power spectrum. If the relativistic effects are ignored, the imprint of interacting dark energy will be incorrectly identified and thus lead to a bias in constraints on interacting dark energy on very large scales.

Quintessence can cluster only on horizon scales. What is the effect on the observed matter distribution? To answer this, we need a relativistic approach that goes beyond the standard Newtonian calculation and deals properly with large scales. Such an approach has recently been developed for the case when dark energy is vacuum energy, which does not cluster at all. We extend this relativistic analysis to deal with dynamical dark energy. Using three quintessence potentials as examples, we compute the angular power spectrum for the case of an HI intensity map survey. Compared to the concordance model with the same small-scale power at z=0, quintessence boosts the angular power by up to ~15% at high redshifts, while power in the two models converges at low redshifts. The difference is mainly due to the background evolution, driven mostly by the normalization of the power spectrum today. The dark energy perturbations make only a small contribution on the largest scales, and a negligible contribution on smaller scales. Ironically, the dark energy perturbations remove the false boost of large-scale power that arises if we impose the (unphysical) assumption that the dark energy is smooth.