Above the equality scale the galaxy bispectrum will be a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. By decomposing the Fourier-space bispectrum into invariant multipoles about the observer's line of sight we examine in detail how the relativistic effects contribute to these. We show how to perform this decomposition analytically, which is significantly faster for subsequent computations. While all multipoles receive a contribution from the relativistic part, odd multipoles arising from the imaginary part of the bispectrum have no Newtonian contribution, making the odd multipoles a smoking gun for a relativistic signature in the bispectrum for single tracers. The dipole and the octopole are significant on equality scales and above where the Newtonian approximation breaks down. This breakdown is further signified by the fact that the even multipoles receive a significant correction on very large scales.

The Fourier-space galaxy bispectrum is complex, with the imaginary part arising from leading-order relativistic corrections, due to Doppler, gravitational redshift and related line-of-sight effects in redshift space. The detection of the imaginary part of the bispectrum is potentially a smoking gun signal of relativistic contributions. We investigate whether next-generation spectroscopic surveys could make such a detection. For a Stage IV spectroscopic $H\alpha$ survey similar to Euclid, we find that the cumulative signal to noise of this relativistic signature is $\mathcal{O}(10)$. Long-mode relativistic effects couple to short-mode Newtonian effects in the galaxy bispectrum, but not in the galaxy power spectrum. This is the basis for detectability of relativistic effects in the bispectrum of a single galaxy survey, whereas the power spectrum requires multiple galaxy surveys to detect the corresponding signal.

We present a systematic study of galaxy bias in the presence of primordial non-Gaussianity in General Relativity (GR) at second order in perturbation theory. The non-linearity of the Poisson equation in GR and primordial non-Gaussianity are consistently included. We show that the inclusion of non-local primordial non-Gaussianity in addition to local non-Gaussianity is important to show the absence of the modulation of small scale clustering by the long-wavelength mode in the single field slow-roll inflation. We study the bispectrum of the relativistic galaxy density in several gauges and identify the effect of primordial non-Gaussianity and GR corrections.

The local galaxy bias formalism relies on the energy constraint equation at the formation time to relate the metric perturbation to the matter density contrast. In the Newtonian approximation, this relationship is linear, which allows us to specify the initial galaxy density as a function of local physical operators. In general relativity however, the relationship is intrinsically nonlinear and a modulation of the short-wavelength mode by the long-wavelength mode might be expected. We describe in detail how to obtain local coordinates where the coupling of the long- to the short-wavelength modes is removed through a change of coordinates (in the absence of the primordial non-Gaussianity). We derive the general-relativistic correction to the galaxy bias expansion at second order. The correction does not come from the modulation of small-scale clustering by the long-wavelength mode; instead, it arises from distortions of the volume element by the long-wavelength mode and it does not lead to new bias parameters.

The bispectrum will play an important role in future galaxy surveys. On large scales it is a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. We show for the first time that the leading relativistic correction from these distortions in the galaxy bispectrum generates a significant dipole, mainly from relativistic redshift space distortions. The amplitude of the dipole can be more than 10% of the monopole even on equality scales. Such a dipole is absent in the Newtonian approximation to the redshift space bispectrum, so it offers a clear signature of relativistic effects on cosmological scales in large scale structure.

We compute the redshift-dependent angular bispectrum of galaxy number counts at tree-level, including nonlinear clustering bias and estimating numerically for the first time the effect of redshift space distortions (RSD). We show that for narrow redshift bins the amplitude of nonlinear RSD is comparable with the matter density perturbations. While our numerical results only include terms relevant on sub-horizon scales, the formalism can readily be extended to the full tree-level bispectrum. Our approach does not rely on the flat-sky approximation and it can be easily generalized to different sources by including the appropriate bias expansion. We test the accuracy of Limber approximation for different z-bins. We highlight the subtle but relevant differences in the angular bispectrum of galaxy number counts with respect to CMB, due to the different scale dependence of perturbations. Our formalism can also be directly applied to the angular HI intensity mapping bispectrum.

The galaxy bispectrum on scales around and above the equality scale receives contributions from relativistic effects. Some of these arise from lightcone deformation effects, which come from local and line-of-sight integrated contributions. Here we calculate the local contributions from the generated vector and tensor background which is formed as scalar modes couple and enter the horizon. We show that these modes are sub-dominant when compared with other relativistic contributions.

Relativistic contributions to the dynamics of structure formation come in a variety of forms, and can potentially give corrections to the standard picture on typical scales of 100 Mpc. These corrections cannot be obtained by Newtonian numerical simulations, so it is important to accurately estimate the magnitude of these relativistic effects. Density fluctuations couple to produce a background of gravitational waves, which is larger than any primordial background. A similar interaction produces a much larger spectrum of vector modes which represent the frame-dragging rotation of spacetime. These can change the metric at the percent level in the concordance model at scales below the equality scale. Vector modes modify the lensing of background galaxies by large-scale structure. This gives in principle the exciting possibility of measuring relativistic frame dragging effects on cosmological scales. The effects of the non-linear tensor and vector modes on the cosmic convergence are computed and compared to first-order lensing contributions from density fluctuations, Doppler lensing, and smaller Sachs-Wolfe effects. The lensing from gravitational waves is negligible so we concentrate on the vector modes. We show the relative importance of this for future surveys such as Euclid and SKA. We find that these non-linear effects only marginally affect the overall weak lensing signal so they can safely be neglected in most analyses, though are still much larger than the linear Sachs-Wolfe terms. The second-order vector contribution can dominate the first-order Doppler lensing term at moderate redshifts and are actually more important for survey geometries like the SKA.

Using the consistency relation in Fourier space, we derive the observed galaxy bispectrum from single-field inflation in the squeezed limit, in which one of the three modes has a wavelength much longer than the other two. This provides a non-trivial check of the full computation of the bispectrum based on second-order cosmological perturbation theory in this limit. We show that gauge modes need to be carefully removed in the second-order cosmological perturbations in order to calculate the observed galaxy bispectrum in the squeezed limit. We then give an estimate of the effective non-Gaussianity due to general relativistic lightcone effects that could mimic a primordial non-Gaussian signal.

The galaxy bispectrum is affected on equality scales and above by relativistic observational effects, at linear and nonlinear order. These lightcone effects include local contributions from Doppler and gravitational potential terms, as well as integrated contributions like lensing, together with all the couplings at nonlinear order. We recently presented the correction to the galaxy bispectrum from all local lightcone effects up to second order in perturbations, using a plane-parallel approximation. Here we update our previous result by including the effects from relativistic nonlinear dynamical evolution. We show that these dynamical effects make a significant contribution to the projection effects.

Neutral hydrogen (HI) will soon be the dark matter tracer observed over the largest volumes of Universe thanks to the 21 cm intensity mapping technique. To unveil cosmological information it is indispensable to understand the HI distribution with respect to dark matter. Using a full one-loop derivation of the power spectrum of HI, we show that higher order corrections change the amplitude and shape of the power spectrum on typical cosmological (linear) scales. These effects go beyond the expected dark matter non-linear corrections and include non-linearities in the way the HI signal traces dark matter. We show that, on linear scales at z = 1, the HI bias drops by up to 15% in both real and redshift space, which results in underpredicting the mass of the halos in which HI lies. Non-linear corrections give rise to a significant scale dependence when redshift space distortions arise, in particular on the scale range of the baryonic acoustic oscillations (BAO). There is a factor of 5 difference between the linear and full HI power spectra over the full BAO scale range, which will modify the ratios between the peaks. This effect will also be seen in other types of survey and it will be essential to take it into account in future experiments in order to match the expectations of precision cosmology.

General relativistic imprints on the galaxy bispectrum arise from observational (or projection) effects. The lightcone projection effects include local contributions from Doppler and gravitational potential terms, as well as lensing and other integrated contributions. We recently presented for the first time, the correction to the galaxy bispectrum from all local lightcone projection effects up to second order in perturbations. Here we provide the details underlying this correction, together with further results and illustrations. For moderately squeezed shapes, the correction to the Newtonian prediction is ~30% on equality scales at z ~ 1. We generalise our recent results to include the contribution, up to second order, of magnification bias (which affects some of the local terms) and evolution bias.

Intensity mapping of the HI brightness temperature provides a unique way of tracing large-scale structures of the Universe up to the largest possible scales. This is achieved by using a low angular resolution radio telescopes to detect emission line from cosmic neutral Hydrogen in the post-reionization Universe. We use general relativistic perturbation theory techniques to derive for the first time the full expression for the HI brightness temperature up to third order in perturbation theory without making any plane-parallel approximation. We use this result and the renormalization prescription for biased tracers to study the impact of nonlinear effects on the power spectrum of HI brightness temperature both in real and redshift space. We show how mode coupling at nonlinear order due to nonlinear bias parameters and redshift space distortion terms modulate the power spectrum on large scales. The large scale modulation may be understood to be due to the effective bias parameter and effective shot noise.

Next-generation galaxy surveys will increasingly rely on the galaxy bispectrum to improve cosmological constraints, especially on primordial non-Gaussianity. A key theoretical requirement that remains to be developed is the analysis of general relativistic effects on the bispectrum, which arise from observing galaxies on the past lightcone, {as well as from relativistic corrections to the dynamics}. {As an initial step towards a fully relativistic analysis of the galaxy bispectrum, we compute for the first time the local relativistic lightcone effects on the bispectrum,} which come from Doppler and gravitational potential contributions. For the galaxy bispectrum, the problem is much more complex than for the power spectrum, since we need the lightcone corrections at second order. Mode-coupling contributions at second order mean that relativistic corrections can be non-negligible at smaller scales than in the case of the power spectrum. In a primordial Gaussian universe, we show that the local lightcone corrections for squeezed shapes at $z\sim1$ mean that the bispectrum can differ from the Newtonian prediction by $\gtrsim 10\%$ when the short modes are $k\lesssim (50\,{\rm Mpc})^{-1}$. These relativistic projection effects, if ignored in the analysis of observations, could be mistaken for primordial non-Gaussianity. For upcoming surveys which probe equality scales and beyond, {all relativistic lightcone effects and relativistic dynamical corrections should be included} for an accurate measurement of primordial non-Gaussianity.

Intensity mapping of the neutral hydrogen brightness temperature promises to provide a three-dimensional view of the universe on very large scales. Nonlinear effects are typically thought to alter only the small-scale power, but we show how they may bias the extraction of cosmological information contained in the power spectrum on ultra-large scales. For linear perturbations to remain valid on large scales, we need to renormalize perturbations at higher order. In the case of intensity mapping, the second-order contribution to clustering from weak lensing dominates the nonlinear contribution at high redshift. Renormalization modifies the mean brightness temperature and therefore the evolution bias. It also introduces a term that mimics white noise. These effects may influence forecasting analysis on ultra-large scales.

We consider the entropy associated with the large-scale structure of the Universe in the linear regime, where the Universe can be described by a perturbed Friedmann-Lema\^itre spacetime. In particular, we compare two different definitions proposed in the literature for the entropy using a spatial averaging prescription. For one definition, the entropy of the large-scale structure for a given comoving volume always grows with time, both for a CDM and a $\Lambda$CDM model. In particular, while it diverges for a CDM model, it saturates to a constant value in the presence of a cosmological constant. The use of a light-cone averaging prescription in the context of the evaluation of the entropy is also discussed.

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.

The success of precision cosmology depends not only on accurate observations, but also on the theoretical model - which must be understood to at least the same level of precision. Subtle relativistic effects can lead to biased measurements if they are neglected. One such effect gives a systematic shift in the distance-redshift relation away from its background value, due to the non- linear relativistic conservation of total photon flux. We also show directly how this shift follows from a fully relativistic analysis of the geodesic deviation equation. We derive the expectation value of the shift using second-order perturbations about a concordance background, and show that the distance to last scattering is increased by 1%. We argue that neglecting this shift could lead to a significant bias in the background cosmological parameters, because it alters the meaning of the background model. A naive adjustment of CMB parameter estimation if this shift is really a correction to the background would raise the H0 value inferred from the CMB by 5%, potentially removing the tension with local measurements of H0. Other CMB parameters which depend on the distance would also be shifted by ~1sigma when combined with local H0 data. While our estimations rely on a simplistic analysis, they nevertheless illustrate that accurately defining the background model in terms of the expectation values of observables is critical when we aim to determine the model parameters at the sub-percent level.

We present a derivation of the cosmological distance-redshift relation up to second order in perturbation theory. In addition, we find the observed redshift and the lensing magnification to second order. We do not require that the density contrast is small, we only that the metric potentials and peculiar velocities are small. Thus our results apply into the nonlinear regime, and can be used for most dark energy models. We present the results in a form which can be readily computed in an N-body simulation. This paper accompanies arXiv:1207.2109, where the key results are summarised in a physically transparent form and applications are discussed.

The next generation of telescopes will usher in an era of precision cosmology, capable of determining the cosmological model to beyond the percent level. For this to be effective, the theoretical model must be understood to at least the same level of precision. A range of subtle relativistic effects remain to be explored theoretically, and offer the potential for probing general relativity in this new regime. We present the distance-redshift relation to second order in cosmological perturbation theory for a general dark energy model. This relation determines the magnification of sources at high precision, as well as redshift space distortions in the mildly non-linear regime. We identify a range of new lensing effects, including: double-integrated and nonlinear integrated Sach-Wolfe contributions, transverse Doppler effects, lensing from the induced vector mode and gravitational wave backgrounds, in addition to lensing from the second-order potential. Modifications to Doppler lensing from redshift-space distortions are identified. Finally, we find a new double-coupling between the density fluctuations integrated along the line of sight, and gradients in the density fluctuations coupled to transverse velocities along the line of sight. These can be large and thus offer important new probes of gravitational lensing and general relativity. This paper accompanies arXiv:1402.1933, where a comprehensive derivation is given.

In this paper, we develop in detail a fully geometrical method for deriving perturbation equations about a spatially homogeneous background. This method relies on the 3+1 splitting of the background space-time and does not use any particular set of coordinates: it is implemented in terms of geometrical quantities only, using the tensor algebra package xTensor in the xAct distribution along with the extension for perturbations xPert. Our algorithm allows one to obtain the perturbation equations for all types of homogeneous cosmologies, up to any order and in all possible gauges. As applications, we recover the well-known perturbed Einstein equations for Friedmann-Lemaitre-Robertson-Walker cosmologies up to second order and for Bianchi I cosmologies at first order. This work paves the way to the study of these models at higher order and to that of any other perturbed Bianchi cosmologies, by circumventing the usually too cumbersome derivation of the perturbed equations.

Light from `point sources' such as supernovae is observed with a beam width of order of the sources' size - typically less than 1 AU. Such a beam probes matter and curvature distributions that are very different from coarse-grained representations in N-body simulations or perturbation theory, which are smoothed on scales much larger than 1 AU. The beam typically travels through unclustered dark matter and hydrogen with a mean density much less than the cosmic mean, and through dark matter halos and hydrogen clouds. Using N-body simulations, as well as a Press-Schechter approach, we quantify the density probability distribution as a function of beam width and show that, even for Gpc-length beams of 500 kpc diameter, most lines of sight are significantly under-dense. From this we argue that modelling the probability distribution for AU-diameter beams is absolutely critical. Standard analyses predict a huge variance for such tiny beam sizes, and nonlinear corrections appear to be non-trivial. It is not even clear whether under-dense regions lead to dimming or brightening of sources, owing to the uncertainty in modelling the expansion rate which we show is the dominant contribution. By considering different reasonable approximations which yield very different cosmologies we argue that modelling ultra-narrow beams accurately remains a critical problem for precision cosmology. This could appear as a discordance between angular diameter and luminosity distances when comparing SN observations to BAO or CMB distances.

Structure occurs over a vast range of scales in the universe. Our large-scale cosmological models are coarse-grained representations of what exists, which have much less structure than there really is. An important problem for cosmology is determining the influence the small-scale structure in the universe has on its large-scale dynamics and observations. Is there a significant, general relativistic, backreaction effect from averaging over structure? One issue is whether the process of smoothing over structure can contribute to an acceleration term and so alter the apparent value of the cosmological constant. If this is not the case, are there other aspects of concordance cosmology that are affected by backreaction effects? Despite much progress, this 'averaging problem' is still unanswered, but it cannot be ignored in an era of precision cosmology.

Smoothing over structures in general relativity leads to a renormalisation of the background, and potentially many other effects which are poorly understood. Observables such as the distance-redshift relation when averaged on the sky do not necessarily yield the same smooth model which arises when performing spatial averages. These issues are thought to be of technical interest only in the standard model of cosmology, giving only tiny corrections. However, when we try to calculate observable quantities such as the all-sky average of the distance-redshift relation, we find that perturbation theory delivers divergent answers in the UV and corrections to the background of order unity. There are further problems. Second-order perturbations are the same size as first-order, and fourth-order at least the same as second, and possibly much larger, owing to the divergences. Much hinges on a coincidental balance of 2 numbers: the primordial power, and the ratio between the comoving Hubble scales at matter-radiation equality and today. Consequently, it is far from obvious that backreaction is irrelevant even in the concordance model, however natural it intuitively seems.

The calculation of the averaged Hubble expansion rate in an averaged perturbed Friedmann-Lemaitre-Robertson-Walker cosmology leads to small corrections to the background value of the expansion rate, which could be important for measuring the Hubble constant from local observations. It also predicts an intrinsic variance associated with the finite scale of any measurement of H_0, the Hubble rate today. Both the mean Hubble rate and its variance depend on both the definition of the Hubble rate and the spatial surface on which the average is performed. We quantitatively study different definitions of the averaged Hubble rate encountered in the literature by consistently calculating the backreaction effect at second order in perturbation theory, and compare the results. We employ for the first time a recently developed gauge-invariant definition of an averaged scalar. We also discuss the variance of the Hubble rate for the different definitions.