Search

Galaxy–galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) data set promise to yield unprecedented combined-probe tests of cosmology and the galaxy–halo connection. In such analyses, it is essential to identify and characterize all relevant statistical and systematic errors. We forecast the covariances of DESI Y1 GGL + clustering measurements and the systematic bias due to redshift evolution in the lens samples. Focusing on the projected clustering and GGL correlations, we compute a Gaussian analytical covariance, using a suite of N -body and lognormal simulations to characterize the effect of the survey footprint. Using the DESI one percent survey data, we measure the evolution of galaxy bias parameters for the DESI luminous red galaxy (LRG) and bright galaxy survey (BGS) samples. We find mild evolution in the LRGs in |$0.4 < z < 0.8$|⁠ , subdominant to the expected statistical errors. For BGS, we find less evolution for brighter absolute magnitude cuts, at the cost of reduced sample size. We find that for a redshift bin width |$\Delta z = 0.1$|⁠ , evolution effects on DESI Y1 GGL is negligible across all scales, all fiducial selection cuts, all fiducial redshift bins. Galaxy clustering is more sensitive to evolution due to the bias squared scaling. Nevertheless the redshift evolution effect is insignificant for clustering above the 1-halo scale of |$0.1h^{-1}$|  Mpc. For studies that wish to reliably access smaller scales, additional treatment of redshift evolution is likely needed. This study serves as a reference for GGL and clustering studies using the DESI Y1 sample. [ABSTRACT FROM AUTHOR]

We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter f NL. Our sample comprises over 12 million targets, covering 14 000 deg2 of the sky, with redshifts in the range 0.2 < z < 1.35. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without f NL and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find |$f_{\rm NL} = 34^{+24(+50)}_{-44(-73)}$| at 68 per cent (95 per cent) confidence. We apply a series of robustness tests (e.g. cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power spectrum and degrades the f NL constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid overcorrection, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of f NL with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the f NL uncertainty. [ABSTRACT FROM AUTHOR]

We present the first comprehensive halo occupation distribution (HOD) analysis of the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey luminous red galaxy (LRG) and Quasi Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32  h −1 Mpc in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of |$f_\mathrm{sat} = 11\pm 1~{y{\ \mathrm{per\,cent}}}$|⁠ , a mean halo mass of |$\log _{10}\overline{M}_h/M_\odot =13.40^{+0.02}_{-0.02}$|⁠ , and a linear bias of |$b_\mathrm{lin} = 1.93_{-0.04}^{+0.06}$|⁠. For LRGs in 0.6 < z < 0.8, we find |$f_\mathrm{sat}=14\pm 1~{{\ \mathrm{per\,cent}}}$|⁠ , |$\log _{10}\overline{M}_h/M_\odot =13.24^{+0.02}_{-0.02}$|⁠ , and |$b_\mathrm{lin}=2.08_{-0.03}^{+0.03}$|⁠. For QSOs, we infer |$f_\mathrm{sat}=3^{+8}_{-2}\mathrm{per\,cent}$|⁠ , |$\log _{10}\overline{M}_h/M_\odot = 12.65^{+0.09}_{-0.04}$|⁠ , and |$b_\mathrm{lin} = 2.63_{-0.26}^{+0.37}$| in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high fidelity galaxy mocks, forming the basis of systematic tests for DESI Y1 cosmological analyses. We also study the redshift-evolution of the DESI LRG sample from z  = 0.4 up to z  = 1.1, revealling significant and interesting trends in mean halo mass, linear bias, and satellite fraction. [ABSTRACT FROM AUTHOR]

We estimate the redshift-dependent, anisotropic clustering signal in the Dark Energy Spectroscopic Instrument (DESI) Year 1 Survey created by tidal alignments of Luminous Red Galaxies (LRGs) and a selection-induced galaxy orientation bias. To this end, we measured the correlation between LRG shapes and the tidal field with DESI's Year 1 redshifts, as traced by LRGs and Emission-Line Galaxies. We also estimate the galaxy orientation bias of LRGs caused by DESI's aperture-based selection, and find it to increase by a factor of seven between redshifts 0.4−1.1 due to redder, fainter galaxies falling closer to DESI's imaging selection cuts. These effects combine to dampen measurements of the quadrupole of the correlation function (ξ2) caused by structure growth on scales of 10–80 h −1 Mpc by about 0.15 per cent for low redshifts (0.4 < z < 0.6) and 0.8 per cent for high (0.8 < z < 1.1), a significant fraction of DESI's error budget. We provide estimates of the ξ2 signal created by intrinsic alignments that can be used to correct this effect, which is necessary to meet DESI's forecasted precision on measuring the growth rate of structure. While imaging quality varies across DESI's footprint, we find no significant difference in this effect between imaging regions in the Legacy Imaging Survey. [ABSTRACT FROM AUTHOR]

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N -body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. |$k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$| or |$s\gt 10\, \mathrm{Mpc}\, h^{-1}$|⁠. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10  per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N -body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies. [ABSTRACT FROM AUTHOR]

We present fitcov an approach for accurate estimation of the covariance of two-point correlation functions that requires fewer mocks than the standard mock-based covariance. This can be achieved by dividing a set of mocks into jackknife regions and fitting the correction term first introduced in Mohammad & Percival (2022), such that the mean of the jackknife covariances corresponds to the one from the mocks. This extends the model beyond the shot-noise limited regime, allowing it to be used for denser samples of galaxies. We test the performance of our fitted jackknife approach, both in terms of accuracy and precision, using lognormal mocks with varying densities and approximate EZmocks mimicking the Dark Energy Spectroscopic Instrument LRG and ELG samples in the redshift range of z  = [0.8, 1.1]. We find that the Mohammad–Percival correction produces a bias in the two-point correlation function covariance matrix that grows with number density and that our fitted jackknife approach does not. We also study the effect of the covariance on the uncertainty of cosmological parameters by performing a full-shape analysis. We demonstrate that our fitted jackknife approach based on 25 mocks can recover unbiased and as precise cosmological parameters as the ones obtained from a covariance matrix based on 1000 or 1500 mocks, while the Mohammad–Percival correction produces uncertainties that are twice as large. The number of mocks required to obtain an accurate estimation of the covariance for the two-point correlation function is therefore reduced by a factor of 40–60. The fitcov code that accompanies this paper is available at this GitHub repository. [ABSTRACT FROM AUTHOR]

We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, V ceil, v smear}-SHAM and {σ, V ceil, f sat}-SHAM. They are designed to reproduce the clustering on 5–30  |$\, {\, h^{-1}\, {\rm Mpc}}$| of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasi-stellar objects (QSOs) at 0.4 < z < 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. V ceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. v smear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction f sat is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy–halo mass relation than LRGs reflected in their weak constraints on σ. LRGs, QSOs, and ELGs show increasing V ceil values, corresponding to the massive galaxy incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. For LRGs, a Gaussian v smear presents a better profile for subsamples at redshift bins than a Lorentzian profile used for other tracers. The impact of the statistical redshift uncertainty on ELG clustering is negligible. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, lower than previous estimations for ELGs. The mean halo mass log10(〈 M vir〉) in |${{\, h^{-1}\, \mbox{M}_\odot }}{}$| for LRGs, ELGs, and QSOs are 13.16 ± 0.01, 11.90 ± 0.06, and 12.66 ± 0.45, respectively. Our generalized SHAM algorithms facilitate the production of multitracer galaxy mocks for cosmological tests. [ABSTRACT FROM AUTHOR]

We present the one-dimensional Ly α forest power spectrum measurement using the first data provided by the Dark Energy Spectroscopic Instrument (DESI). The data sample comprises 26 330 quasar spectra, at redshift z > 2.1, contained in the DESI Early Data Release and the first 2 months of the main survey. We employ a Fast Fourier Transform (FFT) estimator and compare the resulting power spectrum to an alternative likelihood-based method in a companion paper. We investigate methodological and instrumental contaminants associated with the new DESI instrument, applying techniques similar to previous Sloan Digital Sky Survey (SDSS) measurements. We use synthetic data based on lognormal approximation to validate and correct our measurement. We compare our resulting power spectrum with previous SDSS and high-resolution measurements. With relatively small number statistics, we successfully perform the FFT measurement, which is already competitive in terms of the scale range. At the end of the DESI survey, we expect a five times larger Ly α forest sample than SDSS, providing an unprecedented precise one-dimensional power spectrum measurement. [ABSTRACT FROM AUTHOR]

We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel'dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC. [ABSTRACT FROM AUTHOR]

We measure the tidal alignment of the major axes of luminous red galaxies (LRGs) from the Legacy Imaging Survey and use it to infer the artificial redshift-space distortion signature that will arise from an orientation-dependent, surface-brightness selection in the Dark Energy Spectroscopic Instrument (DESI) survey. Using photometric redshifts to downweight the shape–density correlations due to weak lensing, we measure the intrinsic tidal alignment of LRGs. Separately, we estimate the net polarization of LRG orientations from DESI's fibre-magnitude target selection to be of order 10−2 along the line of sight. Using these measurements and a linear tidal model, we forecast a 0.5 per cent fractional decrease on the quadrupole of the two-point correlation function for projected separations of 40–80 h −1 Mpc. We also use a halo catalogue from the Abacussummit cosmological simulation suite to reproduce this false quadrupole. [ABSTRACT FROM AUTHOR]

We investigate using three-point statistics in constraining the galaxy–halo connection. We show that for some galaxy samples, the constraints on the halo occupation distribution parameters are dominated by the three-point function signal (over its two-point counterpart). We demonstrate this on mock catalogues corresponding to the Luminous red galaxies (LRGs), Emission-line galaxies (ELGs), and quasars (QSOs) targeted by the Dark Energy Spectroscopic Instrument (DESI) Survey. The projected three-point function for triangle sides less up to 20  h −1 Mpc measured from a cubic Gpc of data can constrain the characteristic minimum mass of the LRGs with a preci sion of 0.46 per cent. For comparison, similar constraints from the projected two-point function are 1.55 per cent. The improvements for the ELGs and QSOs targets are more modest. In the case of the QSOs, it is caused by the high shot-noise of the sample, and in the case of the ELGs, it is caused by the range of halo masses of the host haloes. The most time-consuming part of our pipeline is the measurement of the three-point functions. We adopt a tabulation method, proposed in earlier works for the two-point function, to significantly reduce the required compute time for the three-point analysis. [ABSTRACT FROM AUTHOR]

Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N -body simulation codes, abacus , gadget-2 , and swift , to investigate the regimes in which their results agree. We run N -body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109  h −1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for |$s\ \gt\ 20\ h^{-1}$|  Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109  h −1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than |$20\ h^{-1}$|  Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10  h  Mpc−1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi- N -body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy–halo relation will be presented in a follow-up study. [ABSTRACT FROM AUTHOR]

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches |$\sim 20\, h^{-3}\, \mathrm{Gpc}^{3}$|⁠. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with |$200\, h^{-3}\, \mathrm{Gpc}^{3}$| (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi- N -body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. |$k \lt 0.3\, h\, \mathrm{Mpc}^{-1}$| for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data. [ABSTRACT FROM AUTHOR]

We characterize the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one and we also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (r mag ≤ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength r 0 and slope γ by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterized and validated with the spectroscopic data of Survey Validation. [ABSTRACT FROM AUTHOR]

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging data set, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 deg2. We have developed new photometric redshift (photo- z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalogue, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N -body catalogues, and present a new technique for finding best-fitting parameters in this situation. Accompanying this paper, we have released the Photometric Redshifts for the Legacy Surveys catalogue of photo- z 's obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 data set. This catalogue provides accurate photometric redshifts for objects with z < 21 over more than 16 000 deg2 of sky. [ABSTRACT FROM AUTHOR]

We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG), and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits |$r_0 = 7.78 \pm 0.26\, h^{-1}$| Mpc, γ  = 1.98 ± 0.02 for LRGs and |$r_0 = 5.45 \pm 0.1\, h^{-1}$| Mpc, γ  = 1.54 ± 0.01 for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale-dependent biases in the weakly non-linear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer dN/dz and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments. [ABSTRACT FROM AUTHOR]

The Emission Line Galaxy survey made by the Dark Energy Spectroscopic Instrument (DESI) survey will be created from five passes of the instrument on the sky. On each pass, the constrained mobility of the ends of the fibres in the DESI focal plane means that the angular distribution of targets that can be observed is limited. Thus, the clustering of samples constructed using a limited number of passes will be strongly affected by missing targets. In two recent papers, we showed how the effect of missing galaxies can be corrected when calculating the correlation function using a weighting scheme for pairs. Using mock galaxy catalogues we now show that this method provides an unbiased estimator of the true correlation function for the DESI survey after any number of passes. We use multiple mocks to determine the expected errors given one to four passes, compared to an idealized survey observing an equivalent number of randomly selected targets. On BAO scales, we find that the error is a factor of 2 worse after one pass, but that after three or more passes, the errors are very similar. Thus, we find that the fibre assignment strategy enforced by the design of DESI will not affect the cosmological measurements to be made by the survey, and can be removed as a potential risk for this experiment. [ABSTRACT FROM AUTHOR]