Showing 1–50 of 90 results for author: Howlett, C

STag II: Classification of Serendipitous Supernovae Observed by Galaxy Redshift Surveys

Authors: W. Davison, D. Parkinson, S. BenZvi, A. Palmese, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, A. de la Macorra, Arjun Dey, P. Doel, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, S. Juneau, T. Kisner, A. Kremin, A. Lambert, M. Landriau, L. Le Guillou, A. Meisner, R. Miquel, J. Moustakas, A. D. Myers, C. Poppett , et al. (10 additional authors not shown)

Abstract: With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic re… ▽ More With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic regression in order to assign 'tags' to spectra based on spectral features. STag II is a continuation of this work, which now makes use of model supernova spectra combined with real DESI spectra in order to train STag to better deal with realistic data. We also make use of the rlap score as a trustworthiness cut, making for a more robust and accurate supernova classifier than before. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 27 pages, 7 figures. The STag code is available at https://github.com/wdavison909/STag

A Trifecta of Modelling Tools: Navigating COMPAS, Shark and Bayesian Inference for Binary Black Hole Model Selection

Authors: Liana Rauf, Cullan Howlett, Simon Stevenson, Jeff Riley

Abstract: Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the L… ▽ More Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the LVK transient catalogue. For our COMPAS models, we calculate the likelihood with and without the dependence on the predicted number of BBH merger events. We find strong correlations between hyper-parameters governing the specific angular momentum (AM) of mass lost during mass transfer, the mass-loss rates of Wolf--Rayet stars via winds and the chemically homogeneous evolution (CHE) formation channel. We conclude that analysing the marginalised and unmarginalised likelihood is a good indicator of whether the population parameters distribution and number of observed events reflect the LVK data. In doing so, we see that the majority of the models preferred in terms of the population-level parameters of the BBHs greatly overpredict the number of events we should have observed to date. Looking at the smaller number of models which perform well with both likelihoods, we find that those with no CHE, AM loss occurring closer to the donor during the first mass-transfer event, and/or higher rates of mass-loss from Wolf--Rayet winds are generally preferred by current data. We find these conclusions to be robust to our choice of selection criteria. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 27 pages, 19 figures, 6 tables

Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators

Authors: M. Pinon, A. de Mattia, P. McDonald, E. Burtin, V. Ruhlmann-Kleider, M. White, D. Bianchi, A. J. Ross, J. Aguilar, S. Ahlen, D. Brooks, E. Chaussidon, T. Claybaugh, S. Cole, A. de la Macorra, B. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, D. Kirkby, T. Kisner, A. Kremin , et al. (27 additional authors not shown)

Abstract: We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We vali… ▽ More We present a method to mitigate the effects of fiber assignment incompleteness in two-point power spectrum and correlation function measurements from galaxy spectroscopic surveys, by truncating small angular scales from estimators. We derive the corresponding modified correlation function and power spectrum windows to account for the small angular scale truncation in the theory prediction. We validate this approach on simulations reproducing the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) with and without fiber assignment. We show that we recover unbiased cosmological constraints using small angular scale truncated estimators from simulations with fiber assignment incompleteness, with respect to standard estimators from complete simulations. Additionally, we present an approach to remove the sensitivity of the fits to high $k$ modes in the theoretical power spectrum, by applying a transformation to the data vector and window matrix. We find that our method efficiently mitigates the effect of fiber assignment incompleteness in two-point correlation function and power spectrum measurements, at low computational cost and with little statistical loss. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 32 pages, 22 figures

Blinding scheme for the scale-dependence bias signature of local primordial non-Gaussianity for DESI 2024

Authors: E. Chaussidon, A. de Mattia, C. Yèche, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, P. Doel, K. Fanning, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, T. Kisner, A. Lambert, L. Le Guillou, M. Manera, A. Meisner, R. Miquel, G. Niz, N. Palanque-Delabrouille, W. J. Percival, F. Prada, A. J. Ross , et al. (10 additional authors not shown)

Abstract: The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics… ▽ More The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics is important relative to the cosmological signal, and a detection of that signal would have significant implications. This is the case for local primordial non-gaussianity, as probed by the scale-dependent bias of the galaxy power spectrum at large scales that are heavily sensitive to the dependence of the target selection on the imaging quality, known as imaging systematics. We propose a blinding method for the scale-dependent bias signature of local primordial non-gaussianity at the density field level which consists in generating a set of weights for the data that replicate the scale-dependent bias. The applied blinding is predictable, and can be straightforwardly combined with other catalog-level blinding procedures that have been designed for the baryon acoustic oscillation and redshift space distortion signals. The procedure is validated through simulations that replicate data from the first year of observation of the Dark Energy Spectroscopic Instrument, but may find applications to other upcoming spectroscopic surveys. △ Less

Submitted 31 May, 2024; originally announced June 2024.

Archetype-Based Redshift Estimation for the Dark Energy Spectroscopic Instrument Survey

Authors: Abhijeet Anand, Julien Guy, Stephen Bailey, John Moustakas, J. Aguilar, S. Ahlen, A. Bolton, A. Brodzeller, D. Brooks, T. Claybaugh, S. Cole, B. Dey, K. Fanning, J. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, L. Le Guillou, G. Gutierrez, K. Honscheid, C. Howlett, S. Juneau, D. Kirkby, T. Kisner, A. Kremin, A. Lambert , et al. (24 additional authors not shown)

Abstract: We present a computationally efficient galaxy archetype-based redshift estimation and spectral classification method for the Dark Energy Survey Instrument (DESI) survey. The DESI survey currently relies on a redshift fitter and spectral classifier using a linear combination of PCA-derived templates, which is very efficient in processing large volumes of DESI spectra within a short time frame. Howe… ▽ More We present a computationally efficient galaxy archetype-based redshift estimation and spectral classification method for the Dark Energy Survey Instrument (DESI) survey. The DESI survey currently relies on a redshift fitter and spectral classifier using a linear combination of PCA-derived templates, which is very efficient in processing large volumes of DESI spectra within a short time frame. However, this method occasionally yields unphysical model fits for galaxies and fails to adequately absorb calibration errors that may still be occasionally visible in the reduced spectra. Our proposed approach improves upon this existing method by refitting the spectra with carefully generated physical galaxy archetypes combined with additional terms designed to absorb data reduction defects and provide more physical models to the DESI spectra. We test our method on an extensive dataset derived from the survey validation (SV) and Year 1 (Y1) data of DESI. Our findings indicate that the new method delivers marginally better redshift success for SV tiles while reducing catastrophic redshift failure by $10-30\%$. At the same time, results from millions of targets from the main survey show that our model has relatively higher redshift success and purity rates ($0.5-0.8\%$ higher) for galaxy targets while having similar success for QSOs. These improvements also demonstrate that the main DESI redshift pipeline is generally robust. Additionally, it reduces the false positive redshift estimation by $5-40\%$ for sky fibers. We also discuss the generic nature of our method and how it can be extended to other large spectroscopic surveys, along with possible future improvements. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: submitted to AAS journals, 29 pages, 13 figures

Candidate strongly-lensed Type Ia supernovae in the Zwicky Transient Facility archive

Authors: A. Townsend, J. Nordin, A. Sagués Carracedo, M. Kowalski, N. Arendse, S. Dhawan, A. Goobar, J. Johansson, E. Mörtsell, S. Schulze, I. Andreoni, E. Fernández, A. G. Kim, P. E. Nugent, F. Prada, M. Rigault, N. Sarin, D. Sharma, E. C. Bellm, M. W. Coughlin, R. Dekany, S. L. Groom, L. Lacroix, R. R. Laher, R. Riddle , et al. (39 additional authors not shown)

Abstract: Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existe… ▽ More Gravitationally lensed Type Ia supernovae (glSNe Ia) are unique astronomical tools for studying cosmological parameters, distributions of dark matter, the astrophysics of the supernovae and the intervening lensing galaxies themselves. Only a few highly magnified glSNe Ia have been discovered by ground-based telescopes, such as the Zwicky Transient Facility (ZTF), but simulations predict the existence of a fainter, undetected population. We present a systematic search in the ZTF archive of alerts from 1 June 2019 to 1 September 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. Candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. While superluminous supernovae (SLSNe) cannot be fully rejected, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 $\pm$ 3 and 34 $\pm$ 1 days for the two events, respectively, which suggests that we have uncovered a population with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, like the Vera C. Rubin Observatory's Legacy Survey of Space and Time. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 15 figures

ELG Spectroscopic Systematics Analysis of the DESI Data Release 1

Authors: Jiaxi Yu, Ashley J. Ross, Antoine Rocher, Otávio Alves, Arnaud de Mattia, Daniel Forero-Sánchez, Jean-Paul Kneib, Alex Krolewski, TingWen Lan, Michael Rashkovetskyi, Jessica Nicole Aguilar, Steven Ahlen, Stephen Bailey, David Brooks, Edmond Chaussidon, Todd Claybaugh, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Kevin Fanning, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid , et al. (36 additional authors not shown)

Abstract: 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… ▽ More 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 ($η_{\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<z<1.33$ in the \textsc{AbacusSummit} ELG mocks. Their $P_\ell$ show non-negligible bias from the uncontaminated mocks. But their influences on the redshift space distortions (RSD) parameters are smaller than $0.2σ$. The redshift uncertainty of \Yone ELGs is 8.5 km/s with a Lorentzian profile. The code for implementing the catastrophics and redshift uncertainty on mocks can be found in https://github.com/Jiaxi-Yu/modelling_spectro_sys. △ Less

Submitted 26 May, 2024; originally announced May 2024.

The Construction of Large-scale Structure Catalogs for the Dark Energy Spectroscopic Instrument

Authors: A. J. Ross, J. Aguilar, S. Ahlen, S. Alam, A. Anand, S. Bailey, D. Bianchi, S. Brieden, D. Brooks, E. Burtin, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, A. de Mattia, Arjun Dey, Biprateep Dey, P. Doel, K. Fanning, S. Ferraro, J. Ereza, A. Font-Ribera, J. E. Forero-Romero , et al. (59 additional authors not shown)

Abstract: 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 produ… ▽ More 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. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: Supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

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

Authors: Hui Kong, Ashley J. Ross, Klaus Honscheid, Dustin Lang, Anna Porredon, Arnaud de Mattia, Mehdi Rezaie, Rongpu Zhou, Edward Schlafly, John Moustakas, Alberto Rosado-Marin, Jessica Nicole Aguilar, Steven Ahlen, David Brooks, Edmond Chaussidon, Todd Claybaugh, Shaun Cole, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Kevin Fanning, Jaime E. Forero-Romero, Enrique Gaztanaga, Satya Gontcho A Gontcho , et al. (28 additional authors not shown)

Abstract: We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate the DES… ▽ More We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate 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 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 power 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 Large-Scale Structure systematic contamination (e.g. Cosmic-Infrared Background) in the Galactic extinction maps is a likely root cause. We additionally observe that the DESI LRGs sample exhibits a complex dependency on a combination of seeing, depth, and intrinsic galaxy brightness, which is not replicated by Obiwan, suggesting discrepancies between the current simulation settings and the actual observations. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRG sample. △ Less

Submitted 25 May, 2024; originally announced May 2024.

Comments: 46 pages, 26 figures

DESI 2024: Constraints on Physics-Focused Aspects of Dark Energy using DESI DR1 BAO Data

Authors: K. Lodha, A. Shafieloo, R. Calderon, E. Linder, W. Sohn, J. L. Cervantes-Cota, A. de Mattia, J. García-Bellido, M. Ishak, W. Matthewson, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, A. de la Macorra, A. Dey, B. Dey, P. Doel, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, S. Juneau, S. Kent, T. Kisner , et al. (25 additional authors not shown)

Abstract: Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range $z=0.1$-$4.2$. We use this data, together with other distance probes, to constrain the cosmic expansion history using some well-motivated physical classes of dark energy. In particular, we explore thre… ▽ More Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range $z=0.1$-$4.2$. We use this data, together with other distance probes, to constrain the cosmic expansion history using some well-motivated physical classes of dark energy. In particular, we explore three physics-focused behaviors of dark energy from the equation of state and energy density perspectives: the thawing class (matching many simple quintessence potentials), emergent class (where dark energy comes into being recently, as in phase transition models), and mirage class (where phenomenologically the distance to CMB last scattering is close to that from a cosmological constant $Λ$ despite dark energy dynamics). All three classes fit the data at least as well as $Λ$CDM, and indeed can improve on it by $Δχ^2\approx -5$ to $-17$ for the combination of DESI BAO with CMB and supernova data, while having one more parameter. The mirage class does essentially as well as $w_0w_a$CDM while having one less parameter. These classes of dynamical behaviors highlight worthwhile avenues for further exploration into the nature of dark energy. △ Less

Submitted 30 May, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: 20 pages, 8 figures. Metadata updated, comments welcome

Galaxy And Mass Assembly (GAMA): Stellar-to-Dynamical Mass Relation II. Peculiar Velocities

Authors: M. Burak Dogruel, Edward Taylor, Michelle Cluver, Matthew Colless, Anna de Graaff, Alessandro Sonnenfeld, John R. Lucey, Francesco D'Eugenio, Cullan Howlett, Khaled Said

Abstract: Empirical correlations connecting starlight to galaxy dynamics (e.g., the fundamental plane (FP) of elliptical/quiescent galaxies and the Tully--Fisher relation of spiral/star-forming galaxies) provide cosmology-independent distance estimation and are central to local Universe cosmology. In this work, we introduce the mass hyperplane (MH), which is the stellar-to-dynamical mass relation… ▽ More Empirical correlations connecting starlight to galaxy dynamics (e.g., the fundamental plane (FP) of elliptical/quiescent galaxies and the Tully--Fisher relation of spiral/star-forming galaxies) provide cosmology-independent distance estimation and are central to local Universe cosmology. In this work, we introduce the mass hyperplane (MH), which is the stellar-to-dynamical mass relation $(M_\star/M_\mathrm{dyn})$ recast as a linear distance indicator. Building on recent FP studies, we show that both star-forming and quiescent galaxies follow the same empirical MH, then use this to measure the peculiar velocities (PVs) for a sample of 2496 galaxies at $z<0.12$ from GAMA. The limiting precision of MH-derived distance/PV estimates is set by the intrinsic scatter in size, which we find to be $\approx$0.1~dex for both quiescent and star-forming galaxies (when modeled independently) and $\approx$0.11~dex when all galaxies are modeled together; showing that the MH is as good as the FP. To empirically validate our framework and distance/PV estimates, we compare the inferred distances to groups as derived using either quiescent or star-forming galaxies. A good agreement is obtained with no discernible bias or offset, having a scatter of $\approx$0.05~dex $\approx$12\% in distance. Further, we compare our PV measurements for the quiescent galaxies to the previous PV measurements of the galaxies in common between GAMA and the Sloan Digital Sky Survey (SDSS), which shows similarly good agreement. Finally, we provide comparisons of PV measurements made with the FP and the MH, then discuss possible improvements in the context of upcoming surveys such as the 4MOST Hemisphere Survey (4HS). △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted: 15th May 2024

DESI 2024: Reconstructing Dark Energy using Crossing Statistics with DESI DR1 BAO data

Authors: R. Calderon, K. Lodha, A. Shafieloo, E. Linder, W. Sohn, A. de Mattia, J. L. Cervantes-Cota, R. Crittenden, T. M. Davis, M. Ishak, A. G. Kim, W. Matthewson, G. Niz, S. Park, J. Aguilar, S. Ahlen, S. Allen, D. Brooks, T. Claybaugh, A. de la Macorra, A. Dey, B. Dey, P. Doel, J. E. Forero-Romero, E. Gaztañaga , et al. (30 additional authors not shown)

Abstract: We implement Crossing Statistics to reconstruct in a model-agnostic manner the expansion history of the universe and properties of dark energy, using DESI Data Release 1 (DR1) BAO data in combination with one of three different supernova compilations (PantheonPlus, Union3, and DES-SN5YR) and Planck CMB observations. Our results hint towards an evolving and emergent dark energy behaviour, with negl… ▽ More We implement Crossing Statistics to reconstruct in a model-agnostic manner the expansion history of the universe and properties of dark energy, using DESI Data Release 1 (DR1) BAO data in combination with one of three different supernova compilations (PantheonPlus, Union3, and DES-SN5YR) and Planck CMB observations. Our results hint towards an evolving and emergent dark energy behaviour, with negligible presence of dark energy at $z\gtrsim 1$, at varying significance depending on data sets combined. In all these reconstructions, the cosmological constant lies outside the 95\% confidence intervals for some redshift ranges. This dark energy behaviour, reconstructed using Crossing Statistics, is in agreement with results from the conventional $w_0$--$w_a$ dark energy equation of state parametrization reported in the DESI Key cosmology paper. Our results add an extensive class of model-agnostic reconstructions with acceptable fits to the data, including models where cosmic acceleration slows down at low redshifts. We also report constraints on \Hord\ from our model-agnostic analysis, independent of the pre-recombination physics. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 24 pages, 10 figures

The MOST Hosts Survey: spectroscopic observation of the host galaxies of ~40,000 transients using DESI

Authors: Maayane T. Soumagnac, Peter Nugent, Robert A. Knop, Anna Y. Q. Ho, William Hohensee, Autumn Awbrey, Alexis Andersen, Greg Aldering, Matan Ventura, Jessica N. Aguilar, Steven Ahlen, Segev Y. Benzvi, David Brooks, Dillon Brout, Todd Claybaugh, Tamara M. Davis, Kyle Dawson, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Kelly A. Douglass, Jaime E. Forero-Romero, Enrique Gaztanaga, Satya Gontcho A Gontcho , et al. (32 additional authors not shown)

Abstract: We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF,… ▽ More We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF, SDSS II, ZTF, DECAT, DESIRT). Scientific questions for which the MOST Hosts survey will be useful include Type Ia supernova cosmology, fundamental plane and peculiar velocity measurements, and the understanding of the correlations between transients and their host galaxy properties. Here, we present the first release of the MOST Hosts survey: 21,931 hosts of 20,235 transients. These numbers represent 36% of the final MOST Hosts sample, consisting of 60,212 potential host galaxies of 38,603 transients (a transient can be assigned multiple potential hosts). Of these galaxies, 40% do not appear in the DESI primary target list and therefore require a specific program like MOST Hosts. Of all the transients in the MOST Hosts list, only 26.7% have existing classifications, and so the survey will provide redshifts (and luminosities) for nearly 30,000 transients. A preliminary Hubble diagram and a transient luminosity-duration diagram are shown as examples of future potential uses of the MOST Hosts survey. The survey will also provide a training sample of spectroscopically observed transients for photometry-only classifiers, as we enter an era when most newly observed transients will lack spectroscopic classification. The MOST Hosts DESI survey data will be released through the Wiserep platform on a rolling cadence and updated to match the DESI releases. Dates of future releases and updates are available through the https://mosthosts.desi.lbl.gov website. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to ApJS

An analysis of parameter compression and full-modeling techniques with Velocileptors for DESI 2024 and beyond

Authors: M. Maus, S. Chen, M. White, J. Aguilar, S. Ahlen, A. Aviles, S. Brieden, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, Arjun Dey, P. Doel, S. Ferraro, N. Findlay, J. E. Forero-Romero, E. Gaztañaga, H. Gil-Marín, S. Gontcho A Gontcho, C. Hahn, K. Honscheid, C. Howlett, M. Ishak, S. Juneau, A. Kremin , et al. (30 additional authors not shown)

Abstract: In anticipation of forthcoming data releases of current and future spectroscopic surveys, we present the validation tests and analysis of systematic effects within \texttt{velocileptors} modeling pipeline when fitting mock data from the \texttt{AbacusSummit} N-body simulations. We compare the constraints obtained from parameter compression methods to the direct fitting (Full-Modeling) approaches o… ▽ More In anticipation of forthcoming data releases of current and future spectroscopic surveys, we present the validation tests and analysis of systematic effects within \texttt{velocileptors} modeling pipeline when fitting mock data from the \texttt{AbacusSummit} N-body simulations. We compare the constraints obtained from parameter compression methods to the direct fitting (Full-Modeling) approaches of modeling the galaxy power spectra, and show that the ShapeFit extension to the traditional template method is consistent with the Full-Modeling method within the standard $Λ$CDM parameter space. We show the dependence on scale cuts when fitting the different redshift bins using the ShapeFit and Full-Modeling methods. We test the ability to jointly fit data from multiple redshift bins as well as joint analysis of the pre-reconstruction power spectrum with the post-reconstruction BAO correlation function signal. We further demonstrate the behavior of the model when opening up the parameter space beyond $Λ$CDM and also when combining likelihoods with external datasets, namely the Planck CMB priors. Finally, we describe different parametrization options for the galaxy bias, counterterm, and stochastic parameters, and employ the halo model in order to physically motivate suitable priors that are necessary to ensure the stability of the perturbation theory. △ Less

Submitted 17 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: 56 pages, 23 figures. Supporting publication of DESI 2024 V: Analysis of the full shape of two-point clustering statistics from galaxies and quasars

A comparison between Shapefit compression and Full-Modelling method with PyBird for DESI 2024 and beyond

Authors: Y. Lai, C. Howlett, M. Maus, H. Gil-Marín, H. E. Noriega, S. Ramírez-Solano, P. Zarrouk, J. Aguilar, S. Ahlen, O. Alves, A. Aviles, D. Brooks, S. Chen, T. Claybaugh, T. M. Davis, K. Dawson, A. de la Macorra, P. Doel, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, K. Honscheid, S. Juneau, M. Landriau, M. Manera , et al. (18 additional authors not shown)

Abstract: DESI aims to provide one of the tightest constraints on cosmological parameters by analysing the clustering of more than thirty million galaxies. However, obtaining such constraints requires special care in validating the methodology and efforts to reduce the computational time required through data compression and emulation techniques. In this work, we perform a rigorous validation of the PyBird… ▽ More DESI aims to provide one of the tightest constraints on cosmological parameters by analysing the clustering of more than thirty million galaxies. However, obtaining such constraints requires special care in validating the methodology and efforts to reduce the computational time required through data compression and emulation techniques. In this work, we perform a rigorous validation of the PyBird power spectrum modelling code with both a traditional emulated Full-Modelling approach and the model-independent ShapeFit compression approach. By using cubic box simulations that accurately reproduce the clustering and precision of the DESI survey, we find that the cosmological constraints from ShapeFit and Full-Modelling are consistent with each other at the $\sim0.3σ$ level for the $Λ$CDM model. Both ShapeFit and Full-Modelling are also consistent with the true $Λ$CDM simulation cosmology down to a scale of $k_{\mathrm{max}} = 0.20 h\mathrm{Mpc}^{-1}$ even after including the hexadecapole. For extended models such as the wCDM and the oCDM models, we find that including the hexadecapole can significantly improve the constraints and reduce the modelling errors with the same $k_{\mathrm{max}}$. While their discrepancies between the constraints from ShapeFit and Full-Modelling are more significant than $Λ$CDM, they remain consistent within $0.7σ$. Lastly, we also show that the constraints on cosmological parameters with the correlation function evaluated from PyBird down to $s_{\mathrm{min}} = 30 h^{-1} \mathrm{Mpc}$ are unbiased and consistent with the constraints from the power spectrum. △ Less

Submitted 26 June, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 V: Analysis of the full shape of two-point clustering statistics from galaxies and quasars (In prep). 50 pages, 21 figures, and 12 tables. Revised version for JCAP

Validating the Galaxy and Quasar Catalog-Level Blinding Scheme for the DESI 2024 analysis

Authors: U. Andrade, J. Mena-Fernández, H. Awan, A. J. Ross, S. Brieden, J. Pan, A. de Mattia, J. Aguilar, S. Ahlen, O. Alves, D. Brooks, E. Buckley-Geer, E. Chaussidon, T. Claybaugh, S. Cole, A. de la Macorra, Arjun Dey, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, H. Gil-Marín, S. Gontcho A Gontcho, J. Guy, C. Hahn , et al. (38 additional authors not shown)

Abstract: In the era of precision cosmology, ensuring the integrity of data analysis through blinding techniques is paramount -- a challenge particularly relevant for the Dark Energy Spectroscopic Instrument (DESI). DESI represents a monumental effort to map the cosmic web, with the goal to measure the redshifts of tens of millions of galaxies and quasars. Given the data volume and the impact of the finding… ▽ More In the era of precision cosmology, ensuring the integrity of data analysis through blinding techniques is paramount -- a challenge particularly relevant for the Dark Energy Spectroscopic Instrument (DESI). DESI represents a monumental effort to map the cosmic web, with the goal to measure the redshifts of tens of millions of galaxies and quasars. Given the data volume and the impact of the findings, the potential for confirmation bias poses a significant challenge. To address this, we implement and validate a comprehensive blind analysis strategy for DESI Data Release 1 (DR1), tailored to the specific observables DESI is most sensitive to: Baryonic Acoustic Oscillations (BAO), Redshift-Space Distortion (RSD) and primordial non-Gaussianities (PNG). We carry out the blinding at the catalog level, implementing shifts in the redshifts of the observed galaxies to blind for BAO and RSD signals and weights to blind for PNG through a scale-dependent bias. We validate the blinding technique on mocks, as well as on data by applying a second blinding layer to perform a battery of sanity checks. We find that the blinding strategy alters the data vector in a controlled way such that the BAO and RSD analysis choices do not need any modification before and after unblinding. The successful validation of the blinding strategy paves the way for the unblinded DESI DR1 analysis, alongside future blind analyses with DESI and other surveys. △ Less

Submitted 15 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Supporting publication of "DESI 2024 II: Sample definitions, characteristics, and two-point clustering statistics", "DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars", and "DESI 2024 V: Analysis of the full shape of two-point clustering statistics from galaxies and quasars". (v2 - update DESI references)

A comparison of effective field theory models of redshift space galaxy power spectra for DESI 2024 and future surveys

Authors: M. Maus, Y. Lai, H. E. Noriega, S. Ramirez-Solano, A. Aviles, S. Chen, S. Fromenteau, H. Gil-Marín, C. Howlett, M. Vargas-Magaña, M. White, P. Zarrouk, J. Aguilar, S. Ahlen, O. Alves, S. Brieden, D. Brooks, E. Burtin, T. Claybaugh, S. Cole, K. Dawson, M. Icaza-Lizaola, A. de la Macorra, A. de Mattia, P. Doel , et al. (32 additional authors not shown)

Abstract: In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in th… ▽ More In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in the pipelines \texttt{velocileptors}, \texttt{PyBird}, and \texttt{Folps$ν$}. While these models have been validated independently, a detailed comparison with consistent choices has not been attempted. After briefly discussing the theoretical differences between the models we describe how to provide a more apples-to-apples comparison between them. We present the results of fitting mock spectra from the \texttt{AbacusSummit} suite of N-body simulations provided in three redshift bins to mimic the types of dark time tracers targeted by the DESI survey. We show that the theories behave similarly and give consistent constraints in both the forward-modeling and ShapeFit compressed fitting approaches. We additionally generate (noiseless) synthetic data from each pipeline to be fit by the others, varying the scale cuts in order to show that the models agree within the range of scales for which we expect 1-loop perturbation theory to be applicable. This work lays the foundation of Full-Shape analysis with DESI Y1 galaxy samples where in the tests we performed, we found no systematic error associated with the modeling of the galaxy redshift space power spectrum for this volume. △ Less

Submitted 6 June, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: 33 pages, 12 figures. Supporting publication of DESI 2024 V: Analysis of the full shape of two-point clustering statistics from galaxies and quasars

Comparing Compressed and Full-modeling Analyses with FOLPS: Implications for DESI 2024 and beyond

Authors: H. E. Noriega, A. Aviles, H. Gil-Marín, S. Ramirez-Solano, S. Fromenteau, M. Vargas-Magaña, J. Aguilar, S. Ahlen, O. Alves, S. Brieden, D. Brooks, J. L. Cervantes-Cota, S. Chen, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, A. de Mattia, P. Doel, N. Findlay, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, K. Honscheid, J. Hou , et al. (29 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) will provide unprecedented information about the large-scale structure of our Universe. In this work, we study the robustness of the theoretical modelling of the power spectrum of FOLPS, a novel effective field theory-based package for evaluating the redshift space power spectrum in the presence of massive neutrinos. We perform this validation by fit… ▽ More The Dark Energy Spectroscopic Instrument (DESI) will provide unprecedented information about the large-scale structure of our Universe. In this work, we study the robustness of the theoretical modelling of the power spectrum of FOLPS, a novel effective field theory-based package for evaluating the redshift space power spectrum in the presence of massive neutrinos. We perform this validation by fitting the AbacusSummit high-accuracy $N$-body simulations for Luminous Red Galaxies, Emission Line Galaxies and Quasar tracers, calibrated to describe DESI observations. We quantify the potential systematic error budget of FOLPS, finding that the modelling errors are fully sub-dominant for the DESI statistical precision within the studied range of scales. Additionally, we study two complementary approaches to fit and analyse the power spectrum data, one based on direct Full-Modelling fits and the other on the ShapeFit compression variables, both resulting in very good agreement in precision and accuracy. In each of these approaches, we study a set of potential systematic errors induced by several assumptions, such as the choice of template cosmology, the effect of prior choice in the nuisance parameters of the model, or the range of scales used in the analysis. Furthermore, we show how opening up the parameter space beyond the vanilla $Λ$CDM model affects the DESI observables. These studies include the addition of massive neutrinos, spatial curvature, and dark energy equation of state. We also examine how relaxing the usual Cosmic Microwave Background and Big Bang Nucleosynthesis priors on the primordial spectral index and the baryonic matter abundance, respectively, impacts the inference on the rest of the parameters of interest. This paper pathways towards performing a robust and reliable analysis of the shape of the power spectrum of DESI galaxy and quasar clustering using FOLPS. △ Less

Submitted 13 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 VII: Cosmological constraints from full-shape analyses of the two-point clustering statistics measurements, in preparation (2024)

Full Modeling and Parameter Compression Methods in configuration space for DESI 2024 and beyond

Authors: S. Ramirez-Solano, M. Icaza-Lizaola, H. E. Noriega, M. Vargas-Magaña, S. Fromenteau, A. Aviles, F. Rodriguez-Martinez, J. Aguilar, S. Ahlen, O. Alves, S. Brieden, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, Arjun Dey, B. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, H. Gil-Marín, S. Gontcho A Gontcho, K. Honscheid, C. Howlett , et al. (27 additional authors not shown)

Abstract: In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (… ▽ More In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full-Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying $n_s$ and allowing for $w_0w_a$CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within $1σ$ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models. △ Less

Submitted 16 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 KP5

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Authors: Z. Ding, A. Variu, S. Alam, Y. Yu, C. Chuang, E. Paillas, C. Garcia-Quintero, X. Chen, J. Mena-Fernández, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, A. de la Macorra, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, C. Hahn, K. Honscheid, C. Howlett, S. Juneau, R. Kehoe , et al. (22 additional authors not shown)

Abstract: Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, $N$-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by… ▽ More Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, $N$-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-$N$-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations $k_2=2k_1=0.2\hMpc$ by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the BAO scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 36 pages, 15 figures, 3 tables, to be submitted to JCAP; comments welcome

HOD-Dependent Systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

Authors: C. Garcia-Quintero, J. Mena-Fernández, A. Rocher, S. Yuan, B. Hadzhiyska, O. Alves, M. Rashkovetskyi, H. Seo, N. Padmanabhan, S. Nadathur, C. Howlett, M. Ishak, L. Medina-Varela, P. McDonald, A. J. Ross, Y. Xie, X. Chen, A. Bera, J. Aguilar, S. Ahlen, U. Andrade, S. BenZvi, D. Brooks, E. Burtin, S. Chen , et al. (51 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we focus on the robustness of the BAO… ▽ More The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we focus on the robustness of the BAO analysis against the Halo Occupation Distribution (HOD) modeling for the Emission Line Galaxy (ELG) tracer. Based on a common dark matter simulation, our analysis relies on HOD mocks tuned to early DESI data, namely the One-Percent survey data. To build the mocks, we use several HOD models for the ELG tracer as well as extensions to the baseline HOD models. Among these extensions, we consider distinct recipes for galactic conformity and assembly bias. We perform two independent analyses in the Fourier space and in the configuration space. We recover the BAO signal from two-point measurements after performing reconstruction on our mocks. Additionally, we also apply the control variates technique to reduce sample variance noise. Our BAO analysis can recover the isotropic BAO parameter $α_\text{iso}$ within 0.1\% and the Alcock Paczynski parameter $α_\text{AP}$ within 0.3\%. Overall, we find that our systematic error due to the HOD dependence is below 0.17\%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis. △ Less

Submitted 12 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

HOD-Dependent Systematics for Luminous Red Galaxies in the DESI 2024 BAO Analysis

Authors: J. Mena-Fernández, C. Garcia-Quintero, S. Yuan, B. Hadzhiyska, O. Alves, M. Rashkovetskyi, H. Seo, N. Padmanabhan, S. Nadathur, C. Howlett, S. Alam, A. Rocher, A. J. Ross, E. Sanchez, M. Ishak, J. Aguilar, S. Ahlen, U. Andrade, S. BenZvi, D. Brooks, E. Burtin, S. Chen, X. Chen, T. Claybaugh, S. Cole , et al. (50 additional authors not shown)

Abstract: In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-pa… ▽ More In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-parameter vanilla model and two extensions to it, that we refer to as baseline and extended models. The baseline model is described by the 5 vanilla HOD parameters, an incompleteness factor and a velocity bias parameter, whereas the extended one also includes a galaxy assembly bias and a satellite profile parameter. We utilize the 25 dark matter simulations available in the AbacusSummit simulation suite at $z=$ 0.8 and generate mock catalogs for our different HOD models. To test the impact of the HOD modeling in the position of the BAO peak, we run BAO fits for all these sets of simulations and compare the best-fit BAO-scaling parameters $α_{\rm iso}$ and $α_{\rm AP}$ between every pair of HOD models. We do this for both Fourier and configuration spaces independently, using post-reconstruction measurements. We find a 3.3$σ$ detection of HOD systematic for $α_{\rm AP}$ in configuration space with an amplitude of 0.19%. For the other cases, we did not find a 3$σ$ detection, and we decided to compute a conservative estimation of the systematic using the ensemble of shifts between all pairs of HOD models. By doing this, we quote a systematic with an amplitude of 0.07% in $α_{\rm iso}$ for both Fourier and configuration spaces; and of 0.09% in $α_{\rm AP}$ for Fourier space. △ Less

Submitted 5 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 36 pages, 9 figures. Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data

Authors: M. Rashkovetskyi, D. Forero-Sánchez, A. de Mattia, D. J. Eisenstein, N. Padmanabhan, H. Seo, A. J. Ross, J. Aguilar, S. Ahlen, O. Alves, U. Andrade, D. Brooks, E. Burtin, T. Claybaugh, S. Cole, A. de la Macorra, Z. Ding, P. Doel, K. Fanning, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, C. Garcia-Quintero, H. Gil-Marín, S. Gontcho A Gontcho , et al. (34 additional authors not shown)

Abstract: We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 20… ▽ More We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation used in the mocks. Importantly, we find a close agreement (<~ 5% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, in particular, relevant to the broad current and future analyses of DESI data. △ Less

Submitted 5 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars (arXiv:2404.03000). 29 pages, 4 figures. Prepared for submission to JCAP. Code available at https://github.com/oliverphilcox/RascalC and https://github.com/misharash/RascalC-scripts/tree/DESI2024. Data points from the plots available at https://zenodo.org/doi/10.5281/zenodo.10895161

Production of Alternate Realizations of DESI Fiber Assignment for Unbiased Clustering Measurement in Data and Simulations

Authors: J. Lasker, A. Carnero Rosell, A. D. Myers, A. J. Ross, D. Bianchi, M. M. S Hanif, R. Kehoe, A. de Mattia, L. Napolitano, W. J. Percival, R. Staten, J. Aguilar, S. Ahlen, L. Bigwood, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, Z. Ding, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez , et al. (30 additional authors not shown)

Abstract: A critical requirement of spectroscopic large scale structure analyses is correcting for selection of which galaxies to observe from an isotropic target list. This selection is often limited by the hardware used to perform the survey which will impose angular constraints of simultaneously observable targets, requiring multiple passes to observe all of them. In SDSS this manifested solely as the co… ▽ More A critical requirement of spectroscopic large scale structure analyses is correcting for selection of which galaxies to observe from an isotropic target list. This selection is often limited by the hardware used to perform the survey which will impose angular constraints of simultaneously observable targets, requiring multiple passes to observe all of them. In SDSS this manifested solely as the collision of physical fibers and plugs placed in plates. In DESI, there is the additional constraint of the robotic positioner which controls each fiber being limited to a finite patrol radius. A number of approximate methods have previously been proposed to correct the galaxy clustering statistics for these effects, but these generally fail on small scales. To accurately correct the clustering we need to upweight pairs of galaxies based on the inverse probability that those pairs would be observed (Bianchi \& Percival 2017). This paper details an implementation of that method to correct the Dark Energy Spectroscopic Instrument (DESI) survey for incompleteness. To calculate the required probabilities, we need a set of alternate realizations of DESI where we vary the relative priority of otherwise identical targets. These realizations take the form of alternate Merged Target Ledgers (AMTL), the files that link DESI observations and targets. We present the method used to generate these alternate realizations and how they are tracked forward in time using the real observational record and hardware status, propagating the survey as though the alternate orderings had been adopted. We detail the first applications of this method to the DESI One-Percent Survey (SV3) and the DESI year 1 data. We include evaluations of the pipeline outputs, estimation of survey completeness from this and other methods, and validation of the method using mock galaxy catalogs. △ Less

Submitted 22 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 30 pages, 21 figures

Optimal Reconstruction of Baryon Acoustic Oscillations for DESI 2024

Authors: E. Paillas, Z. Ding, X. Chen, H. Seo, N. Padmanabhan, A. de Mattia, A. J. Ross, S. Nadathur, C. Howlett, J. Aguilar, S. Ahlen, O. Alves, U. Andrade, D. Brooks, E. Buckley-Geer, E. Burtin, S. Chen, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, Arjun Dey, P. Doel, K. Fanning, S. Ferraro , et al. (51 additional authors not shown)

Abstract: Baryon acoustic oscillations (BAO) provide a robust standard ruler to measure the expansion history of the Universe through galaxy clustering. Density-field reconstruction is now a widely adopted procedure for increasing the precision and accuracy of the BAO detection. With the goal of finding the optimal reconstruction settings to be used in the DESI 2024 galaxy BAO analysis, we assess the sensit… ▽ More Baryon acoustic oscillations (BAO) provide a robust standard ruler to measure the expansion history of the Universe through galaxy clustering. Density-field reconstruction is now a widely adopted procedure for increasing the precision and accuracy of the BAO detection. With the goal of finding the optimal reconstruction settings to be used in the DESI 2024 galaxy BAO analysis, we assess the sensitivity of the post-reconstruction BAO constraints to different choices in our analysis configuration, performing tests on blinded data from the first year of DESI observations (DR1), as well as on mocks that mimic the expected clustering and selection properties of the DESI DR1 target samples. Overall, we find that BAO constraints remain robust against multiple aspects in the reconstruction process, including the choice of smoothing scale, treatment of redshift-space distortions, fiber assignment incompleteness, and parameterizations of the BAO model. We also present a series of tests that DESI followed in order to assess the maturity of the end-to-end galaxy BAO pipeline before the unblinding of the large-scale structure catalogs. △ Less

Submitted 14 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

Validation of the DESI 2024 Ly$α$ forest BAO analysis using synthetic datasets

Authors: Andrei Cuceu, Hiram K. Herrera-Alcantar, Calum Gordon, Paul Martini, Julien Guy, Andreu Font-Ribera, Alma X. Gonzalez-Morales, M. Abdul Karim, J. Aguilar, S. Ahlen, E. Armengaud, A. Bault, D. Brooks, T. Claybaugh, A. de la Macorra, P. Doel, K. Fanning, S. Ferraro, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, K. Honscheid, C. Howlett, N. G. Karaçaylı , et al. (34 additional authors not shown)

Abstract: The first year of data from the Dark Energy Spectroscopic Instrument (DESI) contains the largest set of Lyman-$α$ (Ly$α$) forest spectra ever observed. This data, collected in the DESI Data Release 1 (DR1) sample, has been used to measure the Baryon Acoustic Oscillation (BAO) feature at redshift $z=2.33$. In this work, we use a set of 150 synthetic realizations of DESI DR1 to validate the DESI 202… ▽ More The first year of data from the Dark Energy Spectroscopic Instrument (DESI) contains the largest set of Lyman-$α$ (Ly$α$) forest spectra ever observed. This data, collected in the DESI Data Release 1 (DR1) sample, has been used to measure the Baryon Acoustic Oscillation (BAO) feature at redshift $z=2.33$. In this work, we use a set of 150 synthetic realizations of DESI DR1 to validate the DESI 2024 Ly$α$ forest BAO measurement. The synthetic data sets are based on Gaussian random fields using the log-normal approximation. We produce realistic synthetic DESI spectra that include all major contaminants affecting the Ly$α$ forest. The synthetic data sets span a redshift range $1.8<z<3.8$, and are analysed using the same framework and pipeline used for the DESI 2024 Ly$α$ forest BAO measurement. To measure BAO, we use both the Ly$α$ auto-correlation and its cross-correlation with quasar positions. We use the mean of correlation functions from the set of DESI DR1 realizations to show that our model is able to recover unbiased measurements of the BAO position. We also fit each mock individually and study the population of BAO fits in order to validate BAO uncertainties and test our method for estimating the covariance matrix of the Ly$α$ forest correlation functions. Finally, we discuss the implications of our results and identify the needs for the next generation of Ly$α$ forest synthetic data sets, with the top priority being to simulate the effect of BAO broadening due to non-linear evolution. △ Less

Submitted 5 May, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: Supporting publication of DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest

DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, B. Bahr-Kalus, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, A. Bera, F. Beutler, D. Bianchi, C. Blake, R. Blum , et al. (178 additional authors not shown)

Abstract: We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ 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 s… ▽ More We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ 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<z<4.2$. DESI BAO data alone are consistent with the standard flat $Λ$CDM cosmological model with a matter density $Ω_\mathrm{m}=0.295\pm 0.015$. Paired with a BBN prior and the robustly measured acoustic angular scale from the CMB, DESI requires $H_0=(68.52\pm0.62)$ km/s/Mpc. In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find $Ω_\mathrm{m}=0.307\pm 0.005$ and $H_0=(67.97\pm0.38)$ km/s/Mpc. Extending the baseline model with a constant dark energy equation of state parameter $w$, DESI BAO alone require $w=-0.99^{+0.15}_{-0.13}$. In models with a time-varying dark energy equation of state parametrized by $w_0$ and $w_a$, combinations of DESI with CMB or with SN~Ia individually prefer $w_0>-1$ and $w_a<0$. This preference is 2.6$σ$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $Λ$CDM model at the $2.5σ$, $3.5σ$ or $3.9σ$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $Λ$CDM model with the sum of neutrino mass $\sum m_ν$ free, combining the DESI and CMB data yields an upper limit $\sum m_ν< 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_ν>0$ $(\sum m_ν>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $Λ$CDM. [Abridged.] △ Less

Submitted 24 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 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). Typos corrected and a new figure and discussion added to Appendix A

DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Bautista, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden , et al. (174 additional authors not shown)

Abstract: We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ 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… ▽ More We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ 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. △ Less

Submitted 12 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 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)

DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden, A. Brodzeller , et al. (171 additional authors not shown)

Abstract: 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<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 qu… ▽ More 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<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 quasars with 0.8<z<2.1, over a ~7,500 square degree footprint. The analysis was blinded at the catalog-level to avoid confirmation bias. All fiducial choices of the BAO fitting and reconstruction methodology, as well as the size of the systematic errors, were determined on the basis of the tests with mock catalogs and the blinded data catalogs. We present several improvements to the BAO analysis pipeline, including enhancing the BAO fitting and reconstruction methods in a more physically-motivated direction, and also present results using combinations of tracers. We present a re-analysis of SDSS BOSS and eBOSS results applying the improved DESI methodology and find scatter consistent with the level of the quoted SDSS theoretical systematic uncertainties. With the total effective survey volume of ~ 18 Gpc$^3$, the combined precision of the BAO measurements across the six different redshift bins is ~0.52%, marking a 1.2-fold improvement over the previous state-of-the-art results using only first-year data. We detect the BAO in all of these six redshift bins. The highest significance of BAO detection is $9.1σ$ at the effective redshift of 0.93, with a constraint of 0.86% placed on the BAO scale. We find our measurements are systematically larger than the prediction of Planck-2018 LCDM model at z<0.8. We translate the results into transverse comoving distance and radial Hubble distance measurements, which are used to constrain cosmological models in our companion paper [abridged]. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 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)

Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions

Authors: Z. Brown, R. Demina, A. G. Adame, S. Avila, E. Chaussidon, S. Yuan, V. Gonzalez-Perez, J. García-Bellido, J. Aguilar, S. Ahlen, R. Blum, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, B. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, K. Honscheid, C. Howlett, S. Juneau, R. Kehoe , et al. (25 additional authors not shown)

Abstract: Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by fNL, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spe… ▽ More Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by fNL, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure fNL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter halos. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of fNL. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our model, and forecast the ability to constrainfNL to a precision of sigma(fNL) = 22 with one year of DESI survey data. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Baryon Acoustic Oscillation Theory and Modelling Systematics for the DESI 2024 results

Authors: Shi-Fan Chen, Cullan Howlett, Martin White, Patrick McDonald, Ashley J. Ross, Hee-Jong Seo, Nikhil Padmanabhan, J. Aguilar, S. Ahlen, S. Alam, O. Alves, R. Blum, D. Brooks, X. Chen, S. Cole, T. M. Davis, K. Dawson, A. de la Macorra, Arjun Dey, Z. Ding, P. Doel, S. Ferraro, A. Font-Ribera, D. Forero-Sánchez, J. E. Forero-Romero , et al. (33 additional authors not shown)

Abstract: This paper provides a comprehensive overview of how fitting of Baryon Acoustic Oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument's (DESI) 2024 results using its DR1 dataset, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential bias… ▽ More This paper provides a comprehensive overview of how fitting of Baryon Acoustic Oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument's (DESI) 2024 results using its DR1 dataset, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential biases in measurements of the isotropic ($α_{\mathrm{iso}}$) and anisotropic ($α_{\mathrm{ap}}$) BAO distance scales, and how these can be effectively removed with an appropriate choice of reconstruction algorithm. We then demonstrate how theory leads to a clear choice for how to model the BAO and develop, implement and validate a new model for the remaining smooth-broadband (i.e., without BAO) component of the galaxy clustering. Finally, we explore the impact of all remaining modelling choices on the BAO constraints from DESI using a suite of high-precision simulations, arriving at a set of best-practices for DESI BAO fits, and an associated theory and modelling systematic error. Overall, our results demonstrate the remarkable robustness of the BAO to all our modelling choices and motivate a combined theory and modelling systematic error contribution to the post-reconstruction DESI BAO measurements of no more than $0.1\%$ ($0.2\%$) for its isotropic (anisotropic) distance measurements. We expect the theory and best-practices laid out to here to be applicable to other BAO experiments in the era of DESI and beyond. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 29 pages, 18 figures, 1 table, submitted to MNRAS

An effective description of Laniakea: impact on cosmology and the local determination of the Hubble constant

Authors: L. Giani, C. Howlett, K. Said, T. Davis, S. Vagnozzi

Abstract: We propose an effective model to describe the bias induced on cosmological observables by Laniakea, the gravitational supercluster hosting the Milky Way, which was defined using peculiar velocity data from Cosmicflows-4 (CF4). The structure is well described by an ellipsoidal shape exhibiting triaxial expansion, reasonably approximated by a constant expansion rate along the principal axes. Our bes… ▽ More We propose an effective model to describe the bias induced on cosmological observables by Laniakea, the gravitational supercluster hosting the Milky Way, which was defined using peculiar velocity data from Cosmicflows-4 (CF4). The structure is well described by an ellipsoidal shape exhibiting triaxial expansion, reasonably approximated by a constant expansion rate along the principal axes. Our best fits suggest that the ellipsoid, after subtracting the background expansion, contracts along the two smaller axes and expands along the longest one, predicting an average expansion of $\sim -1.1 ~\rm{km}/\rm{s}/\rm{Mpc}$. The different expansion rates within the region, relative to the mean cosmological expansion, induce line-of-sight-dependent corrections in the computation of luminosity distances. We apply these corrections to two low-redshift datasets: the Pantheon+ catalog of type Ia Supernovae (SN~Ia), and 63 measurements of Surface Brightness Fluctuations (SBF) of early-type massive galaxies from the MASSIVE survey. We find corrections on the distances of order $\sim 2-3\%$, resulting in a shift in the inferred best-fit values of the Hubble constant $H_0$ of order $ΔH_0^{\rm{SN~Ia}}\approx 0.5 ~\rm{km}/\rm{s}/\rm{Mpc}$ and $ΔH_0^{\rm{SBF}}\approx 1.1 ~\rm{km}/\rm{s}/\rm{Mpc}$, seemingly worsening the Hubble tension. △ Less

Submitted 1 February, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: Updated to match the published version. Eighteen pages + references, 13 figures, 1 appendix, title and abstract slightly changed. Comments are welcome!

Journal ref: JCAP 2401 (2024) 071

Ho'oleilana: An Individual Baryon Acoustic Oscillation?

Authors: R. Brent Tully, Cullan Howlett, Daniel Pomarede

Abstract: Theory of the physics of the early hot universe leads to a prediction of baryon acoustic oscillations that has received confirmation from the pair-wise separations of galaxies in samples of hundreds of thousands of objects. Evidence is presented here for the discovery of a remarkably strong individual contribution to the baryon acoustic oscillation (BAO) signal at z=0.068, an entity that is given… ▽ More Theory of the physics of the early hot universe leads to a prediction of baryon acoustic oscillations that has received confirmation from the pair-wise separations of galaxies in samples of hundreds of thousands of objects. Evidence is presented here for the discovery of a remarkably strong individual contribution to the baryon acoustic oscillation (BAO) signal at z=0.068, an entity that is given the name Ho'oleilana. The radius of the 3D structure is 155/h_{75} Mpc. At its core is the Bootes supercluster. The Sloan Great Wall, CfA Great Wall, and Hercules complex all lie within the BAO shell. The interpretation of Ho'oleilana as a BAO structure with our preferred analysis implies a value of the Hubble constant of 76.9+8.2-4.8 km/s/Mpc. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: One video and two interactive models. 14 pages, 14 figures

Journal ref: 2023, Astrophysical Journal, 954, 169

The Galaxy Number Density Profile of Haloes

Authors: Fei Qin, David Parkinson, Adam R. H. Stevens, Cullan Howlett

Abstract: More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for t… ▽ More More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro-Frenk-White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterised galaxy number density profile model of haloes is fit and tested using the DARKSAGE semi-analytic model of galaxy formation. We find that our galaxy number density profile model of haloes can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the DARKSAGE simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this profile model. We use the SDSS DR10 galaxy group catalogue to validate this galaxy number density profile model of haloes. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function. △ Less

Submitted 25 September, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 13 pages. 10 figures. Appear on ApJ

Evaluating bulk flow estimators for CosmicFlows-4 measurements

Authors: A. M. Whitford, C. Howlett, T. M. Davis

Abstract: For over a decade there have been contradictory claims in the literature about whether the local bulk flow motion of galaxies is consistent or in tension with the $Λ$CDM model. While it has become evident that systematics affect bulk flow measurements, systematics in the estimators have not been widely investigated. In this work, we thoroughly evaluate the performance of four estimator variants, i… ▽ More For over a decade there have been contradictory claims in the literature about whether the local bulk flow motion of galaxies is consistent or in tension with the $Λ$CDM model. While it has become evident that systematics affect bulk flow measurements, systematics in the estimators have not been widely investigated. In this work, we thoroughly evaluate the performance of four estimator variants, including the Kaiser maximum likelihood estimator (MLE) and the minimum variance estimator (MVE). We find that these estimators are unbiased, however their precision may be strongly correlated with the survey geometry. Small biases in the estimators can be present leading to underestimated bulk flows, which we suspect are due to the presence of non-linear peculiar velocities. The uncertainty assigned to the bulk flows from these estimators is typically underestimated, which leads to an overestimate of the tension with $Λ$CDM. We estimate the bulk flow for the CosmicFlows-4 data and use mocks to ensure the uncertainties are appropriately accounted for. Using the MLE we find a bulk flow amplitude of $408\pm165 \mathrm{km s}^{-1}$ at a depth of $49\, \mathrm{Mpc} h^{-1}$, in reasonable agreement with $Λ$CDM. However using the MVE which can probe greater effective depths, we find an amplitude of $428\pm108 \mathrm{km s}^{-1}$ at a depth of $173\, \mathrm{Mpc} h^{-1}$, in tension with the model, having only a 0.11% probability of obtaining a larger $χ^2$. These measurements appear directed towards the Great Attractor region where more data may be needed to resolve tensions. △ Less

Submitted 10 September, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

The Early Data Release of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (240 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 15 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 43 pages, 7 figures, 17 tables, submitted to AJ, DESI EDR references added

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (239 additional authors not shown)

Abstract: 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… ▽ More 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<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$. △ Less

Submitted 12 January, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 42 pages, 18 figures, accepted by AJ

Faster cosmological analysis with power spectrum without simulations

Authors: Yan Lai, Cullan Howlett, Tamara M. Davis

Abstract: $… ▽ More $ $Future surveys could obtain tighter constraints on the cosmological parameters with the galaxy power spectrum than with the Cosmic Microwave Background. However, the inclusion of multiple overlapping tracers, redshift bins, and more non-linear scales means that generating the necessary ensemble of simulations for model-fitting presents a computational burden. In this work, we combine full-shape fitting of galaxy power spectra, analytical covariance matrix estimates, the MOPED (Massively Optimised Parameter Estimation and Data compression) method, and the Taylor expansion interpolation of the power spectrum for the first time to constrain the cosmological parameters directly from a state-of-the-art set of galaxy clustering measurements. We find it takes less than a day to compute the analytical covariance while it takes several months to calculate the simulated ones. Combining MOPED with the Taylor expansion interpolation of the power spectrum, we can constrain the cosmological parameters in just a few hours instead of a few days. We also find that even without a priori knowledge of the best-fit cosmological or galaxy bias parameters, the analytical covariance matrix with the MOPED compression still gives consistent cosmological constraints to within 0.1$σ$ after two iterations. Therefore, the pipeline we have developed here can significantly speed up the analysis for future surveys such as DESI and Euclid. △ Less

Submitted 4 April, 2024; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: Revised version for MNRAS submission. 5 figures, 6 tables, and 13 pages. Comments are welcome

Target Selection for the DESI Peculiar Velocity Survey

Authors: Christoph Saulder, Cullan Howlett, Kelly A. Douglass, Khaled Said, Segev BenZvi, Steven Ahlen, Greg Aldering, Stephen Bailey, David Brooks, Tamara Davis, Axel de la Macorra, Arjun Dey, Andreu Font-Ribera, Jaime E. Forero-Romero, Satya Gontcho A Gontcho, Klaus Honscheid, Alex G. Kim, Theodore Kisner, Anthony Kremin, Martin Landriau, Michael E. Levi, John Lucey, Aaron M. Meisner, Ramon Miquel, John Moustakas , et al. (13 additional authors not shown)

Abstract: We describe the target selection and characteristics of the DESI Peculiar Velocity Survey, the largest survey of peculiar velocities (PVs) using both the fundamental plane (FP) and the Tully-Fisher (TF) relationship planned to date. We detail how we identify suitable early-type galaxies (ETGs) for the FP and suitable late-type galaxies (LTGs) for the TF relation using the photometric data provided… ▽ More We describe the target selection and characteristics of the DESI Peculiar Velocity Survey, the largest survey of peculiar velocities (PVs) using both the fundamental plane (FP) and the Tully-Fisher (TF) relationship planned to date. We detail how we identify suitable early-type galaxies (ETGs) for the FP and suitable late-type galaxies (LTGs) for the TF relation using the photometric data provided by the DESI Legacy Imaging Survey DR9. Subsequently, we provide targets for 373 533 ETGs and 118 637 LTGs within the DESI 5-year footprint. We validate these photometric selections using existing morphological classifications. Furthermore, we demonstrate using survey validation data that DESI is able to measure the spectroscopic properties to sufficient precision to obtain PVs for our targets. Based on realistic DESI fiber assignment simulations and spectroscopic success rates, we predict the final DESI PV Survey will obtain $\sim$133 000 FP-based and $\sim$53 000 TF-based PV measurements over an area of 14 000 $\mathrm{deg^{2}}$. We forecast the ability of using these data to measure the clustering of galaxy positions and PVs from the combined DESI PV and Bright Galaxy Surveys (BGS), which allows for cancellation of cosmic variance at low redshifts. With these forecasts, we anticipate a $4\%$ statistical measurement on the growth rate of structure at $z<0.15$. This is over two times better than achievable with redshifts from the BGS alone. The combined DESI PV and BGS will enable the most precise tests to date of the time and scale dependence of large-scale structure growth at $z<0.15$. △ Less

Submitted 28 July, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 20 pages, 15 figures, 14 tables; accepted in MNRAS

Can Einstein (rings) surf Gravitational Waves?

Authors: Leonardo Giani, Cullan Howlett, Tamara M. Davis

Abstract: How does the appearance of a strongly lensed system change if a gravitational wave is produced by the lens? In this work we address this question by considering a supermassive black hole binary at the center of the lens emitting gravitational waves propagating either colinearly or orthogonally to the line of sight. Specializing to an Einstein ring configuration (where the source, the lens and the… ▽ More How does the appearance of a strongly lensed system change if a gravitational wave is produced by the lens? In this work we address this question by considering a supermassive black hole binary at the center of the lens emitting gravitational waves propagating either colinearly or orthogonally to the line of sight. Specializing to an Einstein ring configuration (where the source, the lens and the observer are aligned), we show that the gravitational wave induces changes on the ring's angular size and on the optical path of photons. The changes are the same for a given pair of antipodal points on the ring, but maximally different for any pair separated by $90^{\circ}$. For realistic lenses and binaries, we find that the change in the angular size of the Einstein ring is dozens of orders of magnitude smaller than the precision of current experiments. On the other hand, the difference in the optical path induced on a photon by a gravitational wave propagating \textit{orthogonally} to the line of sight triggers, at peak strain, time delays in the range $\sim 0.01 - 1$ seconds, making the chance of their detection (and thus the use of Einstein rings as gravitational wave detectors) less hopeless. △ Less

Submitted 24 July, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: v2. Version accepted for publication in the Open Journal of Astrophysics. 8 pages, four figures. Comments are welcome!

Exploring binary black hole mergers and host galaxies with Shark and COMPAS

Authors: Liana Rauf, Cullan Howlett, Tamara M. Davis, Claudia D. P. Lagos

Abstract: We explore the connection between the gravitational wave (GW) merger rates of stellar-mass binary black holes (BBH) and galaxy properties. We do this by generating populations of stars using the binary population synthesis code COMPAS and evolving them in galaxies from the semi-analytic galaxy formation model Shark, to determine the number of mergers occurring in each simulation time-step. We find… ▽ More We explore the connection between the gravitational wave (GW) merger rates of stellar-mass binary black holes (BBH) and galaxy properties. We do this by generating populations of stars using the binary population synthesis code COMPAS and evolving them in galaxies from the semi-analytic galaxy formation model Shark, to determine the number of mergers occurring in each simulation time-step. We find that metal-rich and massive galaxies with star formation rate (SFR) greater than $1M_{\odot}/ \rm yr$ are 10 times more likely to have GW events compared to younger, less massive and metal poor galaxies. Our simulation with the default input parameters predicts a higher local merger rate density compared to the third GW transient catalogue (GWTC-3) prediction from LIGO, VIRGO and KAGRA, due to short coalescence times, low metallicities and a high SFR at low redshift in the simulation, which produces more BBHs that merge within the age of the Universe compared to observations. We identify alternate remnant mass models that more accurately reproduce the volumetric rate and provide updated fits to the merger rate as a function of redshift. We then investigate the relative fraction of GW events in our simulation that are in observable host galaxies from upcoming galaxy surveys, determining which of those are ideal for tracing host galaxies with high merger rates. The implications of this work can be utilised for constraining stellar evolution models, better informing follow-up programs, and placing informative priors on host galaxies when measuring cosmological parameters such as the Hubble constant. △ Less

Submitted 12 June, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 20 pages, 16 figures - Accepted by MNRAS

Cross-correlating radial peculiar velocities and CMB lensing convergence

Authors: Leonardo Giani, Cullan Howlett, Rossana Ruggeri, Federico Bianchini, Khaled Said, Tamara M. Davis

Abstract: We study, for the first time, the cross correlation between the angular distribution of radial peculiar velocities (PV) and the lensing convergence of cosmic microwave background (CMB) photons. We derive theoretical expectations for the signal and its covariance and assess its detectability with existing and forthcoming surveys. We find that such cross-correlations are expected to improve constrai… ▽ More We study, for the first time, the cross correlation between the angular distribution of radial peculiar velocities (PV) and the lensing convergence of cosmic microwave background (CMB) photons. We derive theoretical expectations for the signal and its covariance and assess its detectability with existing and forthcoming surveys. We find that such cross-correlations are expected to improve constraints on different gravitational models by partially breaking degeneracies with the matter density. We identify in the distance-scaling dispersion of the peculiar velocities the most relevant source of noise in the cross correlation. For this reason, we also study how the above picture changes assuming a redshift-independent scatter for the PV, obtained for example using a reconstruction technique. Our results show that the cross correlation might be detected in the near future combining PV measurements from DESI and the convergence map from CMB-S4. Using realistic direct PV measurements we predict a cumulative signal-to-noise ratio of approximately $3.8 σ$ using data on angular scales $3 \leq \ell \leq 200$. For an idealized reconstructed peculiar velocity map extending up to redshift $z=0.15$ and a smoothing scale of $4$ Mpc $h^{-1}$ we predict a cumulative signal-to-noise ratio of approximately $ 27 σ$ from angular scales $3 \leq \ell \leq200 $. We conclude that currently reconstructed peculiar velocities have more constraining power than directly observed ones, even though they are more cosmological-model dependent. △ Less

Submitted 24 July, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: v2, minor changes to match the published version. two figures added. Comments are welcome!

WALLABY Pilot Survey: Public release of HI kinematic models for more than 100 galaxies from phase 1 of ASKAP pilot observations

Authors: N. Deg, K. Spekkens, T. Westmeier, T. N. Reynolds, P. Venkataraman, S. Goliath, A. X. Shen, R. Halloran, A. Bosma, B. Catinella, W. J. G. de Blok, H. Dénes, E. M. Di Teodoro, A. Elagali, B. -Q. For, C. Howlett, G. I. G. Józsa, P. Kamphuis, D. Kleiner, B. Koribalski, K. Lee-Waddell, F. Lelli, X. Lin, C. Murugeshan, S. Oh , et al. (7 additional authors not shown)

Abstract: We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I HI kinematic models. This first data release consists of HI observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique HI detections in t… ▽ More We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I HI kinematic models. This first data release consists of HI observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique HI detections in these fields. The modelling method adopted here - which we call the WALLABY Kinematic Analysis Proto-Pipeline (WKAPP) and for which the corresponding scripts are also publicly available - consists of combining results from the homogeneous application of the FAT and 3DBAROLO algorithms to the subset of 209 detections with sufficient resolution and S/N in order to generate optimized model parameters and uncertainties. The 109 models presented here tend to be gas rich detections resolved by at least 3-4 synthesized beams across their major axes, but there is no obvious environmental bias in the modelling. The data release described here is the first step towards the derivation of similar products for thousands of spatially-resolved WALLABY detections via a dedicated kinematic pipeline. Such a large publicly available and homogeneously analyzed dataset will be a powerful legacy product that that will enable a wide range of scientific studies. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to PASA

WALLABY Pilot Survey: Public release of HI data for almost 600 galaxies from phase 1 of ASKAP pilot observations

Authors: T. Westmeier, N. Deg, K. Spekkens, T. N. Reynolds, A. X. Shen, S. Gaudet, S. Goliath, M. T. Huynh, P. Venkataraman, X. Lin, T. O'Beirne, B. Catinella, L. Cortese, H. Dénes, A. Elagali, B. -Q. For, G. I. G. Józsa, C. Howlett, J. M. van der Hulst, R. J. Jurek, P. Kamphuis, V. A. Kilborn, D. Kleiner, B. S. Koribalski, K. Lee-Waddell , et al. (27 additional authors not shown)

Abstract: We present WALLABY pilot data release 1, the first public release of HI pilot survey data from the Wide-field ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. Phase 1 of the WALLABY pilot survey targeted three $60~{\rm deg}^2$ regions on the sky in the direction of the Hydra and Norma galaxy clusters and the NGC 4636 galaxy group, covering the… ▽ More We present WALLABY pilot data release 1, the first public release of HI pilot survey data from the Wide-field ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. Phase 1 of the WALLABY pilot survey targeted three $60~{\rm deg}^2$ regions on the sky in the direction of the Hydra and Norma galaxy clusters and the NGC 4636 galaxy group, covering the redshift range of z < 0.08. The source catalogue, images and spectra of nearly 600 extragalactic HI detections and kinematic models for 109 spatially resolved galaxies are available. As the pilot survey targeted regions containing nearby group and cluster environments, the median redshift of the sample of z ~ 0.014 is relatively low compared to the full WALLABY survey. The median galaxy HI mass is $2.3 \times 10^{9}~M_{\odot}$. The target noise level of 1.6 mJy per $30''$ beam and 18.5 kHz channel translates into a $5σ$ HI mass sensitivity for point sources of about $5.2 \times 10^{8} \, (D_{\rm L} / \mathrm{100~Mpc})^{2} \, M_{\odot}$ across 50 spectral channels (~200 km/s) and a $5σ$ HI column density sensitivity of about $8.6 \times 10^{19} \, (1 + z)^{4}~\mathrm{cm}^{-2}$ across 5 channels (~20 km/s) for emission filling the $30''$ beam. As expected for a pilot survey, several technical issues and artefacts are still affecting the data quality. Most notably, there are systematic flux errors of up to several 10% caused by uncertainties about the exact size and shape of each of the primary beams as well as the presence of sidelobes due to the finite deconvolution threshold. In addition, artefacts such as residual continuum emission and bandpass ripples have affected some of the data. The pilot survey has been highly successful in uncovering such technical problems, most of which are expected to be addressed and rectified before the start of the full WALLABY survey. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: 24 pages, 13 figures, 5 tables, accepted for publication in PASA

WALLABY Pre-Pilot and Pilot Survey: the Tully Fisher Relation in Eridanus, Hydra, Norma and NGC4636 fields

Authors: Hélène M. Courtois, Khaled Said, Jeremy Mould, T. H. Jarrett, Daniel Pomarède, Tobias Westmeier, Lister Staveley-Smith, Alexandra Dupuy, Tao Hong, Daniel Guinet, Cullan Howlett, Nathan Deg, Bi-Qing For, Dane Kleiner, Bärbel Koribalski, Karen Lee-Waddell, Jonghwan Rhee, Kristine Spekkens, Jing Wang, O. I. Wong, Frank Bigiel, Albert Bosma, Matthew Colless, Tamara Davis, Benne Holwerda , et al. (6 additional authors not shown)

Abstract: The WALLABY pilot survey has been conducted using the Australian SKA Pathfinder (ASKAP). The integrated 21-cm HI line spectra are formed in a very different manner compared to usual single-dish spectra Tully-Fisher measurements. It is thus extremely important to ensure that slight differences (e.g. biases due to missing flux) are quantified and understood in order to maximise the use of the large… ▽ More The WALLABY pilot survey has been conducted using the Australian SKA Pathfinder (ASKAP). The integrated 21-cm HI line spectra are formed in a very different manner compared to usual single-dish spectra Tully-Fisher measurements. It is thus extremely important to ensure that slight differences (e.g. biases due to missing flux) are quantified and understood in order to maximise the use of the large amount of data becoming available soon. This article is based on four fields for which the data are scientifically interesting by themselves. The pilot data discussed here consist of 614 galaxy spectra at a rest wavelength of 21cm. Of these spectra, 472 are of high enough quality to be used to potentially derive distances using the Tully-Fisher relation. We further restrict the sample to the 251 galaxies whose inclination is sufficiently close to edge-on. For these, we derive Tully-Fisher distances using the deprojected WALLABY velocity widths combined with infrared (WISE W1) magnitudes. The resulting Tully-Fisher distances for the Eridanus, Hydra, Norma and NGC 4636 clusters are 21.5, 53.5, 69.4 and 23.0 Mpc respectively, with uncertainties of 5--10\%, which are better or equivalent to the ones obtained in studies using data obtained with giant single dish telescopes. The pilot survey data show the benefits of WALLABY over previous giant single-dish telescope surveys. WALLABY is expected to detect around half a million galaxies with a mean redshift of $z = 0.05 (200 Mpc)$. This study suggests that about 200,000 Tully-Fisher distances might result from the survey. △ Less

Submitted 22 October, 2022; originally announced October 2022.

Comments: 21 pages, 12 Figures, 6 Tables, accepted for publication in the MNRAS

Cosmicflows-4

Authors: R. Brent Tully, Ehsan Kourkchi, Hélène M. Courtois, Gagandeep S. Anand, John P. Blakeslee, Dillon Brout, Thomas de Jaeger, Alexandra Dupuy, Daniel Guinet, Cullan Howlett, Joseph B. Jensen, Daniel Pomarède, Luca Rizzi, David Rubin, Khaled Said, Daniel Scolnic, Benjamin E. Stahl

Abstract: With Cosmicflows-4, distances are compiled for 55,877 galaxies gathered into 38,065 groups. Eight methodologies are employed, with the largest numbers coming from the correlations between the photometric and kinematic properties of spiral galaxies (TF) and elliptical galaxies (FP). Supernovae that arise from degenerate progenitors (type Ia Sne) are an important overlapping component. Smaller contr… ▽ More With Cosmicflows-4, distances are compiled for 55,877 galaxies gathered into 38,065 groups. Eight methodologies are employed, with the largest numbers coming from the correlations between the photometric and kinematic properties of spiral galaxies (TF) and elliptical galaxies (FP). Supernovae that arise from degenerate progenitors (type Ia Sne) are an important overlapping component. Smaller contributions come from distance estimates from the surface brightness fluctuations of elliptical galaxies and the luminosities and expansion rates of core collapse supernovae (SNII). Cepheid period-luminosity relation and tip of the red giant branch observations founded on local stellar parallax measurements along with the geometric maser distance to NGC 4258 provide the absolute scaling of distances. The assembly of galaxies into groups is an important feature of the study in facilitating overlaps between methodologies. Merging between multiple contributions within a methodology and between methodologies is carried out with Bayesian Markov chain Monte Carlo procedures. The final assembly of distances is compatible with a value of the Hubble constant of $H_0=74.6$ km s$^{-1}$ Mpc$^{-1}$ with the small statistical error of $\pm 0.8$ km s$^{-1}$ Mpc$^{-1}$ but a large potential systematic error of ~3 km s$^{-1}$ Mpc$^{-1}$. Peculiar velocities can be inferred from the measured distances. The interpretation of the field of peculiar velocities is complex because of large errors on individual components and invites analyses beyond the scope of this study. △ Less

Submitted 28 December, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 38 pages, 24 figures. catalogs available at edd.ifa.hawaii.edu. Revised version, accepted to ApJ

Using peculiar velocity surveys to constrain the growth rate of structure with the wide-angle effect

Authors: Yan Lai, Cullan Howlett, Tamara M. Davis

Abstract: Amongst the most popular explanations for dark energy are modified theories of gravity. The galaxy overdensity and peculiar velocity fields help us to constrain the growth rate of structure and distinguish different models of gravity. We introduce an improved method for constraining the growth rate of structure with the galaxy overdensity and peculiar velocity fields. This method reduces the model… ▽ More Amongst the most popular explanations for dark energy are modified theories of gravity. The galaxy overdensity and peculiar velocity fields help us to constrain the growth rate of structure and distinguish different models of gravity. We introduce an improved method for constraining the growth rate of structure with the galaxy overdensity and peculiar velocity fields. This method reduces the modelling systematic error by accounting for the wide-angle effect and the zero-point calibration uncertainty during the modelling process. We also speed up the posterior sampling by around 30 times by first calculating the likelihood at a small number of fiducial points and then interpolating the likelihood values during MCMC sampling. We test the new method on mocks and we find it is able to recover the fiducial growth rate of structure. We applied our new method to the SDSS PV catalogue, which is the largest single peculiar velocity catalogue to date. Our constraint on the growth rate of structure is \(fσ_8 = 0.405_{-0.071}^{+0.076}\) (stat) \(\pm 0.009\) (sys) at the effective redshift of 0.073. Our constraint is consistent with a Planck 2018 cosmological model, \(fσ_8 = 0.448\), within one standard deviation. Our improved methodology will enable similar analysis on future data, with even larger sample sizes and covering larger angular areas on the sky. △ Less

Submitted 6 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 18 pages, 11 figures and 2 tables. Accepted by MNRAS. You can find the code https://github.com/YanxiangL/Peculiar_velocity_fitting and data https://zenodo.org/record/6640513. Comments are welcome

HI HOD -- I. The Halo Occupation Distribution of HI Galaxies

Authors: Fei Qin, Cullan Howlett, Adam R. H. Stevens, David Parkinson

Abstract: The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to har… ▽ More The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to harness this information to improve our cosmological understanding, and learn about the physics of dark matter and dark energy, we need to accurately model the manner in which galaxies detected in HI trace the underlying matter distribution of the Universe. For this purpose, we develop a new HI-based Halo Occupation Distribution (HOD) model, which makes predictions for the number of galaxies present in dark matter halos conditional on their HI mass. The parameterised HOD model is fit and validated using the Dark Sage semi-analytic model, where we show that the HOD parameters can be modelled by simple linear and quadratic functions of HI mass. However, we also find that the clustering predicted by the HOD depends sensitively on the radial distributions of the HI galaxies within their host dark matter halos, which does not follow the NFW profile in the Dark Sage simulation. As such, this work enables -- for the first time -- a simple prescription for placing galaxies of different HI mass within dark matter halos in a way that is able to reproduce the HI mass-dependent galaxy clustering and HI mass function simultaneously and without requiring knowledge of the optical properties of the galaxies. Further efforts are required to demonstrate that this model can be used to produce large ensembles of mock galaxy catalogues for upcoming surveys. △ Less

Submitted 1 September, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 18 pages, 13 figures, 1 table. Published in ApJ

The DESI Survey Validation: Results from Visual Inspection of Bright Galaxies, Luminous Red Galaxies, and Emission Line Galaxies

Authors: Ting-Wen Lan, R. Tojeiro, E. Armengaud, J. Xavier Prochaska, T. M. Davis, David M. Alexander, A. Raichoor, Rongpu Zhou, Christophe Yeche, C. Balland, S. BenZvi, A. Berti, R. Canning, A. Carr, H. Chittenden, S. Cole, M. -C. Cousinou, K. Dawson, Biprateep Dey, K. Douglass, A. Edge, S. Escoffier, A. Glanville, S. Gontcho A Gontcho, J. Guy , et al. (57 additional authors not shown)

Abstract: 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 (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications.… ▽ More 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 (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications. We then utilize the VI redshift information to characterize the performance of the DESI operation. Based on the 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 for bright galaxies and ELGs and approximately 40 km/s for LRGs. The average redshift accuracy is within 10 km/s 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. △ Less

Submitted 15 January, 2023; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: 22 pages, 14 figures, 3 tables, one of a suite of 8 papers detailing targeting for DESI. ApJ accepted version with minor textual updates

Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument

Authors: B. Abareshi, J. Aguilar, S. Ahlen, Shadab Alam, David M. Alexander, R. Alfarsy, L. Allen, C. Allende Prieto, O. Alves, J. Ameel, E. Armengaud, J. Asorey, Alejandro Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, S. F. Beltran, B. Benavides, S. BenZvi, A. Berti, R. Besuner, Florian Beutler, D. Bianchi , et al. (242 additional authors not shown)

Abstract: 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 modifi… ▽ More 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Å > 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) △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: 78 pages, 32 figures, submitted to AJ