We study the importance of gravitational lensing in the modelling of the number counts of galaxies. We confirm previous results for photometric surveys, showing that lensing cannot be neglected in a survey like LSST since it would infer a significant shift of cosmological parameters. For a spectroscopic survey like SKA2, we find that neglecting lensing in the monopole, quadrupole and hexadecapole of the correlation function also induces an important shift of parameters. For $\Lambda$CDM parameters, the shift is moderate, of the order of $0.6\sigma$ or less. However, for a model-independent analysis, that measures the growth rate of structure in each redshift bin, neglecting lensing introduces a shift of up to $2.3\sigma$ at high redshift. Since the growth rate is directly used to test the theory of gravity, such a strong shift would wrongly be interpreted as the breakdown of General Relativity. This shows the importance of including lensing in the analysis of future surveys. On the other hand, for a survey like DESI, we find that lensing is not important, mainly due to the value of the magnification bias parameter of DESI, $s(z)$, which strongly reduces the lensing contribution at high redshift. We also propose a way of improving the analysis of spectroscopic surveys, by including the cross-correlations between different redshift bins (which is neglected in spectroscopic surveys) from the spectroscopic survey or from a different photometric sample. We show that including the cross-correlations in the SKA2 analysis does not improve the constraints. On the other hand replacing the cross-correlations from SKA2 by cross-correlations measured with LSST improves the constraints by 10 to 20 %. Interestingly, for $\Lambda$CDM parameters, we find that LSST and SKA2 are highly complementary, since they are affected differently by degeneracies between parameters.

We study the impact of relativistic effects in the 3-dimensional cross-correlation between Lyman-$\alpha$ forest and quasars. Apart from the relativistic effects, which are dominated by the Doppler contribution, several systematic effects are also included in our analysis (intervening metals, unidentified high column density systems, transverse proximity effect and effect of the UV fluctuations). We compute the signal-to-noise ratio for the Baryonic Oscillation Spectroscopic Survey (BOSS), the extended Baryonic Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI) surveys, showing that DESI will be able to detect the Doppler contribution in a Large Scale Structure (LSS) survey for the first time, with a S/N $>7$ for $r_{\rm min} > 10$ Mpc$/h$, where r$_{\rm min}$ denotes the minimum comoving separation between sources. We demonstrate that several physical parameters, introduced to provide a full modelling of the cross-correlation function, are affected by the Doppler contribution. By using a Fisher matrix approach, we establish that if the Doppler contribution is neglected in the data analysis, the derived parameters will be shifted by a non-negligible amount for the upcoming surveys.

We investigate for the first time weak lensing in a high-resolution relativistic N-body simulation, which includes all relevant effects from general relativity. We integrate the photon geodesics backwards from the observer to the emitters. Our analysis is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. We solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. We present the angular power spectra and one-point probability distribution functions for the weak-lensing variables, which we find are broadly in agreement with comparable Newtonian simulations. Our geometric approach, however, is more robust and flexible, and can therefore be applied consistently to non-standard cosmologies and modified theories of gravity.

Clustering of a perfect fluid does not lead to the generation of vorticity. It is the collisionless nature of dark matter, inducing velocity dispersion and shell crossing, which is at the origin of cosmological vorticity generation. In this paper we investigate the generation of vorticity during the formation of cosmological large scale structure using the public relativistic N-body code $\textit{gevolution}$. We test several methods to compute the vorticity power spectrum and we study its convergence with respect to the mass and grid resolution of our simulations. We determine the power spectrum, the spectral index on large-scales, the amplitude of the peak position and their time evolution. We also compare the vorticity extracted from our simulations with the vector perturbations of the metric. Our results are accompanied by resolution studies and compared with previous studies in the literature.

We investigate the future perspectives of the detection of the relativistic dipole by cross-correlating the 21 cm emission in Intensity Mapping (IM) and galaxy surveys at low redshift. We model the neutral hydrogen (HI) and the galaxy population by means of the halo model to relate the parameters that affect the dipole signal such as the biases of the two tracers and the Poissonian noise. We investigate the behavior of the signal-to-noise as a function of the galaxy and magnification biases, for two fixed models of the neutral hydrogen. In both cases we found that the signal-to-noise does not grow by increasing the difference between the biases of the two tracers, due to the larger shot-noise yields by highly biased tracers. We also study and provide an optimal luminosity-threshold galaxy catalogue to enhance the signal-to-noise ratio of the relativistic dipole. Interestingly, we show that the maximum magnitude provided by the survey does not lead to the maximum signal-to-noise for detecting relativistic effects and we predict the optimal value for the limiting magnitude. Our work suggests that an optimal analysis could increase the signal-to-noise ratio up to a factor five compared to a standard one.

We study the impact of neglecting lensing magnification in galaxy clustering analyses for future galaxy surveys, considering the LambdaCDM model and two extensions: massive neutrinos and modifications of General Relativity. Our study focuses on the biases on the constraints and on the estimation of the cosmological parameters. We perform a comprehensive investigation of these two effects for the upcoming photometric and spectroscopic galaxy surveys Euclid and SKA for different redshift binning configurations. We also provide a fitting formula for the magnification bias of SKA. Our results show that the information present in the lensing contribution does improve the constraints on the modified gravity parameters whereas the lensing constraining power is negligible for the LambdaCDM parameters. For photometric surveys the estimation is biased for all the parameters if lensing is not taken into account. This effect is particularly significant for the modified gravity parameters. Conversely for spectroscopic surveys the bias is below one sigma for all the parameters. Our findings show the importance of including lensing in galaxy clustering analyses for testing General Relativity and to constrain the parameters which describe its modifications.

We investigate the Alcock Paczynski (AP) test applied to the Baryon Acoustic Oscillation (BAO) feature in the galaxy correlation function. By using a general formalism that includes relativistic effects, we quantify the importance of the linear redshift space distortions and gravitational lensing corrections to the galaxy number density fluctuation. We show that redshift space distortions significantly affect the shape of the correlation function, both in radial and transverse directions, causing different values of galaxy bias to induce offsets up to 1% in the AP test. On the other hand, we find that the lensing correction around the BAO scale modifies the amplitude but not the shape of the correlation function and therefore does not introduce any systematic effect. Furthermore, we investigate in details how the AP test is sensitive to redshift binning: a window function in transverse direction suppresses correlations and shifts the peak position toward smaller angular scales. We determine the correction that should be applied in order to account for this effect, when performing the test with data from three future planned galaxy redshift surveys: Euclid, the Dark Energy Spectroscopic Instrument (DESI) and the Square Kilometer Array (SKA).