Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We use subhalo abundance and age distribution matching to create magnitude-limited mock galaxy catalogs at z ∼ 0.43, 0.52, and 0.63 with z-band and 3.4 μm W1-band absolute magnitudes and r − z and r − W1 colors. From these magnitude-limited mocks, we select mock luminous red galaxy (LRG) samples according to the (r − z)-based (optical) and (r − W1)-based (infrared) selection criteria for the LRG sample of the Dark Energy Spectroscopic Instrument (DESI) survey. Our models reproduce the number densities, luminosity functions, color distributions, and projected clustering of the DESI Legacy Surveys that are the basis for DESI LRG target selection. We predict the halo occupation statistics of both optical and IR DESI LRGs at fixed cosmology and assess the differences between the two LRG samples. We find that IR-based SHAM modeling represents the differences between the optical and IR LRG populations better than using the z band and that age distribution matching overpredicts the clustering of LRGs, implying that galaxy color is uncorrelated with halo age in the LRG regime. Both the optical and IR DESI LRG target selections exclude some of the most luminous galaxies that would appear to be LRGs based on their position on the red sequence in optical color–magnitude space. Both selections also yield populations with a nontrivial LRG–halo connection that does not reach unity for the most massive halos. We find that the IR selection achieves greater completeness (≳90%) than the optical selection across all redshift bins studied.

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a five-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range 0.4 < z ≲ 1.0. Here we present the selection of the DESI LRG sample and assess its spectroscopic performance using data from Survey Validation (SV) and the first two months of the Main Survey. The DESI LRG sample, selected using g, r, z, and W1 photometry from the DESI Legacy Imaging Surveys, is highly robust against imaging systematics. The sample has a target density of 605 deg−2 and a comoving number density of 5 × 10−4 h 3 Mpc−3 in 0.4 < z < 0.8; this is a significantly higher density than previous LRG surveys (such as SDSS, BOSS, and eBOSS) while also extending to z ∼ 1. After applying a bright star veto mask developed for the sample, 98.9% of the observed LRG targets yield confident redshifts (with a catastrophic failure rate of 0.2% in the confident redshifts), and only 0.5% of the LRG targets are stellar contamination. The LRG redshift efficiency varies with source brightness and effective exposure time, and we present a simple model that accurately characterizes this dependence. In the appendices, we describe the extended LRG samples observed during SV.

We analyse the large-scale clustering of the Luminous Red Galaxy (LRG) and Quasar (QSO) sample from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). In particular, we constrain the primordial non-Gaussianity (PNG) parameter f NL loc via the large-scale scale-dependent bias in the power spectrum using 1,631,716 LRGs (0.6 < z < 1.1) and 1,189,129 QSOs (0.8 < z < 3.1). This new measurement takes advantage of the enormous statistical power at large scales of DESI DR1 data, surpassing the latest data release (DR16) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). For the first time in this kind of analysis, we use a blinding procedure to mitigate the risk of confirmation bias in our results. We improve the model of the radial integral constraint proposing an innovative technique allowing the correction through the window matrix convolution. We also carefully test the mitigation of the dependence of the target selection on the photometry qualities by incorporating an angular integral constraint contribution to the window function, and validate our methodology with the blinded data. Finally, combining the two samples, we measure f NL loc = -3.6-9.1 +9.0 at 68% confidence, where we assume the universality relation for the LRG sample and a recent merger model for the QSO sample about the response of bias to primordial non-Gaussianity. Adopting the universality relation for the PNG bias in the QSO analysis leads to f NL loc = 3.5-7.4 +10.7 at 68% confidence. Due to restricted selection in the LRG sample, the inclusion of the LRGs allows for 10% improvement. This measurement is the most precise determination of primordial non-Gaussianity using large-scale structure to date, surpassing the latest result from eBOSS by a factor of 2.3.

The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect DESI spectra of approximately 2500 bright galaxies, 3500 luminous red galaxies (LRGs), and 10,000 emission-line galaxies (ELGs) to obtain robust redshift identifications. We then utilize the visually inspected redshift information to characterize the performance of the DESI operation. Based on the visual inspection (VI) catalogs, our results show that the final survey design yields samples of bright galaxies, LRGs, and ELGs with purity greater than 99%. Moreover, we demonstrate that the precision of the redshift measurements is approximately 10 km s−1 for bright galaxies and ELGs and approximately 40 km s−1 for LRGs. The average redshift accuracy is within 10 km s−1 for the three types of galaxies. The VI process also helps improve the quality of the DESI data by identifying spurious spectral features introduced by the pipeline. Finally, we show examples of unexpected real astronomical objects, such as Lyα emitters and strong lensing candidates, identified by VI. These results demonstrate the importance and utility of visually inspecting data from incoming and upcoming surveys, especially during their early operation phases.

We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as ‘galaxy assembly bias’. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15 ${{\ \rm per\ cent}}$ excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures.

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 &amp;lt; z &amp;lt; 0.6) and 0.8 per cent for high (0.8 &amp;lt; z &amp;lt; 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.

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 &amp;lt; r &amp;lt; 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 &amp;lt; z &amp;lt; 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 &amp;lt; z &amp;lt; 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 &amp;lt; z &amp;lt; 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.

Aims. We use the angular cross-correlation between a luminous red galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey data release DR9 and the Planck cosmic microwave background (CMB) lensing maps to constrain the local primordial non-Gaussianity parameter, fNL, using the scale-dependent galaxy bias effect. The galaxy sample covers approximately 40% of the sky, contains galaxies up to redshift z ∼ 1.4, and is calibrated with the LRG spectra that have been observed for DESI Year 1 (Y1). Methods. We apply a nonlinear imaging systematics treatment based on neural networks to remove observational effects that could potentially bias the fNL measurement. Our measurement is performed without blinding, but the full analysis pipeline is tested with simulations including systematics. Results. Using the two-point angular cross-correlation between LRG and CMB lensing only, we find fNL = 39−38+40 at the 68% confidence level, and our result is robust in terms of systematics and cosmological assumptions. If we combine this information with the autocorrelation of LRG, applying a scale cut to limit the impact of systematics, we find fNL = 24−21+20 at the 68% confidence level. Our results motivate the use of CMB lensing cross-correlations to measure fNL with future datasets, given its stability in terms of observational systematics compared to the angular autocorrelation. Furthermore, performing accurate systematics mitigation is crucially important in order to achieve competitive constraints on fNL from CMB lensing cross-correlation in combination with the tracers’ autocorrelation.

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies, luminous red galaxies, and quasars. In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey and the Milky Way Survey. DESI also observes a selection of “secondary” targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe “supporting” DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target-selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data.

We present the first joint analysis of the power spectrum and bispectrum using the Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI), focusing on Luminous Red Galaxies (LRGs) and quasars (QSOs) across a redshift range of 0.4 ≤ z ≤ 2.1. By combining the two- and three-point statistics, we are able to partially break the degeneracy between the logarithmic growth rate, f(z), and the amplitude of dark matter fluctuations, σ s8(z), which cannot be measured separately in analyses that only involve the power spectrum. In comparison with the (fiducial) Planck ΛCDM cosmology we obtain f/f fid = {0.888-0.089 +0.186,0.977-0.220 +0.182,1.030-0.085 +0.368}, σ s8/σ fid s8 = {1.224-0.133 +0.091,1.071-0.163 +0.278,1.00 0-0.223 +0.088} respectively for the three LRG redshift bins, corresponding to a cumulative 10.1% constraint on f, and of 8.4% on σ s8, including the systematic error budget. Additionally, we obtain constraints for the ShapeFit compressed parameters describing the isotropic scaling parameter, α iso(z), the Alcock-Paczyński parameter, α AP(z), the combined growth of structure parameter fσ s8(z), and the combined shape parameter m(z)+n(z). Their cumulative constraints from our joint power spectrum-bispectrum analysis are respectively σ α iso = 0.9% (9% improvement with respect to our power spectrum-only analysis); σ α AP = 2.3% (no improvement with respect to power spectrum-only analysis, which is expected given that the bispectrum monopole has no significant anisotropic signal); σ fσ s8 = 5.1% (9% improvement); σ m+n = 2.3% (11% improvement). These results are fully consistent with the main DESI power spectrum analysis, demonstrating the robustness of the DESI cosmological constraints, and compatible with Planck ΛCDM cosmology.

In the current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping the dark matter distribution at z ∼ 1. We measure the auto and cross correlation functions of ELGs and LRGs at 0.8 < z ≤ 1.0 from the DESI One-Percent survey. Following Gao et al., we construct the galaxy–halo connections for ELGs and LRGs simultaneously. With the stellar–halo mass relation for the whole galaxy population (i.e., normal galaxies), LRGs can be selected directly by stellar mass, while ELGs can also be selected randomly based on the observed number density of each stellar mass, once the probability P sat of a satellite galaxy becoming an ELG is determined. We demonstrate that the observed small scale clustering prefers a halo mass-dependent P sat model rather than a constant. With this model, we can well reproduce the auto correlations of LRGs and the cross correlations between LRGs and ELGs at r p > 0.1 Mpc h −1. We can also reproduce the auto correlations of ELGs at r p > 0.3 Mpc h −1 (s > 1 Mpc h −1) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduces the observed auto correlations of ELGs in the whole range of 0.8 < z ≤ 1.6, which enables us to generate a lightcone ELG mock for DESI. With the above model, we further derive halo occupation distributions for ELGs, which can be used to produce ELG mocks in coarse simulations without resolving subhalos.

We present a high-significance cross-correlation of CMB lensing maps from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) with luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI) Legacy Survey spectroscopically calibrated by DESI. We detect this cross-correlation at a significance of 38σ; combining our measurement with the Planck Public Release 4 (PR4) lensing map, we detect the cross-correlation at 50σ. Fitting this jointly with the galaxy auto-correlation power spectrum to break the galaxy bias degeneracy with σ 8, we perform a tomographic analysis in four LRG redshift bins spanning 0.4 ≤ z ≤ 1.0 to constrain the amplitude of matter density fluctuations through the parameter combination S 8 × = σ 8 (Ω m / 0.3)0.4. Prior to unblinding, we confirm with extragalactic simulations that foreground biases are negligible and carry out a comprehensive suite of null and consistency tests. Using a hybrid effective field theory (HEFT) model that allows scales as small as k max = 0.6 h/ Mpc, we obtain a 3.3% constraint on S 8 × = σ 8 (Ω m / 0.3)0.4 = 0.792+0.024 -0.028 from ACT data, as well as constraints on S 8 ×(z) that probe structure formation over cosmic time.Our result is consistent with the early-universe extrapolation from primary CMB anisotropies measured by Planck PR4 within 1.2σ. Jointly fitting ACT and Planck lensing cross-correlations we obtain a 2.7% constraint of S 8 × = 0.776+0.019 -0.021, which is consistent with the Planck early-universe extrapolation within 2.1σ, with the lowest redshift bin showing the largest difference in mean. The latter may motivate further CMB lensing tomography analyses at z < 0.6 to assess the impact of potential systematics or the consistency of the ΛCDM model over cosmic time.

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.

Using the Luminous Red Galaxy (LRG) sample of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we investigate the environmental dependence of stellar mass, star formation rate (SFR) and specific star formation rate (SSFR) of LRGs. It is found that stellar mass of LRGs nearly is independent of local environments, and that the environmental dependence of SFR and SSFR in the LRG sample is much weaker than the one in the Main galaxy sample. One possible explanation is that galaxy color and morphology are a pair of galaxy properties most predictive of local environments, while LRGs are a group of galaxies that are likely to be luminous, red and of early types (© 2012 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, Vceil, vsmear}-SHAM and {σ, Vceil, fsat}-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 &amp;lt; z &amp;lt; 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. Vceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. vsmear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction fsat 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 Vceil 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 vsmear 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(〈Mvir〉) 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.