Fiducial Cosmology Research Articles

Alcock–Paczyński effect on void-finding

Under the assumption of statistical isotropy, and in the absence of directional selection effects, a stack of voids is expected to be spherically symmetric, which makes it an excellent object to use for an Alcock–Paczyński (AP) test. This test is commonly carried out using the void-galaxy cross-correlation function (CCF), which has emerged as a competitive probe, especially in combination with the galaxy-galaxy auto-correlation function. Current studies of the AP effect around voids assume that void-centre positions are influenced by the choice of fiducial cosmology in the same way as galaxy positions. We show that this assumption, though prevalent in the literature, is complicated by the response of void-finding algorithms to shifts in tracer positions. Using stretched simulation boxes to emulate the AP effect, we investigate how the void-galaxy CCF changes due to its presence, revealing an additional effect imprinted in the CCF that must be accounted for. The effect originates from the response of void finders to the distorted tracer field – which leads to reduction of the amplitude of the AP signal in the CCF – and thus depends on the specific void-finding algorithm used. We present results for four different void-finding packages, namely REVOLVER, VIDE, voxel, and the spherical void finder in the Pylians3 library, demonstrating how incorrect treatment of the AP effect results in biases in the recovered parameters, regardless of the technique used. Finally, we propose a method to alleviate this issue without resorting to complex and finder-specific modelling of the void-finder response to AP.

The clustering of Lyman Alpha Emitting galaxies at 𝗓=2–3

We measure the clustering of Lyman Alpha Emitting galaxies (LAEs) selected from the One-hundred-square-degree DECam Imaging in Narrowbands (ODIN) survey, with spectroscopic follow-up from Dark Energy Spectroscopic Instrument (DESI). We use DESI spectroscopy to optimize our selection and to constrain the interloper fraction and redshift distribution of our narrow-band selected sources. We select samples of 4000 LAEs at z = 2.45 and 3.1 in 9 sq.deg. centered on the COSMOS field with median Lyα fluxes of ≈ 10-16 erg s-1 cm-2. Covariances and cosmological inferences are obtained from a series of mock catalogs built upon high-resolution N-body simulations that match the footprint, number density, redshift distribution and observed clustering of the sample. We find that both samples have a correlation length of r 0 = 3.0 ± 0.2 h-1 Mpc. Within our fiducial cosmology these correspond to 3D number densities of ≈ 10-3 h3 Mpc-3 and, from our mock catalogs, biases of 1.7 and 2.0 at z = 2.45 and 3.1, respectively. We discuss the implications of these measurements for the use of LAEs as large-scale structure tracers for high-redshift cosmology.

CLASS-OneLoop: accurate and unbiased inference from spectroscopic galaxy surveys

The power spectrum is the most commonly applied summary statistics to extract cosmological information from the observed three-dimensional distribution of galaxies in spectroscopic surveys. We present CLASS-OneLoop, a new numerical tool, fully integrated into the Boltzmann code CLASS, enabling the calculation of the one-loop power spectrum of biased tracers in spectroscopic surveys. Built upon the Eulerian moment expansion framework for redshift-space distortions, the implemented model incorporates a complete set of nonlinear biases, counterterms, and stochastic contributions, and includes the infrared resummation and the Alcock-Paczynski effect. The code features an evaluation of the loops by either direct numerical integration or Fast Fourier Transform, and employs a fast-slow parameter decomposition, which is essential for accelerating MCMC runs. After presenting performance and validation tests, as an illustration of the capabilities of the code, we apply it to fit the measured redshift-space halo power spectrum wedges on a ΛCDM subset of the AbacusSummit simulation suite and considering scales up to kmax = 0.3 h/Mpc. We find that the one-loop model adeptly recovers the fiducial cosmology of the simulation, while a simplified model commonly used in the literature for sensitivity forecasts yields significantly biased results. Furthermore, we conduct Monte Carlo Markov Chain (MCMC) forecasts for a DESI-like survey, considering a model with a dynamical dark energy component. Our results demonstrate the ability to independently constrain cosmological and nuisance parameters, even in the presence of a large parameter space with twenty-nine variables.

Angular correlation function from sample covariance with BOSS and eBOSS LRG

The Baryon Acoustic Oscillations (BAO) are one of the most used probes to understand the accelerated expansion of the Universe. Traditional methods rely on fiducial model information within their statistical analysis, which may be a problem when constraining different families of models. This work aims to provide a method that constrains θBAO\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{BAO}$$\\end{document} through a model-independent approach using the covariance matrix from the galaxy sample from thin redshift bins, later validated with a mock sample covariance matrix. We used widths of δz=0.002\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\delta z = 0.002$$\\end{document} separation for all samples as the basis for a sample covariance matrix weighted by the statistical importance of the redshift bin. Each sample belongs to the Sloan Digital Sky Survey: BOSS1, BOSS2, and eBOSS, with effective redshift zeff\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z_{eff}$$\\end{document}: 0.35, 0.51, 0.71, and different numbers of bins with 50, 100, and 200. To get θBAO\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{BAO}$$\\end{document}, we correct the angular separation from the polynomial fit (θfit\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{fit}$$\\end{document}) by comparing each bin correlation function with the correlation function of the whole set, a parameter named α~\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ ilde{\\alpha }$$\\end{document}. We also tested such correction by choosing the bin at zeff\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z_{eff}$$\\end{document} and found that for eBOSS θBAO\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{BAO}$$\\end{document} is in 1σ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1 \\sigma $$\\end{document} agreement with the Planck 18 model. Finally, we found that the sample covariances are noisy compared to the mocks for lower z samples, something expected due to nonlinear effects. Such noise impact can be seen in the parameter constraints but does not affect the eBOSS covariance sample. It is shown that mocks’ results do tend to its chosen fiducial cosmology θBAO\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{BAO}$$\\end{document}. BOSS1 and BOSS2 showed agreement with Planck 18 and an agreement with Pantheon + SH0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_0$$\\end{document}ES when α~\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ ilde{\\alpha }$$\\end{document} is based on the bin z=zeff\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z=z_{eff}$$\\end{document}.

BAO cosmology in non-spatially flat background geometry from BOSS+eBOSS and lessons for future surveys

We study the impact of the assumption of a non-flat fiducial cosmology on the measurement, analysis and interpretation of BAO distance variables, along and across the line-of-sight. The assumption about cosmology enters in the choice of the base template, as well as on the transformation of tracer's redshifts into distances (the catalog cosmology): here we focus on the curvature assumption, separately and jointly, on both.We employ BOSS and eBOSS publicly available data and show that for the statistical precision of this data set, distance measures and thus cosmological inference are robust to assumptions about curvature both of the template and the catalog. Thus the usual assumptions of flat fiducial cosmologies (but also assumptions of non-flat cosmologies) do not produce any detectable systematic effects.For forthcoming large-volume surveys, however, small but appreciable residual systematic shifts can be generated which may require some care. These are mostly driven by the choice of catalog cosmology if it is significantly different from true cosmology.In particular, the catalog (and template) cosmology should be chosen, possibly iteratively, in such a way that the recovered BAO scaling variables are sufficiently close to unity. At this level of precision, however, other previously overlooked effects become relevant, such as a mismatch between the sound horizon as seen in the BAO and the actual sound horizon in the early Universe. If unaccounted for, such effect may be misinterpreted as cosmological and thus bias the curvature (and cosmology) constraints. We present best practices to avoid this.

Probing cosmic homogeneity in the Local Universe

ABSTRACT We investigate the transition scale to homogeneity, RH, using as cosmic tracer the spectroscopic sample of blue galaxies from the Sloan Digital Sky Survey (SDSS). Considering the spatial distribution of the galaxy sample, we compute the two point correlation function ξ(r), the scaled counts in spheres $\mathcal {N}(\lt r)$, and the fractal dimension $\mathcal {D}_2(r)$ to quantify the homogeneity scale in the Local Universe (0.04 < z < 0.20). The sample in analysis is compared with random and mock catalogues with the same geometry, and the same number of synthetic cosmic objects as the data set, to calculate the covariance matrix for the errors determination. The criteria adopted for the transition-to-homogeneity follows the literature, it is attained when $\mathcal {D}_2(r)$ reaches the 1 per cent level of the limit value 3 (i.e. where it reaches 2.97) as the scale increases. We obtain RH = 70.33 ± 10.74 Mpc h−1, at the effective redshift zeff = 0.128, for a sample containing $150\, 302$ SDSS blue galaxies with 0.04 < z < 0.20. Additionally, we perform robustness tests by analysing the homogeneity scale in sub-volumes of the original one, obtaining coherent results; we also check for a possible artefact in our procedure examining a homogeneous synthetic data set as a pseudo-data, verifying that such systematic is absent. Because our analyses concentrate in data at low redshifts, z < 0.20, we find interesting to use cosmography to calculate the radial comoving distances; therefore in this subject our analyses do not use fiducial cosmological model. For completeness, we evaluate the difference of the comoving distances estimation using cosmography and fiducial cosmology.

Model-agnostic cosmological constraints from the baryon acoustic oscillation feature in redshift space

ABSTRACT We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, without relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point rLP, and substantially extends it to simultaneously model the multipoles ℓ = 0, 2, 4 of the anisotropic galaxy 2-point correlation function (2pcf). The approach is ‘model-agnostic’: it assumes that the nonlinear growth of structure smears the BAO feature by an approximately Gaussian kernel with a smearing scale σv, but does not assume any fiducial cosmology for describing the shape of the feature itself. Using mock observations for two realistic survey configurations assuming Λ cold dark matter (ΛCDM), combined with Bayesian parameter inference, we show that the linear point rLP and smearing scale σv can be accurately recovered by our method in both existing and upcoming surveys. The precision of the recovery of rLP is always better than $1{{\ \rm per\ cent}}$, while σv can be recovered with $\lesssim 10{{\ \rm per\ cent}}$ uncertainty provided the linear galaxy bias b is separately constrained, e.g. using weak lensing observations. Our method is also sensitive to the linear growth rate f, albeit with larger uncertainties and systematic errors, especially for upcoming surveys such as DESI. We discuss how our model can be modified to improve the recovery of f, such that the resulting constraints on {f, σv, rLP} can potentially be used as a test of cosmological models including and beyond ΛCDM.

Synthetic light-cone catalogues of modern redshift and weak lensing surveys with abacussummit

ABSTRACT The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy, and dark matter. However, the development of high-fidelity theoretical models is a necessary stepping stone. Here, we present public high-resolution weak lensing maps on the light-cone, generated using the N-body simulation suite abacussummit, and accompanying weak lensing mock catalogues, tuned to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument. Available in this release are maps of the cosmic shear, deflection angle, and convergence fields at source redshifts ranging from z = 0.15 to 2.45 as well as cosmic microwave background convergence maps for each of the 25 base-resolution simulations ($L_{\rm box} = 2000\, h^{-1}\, {\rm Mpc}$ and Npart = 69123) as well as for the two huge simulations ($L_{\rm box} = 7500\, h^{-1}\, {\rm Mpc}$ and Npart = 86403) at the fiducial abacussummit cosmology. The pixel resolution of each map is 0.21 arcmin, corresponding to a healpix Nside of 16 384. The sky coverage of the base simulations is an octant until z ≈ 0.8 (decreasing to about 1800 deg2 at z ≈ 2.4), whereas the huge simulations offer full-sky coverage until z ≈ 2.2. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper Suprime-Cam data sets. The mock catalogues are validated against theoretical predictions for various clustering and lensing statistics, such as correlation multipoles, galaxy–shear, and shear–shear, showing excellent agreement. All products can be downloaded via a Globus endpoint (see Data Availability section).

Planting a Lyman alpha forest on AbacusSummit

ABSTRACT The full-shape correlations of the Lyman alpha (Ly α) forest contain a wealth of cosmological information through the Alcock–Paczyński effect. However, these measurements are challenging to model without robustly testing and verifying the theoretical framework used for analysing them. Here, we leverage the accuracy and volume of the N-body simulation suite AbacusSummit to generate high-resolution Ly α skewers and quasi-stellar object (QSO) catalogues. One of the main goals of our mocks is to aid in the full-shape Ly α analysis planned by the Dark Energy Spectroscopic Instrument (DESI) team. We provide optical depth skewers for six of the fiducial cosmology base-resolution simulations ($L_{\rm box} = 2\, h^{-1}\, {\rm Gpc}$, N = 69123) at z = 2.5. We adopt a simple recipe based on the Fluctuating Gunn–Peterson Approximation (FGPA) for constructing these skewers from the matter density in an N-body simulation and calibrate it against the 1D and 3D Ly α power spectra extracted from the hydrodynamical simulation IllustrisTNG (TNG; $L_{\rm box} = 205\, h^{-1}\, {\rm Mpc}$, N = 25003). As an important application, we study the non-linear broadening of the baryon acoustic oscillation (BAO) peak and show the cross-correlation between DESI-like QSOs and our Ly α forest skewers. We find differences on small scales between the Kaiser approximation prediction and our mock measurements of the Ly α × QSO cross-correlation, which would be important to account for in upcoming analyses. The AbacusSummit Ly α forest mocks open up the possibility for improved modelling of cross-correlations between Ly α and cosmic microwave background (CMB) lensing and Ly α and QSOs, and for forecasts of the 3-point Ly α correlation function. Our catalogues and skewers are publicly available on Globus via the National Energy Research Scientific Computing Center (NERSC) (full link under the section ‘Data Availability’).

Cosmological forecasts for future galaxy surveys with the linear point standard ruler: Toward consistent BAO analyses far from a fiducial cosmology

The linear point (LP) standard ruler was identified as the basis of a purely geometric method for exploiting the baryon acoustic oscillations (BAO). The LP exploits the BAO feature imprinted in the galaxy two-point correlation function to measure cosmological distances independent of any specific cosmological model. We forecast the expected precision of future and ongoing spectroscopic galaxy surveys to measure distances leveraging the linear point. We investigate the cosmological implications of our forecasted results. We focus in particular on a relevant working example: the detection of the late-time cosmic acceleration independent of other cosmological probes. Our findings show that, even within the $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ standard cosmological paradigm, estimated distances need to be reliable over a very wide parameter range in order to realize their maximum utility. This is particularly relevant if we aim to properly characterize cosmological tensions. The LP is a promising candidate approach to achieve this reliability. In contrast, widely employed procedures in BAO analysis estimate distances keeping fixed cosmological parameters to fiducial values close to cosmic-microwave-background constraints in flat-$\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$. It is unclear whether they are purely geometric methods. Moreover, they rely on untested extrapolations to explore the parameter space away from those fiducial flat-$\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ values. We recommend that all BAO methodologies be validated across the full range of models and parameters over which their results are quoted, first by means of linear predictions and then $N$-body simulations.

The XXL Survey

Context.X-ray observations of galaxy clusters are impacted by the presence of active galactic nuclei (AGNs) in a manner that is challenging to quantify, leading to biases in the detection and measurement of cluster properties for both astrophysics and cosmological applications.Aims.We detect and characterise clusters contaminated by central AGNs within the XXL survey footprint and provide a systematic assessment of the cosmological impact of such systems in X-ray cluster samples.Methods.We introduce a new automated class for AGN-contaminated (AC) clusters in the XXL source detection pipeline. The majority of these systems are otherwise missed by current X-ray cluster-detection methods. The AC selection is also effective in distinguishing AGN and cool-core presence using supplementary optical and infrared information.Results.We present 33 AC objects, including 25 clusters in the redshift range, 0.14 ≤ z ≤ 1.03, and eight other sources with significantly peaked central profiles based on X-ray observations. Six of these are new confirmed clusters. We computed the missed fraction of the XXL survey, which is defined as the fraction of genuine clusters that are undetected due to their centrally peaked X-ray profiles. We report seven undetected AC clusters abovez > 0.6, in the range where X-ray cluster detection efficiency drops significantly. The missed fraction is estimated to be at the level of 5% for the 50 square-degree XXL area. The impact on cosmological estimates from missed clusters is negligible for XXL, but it produces a tension of ∼3σwith the fiducial cosmology when considering larger survey areas.Conclusions.This work demonstrates the first systematic attempt to quantify the percentage of missed clusters in X-ray surveys as a result of central AGN contamination. Looking towards surveys such as eROSITA and Athena, larger areas and increased sensitivity will significantly enhance cluster detection, and therefore robust methods for characterising AGN contamination will be crucial for precise cluster cosmology, particularly in the redshiftz > 1 regime.

What It Takes to Solve the Hubble Tension through Modifications of Cosmological Recombination.

We construct data-driven solutions to the Hubble tension which are perturbative modifications to the fiducial ΛCDM cosmology, using the Fisher bias formalism. Taking as proof of principle the case of a time-varying electron mass and fine structure constant, and focusing first on Planck CMB data, we demonstrate that a modified recombination can solve the Hubble tension and lower S_{8} to match weak lensing measurements. Once baryonic acoustic oscillation and uncalibrated supernovae data are included, however, it is not possible to fully solve the tension with perturbative modifications to recombination.

Guess the cheese flavour by the size of its holes: a cosmological test using the abundance of popcorn voids

ABSTRACT We present a new definition of cosmic void and a publicly available code with the algorithm that implements it. Underdense regions are defined as free-form objects, called popcorn voids, made from the union of spheres of maximum volume with a given joint integrated underdensity contrast. The method is inspired by the excursion-set theory and consequently no rescaling processing is needed, the removal of overlapping voids and objects with sizes below the shot noise threshold is inherent in the algorithm. The abundance of popcorn voids in the matter field can be fitted using the excursion-set theory provided the relationship between the linear density contrast of the barrier and the threshold used in void identification is modified relative to the spherical evolution model. We also analysed the abundance of voids in biased tracer samples in redshift space. We show how the void abundance can be used to measure the geometric distortions due to the assumed fiducial cosmology, in a test similar to an Alcock–Paczyński test. Using the formalism derived from previous works, we show how to correct the abundance of popcorn voids for redshift-space distortion effects. Using this treatment, in combination with the excursion-set theory, we demonstrate the feasibility of void abundance measurements as cosmological probes. We obtain unbiased estimates of the target parameters, albeit with large degeneracies in the parameter space. Therefore, we conclude that the proposed test in combination with other cosmological probes has potential to improve current cosmological parameter constraints.

CosmoGridV1: a simulated \U0001d5d0CDM theory prediction for map-level cosmological inference

We present CosmoGridV1: a large set of lightcone simulations for map-level cosmological inference with probes of large scale structure. It is designed for cosmological parameter measurement based on Stage-III photometric surveys with non-Gaussian statistics and machine learning. CosmoGridV1 spans the wCDM model by varying Ωm, σ 8, w 0, H 0, n s, Ω b , and assumes three degenerate neutrinos with fixed ∑ mν = 0.06 eV. This space is covered by 2500 grid points on a Sobol sequence. At each grid point, we run 7 simulations with PkdGrav3 and store 69 particle maps at nside = 2048 up to z = 3.5, as well as halo catalog snapshots. The fiducial cosmology has 200 independent simulations, along with their stencil derivatives. An important part of CosmoGridV1 is the benchmark set of 28 simulations, which include larger boxes, higher particle counts, and higher redshift resolution of shells. They allow for testing if new types of analyses are sensitive to choices made in CosmoGridV1. We add baryon feedback effects on the map level, using shell-based baryon correction model. The shells are used to create maps of weak gravitational lensing, intrinsic alignment, and galaxy clustering, using the UFalcon code. The main part of CosmoGridV1 are the raw particle count shells that can be used to create full-sky maps for a given n(z). We also release projected maps for a Stage-III forecast, as well as maps used previously in KiDS-1000 deep learning constraints with CosmoGridV1. The data is available at http://www.cosmogrid.ai/.

COMET: Clustering observables modelled by emulated perturbation theory

ABSTRACT In this paper, we present COMET, a Gaussian process emulator of the galaxy power spectrum multipoles in redshift space. The model predictions are based on one-loop perturbation theory and we consider two alternative descriptions of redshift-space distortions: one that performs a full expansion of the real- to redshift-space mapping, as in recent effective field theory models, and another that preserves the non-perturbative impact of small-scale velocities by means of an effective damping function. The outputs of COMET can be obtained at arbitrary redshifts, for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters ωc, ωb, and ns, as well as the evolution parameters h, As, ΩK, w0, and wa. This flexibility does not impair COMET’s accuracy, since we exploit an exact degeneracy between the evolution parameters that allows us to train the emulator on a significantly reduced parameter space. While the predictions are sped up by two orders of magnitude, validation tests reveal an accuracy of $0.1\, {{\ \rm per\ cent}}$ for the monopole and quadrupole ($0.3\, {{\ \rm per\ cent}}$ for the hexadecapole), or alternatively, better than $0.25\, \sigma$ for all three multipoles in comparison to statistical uncertainties expected for the Euclid survey with a tenfold increase in volume. We show that these differences translate into shifts in mean posterior values that are at most of the same size, meaning that COMET can be used with the same confidence as the exact underlying models. COMET is a publicly available python package that also provides the tree-level bispectrum multipoles and Gaussian covariance matrices.

Accelerating BAO scale fitting using Taylor series

ABSTRACT The Universe is currently undergoing accelerated expansion driven by dark energy. Dark energy’s essential nature remains mysterious: one means of revealing it is by measuring the Universe’s size at different redshifts. This may be done using the baryon acoustic oscillation (BAO) feature, a standard ruler in the galaxy two-point correlation function (2PCF). In order to measure the distance scale, one dilates and contracts a template for the 2PCF in a fiducial cosmology, using a scaling factor α. The standard method for finding the best-fitting α is to compute the likelihood over a grid of roughly 100 values of it. This approach is slow; in this work, we propose a significantly faster way. Our method writes the 2PCF as a polynomial in α by Taylor-expanding it about α = 1, exploiting that we know the fiducial cosmology sufficiently well that α is within a few per cent of unity. The likelihood resulting from this expansion may then be analytically solved for the best-fitting α. Our method is 48–85× faster than a directly comparable approach in which we numerically minimize α, and ∼12 000× faster than the standard iterative method. Our work will be highly enabling for upcoming large-scale structure redshift surveys such as that by Dark Energy Spectroscopic Instrument.

Intermediate redshift calibration of gamma-ray bursts and cosmic constraints in non-flat cosmology

ABSTRACT We propose a new method to calibrate gamma-ray burst (GRB) correlations employing intermediate redshift data sets, instead of limiting to z ≃ 0 catalogues, and applied it to the well-consolidated Amati correlation. This model-independent calibration technique is based on the Bézier polynomial interpolation of the most updated observational Hubble data and baryonic acoustic oscillations (BAO) and alleviates de facto the well-known circularity problem affecting GRB correlations. In doing so, we also investigate the influence of the BAO scales rs, got from Planck results, and $r_{\rm s}^{\rm fid}$, got from the considered fiducial cosmology, by considering the cases $(r_{\rm s}/r_{\rm s}^{\rm fid})=1$ and $(r_{\rm s}/r_{\rm s}^{\rm fid})\ne 1$. We get constraints on the cosmic parameters, using Markov chain–Monte Carlo simulations, first fixing and then leaving free the spatial curvature density parameter Ωk. Compared to previous literature, we obtain tighter constraints on the mass density parameter Ωm. In particular, our findings turn out to be highly more compatible with those got from standard candle indicators within the Lambda cold dark matter paradigm. Finally, we critically re-examine the recent H0 tension in view of our outcomes.

Observational constraints on f(T) gravity from model-independent data

We establish new constraints on f(T) gravity models by using cosmological data. In particular, we investigate the restrictions given by the gas mass fraction measurements of galaxy clusters and transversal BAO data. Both data sets are regarded as weakly dependent on a fiducial cosmology. In addition, we also include a CMB measurement of the temperature power spectrum first peak, along with H(z) values from cosmic chronometers and supernovae data from the Pantheon data set. We also perform a forecast for future constraints on the deviation of f(T) models from the Lambda CDM scenario by following the specifications of the J-PAS and Euclid surveys and find significant improvements on the constraints of the b-parameter, when compared to the results of the statistical analysis.

Combined full shape analysis of BOSS galaxies and eBOSS quasars using an iterative emulator

ABSTRACT Standard full-shape clustering analyses in Fourier space rely on a fixed power spectrum template, defined at the fiducial cosmology used to convert redshifts into distances, and compress the cosmological information into the Alcock–Paczynski parameters and the linear growth rate of structure. In this paper, we propose an analysis method that operates directly in the cosmology parameter space and varies the power spectrum template accordingly at each tested point. Predictions for the power spectrum multipoles from the TNS model are computed at different cosmologies in the framework of $\Lambda \rm {CDM}$. Applied to the final eBOSS QSO and LRG samples together with the low-z DR12 BOSS galaxy sample, our analysis results in a set of constraints on the cosmological parameters Ωcdm, H0, σ8, Ωb, and ns. To reduce the number of computed models, we construct an iterative process to sample the likelihood surface, where each iteration consists of a Gaussian process regression. This method is validated with mocks from N-body simulations. From the combined analysis of the (e)BOSS data, we obtain the following constraints: σ8 = 0.877 ± 0.049 and $\Omega _{\rm m}=0.304^{+0.016}_{-0.010}$ without any external prior. The eBOSS quasar sample alone shows a 3.1σ discrepancy compared to the Planck prediction.

Impact of lensing magnification on the analysis of galaxy clustering in redshift space

We study the impact of lensing magnification on the observed three-dimensional galaxy clustering in redshift space. We used the RayGal suite of N-body simulations, from which we extracted samples of dark matter particles and haloes in the redshift regime of interest for future large redshift surveys. Several magnitude-limited samples were built that reproduce various levels of magnification bias ranging from s = 0 to s = 1.2, where s is the logarithmic slope of the cumulative magnitude number counts, in three redshift intervals within 1 < z < 1.95. We studied the two-point correlation function multipole moments in the different cases in the same way as would be applied to real data, and investigated how well the growth rate of structure parameter could be recovered. In the analysis, we used an hybrid model that combines non-linear redshift-space distortions and linear curved-sky lensing magnification. We find that the growth rate is underestimated when magnification bias is not accounted for in the modelling. This bias becomes non-negligible for z > 1.3 and can reach 10% at z ≃ 1.8, depending on the properties of the target sample. In our data, adding the lensing linear correction allowed us to recover an unbiased estimate of the growth rate in most cases when the correction was small, even when the fiducial cosmology was different from that of the data. For larger corrections (high redshifts, low bias, and high s value), we find that the weak-lensing limit has to be treated with caution as it may no longer be a good approximation. Our results also show the importance of knowing s in advance instead of letting this parameter free with flat priors because in this case, the error bars increase significantly.