Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies

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.

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.

We study the cross-correlation between 212 MgII quasar absorption systems and \~20,000 Luminous Red Galaxies (LRGs) selected from the Sloan Digital Sky Survey Data Release 1 in the redshift range 0.4 =1.0 Angstrom and identifications confirmed by the FeII 2600 or MgI 2852 lines. Over comoving scales 0.05--13 h^-1 Mpc, the MgII--LRG cross-correlation has an amplitude 0.69+/-0.09 times that of the LRG--LRG auto-correlation. Since LRGs have halo-masses greater than 3.5 x 10^12 solar masses for M_R =0.3 Angstrom are associated with ~0.7 L^*_B galaxies.

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.

We present the first scientific results from the luminous red galaxy (LRG) sample of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) combined with the high-redshift galaxies of the previous BOSS sample. We measure the small- and intermediate-scale clustering from a sample of more than 97,000 galaxies in the redshift range . We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of this sample of LRGs is 2.30 ± 0.03, with a satellite fraction of 13% ± 3% and a mean halo mass of . These results are consistent with expectations, demonstrating that these LRGs will be reliable tracers of large-scale structure at . The galaxy bias implies a scatter of luminosity at fixed halo mass, , of 0.19 dex. Using the clustering of massive galaxies from BOSS CMASS, BOSS LOWZ, and SDSS, we find that is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows the investigation of the evolution of massive galaxies over the past ∼7 Gyr.

We model the apparent clustering anisotropy of Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey using subhalos identified in cosmological $N$-body simulations. We first conduct a Markov-chain Monte Carlo analysis on the parameters characterizing subhalos hosting LRGs assuming a specific $\Lambda$CDM cosmology on which we run the simulations. We show that simple models with central and satellite subhalos can explain the observed multipole moments of the power spectrum up to hexadecapole on large scales ($k\lesssim0.3~h\mathrm{Mpc}^{-1}$). A satellite fraction of $20$ to $30$ per cent is favored weakly depending on the detail of the model. The fraction is shown to be robust when we adopt a more refined model based on the halo occupation number from the literature. We then vary cosmological parameters controlling the anisotropy in redshift-space effectively by deforming the simulation box (the Alcock-Paczynski effect) and changing the amplitude of the velocities (the redshift-space distortions). We demonstrate that we can constrain the geometry of the universe, the structure growth rate, and the parameters characterizing LRGs simultaneously. This is a step toward cosmological analysis with realistic bias description beyond empirical bias functions with nuisance parameters.

A photometric redshift sample of luminous red galaxies (LRGs) obtained from the DECam Legacy Survey (DECaLS) is analyzed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in previous work (Sridhar & Song 2019) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analyzing the LRG sample from DECaLS, which will provide the optical imaging for targeting two-thirds of the Dark Energy Spectroscopic Instrument footprint and measure the angular diameter distances at z = 0.69 and z = 0.87 to be and with a fractional error of 4.77% and 6.09%, respectively. We obtain a value of H 0 = 66.58 ± 5.31 km s−1 Mpc−1, which supports the H 0 measured by all other baryon acoustic oscillation results and is consistent with the ΛCDM model.

We usethe forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

We apply a new model for the spherically averaged correlation function at large pair separations to the measurement of the clustering of luminous red galaxies (LRGs) made from the SDSS by Cabre and Gaztanaga(2009). Our model takes into account the form of the BAO peak and the large scale shape of the correlation function. We perform a Monte Carlo Markov chain analysis for different combinations of datasets and for different parameter sets. When used in combination with a compilation of the latest CMB measurements, the LRG clustering and the latest supernovae results give constraints on cosmological parameters which are comparable and in remarkably good agreement, resolving the tension reported in some studies. The best fitting model in the context of a flat, Lambda-CDM cosmology is specified by Omega_m=0.261+-0.013, Omega_b=0.044+-0.001, n_s=0.96+-0.01, H_0=71.6+-1.2 km/s/Mpc and sigma_8=0.80+-0.02. If we allow the time-independent dark energy equation of state parameter to vary, we find results consistent with a cosmological constant at the 5% level using all data sets: w_DE=-0.97+-0.05. The large scale structure measurements by themselves can constrain the dark energy equation of state parameter to w_DE=-1.05+-0.15, independently of CMB or supernovae data. We do not find convincing evidence for an evolving equation of state. We provide a set of "extended distance priors" that contain the most relevant information from the CMB power spectrum and the shape of the LRG correlation function which can be used to constrain dark energy models and spatial curvature. Our model should provide an accurate description of the clustering even in much larger, forthcoming surveys, such as those planned with NASA's JDEM or ESA's Euclid mission.

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.

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 present the 3-point function \xi_3 and Q_3=\xi_3/\xi_2^2 for a spectroscopic sample of luminous red galaxies (LRG) from SDSS DR6 & DR7. We find a strong (S/N$>$6) detection of $Q_3$ on scales of 55-125 Mpc/h, with a well defined peak around 105 Mpc/h in all \xi_2, \xi_3 and Q_3, in excellent agreement with the predicted shape and location of the imprint of the baryon acoustic oscillations (BAO). We use very large simulations to asses and test the significance of our measurement. Models without the BAO peak are ruled out by the $Q_3$ data with 99% confidence. Our measurements show the expected shape for Q_3 as a function of the triangular configuration. This provides a first direct measurement of the non-linear mode coupling coefficients of density and velocity fluctuations which, on these large scales, are independent of cosmic time, the amplitude of fluctuations or cosmological parameters. The location of the BAO peak in the data indicates \Omega_m =0.28 \pm 0.05 and \Omega_B=0.079 \pm 0.025 (h=0.70) after marginalization over spectral index (n_s=0.8-1.2) linear b_1 and quadratic c_2 bias,which are found to be in the range: b_1=1.7-2.2 and c_2=0.75-3.55. The data allows a hierarchical contribution from primordial non-Gaussianities in the range Q_3=0.55-3.35. These constraints are independent and complementary to the ones that can be obtained using the 2-point function, which are presented in a separate paper. This is the first detection of the shape of $Q_3$ on BAO scales, but our errors are shot-noise dominated and the SDSS volume is still relatively small, so there is ample room for future improvement in this type of measurements.

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.

The use of standard rulers, such as the scale of the Baryonic Acoustic oscillations (BAO), has become one of the more powerful techniques employed in cosmology to probe the entity driving the accelerating expansion of the Universe. In this paper, the topology of large scale structure (LSS) is used as one such standard ruler to study this mysterious `dark energy'. By following the redshift evolution of the clustering of luminous red galaxies (LRGs) as measured by their 3D topology (counting structures in the cosmic web), we can chart the expansion rate and extract information about the equation of state of dark energy. Using the technique first introduced in (Park & Kim, 2009), we evaluate the constraints that can be achieved using 3D topology measurements from next-generation LSS surveys such as the Baryonic Oscillation Spectroscopic Survey (BOSS). In conjunction with the information that will be available from the Planck satellite, we find a single topology measurement on 3 different scales is capable of constraining a single dark energy parameter to within 5% and 10% when dynamics are permitted. This offers an alternative use of the data available from redshift surveys and serves as a cross-check for BAO studies.

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.