We show that correlations between the phases of the galaxy density field in redshift space provide additional information about the growth rate of large-scale structure that is complementary to the power spectrum multipoles. In particular, we consider the multipoles of the line correlation function (LCF), which correlates phases between three collinear points, and use the Fisher forecasting method to show that the LCF multipoles can break the degeneracy between the measurement of the growth rate of structure $f$ and the amplitude of perturbations $\sigma_8$ that is present in the power spectrum multipoles at large scales. This leads to an improvement in the measurement of $f$ and $\sigma_8$ by up to 220 per cent for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ and up to 50 per cent for $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$ at redshift $z=0.25$, with respect to power spectrum measurements alone for the upcoming generation of galaxy surveys like DESI and Euclid. The average improvements in the constraints on $f$ and $\sigma_8$ for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ are $\sim 90$ per cent for the DESI BGS sample with mean redshift $\overline{z}=0.25$, $\sim 40$ per cent for the DESI ELG sample with $\overline{z}=1.25$, and $\sim 40$ per cent for the Euclid H$\alpha$ galaxies with $\overline{z}=1.3$. For $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$, the average improvements are $\sim 40$ per cent for the DESI BGS sample and $\sim 20$ per cent for both the DESI ELG and Euclid H$\alpha$ galaxies.

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.

The weak equivalence principle is one of the cornerstone of general relativity. Its validity has been tested with impressive precision in the Solar System, with experiments involving baryonic matter and light. However, on cosmological scales and when dark matter is concerned, the validity of this principle is still unknown. In this paper we construct a null test that probes the validity of the equivalence principle for dark matter. Our test has the strong advantage that it can be applied on data without relying on any modelling of the theory of gravity. It involves a combination of redshift-space distortions and relativistic effects in the galaxy number-count fluctuation, that vanishes if and only if the equivalence principle holds. We show that the null test is very insensitive to typical uncertainties in other cosmological parameters, including the magnification bias parameter, and to non-linear effects, making this a robust null test for modified gravity.

We investigate the ability of the Laser Interferometer Space Antenna (LISA) to measure the center of mass acceleration of stellar-origin black hole binaries emitting gravitational waves. Our analysis is based on the idea that the acceleration of the center of mass induces a time variation in the redshift of the gravitational wave, which in turn modifies its waveform. We confirm that while the cosmological acceleration is too small to leave a detectable imprint on the gravitational waveforms observable by LISA, larger peculiar accelerations may be measurable for sufficiently long lived sources. We focus on stellar mass black hole binaries, which will be detectable at low frequencies by LISA and near coalescence by ground based detectors. These sources may have large peculiar accelerations, for instance, if they form in nuclear star clusters or in AGN accretion disks. If that is the case, we find that in an astrophysical population calibrated to the LIGO-Virgo observed merger rate, LISA will be able to measure the peculiar acceleration of a small but significant fraction of the events if the mission lifetime is extended beyond the nominal duration of 4 years. In this scenario LISA will be able to assess whether black hole binaries form close to galactic centers, particularly in AGN disks, and will thus help discriminate between different formation mechanisms. Although for a nominal 4 years LISA mission the peculiar acceleration effect cannot be measured, a consistent fraction of events may be biased by strong peculiar accelerations which, if present, may imprint large systematic errors on some waveform parameters. In particular, estimates of the luminosity distance could be strongly biased and consequently induce large systematic errors on LISA measurements of the Hubble constant with stellar mass black hole binaries.

The main science driver for the coming generation of cosmological surveys is understanding dark energy which relies on testing General Relativity on the largest scales. Once we move beyond the simplest explanation for dark energy of a cosmological constant, the space of possible theories becomes both vast and extremely hard to compute realistic observables. A key discriminator of a cosmological constant, however, is that the growth of structure is scale-invariant on large scales. By carefully weighting observables derived from distributions of galaxies and a dipole pattern in their apparent sizes, we construct a null test which vanishes for any model of gravity or dark energy where the growth of structure is scale-independent. It relies only on very few assumptions about cosmology, and does not require any modelling of the growth of structure. We show that with a survey like DESI a scale-dependence of the order of 10-20 percent can be detected at 3 sigma with the null test, which will drop by a factor of 2 for a survey like the Square Kilometre Array. We also show that the null test is very insensitive to typical uncertainties in other cosmological parameters including massive neutrinos and scale-dependent bias, making this a key null test for dark energy.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We propose a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$\times$galaxy $-$ galaxy$\times$IM; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As a pure cross-correlation estimator, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. For cosmic variance dominated surveys the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of DES cross-correlated with HIRAX, the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

The apparent sizes and brightnesses of galaxies are correlated in a dipolar pattern around matter overdensities in redshift space, appearing larger on their near side and smaller on their far side. The opposite effect occurs for galaxies around an underdense region. These patterns of apparent magnification induce dipole and higher multipole terms in the cross-correlation of galaxy number density fluctuations with galaxy size/brightness (which is sensitive to the convergence field). This provides a means of directly measuring peculiar velocity statistics at low and intermediate redshift, with several advantages for performing cosmological tests of GR. In particular, it does not depend on empirically-calibrated scaling relations like the Tully-Fisher and Fundamental Plane methods. We show that the next generation of spectroscopic galaxy redshift surveys will be able to measure the Doppler magnification effect with sufficient signal-to-noise to test GR on large scales. We illustrate this with forecasts for the constraints that can be achieved on parametrised deviations from GR for forthcoming low-redshift galaxy surveys with DESI and SKA2. Although the cross-correlation statistic considered has a lower signal to noise than RSD, it will be a useful probe of GR since it is sensitive to different systematics.

Redshift-space distortions are a sensitive probe of the growth of large-scale structure. In the linear regime, redshift-space distortions are fully described by the multipoles of the two-point correlation function. In the nonlinear regime, however, higher-order statistics are needed to capture the full information of the galaxy density field. In this paper, we show that the redshift-space line correlation function--which is a measure of Fourier phase correlations--is sensitive to the nonlinear growth of the density and velocity fields and to the nonlinear mapping between real and redshift space. We expand the line correlation function in multipoles, and we show that almost all of the information is encoded in the monopole, quadrupole, and hexadecapole. We argue that these multipoles are highly complementary to the multipoles of the two-point correlation function: first, because they are directly sensitive to the difference between the density and the velocity coupling kernels, which is a purely nonlinear quantity; and second, because the multipoles are proportional to different combinations of $f$ and $\sigma_8$. Measured in conjunction with the two-point correlation function and the bispectrum, the multipoles of the line correlation function could therefore allow us to disentangle efficiently these two quantities and to test modified theories of gravity.

We present a full-sky derivation of weak lensing observables in the Post-Friedmann (PF) formalism. Weak lensing has the characteristic of mixing small scales and large scales since it is affected by inhomogeneities integrated along the photon trajectory. With the PF formalism, we develop a modelling of lensing observables which encompasses both leading order relativistic effects and effects that are due to the fully non-linear matter distribution at small scales. We derive the reduced shear, convergence and rotation up to order $1/c^4$ in the PF approximation, accounting for scalar, vector and tensor perturbations, as well as galaxies' peculiar velocities. We discuss the various contributions that break the Kaiser-Squires relation between the shear and the convergence at different orders. We pay particular attention to the impact of the frame-dragging vector potential on lensing observables and we discuss potential ways to measure this effect in future lensing surveys.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We study the impact on the galaxy correlation function of the presence of a vector component in the tracers' peculiar velocities, in the case in which statistical isotropy is violated. We present a general framework - based on the bipolar spherical harmonics expansion - to study this effect in a model independent way, without any hypothesis on the origin or the properties of these vector modes. We construct six new observables, that can be directly measured in galaxy catalogs in addition to the standard monopole, quadrupole and hexadecapole, and we show that they completely describe any deviations from isotropy. We then perform a Fisher analysis in order to quantify the constraining power of future galaxy surveys. As an example, we show that the SKA2 would be able to detect anisotropic rotational velocities with amplitudes as low as 1% of that of the vorticity generated during shell-crossing in standard dark matter scenarios.

We present a public version of the code COFFE (COrrelation Function Full-sky Estimator) available at https://github.com/JCGoran/coffe. The code computes the galaxy two-point correlation function and its multipoles in linear perturbation theory, including all relativistic and wide angle corrections. COFFE also calculates the covariance matrix for two physically relevant estimators of the correlation function multipoles. We illustrate the usefulness of our code by a simple but relevant example: a forecast of the detectability of the lensing signal in the multipoles of the two-point function. In particular, we show that lensing should be detectable in the multipoles of the two-point function, with a signal-to-noise larger than 10, in future surveys like Euclid or the SKA.

Most statistical inference from cosmic large-scale structure relies on two-point statistics, i.e.\ on the galaxy-galaxy correlation function (2PCF) or the power spectrum. These statistics capture the full information encoded in the Fourier amplitudes of the galaxy density field but do not describe the Fourier phases of the field. Here, we quantify the information contained in the line correlation function (LCF), a three-point Fourier phase correlation function. Using cosmological simulations, we estimate the Fisher information (at redshift $z=0$) of the 2PCF, LCF and their combination, regarding the cosmological parameters of the standard $\Lambda$CDM model, as well as a Warm Dark Matter (WDM) model and the $f(R)$ and Symmetron modified gravity models. The galaxy bias is accounted for at the level of a linear bias. The relative information of the 2PCF and the LCF depends on the survey volume, sampling density (shot noise) and the bias uncertainty. For a volume of $1h^{-3}\rm Gpc^3$, sampled with points of mean density $\bar{n} = 2\times10^{-3} h^{3}\ \rm Mpc^{-3}$ and a bias uncertainty of 13\%, the LCF improves the parameter constraints by about 20\% in the $\Lambda$CDM cosmology and potentially even more in alternative models. Finally, since a linear bias only affects the Fourier amplitudes (2PCF), but not the phases (LCF), the combination of the 2PCF and the LCF can be used to break the degeneracy between the linear bias and $\sigma_8$, present in 2-point statistics.

The equivalence principle, that is one of the main pillars of general relativity, is very well tested in the Solar system; however, its validity is more uncertain on cosmological scales, or when dark matter is concerned. This article shows that relativistic effects in the large-scale structure can be used to directly test whether dark matter satisfies Euler's equation, i.e. whether its free fall is characterised by geodesic motion, just like baryons and light. After having proposed a general parametrisation for deviations from Euler's equation, we perform Fisher-matrix forecasts for future surveys like DESI and the SKA, and show that such deviations can be constrained with a precision of order 10%. Deviations from Euler's equation cannot be tested directly with standard methods like redshift-space distortions and gravitational lensing, since these observables are not sensitive to the time component of the metric. Our analysis shows therefore that relativistic effects bring new and complementary constraints to alternative theories of gravity.

We derive an exact expression for the correlation function in redshift shells including all the relativistic contributions. This expression, which does not rely on the distant-observer or flat-sky approximation, is valid at all scales and includes both local relativistic corrections and integrated contributions, like gravitational lensing. We present two methods to calculate this correlation function, one which makes use of the angular power spectrum C_ell(z1,z2) and a second method which evades the costly calculations of the angular power spectra. The correlation function is then used to define the power spectrum as its Fourier transform. In this work theoretical aspects of this procedure are presented, together with quantitative examples. In particular, we show that gravitational lensing modifies the multipoles of the correlation function and of the power spectrum by a few percent at redshift z=1 and by up to 30% and more at z=2. We also point out that large-scale relativistic effects and wide-angle corrections generate contributions of the same order of magnitude and have consequently to be treated in conjunction. These corrections are particularly important at small redshift, z=0.1, where they can reach 10%. This means in particular that a flat-sky treatment of relativistic effects, using for example the power spectrum, is not consistent.

We investigate the feasibility of measuring the effects of peculiar velocities in large-scale structure using the dipole of the redshift-space cross-correlation function. We combine number counts of galaxies with brightness-temperature fluctuations from 21cm intensity mapping, demonstrating that the dipole may be measured at modest significance ($\lesssim 2\sigma$) by combining the upcoming radio survey CHIME with the future redshift surveys of DESI and Euclid. More significant measurements ($\lesssim~10\sigma$) will be possible by combining intensity maps from the SKA with these of DESI or Euclid, and an even higher significance measurement ($\lesssim 100\sigma$) may be made by combining observables completely internally to the SKA. We account for effects such as contamination by wide-angle terms, interferometer noise and beams in the intensity maps, non-linear enhancements to the power spectrum, stacking multiple populations, sensitivity to the magnification slope, and the possibility that number counts and intensity maps probe the same tracers. We also derive a new expression for the covariance matrix of multi-tracer redshift-space correlation function estimators with arbitrary orientation weights, which may be useful for upcoming surveys aiming at measuring redshift-space clustering with multiple tracers.

We compute a general expression for the contribution of vector perturbations to the redshift-space distortion of galaxy surveys. We show that they contribute to the same multipoles of the correlation function as scalar perturbations and should thus in principle be taken into account in data analysis. We derive constraints for next-generation surveys on the amplitude of two sources of vector perturbations, namely non-linear clustering and topological defects. While topological defects leave a very small imprint on redshift-space distortions, we show that the multipoles of the correlation function are sensitive to vorticity induced by non-linear clustering. Therefore future redshift surveys such as DESI or the SKA should be capable of measuring such vector modes, especially with the hexadecapole which appears to be the most sensitive to the presence of vorticity.

Objects falling into an overdensity appear larger on its near side and smaller on its far side than other objects at the same redshift. This produces a dipolar pattern of magnification, primarily as a consequence of the Doppler effect. At low redshift this Doppler magnification completely dominates the usual integrated gravitational lensing contribution to the lensing magnification. We show that one can optimally observe this pattern by extracting the dipole in the cross-correlation of number counts and galaxy sizes. This dipole allows us to almost completely remove the contribution from gravitational lensing up to redshift 0.5, and even at high redshift z~1 the dipole picks up the Doppler magnification predominantly. Doppler magnification should be easily detectable in current and upcoming optical and radio surveys; by forecasting for telescopes such as the SKA, we show that this technique is competitive with using peculiar velocities via redshift-space distortions to constrain dark energy. It produces similar yet complementary constraints on the cosmological model to those found using measurements of the cosmic shear.

Third generation ground-based interferometers as well as the planned space-based interferometer LISA are expected to detect a plethora of gravitational wave signals from coalescing binaries at cosmological distance. The emitted gravitational waves propagate in the expanding universe through the inhomogeneous distribution of matter. Here we show that the acceleration of the universe and the peculiar acceleration of the binary with respect to the observer distort the gravitational chirp signal from the simplest General Relativity prediction beyond a mere time independent rescaling of the chirp mass, affecting intrinsic parameter estimations for the binaries visible by LISA. We find that the effect due to the peculiar acceleration can be much larger than the one due to the universe acceleration. Moreover, peculiar accelerations can introduce a bias in the estimation of parameters such as the time of coalescence and the individual masses of the binary. An error in the estimation of the arrival time will have an impact in the case of sources visible first by LISA and later by ground based interferometers.

It has been shown recently that relativistic distortions generate a dipolar modulation in the two-point correlation function of galaxies. To measure this relativistic dipole it is necessary to cross-correlate different populations of galaxies with for example different luminosities or colours. In this paper, we construct an optimal estimator to measure the dipole with multiple populations. We show that this estimator increases the signal-to-noise of the dipole by up to 35 percent. Using 6 populations of galaxies, in a survey with halos and number densities similar to those of the millennium simulation, we forecast a cumulative signal-to-noise of 4.4. For the main galaxy sample of SDSS at low redshift z<0.2 our optimal estimator predicts a cumulative signal-to-noise of 2.4. Finally we forecast a cumulative signal-to-noise of 7.4 in the upcoming DESI survey. These forecasts indicate that with the appropriate choice of estimator the relativistic dipole should be detectable in current and future surveys.

It is usually assumed that in the linear regime the two-point correlation function of galaxies contains only a monopole, quadrupole and hexadecapole. Looking at cross-correlations between different populations of galaxies, this turns out not to be the case. In particular, the cross-correlations between a bright and a faint population of galaxies contain also a dipole. In this paper we present the first attempt to measure this dipole. We discuss the four types of effects that contribute to the dipole: relativistic distortions, evolution effect, wide-angle effect and large-angle effect. We show that the first three contributions are intrinsic anti-symmetric contributions that do not depend on the choice of angle used to measure the dipole. On the other hand the large-angle effect appears only if the angle chosen to extract the dipole breaks the symmetry of the problem. We show that the relativistic distortions, the evolution effect and the wide-angle effect are too small to be detected in the LOWz and CMASS sample of the BOSS survey. On the other hand with a specific combination of angles we are able to measure the large-angle effect with high significance. We emphasise that this large-angle dipole does not contain new physical information, since it is just a geometrical combination of the monopole and the quadrupole. However this measurement, which is in excellent agreement with theoretical predictions, validates our method for extracting the dipole from the two-point correlation function and it opens the way to the detection of relativistic effects in future surveys like e.g. DESI.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We show that at second order, ensemble averages of observables and directional averages do not commute due to gravitational lensing -- observing the same thing in many directions over the sky is not the same as taking an ensemble average. In principle this non-commutativity is significant for a variety of quantities that we often use as observables and can lead to a bias in parameter estimation. We derive the relation between the ensemble average and the directional average of an observable, at second order in perturbation theory. We discuss the relevance of these two types of averages for making predictions of cosmological observables, focusing on observables related to distances and magnitudes. In particular, we show that the ensemble average of the distance in a given observed direction is increased by gravitational lensing, whereas the directional average of the distance is decreased. For a generic observable, there exists a particular function of the observable that is not affected by second-order lensing perturbations. We also show that standard areas have an advantage over standard rulers, and we discuss the subtleties involved in averaging in the case of supernova observations.

It has recently been shown that second-order corrections to the background distance-redshift relation can build up significantly at large redshifts, due to an aggregation of gravitational lensing events. This shifts the expectation value of the distance to the CMB by 1%. In this paper we show that this shift is already properly accounted for in standard CMB analyses. We clarify the role that the area distance to the CMB plays in the presence of second-order lensing corrections.

We derive an analytical expression for a novel large-scale structure observable: the line correlation function. The line correlation function, which is constructed from the three-point correlation function of the phase of the density field, is a robust statistical measure allowing the extraction of information in the non-linear and non-Gaussian regime. We show that, in perturbation theory, the line correlation is sensitive to the coupling kernel F_2, which governs the non-linear gravitational evolution of the density field. We compare our analytical expression with results from numerical simulations and find a 1-sigma agreement for separations r<30 Mpc/h. Fitting formulae for the power spectrum and the non-linear coupling kernel at small scales allow us to extend our prediction into the strongly non-linear regime where we find a 1-sigma agreement with the simulations for r<2 Mpc/h. We discuss the advantages of the line correlation relative to standard statistical measures like the bispectrum. Unlike the latter, the line correlation is independent of the bias, in the regime where the bias is local and linear. Furthermore, the variance of the line correlation is independent of the Gaussian variance on the modulus of the density field. This suggests that the line correlation can probe more precisely the non-linear regime of gravity, with less contamination from the power spectrum variance.

We present a fully relativistic calculation of the observed galaxy number counts in the linear regime. We show that besides the density fluctuations and redshift-space distortions, various relativistic effects contribute to observations at large scales. These effects all have the same physical origin: they result from the fact that our coordinate system, namely the galaxy redshift and the incoming photons' direction, is distorted by inhomogeneities in our universe. We then discuss the impact of the relativistic effects on the angular power spectrum and on the two-point correlation function in configuration space. We show that the latter is very well adapted to isolate the relativistic effects since it naturally makes use of the symmetries of the different contributions. In particular, we discuss how the Doppler effect and the gravitational redshift distortions can be isolated by looking for a dipole in the cross-correlation function between a bright and a faint population of galaxies.

We study the two-point cross-correlation function between two populations of galaxies: for instance a bright population and a faint population. We show that this cross-correlation is asymmetric under the exchange of the line-of-sight coordinate of the galaxies, i.e. that the correlation is different if the bright galaxy is in front of, or behind, the faint galaxy. We give an intuitive, quasi-Newtonian derivation of all the effects that contribute to such an asymmetry in large-scale structure: gravitational redshift, Doppler shift, lensing, light-cone, evolution and Alcock-Paczynski effects - interestingly, the gravitational redshift term is exactly canceled by some of the others, assuming geodesic motion. Most of these effects are captured by previous calculations of general relativistic corrections to the observed galaxy density fluctuation; the asymmetry arises from terms that are suppressed by the ratio H/k - H is the Hubble constant and k is the wavenumber - which are more readily observable than the terms suppressed by (H/k)^2. Some of the contributions to the asymmetry, however, arise from terms that are generally considered 'Newtonian' - the lensing and evolution - and thus represent a contaminant in the search for general relativistic corrections. We propose methods to disentangle these different contributions. A simple method reduces the contamination to a level of < 10% for redshifts z<1. We also clarify the relation to recent work on measuring gravitational redshifts by stacking clusters.

If the B-mode signal in the CMB polarization seen by the BICEP2 experiment is confirmed, it has dramatic implications for models of inflation. The result is also in tension with Planck limits on standard inflationary models. It is therefore important to investigate whether this signal can arise from alternative sources. If so, this could lessen the pressure on inflationary models and the tension with Planck data. We investigate whether vector and tensor modes from primordial magnetic fields can explain the signal. We find that in principle, magnetic fields generated during inflation can indeed produce the required B-mode, for a suitable range of energy scales of inflation. In this case, the primordial gravitational wave amplitude is negligible, so that there is no tension with Planck and no problems posed for current inflationary models. However, the simplest magnetic model is in tension with Planck limits on non-Gaussianity in the trispectrum. It may be possible to fine-tune the magnetogenesis model so that this non-Gaussianity is suppressed. Alternatively, a weaker magnetic field can pass the non-Gaussianity constraints and allow the primordial tensor mode to be reduced to r~0.09, thus removing the tension with Planck data and alleviating the problems with simple inflationary models.

We study the effect of peculiar motion in weak gravitational lensing. We derive a fully relativistic formula for the cosmic shear and the convergence in a perturbed Friedmann Universe. We find a new contribution related to galaxies peculiar velocity. This contribution does not affect cosmic shear in a measurable way, since it is of second order in the velocity. However, its effect on the convergence (and consequently on the magnification, which is a measurable quantity) is important, especially for redshifts z < 1. As a consequence, peculiar motion modifies also the relation between the shear and the convergence.

The cosmic microwave background (CMB) temperature anisotropies from primordial magnetic fields are studied. In addition to the known passive and the compensated mode we discuss an inflationary magnetic mode in the curvature perturbation, present when magnetic fields are generated during inflation. This mode is absent if the generation mechanism is causal, e.g. a phase transition. We compute and discuss the effect of this mode on the observed CMB anisotropy spectrum, in comparison with the passive and compensated ones. We find that it dominates the CMB anisotropy, and consequently leads to stronger constraints on the amplitude $B_{1\,{\rm Mpc}}$ and spectral index $n_B$ of the magnetic field than what usually found in CMB analyses from the compensated mode. This happens in particular for spectral indexes $n_B>-3$: the inflationary magnetic mode is always scale invariant, therefore through this mode even a magnetic field with a spectrum which is not scale invariant can leave a detectable signal in the CMB at large scales.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We investigate the prospects for constraining alternative theories of gravity with a typical near-term low-budget 21 cm intensity mapping experiment. We derive the 21 cm brightness temperature perturbation consistently in linear theory including all line-of-sight and relativistic effects. We uncover new terms that are a small correction on large scales, analogous to those recently found in the context of galaxy surveys. We then perform a Fisher matrix analysis of the B_0 parametrization of f(R) gravity, where B_0 is proportional to the square of Compton wavelength of the scalaron. We find that our 21 cm survey, in combination with CMB information from Planck, will be able to place a 95% upper limit of 7 x 10^{-5} on B_0 in flat models with a LCDM expansion history, improving on current cosmological constraints by several orders of magnitude. We argue that this constraint is limited by our ability to model the mildly non-linear regime of structure formation in modified gravity. We also perform a model-independent principal component analysis on the free functions introduced into the field equations by modified gravity, mu and Sigma. We find that 20--30 modes of the free functions will be `well-constrained' by our combination of observables, the lower and upper limits dependent on the criteria used to define the `goodness' of the constraint. These constraints are found to be robust to uncertainties in the time-dependence of the bias. Our analysis reveals that our observables are sensitive primarily to temporal variations in Sigma and scale variations in mu. We argue that the inclusion of 21 cm intensity maps will significantly improve constraints on any cosmological deviations from General Relativity in large-scale structure in a very cost-effective manner.

We discover that the mass of dark matter particles mDM is imprinted in phase-correlations of the cosmic density field more significantly than in the 2-point correlation. In particular, phase-correlations trace mDM out to scales about five times larger than the 2-point correlation. This result relies on a new estimator l(r) of pure phase-information in Fourier space, which can be interpreted as a parameter-free and scale-invariant tracer of filament-like structure. Based on simulated density fields we show how mDM can, in principle, be measured using l(r), given a suitably reconstructed density field.

We compute the contribution to the scalar metric perturbations from large-scale magnetic fields which are generated during inflation. We show that apart from the usual passive and compensated modes, the magnetic fields also contribute to the constant mode from inflation. This is different from the causal (post inflationary) generation of magnetic fields where such a mode is absent and it might lead to significant, non-Gaussian CMB anisotropies.

In this paper we compute the quantity which is truly measured in a large galaxy survey. We take into account the effects coming from the fact that we actually observe galaxy redshifts and sky positions and not true spatial positions. Our calculations are done within linear perturbation theory for both the metric and the observer velocities but they can be used for non-linear matter power spectra. We shall see that the complications due to the fact that we only observe on our background lightcone and that we do not truly know the distance of the observed galaxy, but only its redshift is not only an additional difficulty, but even more a new opportunity for future galaxy surveys.

This paper has been withdrawn by the authors due to a mistake in the calculations leading to an erroneous conclusion.

Future lensing surveys will be nearly full-sky and reach an unprecedented depth, probing scales closer and closer to the Hubble radius. This motivates the study of the cosmic shear beyond the small-angle approximation and including general relativistic corrections that are usually suppressed on sub-Hubble scales. The complete expression of the reduced cosmic shear at second order including all relativistic effects was derived in [1]. In the present paper we compute the resulting cosmic shear bispectrum when all these effects are properly taken into account and we compare it to primordial non-Gaussianity of the local type. The new general relativistic effects are generically smaller than the standard non-linear couplings. However, their relative importance increases at small multipoles and for small redshifts of the sources. The dominant effect among these non standard corrections is due to the inhomogeneity of the source redshift. In the squeezed limit, its amplitude can become of the order of the standard couplings when the redshift of the sources is below 0.5. Moreover, while the standard non-linear couplings depend on the angle between the short and long mode, the relativistic corrections do not and overlap almost totally with local type non-Gaussianity. We find that they can contaminate the search for a primordial local signal by f_NL>10.

We derive an expression for the luminosity distance in a perturbed Friedmann universe. We define the correlation function and the power spectrum of the luminosity distance fluctuations and express them in terms of the initial spectrum of the Bardeen potential. We present semi-analytical results for the case of a pure CDM universe. We argue that the luminosity distance power spectrum represents a new observational tool which can be used to determine cosmological parameters. In addition, our results shed some light into the debate whether second order small scale fluctuations can mimic an accelerating universe.

We present an analytical derivation of the Sachs Wolfe effect sourced by a primordial magnetic field. In order to consistently specify the initial conditions, we assume that the magnetic field is generated by a causal process, namely a first order phase transition in the early universe. As for the topological defects case, we apply the general relativistic junction conditions to match the perturbation variables before and after the phase transition which generates the magnetic field, in such a way that the total energy momentum tensor is conserved across the transition and Einstein's equations are satisfied. We further solve the evolution equations for the metric and fluid perturbations at large scales analytically including neutrinos, and derive the magnetic Sachs Wolfe effect. We find that the relevant contribution to the magnetic Sachs Wolfe effect comes from the metric perturbations at next-to-leading order in the large scale limit. The leading order term is in fact strongly suppressed due to the presence of free-streaming neutrinos. We derive the neutrino compensation effect dynamically and confirm that the magnetic Sachs Wolfe spectrum from a causal magnetic field behaves as l(l+1)C_l^B \propto l^2 as found in the latest numerical analyses.

We present an analytical derivation of the Sachs Wolfe effect sourced by a primordial magnetic field, generated by a causal process, such as a first order phase transition in the early universe. As for the topological defects case, we apply the general relativistic junction conditions to match the perturbation variables before and after the phase transition, in such a way that the total energy momentum tensor is conserved across the transition. We find that the relevant contribution to the magnetic Sachs Wolfe effect comes from the metric perturbations at next-to-leading order in the large scale limit. The leading order term is strongly suppressed due to the presence of free-streaming neutrinos. We derive the neutrino compensation effect and confirm that the magnetic Sachs Wolfe spectrum from a causal magnetic field behaves as l(l+1)C_l^B ~ l^2 as found in the latest numerical analyses.

We compute the reduced cosmic shear up to second order in the gravitational potential without relying on the small angle or thin-lens approximation. This is obtained by solving the Sachs equation which describes the deformation of the infinitesimal cross-section of light bundle in the optical limit, and maps galaxy intrinsic shapes into their angular images. The calculation is done in the Poisson gauge without a specific matter content, including vector and tensor perturbations generated at second order and taking account of the inhomogeneities of a fixed redshift source plane. Our final result is expressed in terms of spin-2 operators on the sphere and is valid on the full sky. Beside the well known lens-lens and Born corrections that dominate on small angular scales, we find new non-linear couplings. These are a purely general relativistic intrinsic contribution, a coupling between the gravitational potential at the source with the lens, couplings between the time delay with the lens, couplings between two photon deflections, as well as non-linear couplings due to the second-order vector and tensor components. The inhomogeneity in the redshift of the source induces a coupling between the photon redshift with the lens. All these corrections become important on large angular scales and should thus be included when computing higher-order observables such as the bispectrum, in full or partially full-sky surveys.

It has been shown that generalized Einstein-Aether theories may lead to significant modifications to the non-relativistic limit of the Einstein equations. In this paper we study the effect of a general class of such theories on the Solar System. We consider corrections to the gravitational potential in negative and positive powers of distance from the source. Using measurements of the perihelion shift of Mercury and time delay of radar signals to Cassini, we place constraints on these corrections. We find that a subclass of generalized Einstein-Aether theories are compatible with these constraints.

We discuss some properties of Lorentz invariant theories which allow for superluminal motion. We show that, if signals are always sent forward in time, closed curves along which signals propagate can be formed. This leads to problems with causality and with the second law of thermodynamics. Only if one singles out one frame with respect to which all signals travel forward in time, the formation of 'closed signal curves' can be prevented. However, the price to pay is that in some reference frames perturbations propagate towards the past or towards the future, depending on the direction of emission.

We demonstrate that if k-essence can solve the coincidence problem and play the role of dark energy in the universe, the fluctuations of the field have to propagate superluminally at some stage. We argue that this implies that successful k-essence models violate causality. It is not possible to define a time ordered succession of events in a Lorentz invariant way. Therefore, k-essence cannot arise as low energy effective field theory of a causal, consistent high energy theory.

We show that the dipole of the luminosity distance is a useful observational tool which allows us to determine the Hubble parameter as a function of the redshift H(z). We determine the number of supernovae needed to achieve a given precision for H(z) and to distinguish between different models for dark energy. We analyse a sample of nearby supernovae and find a dipole consistent with the cosmic microwave background at a significance of more than 2 sigma.