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.

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 point out the existence of a second feature in the matter two-point function, besides the acoustic peak, due to the baryon-baryon correlation in the early universe and positioned at twice the distance of the peak. We discuss how the existence of this feature is implied by the well-known heuristic argument that explains the baryon bump in the correlation function. A standard $\chi^2$ analysis to estimate the detection significance of the second feature is mimicked. We conclude that for realistic values of the baryon density, an SKA-like galaxy survey will not be able to detect this feature with standard correlation function analysis.

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.

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.

Primordial inflation may represent the most powerful collider to test high-energy physics models. In this paper we study the impact on the inflationary power spectrum of the comoving curvature perturbation in the specific model where massive higher spin fields are rendered effectively massless during a de Sitter epoch through suitable couplings to the inflaton field. In particular, we show that such fields with spin $s$ induce a distinctive statistical anisotropic signal on the power spectrum, in such a way that not only the usual $g_{2M}$-statistical anisotropy coefficients, but also higher-order ones (i.e., $g_{4M}$, $g_{6M}$, $\cdots$, $g_{(2s-2)M}$ and $g_{(2s) M}$) are nonvanishing. We examine their imprints in the cosmic microwave background and galaxy power spectra. Our Fisher matrix forecasts indicate that the detectability of $g_{LM}$ depends very weakly on $L$: all coefficients could be detected in near future if their magnitudes are bigger than about $10^{-3}$.

We compute the generation of vorticity from velocity dispersion in the dark matter fluid. For dark matter at zero temperature Helmholtz's theorem dictates that no vorticity is generated and we therefore allow the dark matter fluid to have a non-vanishing velocity dispersion. This implies a modification to the usual hydrodynamical system (continuity and Euler equations): we have to consider the Boltzmann hierarchy up to the second moment. This means that the Euler equation is modified with a source term that describes the effect of non-zero velocity dispersion. We write an equation for the Eulerian vorticity in Lagrangian coordinates and show that it has a growing mode already at second order in perturbation theory. We compute the power spectrum of the vorticity and the rotational velocity at second order in perturbation theory.

We derive the contribution to relativistic galaxy number count fluctuations from vector and tensor perturbations within linear perturbation theory. Our result is consistent with the the relativistic corrections to number counts due to scalar perturbation, where the Bardeen potentials are replaced with line-of-sight projection of vector and tensor quantities. Since vector and tensor perturbations do not lead to density fluctuations the standard density term in the number counts is absent. We apply our results to vector perturbations which are induced from scalar perturbations at second order and give numerical estimates of their contributions to the power spectrum of relativistic galaxy number counts.

We compute the angular power spectra of the E-type and B-type lensing potentials for gravitational waves from inflation and for tensor perturbations induced by scalar perturbations. We derive the tensor-lensed CMB power spectra for both cases. We also apply our formalism to determine the linear lensing potential for a Bianchi I spacetime with small anisotropy.