We study cosmological observables on the past light cone of a fixed observer in the context of clustering dark energy. We focus on observables that probe the gravitational field directly, namely the integrated Sachs-Wolfe and non-linear Rees-Sciama effect (ISW-RS), weak gravitational lensing, gravitational redshift and Shapiro time delay. With our purpose-built $N$-body code "$k$-evolution" that tracks the coupled evolution of dark matter particles and the dark energy field, we are able to study the regime of low speed of sound $c_s$ where dark energy perturbations can become quite large. Using ray tracing we produce two-dimensional sky maps for each effect and we compute their angular power spectra. It turns out that the ISW-RS signal is the most promising probe to constrain clustering dark energy properties coded in $w-c_s^2$, as the $\textit{linear}$ clustering of dark energy would change the angular power spectrum by $\sim 30\%$ at low $\ell$ when comparing two different speeds of sound for dark energy. Weak gravitational lensing, Shapiro time-delay and gravitational redshift are less sensitive probes of clustering dark energy, showing variations of a few percent only. The effect of dark energy $\textit{non-linearities}$ in all the power spectra is negligible at low $\ell$, but reaches about $2\%$ and $3\%$, respectively, in the convergence and ISW-RS angular power spectra at multipoles of a few hundred when observed at redshift $\sim 0.85$. Future cosmological surveys achieving percent precision measurements will allow to probe the clustering of dark energy to a high degree of confidence. Clear evidence for clustering dark energy at any scale would rule out a cosmological constant as the leading contender for dark energy, and may help to distinguish between different dark energy models.

Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $\Lambda$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $\Lambda$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

We introduce $k$-evolution, a relativistic $N$-body code based on $\textit{gevolution}$, which includes clustering dark energy among its cosmological components. To describe dark energy, we use the effective field theory approach. In particular, we focus on $k$-essence with a speed of sound much smaller than unity but we lay down the basis to extend the code to other dark energy and modified gravity models. We develop the formalism including dark energy non-linearities but, as a first step, we implement the equations in the code after dropping non-linear self-coupling in the $k$-essence field. In this simplified setup, we compare $k$-evolution simulations with those of $\texttt{CLASS}$ and $\textit{gevolution}$ 1.2, showing the effect of dark matter and gravitational non-linearities on the power spectrum of dark matter, of dark energy and of the gravitational potential. Moreover, we compare $k$-evolution to Newtonian $N$-body simulations with back-scaled initial conditions and study how dark energy clustering affects massive halos.

In this paper, we quantify the non-linear effects from $k$-essence dark energy through an effective parameter $\mu$ that encodes the additional contribution of a dark energy fluid or a modification of gravity to the Poisson equation. This is a first step toward quantifying non-linear effects of dark energy/modified gravity models in a more general approach. We compare our $N$-body simulation results from $k$-evolution with predictions from the linear Boltzmann code $\texttt{CLASS}$, and we show that for the $k$-essence model one can safely neglect the difference between the two potentials, $ \Phi -\Psi$, and short wave corrections appearing as higher order terms in the Poisson equation, which allows us to use single parameter $\mu$ for characterizing this model. We also show that for a large $k$-essence speed of sound the $\texttt{CLASS}$ results are sufficiently accurate, while for a low speed of sound non-linearities in matter and in the $k$-essence field are non-negligible. We propose a $\tanh$-based parameterisation for $\mu$, motivated by the results for two cases with low ($c_s^2=10^{-7}$) and high ($c_s^2=10^{-4}$) speed of sound, to include the non-linear effects based on the simulation results. This parametric form of $\mu$ can be used to improve Fisher forecasts or Newtonian $N$-body simulations for $k$-essence models.

Cosmological N-body simulations are done on massively parallel computers. This necessitates the use of simple time integrators, and, additionally, of mesh-grid approximations of the potentials. Recently, Adamek et al. (2015); Barrera-Hinojosa et al. (2019) have developed general relativistic N-body simulations to capture relativistic effects mainly for cosmological purposes. We therefore ask whether, with the available technology, relativistic effects like perihelion advance can be detected numerically to a relevant precision. We first study the spurious perihelion shift in the Kepler problem, as a function of the integration method used, and then as a function of an additional interpolation of forces on a 2-dimensional lattice. This is done for several choices of eccentricities and semi-major axes. Using these results, we can predict which precisions and lattice constants allow for a detection of the relativistic perihelion advance in N-body simulation. We find that there are only small windows of parameters---such as eccentricity, distance from the central object and the Schwarzschild radius---for which the corrections can be detected in the numerics.

We present $\texttt{MG-evolution}$, an $N$-body code simulating the cosmological structure formation for parametrised modifications of gravity. It is built from the combination of parametrised linear theory with a parametrisation of the deeply nonlinear cosmological regime extrapolated from modified spherical collapse computations that cover the range of known screening mechanisms. We test $\texttt {MG-evolution}$, which runs at the speed of conventional $\Lambda$CDM simulations, against a suit of existing exact model-specific codes, encompassing linearised and chameleon $f(R)$ gravity as well as the normal branch of the Dvali-Gabadadz-Porrati braneworld model, hence covering both large-field value and large-derivative screening effects. We compare the nonlinear power spectra produced by the parametrised and model-specific approaches over the full range of scales set by the box size and resolution of our simulations, $k=(0.05-2.5)$~h/Mpc, and for two redshift slices, $z=0$ and $z=1$. We find sub-percent to one-percent level recovery of all the power spectra generated with the model-specific codes for the full range of scales. $\texttt {MG-evolution}$ can be used for generalised and accurate tests of gravity and dark energy with the increasing wealth of high-precision cosmological survey data becoming available over the next decade.

The Cosmic Microwave Background (CMB) radiation lensing is a promising tool to study the physics of early universe. In this work we probe the imprints of deviations from isotropy and scale invariance of primordial curvature perturbation power spectrum on CMB lensing potential and convergence. Specifically, we consider a scale-dependent hemispherical asymmetry in primordial power spectrum. We show that the CMB lensing potential and convergence and also the cross-correlation of the CMB lensing and late time galaxy convergence can probe the amplitude and the scale dependence of the dipole modulation. As another example, we consider a primordial power spectrum with local feature. We show that the CMB lensing and the cross-correlation of the CMB lensing and galaxy lensing can probe the amplitude and the shape of the local feature. We show that the cross correlation of CMB lensing convergence and galaxy lensing is capable to probe the effects of local features in power spectrum on smaller scales than the CMB lensing. Finally we showed that the current data can constrain the amplitude and moment dependence of dipole asymmetry.