Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation and the possibility of constraining novel gravitational phenomena beyond general relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum are insufficient, and we must consider a fully relativistic and gauge-independent treatment of observables such as galaxy number counts in order to avoid subtle biases, e.g. in the determination of the fNL parameter.In this work, we present an initial application of an analysis pipeline capable of accurately modelling and recovering relativistic spectra and correlation functions. As a proof of concept, we focus on the non-zero dipole of the redshift-space power spectrum that arises in the cross-correlation of different mass bins of dark matter haloes, using strictly gauge-independent observable quantities evaluated on the past light cone of a fully relativistic N-body simulation in a redshift bin 1.7 ≤ z ≤ 2.9. We pay particular attention to the correct estimation of power spectrum multipoles, comparing different methods of accounting for complications such as the survey geometry (window function) and evolution/bias effects on the past light cone, and discuss how our results compare with previous attempts at extracting novel GR signatures from relativistic simulations.

The standard cosmological model is inherently relativistic, and yet a wide range of cosmological observations can be predicted accurately from essentially Newtonian theory. This is not the case on ‘ultra-large’ distance scales, around the cosmic horizon size, however, where relativistic effects can no longer be neglected. In this paper, we present a novel suite of 53 fully relativistic simulations generated using the gevolution code, each covering the full sky out to z ≈ 0.85, and approximately 1930 square degrees out to z ≈ 3.55. These include a relativistic treatment of massive neutrinos, as well as the gravitational potential that can be used to exactly calculate observables on the past light cone. The simulations are divided into two sets, the first being a set of 39 simulations of the same fiducial cosmology (based on the Euclid Flagship 2 cosmology) with different realisations of the initial conditions, and the second which fixes the initial conditions, but varies each of seven cosmological parameters in turn. Taken together, these simulations allow us to perform statistical studies and calculate derivatives of any relativistic observable with respect to cosmological parameters. As an example application, we compute the cross-correlation between the Doppler magnification term in the convergence, κv, and the CDM+baryon density contrast, δcb, which arises only in a (special) relativistic treatment. We are able to accurately recover this term as predicted by relativistic perturbation theory, and study its sample variance and derivatives with respect to cosmological parameters.

No abstract

The Cosmological Principle (CP) is part of the foundation that underpins the standard model of the Universe. In the era of precision cosmology, when stress tests of the standard model are uncovering various tensions and possible anomalies, it is critical to check the viability of this principle. A key test is the consistency between the kinematic dipoles of the cosmic microwave background and of the large-scale matter distribution. Results using radio continuum and quasar samples indicate a rough agreement in the directions of the two dipoles, but a larger than expected amplitude of the matter dipole. The resulting tension with the radiation dipole has been estimated at ∼5σ for some cases, suggesting a potential new cosmological tension and a possible violation of the CP. However, the standard formalism for predicting the dipole in the two-dimensional projection of sources overlooks possible evolution effects in the luminosity function. In fact, radial information from the luminosity function is necessary for a correct projection of the three-dimensional source distribution. Using a variety of current models of the quasar luminosity function, we show that neglecting redshift evolution can significantly overestimate the relative velocity amplitude. While the models we investigate are consistent with each other and with current data, the dipole derived from these, which depends on derivatives of the luminosity function, can disagree by more than 3σ. This theoretical systematic bias needs to be resolved before robust conclusions can be made about a new cosmic tension.

Early Universe physics leaves distinct imprints on the Cosmic Microwave Background (CMB) and Large-Scale Structure (LSS). The current cosmological paradigm to explain the origin of the structures we see in the Universe today (CMB and LSS), named Inflation, says that the Universe went through a period of accelerated expansion. Density fluctuations that eventually have grown into the temperature fluctuations of the CMB and the galaxies and other structures we see in the LSS come from the quantization of the scalar field (inflaton) which provokes the accelerated expansion. The most simple inflationary model, which contains only one slowly-rolling scalar field with canonical kinetic term in the action, produces a power-spectrum (Fourier transform of the two-point correlation function) approximately scale invariant and an almost null bispectrum (Fourier transform of the three-point correlation function). This characteristic is called Gaussianity, once random fields that follow a normal distribution have all the odd moments null. Yet, more complex inflationary models (with more scalar fields and/or non-trivial kinetic terms in the action, etc) and possible alternatives to inflation have a non-vanishing bispectrum which can be parametrized by a non-linearity parameter f NL , whose value differs from model to model. In this work we studied the basic ingredients to understand such statements and focused on the observational evidences of this parameters and how the current and upcoming galaxy surveys are able to impose constraints to the value of f NL with a better accuracy, through the multi-tracer technique, than those obtained by means of CMB measurements.