Search

The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1., Comment: 29 pages, 17 figures, submitted for publication in The Open Journal of Astrophysics, comments welcome

We use the angular cross-correlation between a Luminous Red Galaxy (LRG) sample from the DR9 DESI Legacy Survey and the $Planck$ PR4 CMB lensing maps to constrain the local primordial non-Gaussianity parameter $f_{\rm NL}$ using the scale-dependent galaxy bias effect. The galaxy sample covers $\sim$ 40% of the sky and contains galaxies up to $z \sim 1.4$, and is calibrated with the LRG spectra that have been observed for the DESI Survey Validation. We apply a nonlinear imaging systematics treatment based on neural networks to remove observational effects that could potentially bias the $f_{\rm NL}$ measurement. Our measurement is performed without blinding, but the full analysis pipeline is tested with simulations including systematics. Using the two-point angular cross-correlation between LRG and CMB lensing only ($C_\ell^{\kappa G}$) we find $f_{\rm NL} = 39_{-38}^{+40}$ at 68% confidence level, and our result is robust in terms of systematics and cosmology assumptions. If we combine this information with the autocorrelation of LRG ($C_\ell^{GG}$) applying a $\ell_{\rm min}$ scale cut to limit the impact of systematics, we find $f_{\rm NL} = 24_{-21}^{+20}$ at 68% confidence level. Our results motivate the use of CMB lensing cross-correlations for measuring $f_{\rm NL}$ with future datasets given its stability in terms of observational systematics compared to the angular auto-correlation., Comment: 13 pages, 11 figures, to be submitted to A&A

Super sample covariance (SSC) is important when estimating covariance matrices using a set of mock catalogues for galaxy surveys. If the underlying cosmological simulations do not include the variation in background parameters appropriate for the simulation sizes, then the scatter between mocks will be missing the SSC component. The coupling between large and small modes due to non-linear structure growth makes this pernicious on small scales. We compare different methods for generating ensembles of mocks with SSC built in to the covariance, and contrast against methods where the SSC component is computed and added to the covariance separately. We find that several perturbative expansions, developed to derive background fluctuations, give similar results. We then consider scaling covariance matrices calculated for simulations of different volumes to improve the accuracy of covariance matrix estimation for a given amount of computational time. On large scales, we find that the primary limitation is from the discrete number of modes contributing to the measured power spectrum, and we propose a new method for correcting this effect. Correct implementation of SSC and the effect of discrete mode numbers allows covariance matrices created from mocks to be scaled between volumes, potentially leading to a significant saving on computational resources when producing covariance matrices. We argue that a sub-percent match is difficult to achieve because of the effects of modes on scales between the box sizes, which cannot be easily included. Even so, a 3% match is achievable on scales of interest for current surveys scaling the simulation volume by 512x, costing a small fraction of the computational time of running full-sized simulations. This is comparable to the agreement between analytic and mock-based covariance estimates to be used with DESI Y1 results., Comment: 26 pages, 9 figures. For submission to JCAP

The physical size of the sound horizon at recombination is a powerful source of information for early-time measurements of the Hubble constant $H_0$, and many proposed solutions to the Hubble tension therefore involve modifications to this scale. In light of this, there has been growing interest in measuring $H_0$ independently of the sound horizon. We present the first such measurement to use data from the Dark Energy Spectroscopic Instrument (DESI), jointly analyzing the full-shape galaxy power spectra of DESI luminous red galaxies, emission line galaxies, quasars, and the bright galaxy sample, in a total of six redshift bins. Information from the sound horizon scale is removed from our constraints via a rescaling procedure at the power spectrum level, with our sound horizon-marginalized measurement being driven instead primarily by the matter-radiation equality scale. This measurement is then combined with additional sound horizon-free information from Planck+ACT CMB lensing, uncalibrated type Ia supernovae, and the DESI Lyman-$\alpha$ forest. We agnostically combine with the DESY5, Pantheon+, and Union3 supernova datasets, with our tightest respective constraints being $H_0=66.7^{+1.7}_{-1.9},~67.9^{+1.9}_{-2.1},$ and $67.8^{+2.0}_{-2.2}$ km s-1 Mpc-1. This corresponds to a sub-3% sound horizon-free constraint of the Hubble constant, and is the most precise measurement of its kind to date. Even without including information from the sound horizon, our measurement is still in 2.2-3.0$\sigma$ tension with SH0ES. Additionally, the consistency between our result and other measurements that do rely on the sound horizon scale provides no evidence for new early-Universe physics (e.g. early dark energy). Future DESI data releases will allow unprecedented measurements of $H_0$ and place strong constraints on models that use beyond-$\Lambda$CDM physics to ameliorate the Hubble tension., Comment: 17+7 pages; 5 figures. Submitted to JCAP

Emission Line Galaxies (ELGs) are one of the main tracers that the Dark Energy Spectroscopic Instrument (DESI) uses to probe the universe. However, they are afflicted by strong spurious correlations between target density and observing conditions known as imaging systematics. We present the imaging systematics mitigation applied to the DESI Data Release 1 (DR1) large-scale structure catalogs used in the DESI 2024 cosmological analyses. We also explore extensions of the fiducial treatment. This includes a combined approach, through forward image simulations in conjunction with neural network-based regression, to obtain an angular selection function that mitigates the imaging systematics observed in the DESI DR1 ELGs target density. We further derive a line-of-sight selection function from the forward model that removes the strong redshift dependence between imaging systematics and low redshift ELGs. Combining both angular and redshift-dependent systematics, we construct a 3D selection function and assess the impact of all selection functions on clustering statistics. We quantify differences between these extended treatments and the fiducial treatment in terms of the measured 2-point statistics. We find that the results are generally consistent with the fiducial treatment and conclude that the differences are far less than the imaging systematics uncertainty included in DESI 2024 full-shape measurements. We extend our investigation to the ELGs at $0.6

We present cosmological results from the measurement of clustering of galaxy, quasar and Lyman-$\alpha$ forest tracers from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). We adopt the full-shape (FS) modeling of the power spectrum, including the effects of redshift-space distortions, in an analysis which has been validated in a series of supporting papers. In the flat $\Lambda$CDM cosmological model, DESI (FS+BAO), combined with a baryon density prior from Big Bang Nucleosynthesis and a weak prior on the scalar spectral index, determines matter density to $\Omega_\mathrm{m}=0.2962\pm 0.0095$, and the amplitude of mass fluctuations to $\sigma_8=0.842\pm 0.034$. The addition of the cosmic microwave background (CMB) data tightens these constraints to $\Omega_\mathrm{m}=0.3056\pm 0.0049$ and $\sigma_8=0.8121\pm 0.0053$, while further addition of the the joint clustering and lensing analysis from the Dark Energy Survey Year-3 (DESY3) data leads to a 0.4% determination of the Hubble constant, $H_0 = (68.40\pm 0.27)\,{\rm km\,s^{-1}\,Mpc^{-1}}$. In models with a time-varying dark energy equation of state, combinations of DESI (FS+BAO) with CMB and type Ia supernovae continue to show the preference, previously found in the DESI DR1 BAO analysis, for $w_0>-1$ and $w_a<0$ with similar levels of significance. DESI data, in combination with the CMB, impose the upper limits on the sum of the neutrino masses of $\sum m_\nu < 0.071\,{\rm eV}$ at 95% confidence. DESI data alone measure the modified-gravity parameter that controls the clustering of massive particles, $\mu_0=0.11^{+0.45}_{-0.54}$, while the combination of DESI with the CMB and the clustering and lensing analysis from DESY3 constrains both modified-gravity parameters, giving $\mu_0 = 0.04\pm 0.22$ and $\Sigma_0 = 0.044\pm 0.047$, in agreement with general relativity. [Abridged.], Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers/). 55 pages, 10 figures

We present the samples of galaxies and quasars used for DESI 2024 cosmological analyses, drawn from the DESI Data Release 1 (DR1). We describe the construction of large-scale structure (LSS) catalogs from these samples, which include matched sets of synthetic reference `randoms' and weights that account for variations in the observed density of the samples due to experimental design and varying instrument performance. We detail how we correct for variations in observational completeness, the input `target' densities due to imaging systematics, and the ability to confidently measure redshifts from DESI spectra. We then summarize how remaining uncertainties in the corrections can be translated to systematic uncertainties for particular analyses. We describe the weights added to maximize the signal-to-noise of DESI DR1 2-point clustering measurements. We detail measurement pipelines applied to the LSS catalogs that obtain 2-point clustering measurements in configuration and Fourier space. The resulting 2-point measurements depend on window functions and normalization constraints particular to each sample, and we present the corrections required to match models to the data. We compare the configuration- and Fourier-space 2-point clustering of the data samples to that recovered from simulations of DESI DR1 and find they are, generally, in statistical agreement to within 2\% in the inferred real-space over-density field. The LSS catalogs, 2-point measurements, and their covariance matrices will be released publicly with DESI DR1., Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers/)

We present the measurements and cosmological implications of the galaxy two-point clustering using over 4.7 million unique galaxy and quasar redshifts in the range $0.1

This study presents the modeling of the gravitational wave (GW) bias parameter by bridging a connection between simulated GW sources and galaxies in low redshift galaxy surveys 2MPZ and WISExSCOS (WISC). We study this connection by creating a mock GW catalog, populating galaxy surveys with binary black holes (BBHs) for different scenarios of the GW occupation fraction (or selection function) as a function of the galaxy stellar mass. We probe the observable consequences of this connection by exploring the spatial clustering of the GW sources in terms of the GW bias parameter. We consider a phenomenological broken power law model for the selection function, with a potential turnover $M_{K}$ at high stellar mass ($10^{11}$ $M_{\odot}$ in the fiducial model) where the star formation efficiency begins to drop. We vary the parameters of the selection function and find that generically the GW bias increases as $M_{K}$ increases (and gets suppressed as $M_{K}$ decreases). The change in the GW bias parameter shows a maximum change of about $30\%$ for different scenarios explored in this work in comparison to the galaxy bias. Future measurements of the GW bias can help constrain $M_{K}$ and the slopes of the selection function and thus offer insights into the underlying astrophysical processes., Comment: 29 pages (+15 pages in appendices), 15 figures (+9 figures in appendices), to be submitted to JCAP

Light cone selection effects on cosmic observables must be precisely accounted for in the next generation of surveys, including the Dark Energy Spectroscopic Instrument (DESI) survey. This will allow us to correctly model the data and extract subtle shifts from general-relativistic effects. We examine the effects of peculiar velocities on color selection in spectroscopic galaxy surveys, with a focus on their implications for the galaxy clustering dipole $P_1(k)$. Using DESI Emission Line Galaxy (ELG) targets, we show that peculiar velocities can shift spectral emission features into or out of filter bands, modifying galaxy colors and thereby changing galaxy selection. This phenomenon mimics the effect of evolution bias, and we refer to it as the Doppler bias, $b_D$. The Doppler bias is of comparable size to the evolution bias at $0.8 < z < 1$, where it is largest. This enhances the ELG-LRG (Luminous Red Galaxy) cross-correlation dipole by 15-30%. This could be detectable at the $\sim$4$\sigma$ level for the full DESI survey. Additionally, we found that our $b_D$ estimate is impacted by the incompleteness of the parent ELG sample. Therefore, this work highlights the essential need for careful consideration of spectral-dependent biases caused by peculiar velocities during the selection phase of galaxy surveys, to enable unbiased analyses., Comment: 31 pages, 10 figures

We present a joint analysis of the CMB lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrised by the best constrained parameter combination $S_8^{\rm 3x2pt}\equiv\sigma_8 (\Omega_m/0.3)^{0.4}=0.815\pm0.012$. The commonly used $S_8\equiv\sigma_8 (\Omega_m/0.3)^{0.5}$ parameter is constrained to $S_8=0.816\pm0.015$. In combination with baryon acoustic oscillation (BAO) measurements we find $\sigma_8=0.815\pm 0.012$. We also present sound-horizon-independent estimates of the present day Hubble rate of $H_0=66.4^{+3.2}_{-3.7} \,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ from our large scale structure data alone and $H_0=64.3^{+2.1}_{-2.4}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and auto-correlations, we derive a 3.3% constraint on the integrated matter density fluctuations above $z=2.4$, one of the tightest constraints in this redshift range and fully consistent with a $\Lambda$CDM model fit to the primary CMB from Planck. Combining with primary CMB observations and using the extended low redshift coverage of these combined data sets we derive constraints on a variety of extensions to the $\Lambda$CDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat $\Lambda$CDM $\sum m_\nu<0.12$ eV at 95% confidence using the LSS data, BAO measurements from SDSS and primary CMB observations., Comment: 30 pages, 15 figures, to be submitted to PRD, comments welcome

Measuring the growth of structure is a powerful probe for studying the dark sector, especially in light of the $\sigma_8$ tension between primary CMB anisotropy and low-redshift surveys. This paper provides a new measurement of the amplitude of the matter power spectrum, $\sigma_8$, using galaxy-galaxy and galaxy-CMB lensing power spectra of Dark Energy Spectroscopic Instrument Legacy Imaging Surveys Emission-Line Galaxies and the $\textit{Planck}$ 2018 CMB lensing map. We create an ELG catalog composed of $27$ million galaxies and with a purity of $85\%$, covering a redshift range $0 < z < 3$, with $z_{\rm mean} = 1.09$. We implement several novel systematic corrections, such as jointly modeling the contribution of imaging systematics and photometric redshift uncertainties to the covariance matrix. We also study the impacts of various dust maps on cosmological parameter inference. We measure the cross-power spectra over $f_{\rm sky} = 0.25$ with a signal-to-background ratio of up to $ 30\sigma$. We find that the choice of dust maps to account for imaging systematics in estimating the ELG overdensity field has a significant impact on the final estimated values of $\sigma_8$ and $\Omega_{\rm M}$, with far-infrared emission-based dust maps preferring $\sigma_8$ to be as low as $0.702 \pm 0.030$, and stellar-reddening-based dust maps preferring as high as $0.719 \pm 0.030$. The highest preferred value is at $\sim 3 \sigma$ tension with the $\textit{Planck}$ primary anisotropy results. These findings indicate a need for tomographic analyses at high redshifts and joint modeling of systematics., Comment: 50 pages, 24 figures (figure data can be obtained at https://zenodo.org/records/13381499)

Recent studies have found evidence for parity violation in the BOSS spectroscopic galaxy survey, with statistical significance as high as $7\sigma$. These analyses assess the significance of the parity-odd four-point correlation function (4PCF) with a statistic called $\chi^2$. This statistic is biased if the parity-even eight-point correlation function (8PCF) of the data differs from the mock catalogs. We construct new statistics $\chi^2_\times$, $\chi^2_{\mathrm{null}}$ that separate the parity violation signal from the 8PCF bias term, allowing them to be jointly constrained. Applying these statistics to BOSS, we find that the parity violation signal ranges from $0$ to $2.5\sigma$ depending on analysis choices, whereas the 8PCF bias term is $\sim 6\sigma$. We conclude that there is no compelling evidence for parity violation in BOSS. Our new statistics can be used to search for parity violation in future surveys, such as DESI, without 8PCF biases., Comment: 46 pages, 11 figures; comments welcome!

The large scale structure catalogs within DESI Data Release 1 (DR1) use nearly 6 million galaxies and quasars as tracers of the large-scale structure of the universe to measure the expansion history with baryon acoustic oscillations and the growth of structure with redshift-space distortions. In order to take advantage of DESI's unprecedented statistical power, we must ensure that the galaxy clustering measurements are unaffected by non-cosmological density fluctuations. One source of spurious fluctuations comes from variation in galaxy density with spectroscopic observing conditions, lowering the redshift efficiency (and thus galaxy density) in certain areas of the sky. We measure the uniformity of the redshift success rate for DESI luminous red galaxies (LRG), bright galaxies (BGS) and quasars (QSO), complementing the detailed discussion of emission line galaxy (ELG) systematics in a companion paper (Yu et al., 2024). We find small but significant fluctuations of up to 3% in redshift success rate with the effective spectroscopic signal-to-noise, and create and describe weights that remove these fluctuations. We also describe the process to identify and remove data from certain poorly performing fibers from DESI DR1, and measure the stability of the redshift success rate with time. Finally, we find small but significant correlations of redshift success rate with position on the focal plane, survey speed, and number of exposures required, and show the impact of weights correcting these trends on the power spectrum multipoles and on cosmological parameters from BAO and RSD fits. These corrections change the best-fit parameters by $<15\%$ of their statistical errors, and thus contribute negligibly to the overall DESI error budget., Comment: 53 pages, 41 figures. Supporting paper for DESI DR1 cosmological measurements

Dark Energy Spectroscopic Instrument (DESI) uses more than 2.4 million Emission Line Galaxies (ELGs) for 3D large-scale structure (LSS) analyses in its Data Release 1 (DR1). Such large statistics enable thorough research on systematic uncertainties. In this study, we focus on spectroscopic systematics of ELGs. The redshift success rate ($f_{\rm goodz}$) is the relative fraction of secure redshifts among all measurements. It depends on observing conditions, thus introduces non-cosmological variations to the LSS. We, therefore, develop the redshift failure weight ($w_{\rm zfail}$) and a per-fibre correction ($\eta_{\rm zfail}$) to mitigate these dependences. They have minor influences on the galaxy clustering. For ELGs with a secure redshift, there are two subtypes of systematics: 1) catastrophics (large) that only occur in a few samples; 2) redshift uncertainty (small) that exists for all samples. The catastrophics represent 0.26\% of the total DR1 ELGs, composed of the confusion between O\,\textsc{ii} and sky residuals, double objects, total catastrophics and others. We simulate the realistic 0.26\% catastrophics of DR1 ELGs, the hypothetical 1\% catastrophics, and the truncation of the contaminated $1.31

We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produce catalogs that are weighted subsamples of the observed data, each matched to a weighted `random' catalog that forms an unclustered sampling of the probability density that DESI could have observed those data at each location. Precise knowledge of the DESI observing history and associated hardware performance allows for a determination of the DESI footprint and the number of times DESI has covered it at sub-arcsecond level precision. This enables the completeness of any DESI sample to be modeled at this same resolution. The pipeline developed to create LSS catalogs has been designed to easily allow robustness tests and enable future improvements. We describe how it allows ongoing work improving the match between galaxy and random catalogs, such as including further information when assigning redshifts to randoms, accounting for fluctuations in target density, accounting for variation in the redshift success rate, and accommodating blinding schemes., Comment: Accepted (by JCAP) version of supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

The Dark Energy Spectroscopic Instrument (DESI) survey will measure spectroscopic redshifts for millions of galaxies across roughly $14,000 \, \mathrm{deg}^2$ of the sky. Cross-correlating targets in the DESI survey with complementary imaging surveys allows us to measure and analyze shear distortions caused by gravitational lensing in unprecedented detail. In this work, we analyze a series of mock catalogs with ray-traced gravitational lensing and increasing sophistication to estimate systematic effects on galaxy-galaxy lensing estimators such as the tangential shear $\gamma_{\mathrm{t}}$ and the excess surface density $\Delta\Sigma$. We employ mock catalogs tailored to the specific imaging surveys overlapping with the DESI survey: the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS). Among others, we find that fiber incompleteness can have significant effects on galaxy-galaxy lensing estimators but can be corrected effectively by up-weighting DESI targets with fibers by the inverse of the fiber assignment probability. Similarly, we show that intrinsic alignment and lens magnification are expected to be statistically significant given the precision forecasted for the DESI year-1 data set. Our study informs several analysis choices for upcoming cross-correlation studies of DESI with DES, HSC, and KiDS., Comment: 18 pages, 13 figures, version accepted for publication

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$\alpha$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range $0.1-1$ and $w_a<0$. This preference is 2.6$\sigma$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $\Lambda$CDM model at the $2.5\sigma$, $3.5\sigma$ or $3.9\sigma$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $\Lambda$CDM model with the sum of neutrino mass $\sum m_\nu$ free, combining the DESI and CMB data yields an upper limit $\sum m_\nu < 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_\nu>0$ $(\sum m_\nu>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $\Lambda$CDM. [Abridged.], Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). 68 pages, 15 figures. Version accepted for publication in JCAP

We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$\alpha$ (Ly$\alpha$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$\alpha$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a new analysis methodology on a blinded dataset. We conducted rigorous tests using synthetic data to ensure the reliability of our methodology and findings before unblinding. Additionally, we conducted multiple data splits to assess the consistency of the results and scrutinized various analysis approaches to confirm their robustness. For a given value of the sound horizon ($r_d$), we measure the expansion at $z_{\rm eff}=2.33$ with 2\% precision, $H(z_{\rm eff}) = (239.2 \pm 4.8) (147.09~{\rm Mpc} /r_d)$ km/s/Mpc. Similarly, we present a 2.4\% measurement of the transverse comoving distance to the same redshift, $D_M(z_{\rm eff}) = (5.84 \pm 0.14) (r_d/147.09~{\rm Mpc})$ Gpc. Together with other DESI BAO measurements at lower redshifts, these results are used in a companion paper to constrain cosmological parameters., Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). Minor changes in v4, version accepted for publication in JCAP

We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1

The amplitude of the baryon signature in galaxy clustering depends on the cosmological baryon fraction. We consider two ways to isolate this signal in galaxy redshift surveys. First, we extend standard template-based Baryon Acoustic Oscillation (BAO) models to include the amplitude of the baryonic signature by splitting the transfer function into baryon and cold dark matter components with freely varying proportions. Second, we include the amplitude of the split as an extra parameter in Effective Field Theory of Large Scale Structure (EFT) models of the full galaxy clustering signal. We find similar results from both approaches. For the Baryon Oscillation Spectroscopic Survey (BOSS) data we find $f_b\equiv\Omega_b/\Omega_m=0.173\pm0.027$ for template fits post-reconstruction, $f_b=0.153\pm0.029$ for template fits pre-reconstruction, and $f_b=0.154\pm0.022$ for EFT fits, with an estimated systematic error of 0.013 for all three methods. Using reconstruction only produces a marginal improvement for these measurements. Although significantly weaker than constraints on $f_b$ from the Cosmic Microwave Background, these measurements rely on very simple physics and, in particular, are independent of the sound horizon. In a companion paper we show how they can be used, together with Big Bang Nucleosynthesis measurements of the physical baryon density and geometrical measurements of the matter density from the Alcock-Paczynski effect, to constrain the Hubble parameter. While the constraints on $H_0$ based on density measurements from BOSS are relatively weak, measurements from DESI and Euclid will lead to errors on $H_0$ that are competitive with those from local distance ladder measurements., Comment: 39 pages, 20 figures, v2 updated to match accepted version in PRD

We introduce a new method for measuring the Hubble parameter from low-redshift large-scale observations that is independent of the comoving sound horizon. The method uses the baryon-to-photon ratio determined by the primordial deuterium abundance, together with Big Bang Nucleosynthesis (BBN) calculations and the present-day CMB temperature to determine the physical baryon density $\Omega_b h^2$. The baryon fraction $\Omega_b/\Omega_m$ is measured using the relative amplitude of the baryonic signature in galaxy clustering measured by the Baryon Oscillation Spectroscopic Survey, scaling the physical baryon density to the physical matter density. The physical density $\Omega_mh^2$ is then compared with the geometrical density $\Omega_m$ from Alcock-Paczynski measurements from Baryon Acoustic Oscillations (BAO) and voids, to give $H_0$. Including type Ia supernovae and uncalibrated BAO, we measure $H_0 = 67.1^{+6.3}_{-5.3}$ km s$^{-1}$ Mpc$^{-1}$. We find similar results when varying analysis choices, such as measuring the baryon signature from the reconstructed correlation function, or excluding supernovae or voids. This measurement is currently consistent with both the distance-ladder and CMB $H_0$ determinations, but near-future large-scale structure surveys will obtain 3--4$\times$ tighter constraints., Comment: 8 pages, 2 figures, comments welcome!

When fitting cosmological models to data, a Bayesian framework is commonly used, requiring assumptions on the form of the likelihood and model prior. In light of current tensions between different data, it is interesting to investigate the robustness of cosmological measurements to statistical assumptions about the likelihood distribution from which the data was drawn. We consider the impact of changes to the likelihood caused by uncertainties due to the finite number of mock catalogs used to estimate the covariance matrix, leading to the replacement of the standard Gaussian likelihood with a multivariate $t$-distribution. These changes to the likelihood have a negligible impact on recent cosmic microwave background (CMB) lensing and baryon acoustic oscillation (BAO) measurements. We then extend our analysis to perform a sensitivity test on the Gaussian likelihoods typically adopted, considering how increasing the size of the tails of the likelihood (again using a $t$-distribution) affects cosmological inferences. For an open $\Lambda$CDM model constrained by BAO alone, we find that increasing the weight in the tails shifts and broadens the resulting posterior on the parameters, with a $\sim$0.2-0.4$\sigma$ effect on $\Omega_\Lambda$ and $\Omega_{\rm k}$. In contrast, the CMB temperature and polarization constraints in $\Lambda$CDM showed less than 0.03$\sigma$ changes in the parameters, except for $\{\tau$, ln($10^{10}A_{\rm s})$, $\sigma_8$, $S_8$, $\sigma_8\Omega_{\rm m}^{0.25}$, $z_{\rm re}$, $10^9A_{\rm s}e^{-2\tau}\}$ which shifted by around 0.1-0.2$\sigma$. If we use solely $\ell < 30$ data, the amplitude $A_{\rm s} e^{-2\tau}$ varies in the posterior mean by 0.7$\sigma$ and the error bars increase by 6%. We conclude, at least for current-generation CMB and BAO measurements, that uncertainties in the shape and tails of the likelihood do not contribute to current tensions., Comment: All comments are welcome, 22 pages, 4 figures

Galaxy-galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) dataset promise to yield unprecedented combined-probe tests of cosmology and the galaxy-halo connection. In such analyses, it is essential to identify and characterise all relevant statistical and systematic errors. In this paper, we forecast the covariances of DESI Y1 GGL+clustering measurements and characterise the systematic bias due to redshift evolution in the lens samples. Focusing on the projected clustering and galaxy-galaxy lensing correlations, we compute a Gaussian analytical covariance, using a suite of N-body and log-normal simulations to characterise 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 compared to the expected statistical errors. For BGS, we find less evolution effects for brighter absolute magnitude cuts, at the cost of reduced sample size. We find that with a fiducial redshift bin width delta z = 0.1, evolution effects on GGL is negligible across all scales, all fiducial selection cuts, all fiducial redshift bins, given DESI Y1 sample size. 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.1Mpc/h. 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, Comment: 19 pages, 15 figures, submitted to MNRAS

We present a first measurement of the galaxy-galaxy-CMB lensing bispectrum. The signal is detected at $26\sigma$ and $22\sigma$ significance using two samples from the unWISE galaxy catalog at mean redshifts $\bar{z}=0.6$ and $1.1$ and lensing reconstructions from Planck PR4. We employ a compressed bispectrum estimator based on the cross-correlation between the square of the galaxy overdensity field and CMB lensing reconstructions. We present a series of consistency tests to ensure the cosmological origin of our signal and rule out potential foreground contamination. We compare our results to model predictions from a halo model previously fit to only two-point spectra, finding reasonable agreement when restricting our analysis to large scales. Such measurements of the CMB lensing galaxy bispectrum will have several important cosmological applications, including constraining the uncertain higher-order bias parameters that currently limit lensing cross-correlation analyses., Comment: 10 pages, 4 figures, 1 table; v2: replaced HOD model from arXiv:2203.12583 with updated/improved model from arXiv:2303.08121, conclusions unchanged; to be submitted to PRL

We present tomographic measurements of structure growth using cross-correlations of Atacama Cosmology Telescope (ACT) DR6 and Planck CMB lensing maps with the unWISE Blue and Green galaxy samples, which span the redshift ranges $0.2 \lesssim z \lesssim 1.1$ and $0.3 \lesssim z \lesssim 1.8$, respectively. We improve on prior unWISE cross-correlations not just by making use of the new, high-precision ACT DR6 lensing maps, but also by including additional spectroscopic data for redshift calibration and by analysing our measurements with a more flexible theoretical model. An extensive suite of systematic and null tests within a blind analysis framework ensures that our results are robust. We determine the amplitude of matter fluctuations at low redshifts ($z\simeq 0.2-1.6$), finding $S_8 \equiv \sigma_8 (\Omega_m / 0.3)^{0.5} = 0.813 \pm 0.021$ using the ACT cross-correlation alone and $S_8 = 0.810 \pm 0.015$ with a combination of Planck and ACT cross-correlations; these measurements are fully consistent with the predictions from primary CMB measurements assuming standard structure growth. The addition of Baryon Acoustic Oscillation data breaks the degeneracy between $\sigma_8$ and $\Omega_m$, allowing us to measure $\sigma_8 = 0.813 \pm 0.020$ from the cross-correlation of unWISE with ACT and $\sigma_8 = 0.813\pm 0.015$ from the combination of cross-correlations with ACT and Planck. These results also agree with the expectations from primary CMB extrapolations in $\Lambda$CDM cosmology; the consistency of $\sigma_8$ derived from our two redshift samples at $z \sim 0.6$ and $1.1$ provides a further check of our cosmological model. Our results suggest that structure formation on linear scales is well described by $\Lambda$CDM even down to low redshifts $z\lesssim 1$., Comment: 73 pages (incl. 31 pages of appendices), 52 figures, 16 tables, published in ApJ. Watch G. S. Farren and A. Krolewski discuss the analysis and results under https://cosmologytalks.com/2023/09/11/act-unwise

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 $\fnl$. Our sample comprises over 12 million targets, covering 14,000 square degrees 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 $\fnl$ 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 $\fnl = 34^{+24(+50)}_{-44(-73)}$ at 68\%(95\%) 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 $\fnl$ 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 over-correction, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of $\fnl$ with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the $\fnl$ uncertainty., Comment: 21 pages, 17 figures, 7 tables (Appendix excluded). Published in MNRAS

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra., Comment: 43 pages, 7 figures, 17 tables, accepted for publication in the Astronomical Journal

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1

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 that cover a wide range of scales and redshifts 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 in the Born approximation, and accompanying weak lensing mock catalogues, tuned via fits to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument (DESI). 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 with $\Delta z = 0.05$ as well as CMB convergence maps ($z \approx 1090$) for each of the 25 ${\tt base}$-resolution simulations ($L_{\rm box} = 2000\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 6912^3$) as well as for the two ${\tt huge}$ simulations ($L_{\rm box} = 7500\,h^{-1}{\rm Mpc}$, $N_{\rm part} = 8640^3$) at the fiducial AbacusSummit cosmology ($Planck$ 2018). The pixel resolution of each map is 0.21 arcmin, corresponding to a HEALPiX $N_{\rm side}$ of 16384. The sky coverage of the ${\tt base}$ simulations is an octant until $z \approx 0.8$ (decreasing to about 1800 deg$^2$ at $z \approx 2.4$), whereas the ${\tt huge}$ simulations offer full-sky coverage until $z \approx 2.2$. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper-Suprime Cam weak lensing datasets. The produced mock catalogues are validated against theoretical predictions for various clustering and lensing statistics such as galaxy clustering multipoles, galaxy-shear and shear-shear, showing excellent agreement., Comment: 21 pages, 13 figures, data available at https://app.globus.org/file-manager?origin_id=3dd6566c-eed2-11ed-ba43-09d6a6f08166&path=%2F

We detect the cross-correlation between 2.7 million DESI quasar targets across 14,700 deg$^2$ (180 quasars deg$^{-2}$) and Planck 2018 CMB lensing at $\sim$30$\sigma$. We use the cross-correlation on very large scales to constrain local primordial non-Gaussianity via the scale dependence of quasar bias. The DESI quasar targets lie at an effective redshift of 1.51 and are separated into four imaging regions of varying depth and image quality. We select quasar targets from Legacy Survey DR9 imaging, apply additional flux and photometric redshift cuts to improve the purity and reduce the fraction of unclassified redshifts, and use early DESI spectroscopy of 194,000 quasar targets to determine their redshift distribution and stellar contamination fraction (2.6%). Due to significant excess large-scale power in the quasar autocorrelation, we apply weights to mitigate contamination from imaging systematics such as depth, extinction, and stellar density. We use realistic contaminated mocks to determine the greatest number of systematic modes that we can fit, before we are biased by overfitting and spuriously remove real power. We find that linear regression with one to seven imaging templates removed per region accurately recovers the input cross-power, $f_{\textrm{NL}}$ and linear bias. As in previous analyses, our $f_{\textrm{NL}}$ constraint depends on the linear primordial non-Gaussianity bias parameter, $b_{\phi} = 2(b - p)\delta_c$ assuming universality of the halo mass function. We measure $f_{\textrm{NL}} = -26^{+45}_{-40}$ with $p=1.6$ $(f_{\textrm{NL}} = -18^{+29}_{-27}$ with $p=1.0$), and find that this result is robust under several systematics tests. Future spectroscopic quasar cross-correlations with Planck lensing lensing can tighten the $f_{\textrm{NL}}$ constraints by a factor of 2 if they can remove the excess power on large scales in the quasar auto power spectrum., Comment: 57 pages, 25 figures, updated to match accepted version

Galaxy formation theories predict that galaxy shapes and angular momenta have non-random alignments with the cosmic web. This leads to so-called intrinsic alignment between pairs of galaxies, which is important to quantify as a nuisance parameter for weak lensing. We study galaxy-cosmic web alignment in the IllustrisTNG suite of hydrodynamical simulations at redshifts 1 and 2, finding that alignment trends are consistent with previous studies. However, we find that the magnitude of the spin alignment signal is $\sim 2.4 \times$ weaker than seen in previous studies of the Horizon-AGN simulation, suggesting that this signal may have significant dependence on subgrid physics. Based on IllustrisTNG, we then construct mock observational spectroscopic surveys that can probe shape-cosmic web alignment at $z \sim 1-2$, modeled on the low-$z$ galaxy redshift and IGM tomography surveys on the upcoming Subaru Prime Focus Spectrograph Galaxy Evolution (PFS GE) survey. However, even over box sizes of $L=205\,h^{-1}\,\mathrm{Mpc}$, we find that global anisotropies induce a sample variance in the 2D projected alignment signal that depend on the projected direction -- this induces significant errors in the observed alignment. We predict a $5.3\sigma$ detection of IllustrisTNG's shape alignment signal at $z \sim 1$ from Subaru PFS GE, although a detection would be challenging at $z \sim 2$. However, a rough rescaling of the relative alignment signal strengths between the TNG and HorizonAGN simulations suggests that PFS GE should be able to more easily constrain the latter's stronger signal., Comment: 28 pages, 10 figures. Accepted by The Astrophysical Journal

The kinetic Sunyaev-Zel'dovich (kSZ) effect, i.e., the Doppler boost of cosmic microwave background (CMB) photons caused by their scattering off free electrons in galaxy clusters and groups with non-zero bulk velocity, is a powerful window on baryons in the universe. We present the first halo-model computation of the cross-power spectrum of the ``projected-field'' kSZ signal with large-scale structure (LSS) tracers. We compare and validate our calculations against previous studies, which relied on $N$-body-calibrated effective formulas rather than the halo model. We forecast results for CMB maps from the Atacama Cosmology Telescope (AdvACT), Simons Observatory (SO), and CMB-S4, and LSS survey data from the Dark Energy Survey, the Vera C.~Rubin Observatory (VRO), and \textit{Euclid}. In cross-correlation with galaxy number density, for AdvACT $\times$ \textit{unWISE} we forecast an 18$\sigma$ projected-field kSZ detection using data already in hand. Combining SO CMB maps and \textit{unWISE} galaxy catalogs, we expect a $62\sigma$ detection, yielding precise measurements of the gas density profile radial slopes. Additionally, we forecast first detections of the kSZ -- galaxy weak lensing cross-correlation with AdvACT $\times$ VRO/\textit{Euclid} (at 6$\sigma$) and of the kSZ -- CMB weak lensing cross-correlation with SO (at 16$\sigma$). Finally, $\approx 10-20$\% precision measurements of the shape of the gas density profile should be possible with CMB-S4 kSZ -- CMB lensing cross-correlation without using any external datasets.

Combining different observational probes, such as galaxy clustering and weak lensing, is a promising technique for unveiling the physics of the Universe with upcoming dark energy experiments. The galaxy redshift sample from the Dark Energy Spectroscopic Instrument (DESI) will have a significant overlap with major ongoing imaging surveys specifically designed for weak lensing measurements: the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES) and the Hyper Suprime-Cam (HSC) survey. In this work we analyse simulated redshift and lensing catalogues to establish a new strategy for combining high-quality cosmological imaging and spectroscopic data, in view of the first-year data assembly analysis of DESI. In a test case fitting for a reduced parameter set, we employ an optimal data compression scheme able to identify those aspects of the data that are most sensitive to the cosmological information, and amplify them with respect to other aspects of the data. We find this optimal compression approach is able to preserve all the information related to the growth of structure; we also extend this scheme to derive weights to be applied to individual galaxies, and show that these produce near-optimal results., Comment: 14 pages, 12 Figures, DESI collaboration article

The Dark Energy Spectroscopic Instrument (DESI) has embarked on an ambitious five-year survey to explore the nature of dark energy with spectroscopy of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the Baryon Acoustic Oscillation method to measure distances from the nearby universe to z > 3.5, as well as measure the growth of structure and probe potential modifications to general relativity. In this paper we describe the significant instrumentation we developed for the DESI survey. The new instrumentation includes a wide-field, 3.2-deg diameter prime-focus corrector that focuses the light onto 5020 robotic fiber positioners on the 0.812 m diameter, aspheric focal surface. The positioners and their fibers are divided among ten wedge-shaped petals. Each petal is connected to one of ten spectrographs via a contiguous, high-efficiency, nearly 50 m fiber cable bundle. The ten spectrographs each use a pair of dichroics to split the light into three channels that together record the light from 360 - 980 nm with a resolution of 2000 to 5000. We describe the science requirements, technical requirements on the instrumentation, and management of the project. DESI was installed at the 4-m Mayall telescope at Kitt Peak, and we also describe the facility upgrades to prepare for DESI and the installation and functional verification process. DESI has achieved all of its performance goals, and the DESI survey began in May 2021. Some performance highlights include RMS positioner accuracy better than 0.1", SNR per \sqrt{\AA} > 0.5 for a z > 2 quasar with flux 0.28e-17 erg/s/cm^2/A at 380 nm in 4000s, and median SNR = 7 of the [OII] doublet at 8e-17 erg/s/cm^2 in a 1000s exposure for emission line galaxies at z = 1.4 - 1.6. We conclude with highlights from the on-sky validation and commissioning of the instrument, key successes, and lessons learned. (abridged), Comment: 78 pages, 32 figures, submitted to AJ

We present the first detailed analysis of the connection between galaxies and their dark matter halos for the unWISE galaxy catalog -- a full-sky, infrared-selected sample built from WISE data, containing over 500 million galaxies. Using unWISE galaxy-galaxy auto-correlation and Planck CMB lensing-galaxy cross-correlation measurements down to 10 arcmin angular scales, we constrain the halo occupation distribution (HOD), a model describing how central and satellite galaxies are distributed within dark matter halos, for three unWISE} galaxy samples at mean redshifts $\bar{z} \approx 0.6$, $1.1$, and $1.5$. We constrain the characteristic minimum halo mass to host a central galaxy, $M_\mathrm{min}^\mathrm{HOD} = 1.83^{+0.41}_{-1.63} \times 10^{12} M_\odot/h$, $5.22^{+0.34}_{-4.80} \times 10^{12} M_\odot/h$, $6.60 ^{+0.30}_{-1.11} \times 10^{13} M_\odot/h$ for the unWISE samples at $\bar{z}\approx 0.6$, $1.1$, and $1.5$, respectively. We find that all three samples are dominated by central galaxies, rather than satellites. Using our constrained HOD models, we infer the effective linear galaxy bias for each unWISE sample, and find that it does not evolve as steeply with redshift as found in previous perturbation-theory-based analyses of these galaxies. We discuss possible sources of systematic uncertainty in our results, the most significant of which is the uncertainty on the galaxy redshift distribution. Our HOD constraints provide a detailed, quantitative understanding of how the unWISE galaxies populate the underlying dark matter halo distribution. These constraints will have a direct impact on future studies employing the unWISE galaxies as a cosmological and astrophysical probe, including measurements of ionized gas thermodynamics and dark matter profiles via Sunyaev-Zel'dovich and lensing cross-correlations.

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan., Comment: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

We apply the cross-correlation technique to infer the Hubble constant ($H_0$) of the Universe using gravitational wave (GW) sources without electromagnetic counterparts (dark sirens) from the third GW Transient Catalog (GWTC-3) and the photometric galaxy surveys 2MPZ and WISE-SuperCOSMOS, and combine these with the bright siren measurement of $H_0$ from GW170817. The posterior on $H_0$ with only dark sirens is uninformative due to the small number of well-localised GW sources. Using the eight well-localized dark sirens and the binary neutron star GW170817 with EM counterpart, we obtain a value of the Hubble constant $H_0= 75.4_{-6}^{+11}$ km/s/Mpc (median and 68.3$\%$ equal-tailed interval (ETI)) after marginalizing over the matter density and the GW bias parameters. This measurement is mostly driven by the bright siren measurement and any constraint from dark sirens is not statistically significant. In the future, with more well-localized GW events, the constraints on expansion history will improve., Comment: 13 pages, 9 figures. Published in The Astrophysical Journal