This paper focuses on two aspects of the statistics of cosmological observables that are important for the next stages of precision cosmology. First, we note that the theory of reduced angular $N$-point spectra has only been developed in detail up to the trispectrum case and in a fashion that makes it difficult to go beyond. To fill this gap, here we present a constructive approach that provides a systematic description of reduced angular $N$-point spectra and their covariance matrices, for arbitrary $N$. Second, we focus on the common practice in the literature on cosmological observables, which consists in simply discarding a part of the expression, namely, the terms containing fields evaluated at the observer position. We point out that this is not justified beyond linear order in perturbation theory, as these terms contribute to all the multipoles of the corresponding spectra and with a magnitude that is of the same order as the rest of the non-linear corrections. We consider the possibility that the reason for neglecting these terms is a conceptual discomfort when using ensemble averages, which originates in an apparent tension between the ergodic hypothesis and the privileged position of the observer on the light-cone. We clarify this subtle issue by performing a careful derivation of the relation between the theoretical statistical predictions and the observational estimators for all $N$. We conclude that there is no inconsistency whatsoever in ensemble-averaging fields at and near the observer position, thus clearing the way for consistent and robust high-precision calculations.

The standard description of cosmological observables is incomplete, because it does not take into account the correct angular parametrization of the sky, i.e. the one determined by the observer frame. The corresponding corrections must be taken into account for reliable results at non-linear orders. This can be accomplished by introducing an orthonormal basis, or "tetrad", at the observer point, representing the frame with respect to which observations are performed. In this work we consider the tetrad formulation of General Relativity, thus associating tetrads to sources as well, and develop a new formalism for describing cosmological observables associated with localized sources. It is based on a manifold which we call the "observer space-time", whose coordinates are the proper time, redshift and angles an observer uses to parametrize measurements, and on which the rest of the observables are defined. This manifold does not have to be diffeomorphic to the true space-time and allows us to resolve caustics in the latter, in contrast to similar coordinate-based formalisms. As a concrete example, we work out the definitions and equations for the angular diameter distance, weak lensing and number count observables. As for the observables associated to the CMB, they lie inside the phase space distribution of the photon fluid, so we also revisit the construction of general-relativistic matrix kinetic theory from the tetrad formalism viewpoint. Here too the latter appears as the natural approach for relating the macroscopic dynamics to the microscopic quantum field theory, and therefore for constructing the matrix Boltzmann equations, without any approximation on the gravitational side. (... more in the manuscript)

This is the first in a series of papers devoted to fully general-relativistic $N$-body simulations applied to late-time cosmology. The purpose of this paper is to present the combination of a numerical relativity scheme, discretization method and time-integration algorithm that provides satisfyingly stable evolution. More precisely, we show that it is able to pass a robustness test and to follow scalar linear modes around an expanding homogeneous and isotropic space-time. Most importantly, it is able to evolve typical cosmological initial conditions on comoving scales down to tenths of megaparsecs with controlled constraint and energy-momentum conservation violations all the way down to the regime of strong inhomogeneity.

The background photon temperature $\bar T$ is one of the fundamental cosmological parameters. Despite its significance, $\bar T$ has never been allowed to vary in the data analysis, owing to the precise measurement of the comic microwave background (CMB) temperature by COBE FIRAS. However, even in future CMB experiments, $\bar T$ will remain unknown due to the unknown monopole contribution $\Theta_0$ at our position to the observed (angle-averaged) temperature $\langle T\rangle^{\rm obs}$. By fixing $\bar T\equiv\langle T\rangle^{\rm obs}$, the standard analysis underestimates the error bars on cosmological parameters, and the best-fit parameters obtained in the analysis are biased in proportion to the unknown amplitude of $\Theta_0$. Using the Fisher formalism, we find that these systematic errors are smaller than the error bars from the $Planck$ satellite. However, with $\bar T\equiv\langle T\rangle^{\rm obs}$, these systematic errors will always be present and irreducible, and future cosmological surveys might misinterpret the measurements.

We develop a theoretical framework to describe the cosmological observables on the past light cone such as the luminosity distance, weak lensing, galaxy clustering, and the cosmic microwave background anisotropies. We consider that all the cosmological observables include not only the background quantity, but also the perturbation quantity, and they are subject to cosmic variance, which sets the fundamental limits on the cosmological information that can be derived from such observables, even in an idealized survey with an infinite number of observations. To quantify the maximum cosmological information content, we apply the Fisher information matrix formalism and spherical harmonic analysis to cosmological observations, in which the angular and the radial positions of the observables on the light cone carry different information. We discuss the maximum cosmological information that can be derived from five different observables: (1) type Ia supernovae, (2) cosmic microwave background anisotropies, (3) weak gravitational lensing, (4) local baryon density, and (5) galaxy clustering. We compare our results with the cosmic variance obtained in the standard approaches, which treat the light cone volume as a cubic box of simultaneity. We discuss implications of our formalism and ways to overcome the fundamental limit.

We present the gauge-invariant formalism of cosmological weak lensing, accounting for all the relativistic effects due to the scalar, vector, and tensor perturbations at the linear order. While the light propagation is fully described by the geodesic equation, the relation of the photon wavevector to the physical quantities requires the specification of the frames, where they are defined. By constructing the local tetrad bases at the observer and the source positions, we clarify the relation of the weak lensing observables such as the convergence, the shear, and the rotation to the physical size and shape defined in the source rest-frame and the observed angle and redshift measured in the observer rest-frame. Compared to the standard lensing formalism, additional relativistic effects contribute to all the lensing observables. We explicitly verify the gauge-invariance of the lensing observables and compare our results to previous work. In particular, we demonstrate that even in the presence of the vector and tensor perturbations, the physical rotation of the lensing observables vanishes at the linear order, while the tetrad basis rotates along the light propagation compared to a FRW coordinate. Though the latter is often used as a probe of primordial gravitational waves, the rotation of the tetrad basis is indeed not a physical observable. We further clarify its relation to the E-B decomposition in weak lensing. Our formalism provides a transparent and comprehensive perspective of cosmological weak lensing.

Cosmological simulations involving the fully covariant gravitational dynamics may prove relevant in understanding relativistic/non-linear features and, therefore, in taking better advantage of the upcoming large scale structure survey data. We propose a new 3+1 integration scheme for General Relativity in the case where the matter sector contains a minimally-coupled perfect fluid field. The original feature is that we completely eliminate the fluid components through the constraint equations, thus remaining with a set of unconstrained evolution equations for the rest of the fields. This procedure does not constrain the lapse function and shift vector, so it holds in arbitrary gauge and also works for arbitrary equation of state. An important advantage of this scheme is that it allows one to define and pass an adaptation of the robustness test to the cosmological context, at least in the case of pressureless perfect fluid matter, which is the relevant one for late-time cosmology.

We perform a detailed study of the cosmological dynamics of a recently proposed infrared modification of the Einstein equations, based on the introduction of a non-local term constructed with $m^2g_{\mu\nu}\Box^{-1} R$, where $m$ is a mass parameter. The theory generates automatically a dynamical dark energy component, that can reproduce the observed value of the dark energy density without introducing a cosmological constant. Fixing $m$ so to reproduce the observed value $\Omega_{\rm DE}\simeq 0.68$, and writing $w(a)=w_0+(1-a) w_a$, the model provides a neat prediction for the equation of state parameters of dark energy, $w_0\simeq -1.042$ and $w_a\simeq -0.020$. We show that, because of some freedom in the definition of $\Box^{-1}$, one can extend the construction so to define a more general family of non-local models. However, in a first approximation this turns out to be equivalent to adding an explicit cosmological constant term on top of the dynamical dark energy component. This leads to an extended model with two parameters, $\Omega_{\Lambda}$ and $m$. Even in this extension the EOS parameter $w_0$ is always on the phantom side, in the range $-1.33 < w_0\leq -1$, and there is a prediction for the relation between $w_0$ and $w_a$.

Recent work has shown that non-local modifications of the Einstein equations can have interesting cosmological consequences and can provide a dynamical origin for dark energy, consistent with existing data. At first sight these theories are plagued by ghosts. We show that these apparent ghost-like instabilities do not describe actual propagating degrees of freedom, and there is no issue of ghost-induced quantum vacuum decay.

We construct a fully covariant theory of massive gravity which does not require the introduction of an external reference metric, and overcomes the usual problems of massive gravity theories (fatal ghosts instabilities, acausality and/or vDVZ discontinuity). The equations of motion of the theory are non-local, but respect causality. The starting point is the quadratic action proposed in the context of the degravitation idea. We show that it is possible to extended it to a fully non-linear covariant theory. This theory describes the five degrees of freedom of a massive graviton plus a scalar ghost. However, contrary to generic non-linear extensions of Fierz-Pauli massive gravity, the ghost has the same mass m as the massive graviton, independently of the background, and smoothly goes into a non-radiative degree of freedom for m-> 0. As a consequence, for $m\sim H_0$ the vacuum instability induced by the ghost is irrelevant even over cosmological time-scales. We finally show that an extension of the model degravitates a vacuum energy density of order $M_{Planck}^4$ down to a value of order $M_{Planck}^2 m^2$, which for $m \sim H_0$ is of order of the observed value of the vacuum energy density.

We re-examine the classic problem of the renormalization of zero-point quantum fluctuations in a Friedmann-Robertson-Walker background. We discuss a number of issues that arise when regularizing the theory with a momentum-space cutoff, and show explicitly how introducing non-covariant counter-terms allows to obtain covariant results for the renormalized vacuum energy-momentum tensor. We clarify some confusion in the literature concerning the equation of state of vacuum fluctuations. Further, we point out that the general structure of the effective action becomes richer if the theory contains a scalar field phi with mass m smaller than the Hubble parameter H(t). Such an ultra-light particle cannot be integrated out completely to get the effective action. Apart from the volume term and the Einstein-Hilbert term, that are reabsorbed into renormalizations of the cosmological constant and Newton's constant, the effective action in general also has a term proportional to F(phi)R, for some function F(phi). As a result, vacuum fluctuations of ultra-light scalar fields naturally lead to models where the dark energy density has the form rho_{DE}(t)=rho_X(t)+rho_Z(t), where rho_X is the component that accelerates the Hubble expansion at late times and rho_Z(t) is an extra contribution proportional to H^2(t). We perform a detailed comparison of such models with CMB, SNIa and BAO data.

We examine a cosmological model with a dark energy density of the form $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$, where $\rho_X$ is the component that accelerates the Hubble expansion at late times and $\rho_Z(t)$ is an extra contribution proportional to $H^2(t)$. This form of $\rho_Z(t)$ follows from the recent proposal that the contribution of zero-point fluctuations of quantum fields to the total energy density should be computed by subtracting the Minkowski-space result from that computed in the FRW space-time. We discuss theoretical arguments that support this subtraction. By definition, this eliminates the quartic divergence in the vacuum energy density responsible for the cosmological constant problem. We show that the remaining quadratic divergence can be reabsorbed into a redefinition of Newton's constant only under the assumption that the energy-momentum tensor of vacuum fluctuations is conserved in isolation. However, in the presence of an ultra-light scalar field $X$ with $m_X<H_0$, as typical of some dark energy models, the gravity effective action depends both on the gravitational field and on the $X$ field. In this case general covariance only requires the conservation of the total energy-momentum tensor, including both the classical term $T^X_{\mu\nu}$ and the vacuum expectation value of T_{\mu\nu}. If there is an exchange of energy between these two terms, there are potentially observable consequences. We construct an explicit model with an interaction between $\rho_X$ and $\rho_Z$ and we show that the total dark energy density $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$ always remains a finite fraction of the critical density at any time, providing a specific model of early dark energy. We discuss the implication of this result for the coincidence problem and we estimate the model parameters by means of a full likelihood analysis using current CMB, SNe Ia and BAO data.