Engineering topics

Dey, A. (ORCID:0000000249284003)

Publications and source records attributed to Dey, A. (ORCID:0000000249284003).

At least 19 records

DOE OSTI2026

Probing the Environment around GW170817 with DESI: Insights on Galaxy Group Peculiar Velocities for Standard Siren Measurements

We present a new measurement of the Hubble constant, H 0 , following the gravitational-wave event GW170817 and Dark Energy Spectroscopic Instrument (DESI) observations. A standard siren measurement with a nearby (luminosity distance ∼40 Mpc) event such as GW170817 is typically sensitive to the peculiar motion of the host galaxy owing to local dynamics. Previous measurements from this event have taken advantage of peculiar velocity measurements of nearby galaxies, including a handful of objects in the galaxy group that the host of the event, NGC 4993, has been associated with. Still, the group’s properties and NGC 4993’s membership were debated. We present DESI observations of thousands of galaxies in the vicinity of NGC 4993, resulting in 39 group galaxies and a fivefold increase in galaxies compared to previous observations, with many contributing to a peculiar velocity measurement. Examining the local dynamics, our observations support the presence of a galaxy group of which NGC 4993 is a part with a halo mass of order ∼10 13 M ⊙ . Using peculiar velocity measurements from our fundamental plane galaxy observations, we find $H_0 = 70.9^{+6.4}_{-8.5}$ km s −1 Mpc −1 . In addition, using a peculiar velocity measurement for NGC 4993 from surface brightness fluctuations in Cosmicflows-4, we find $H_0 = 73.4^{+3.3}_{-3.9}$ km s −1 Mpc −1 . We study the impact of different galaxy selection criteria on the determination of the peculiar velocity and, in turn, on the H 0 measurement. Our results demonstrate the value of multiplexed spectroscopic observations for probing the local environments of gravitational-wave events used in standard siren measurements.

Amsellem, A. J. [Carnegie Mellon University, Pitts↗

DOE OSTI2026

DESI Strong Lens Foundry. I. HST Observations and Modeling with GIGA-Lens

We present the Dark Energy Spectroscopic Instrument (DESI) Strong Lens Foundry. We discovered ∼3500 new strong gravitational lens candidates in the DESI Legacy Imaging Surveys using residual neural networks (ResNet). We observed a subset (51) of our candidates using the Hubble Space Telescope (HST). Except for one ambiguous case, we have confirmed 50 of the 51 candidates to be strong lenses. We also briefly describe spectroscopic follow-up observations by DESI and Keck NIRES programs. From this very rich data set, a number of studies will be carried out, including evaluating the quality of the ResNet search candidates and lens modeling. In this paper, we present our initial effort in these directions. In particular, as a demonstration, we present the lens model for DESI-165.4754−06.0423, with imaging data from HST, and lens and source redshifts from DESI and Keck NIRES, respectively. In this effort, we have applied a fully forward-modeling Bayesian approach (GIGA-Lens), using multiple GPUs, to a strong lens with HST data, and achieved statistical convergence.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A Comprehensive Characterization of Galaxy-cool CGM Connections at z < 0.4 with DESI Year 1 Data

We investigate the relationships between the cool circumgalactic medium (CGM), traced by Ca II absorption lines, and galaxy properties at z < 0.4 using ∼900,000 galaxy–quasar pairs within 200 kpc from the Year 1 data of the Dark Energy Spectroscopic Instrument (DESI). This large data set enables us to obtain composite spectra with sensitivity reaching to the mÅ level and to explore the Ca II absorption as a function of stellar mass, star formation rate (SFR), redshift, and galaxy types, including active galactic nuclei (AGNs). Our results show a positive correlation between the absorption strength and stellar mass of star-forming galaxies with $\langle$$W$$^{Ca II}_{0}$$\rangle$ α $M$$^{0.5}_{*}$ over 3 orders of magnitude in stellar mass from ∼10 8 to 10 11 M ⊙ , while such a mass dependence is weaker for quiescent galaxies. At a fixed mass, Ca II absorption is stronger around star-forming galaxies than quiescent ones especially within impact parameters <30 kpc. Among star-forming galaxies, the Ca II absorption further correlates with SFR, following ∝SFR 0.3 . However, in contrast to the results at higher redshifts, stronger absorption is not preferentially observed along the minor axis of star-forming galaxies, indicating a possible redshift evolution of CGM dynamics resulting from galactic feedback. Moreover, no significant difference between the properties of the cool gas around AGNs and galaxies is detected. Finally, we measure the absorption profiles with respect to the virial radius of dark matter halos and show that the total Ca II mass in the CGM is comparable to the Ca mass in the ISM of galaxies.

Circumgalactic medium↗

DOE OSTI2025

DESI DR1 Lyα 1D power spectrum: the optimal estimator measurement

The one-dimensional power spectrum P 1D of Lyα forest offers rich insights into cosmological and astrophysical parameters, including constraints on the sum of neutrino masses, warm dark matter models, and the thermal state of the intergalactic medium. We present the measurement of P 1D using the optimal quadratic maximum likelihood estimator applied to over 300,000 Lyα quasars from Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) survey. This sample represents the largest to date for P 1D measurements and is larger than the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) by a factor of 1.7. We conduct a meticulous investigation of instrumental and analysis systematics and quantify their impact on P 1D . This includes the development of a cross-exposure estimator that eliminates the need to model the pipeline noise and has strong potential for future P 1D measurements. We also present new insights into metal contamination through the 1D correlation function. Using a fitting function we measure the evolution of the Lyα forest bias with high precision: b F (z) = (-0.218 ± 0.002) × ((1 + z)/4) 2.96±0.06 . In a companion validation paper, we substantially extend our previous suite of CCD image simulations to quantify the pipeline's exquisite performance accurately. In another companion paper, we present DR1 P 1D measurements using the Fast Fourier Transform (FFT) approach to power spectrum estimation. These two measurements produce a forest bias parameter that differs by 2.2 sigma. However, our model is simplistic, so this disagreement will be investigated in future work.

Lyman alpha forest↗

DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 V: Full-Shape galaxy clustering from galaxies and quasars

We present the measurements and cosmological implications of the galaxy two-point clustering using over 4.7 million unique galaxy and quasar redshifts in the range 0.1 < z < 2.1 divided into six redshift bins over a ∼ 7,500 square degree footprint, from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). By fitting the full power spectrum, we extend previous DESI DR1 baryon acoustic oscillation (BAO) measurements to include redshift-space distortions and signals from the matter-radiation equality scale. For the first time, this Full-Shape analysis is blinded at the catalogue-level to avoid confirmation bias and the systematic errors are accounted for at the two-point clustering level, which automatically propagates them into any cosmological parameter. When analyzing the data in terms of compressed model-agnostic variables, we obtain a combined precision of 4.7% on the amplitude of the redshift space distortion (RSD) signal reaching a similar precision with just one year of DESI data than with twenty years of observation from the previous generation survey. We also analyze the data to directly constrain the cosmological parameters within the ΛCDM model using perturbation theory and combine this information with the reconstructed DESI DR1 galaxy BAO. Using a Big Bang Nucleosynthesis Gaussian prior on the baryon density parameter, ω b , and a weak Gaussian prior on the spectral index, n s , we constrain the matter density is Ω m = 0.296±0.010 and the Hubble constant H 0 = (68.63 ± 0.79)[km s -1 Mpc -1 ]. Additionally, we measure the amplitude of clustering σ 8 = 0.841±0.034. The DESI DR1 galaxy clustering results are in agreement with the ΛCDM model based on general relativity with parameters consistent with those from Planck. The cosmological interpretation of these results in combination with DESI DR1 Ly-α forest data and external datasets are presented in the companion paper [1].

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 VII: cosmological constraints from the full-shape modeling of clustering measurements

We present cosmological results from the measurement of clustering of galaxy, quasar and Lyman-α forest tracers from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). We adopt the full-shape (FS) modeling of the power spectrum, including the effects of redshift-space distortions, in an analysis which has been thoroughly validated in a series of supporting papers as summarised in [1]. We combine the full-shape information with DESI's DR1 constraints from the baryon acoustic oscillations (BAO) of these tracers. In the flat ΛCDM cosmological model, DESI (FS+BAO), combined with a baryon density prior from Big Bang Nucleosynthesis and a weak prior on the scalar spectral index, determines matter density to Ω m = 0.2962 ± 0.0095, and the amplitude of mass fluctuations to σ 8 = 0.842 ± 0.034. The addition of the cosmic microwave background (CMB) data tightens these constraints to Ω m = 0.3056 ± 0.0049 and σ 8 = 0.8121 ± 0.0053, while further addition of the joint clustering and lensing analysis from the Dark Energy Survey Year-3 (DESY3) data further improves these measurements, and leads to a 0.4% determination of the Hubble constant, H 0 = (68.40 ± 0.27) km s -1 Mpc -1 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI (FS+BAO) with CMB and type Ia supernovae continue to show the preference, previously found in the DESI DR1 BAO analysis, for w 0 > -1 and w a < 0 with similar levels of significance. DESI data, in combination with the CMB, improve the upper limits on the sum of the neutrino masses relative to the case when only the DR1 BAO was available, giving ∑m ν < 0.071 eV at 95% confidence. We finally constrain deviations from general relativity represented by two modified gravity parameters. DESI (FS+BAO) data alone measure the parameter that controls the clustering of massive particles, μ 0 = 0.11 +0.45 -0.54 , in agreement with the zero value predicted by general relativity. The combination of DESI with the CMB and the clustering and lensing analysis from DESY3 constrains both modified-gravity parameters, giving μ 0 = 0.04 ± 0.22 and Σ 0 = 0.044 ± 0.047, again in agreement with general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 II: sample definitions, characteristics, and two-point clustering statistics

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

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Mitigating imaging systematics for DESI 2024 emission Line Galaxies and beyond

Emission Line Galaxies (ELGs) are one of the main tracers that the Dark Energy Spectroscopic Instrument (DESI) uses to probe the universe. However, they are afflicted by strong spurious correlations between target density and observing conditions known as imaging systematics. In this paper, we present the imaging systematics mitigation applied to the DESI Data Release 1 (DR1) large-scale structure catalogs used in the DESI 2024 cosmological analyses. We also explore extensions of the fiducial treatment. This includes a combined approach, through forward image simulations (Obiwan) in conjunction with neural network-based regression, to obtain an angular selection function that mitigates the imaging systematics observed in the DESI DR1 ELGs target density. We further derive a line of sight selection function from the forward model that removes the strong redshift dependence between imaging systematics and low redshift ELGs. Combining both angular and redshift-dependent systematics, we construct a three-dimensional selection function and assess the impact of all selection functions on clustering statistics. We quantify differences between these extended treatments and the fiducial treatment in terms of the measured 2-point statistics. We find that the results are generally consistent with the fiducial treatment and conclude that the differences are far less than the imaging systematics uncertainty included in DESI 2024 full-shape measurements. We extend our investigation to the ELGs at 0.6 < z < 0.8, i.e., beyond the redshift range (0.8 < z < 1.6) adopted for the DESI clustering catalog, and demonstrate that determining the full three-dimensional selection function is necessary in this redshift range. Our tests showed that all changes are consistent with statistical noise for BAO analyses indicating they are robust to even severe imaging systematics. Specific tests for the full-shape analysis will be presented in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Full-Shape analysis of the power spectrum and bispectrum of DESI DR1 LRG and QSO samples

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

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 III: baryon acoustic oscillations from galaxies and quasars

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 employ a unified BAO analysis method across all 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 Gpc3, 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 that our observed BAO scales are systematically larger than the prediction of the Planck 2018-ΛCDM 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 VI: cosmological constraints from the measurements of baryon acoustic oscillations

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. To mitigate confirmation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data alone are consistent with the standard flat ΛCDM cosmological model with a matter density Ω m =0.295±0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and the robustly measured acoustic angular scale from the cosmic microwave background (CMB), DESI requires H 0 =(68.52±0.62) km s -1 Mpc -1 . In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find Ω m =0.307± 0.005 and H 0 =(67.97±0.38) km s -1 Mpc -1 . Extending the baseline model with a constant dark energy equation of state parameter w, DESI BAO alone requirew=-0.99 +0.15 -0.13 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI with CMB or with type Ia supernovae (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 the Pantheon+, Union3, or DES-SN5YR supernova datasets respectively. For the flat ΛCDM model with the sum of neutrino mass ∑ m ν free, combining the DESI and CMB data yields an upper limit ∑ m ν < 0.072 (0.113) eV at 95% confidence for a ∑ m ν > 0 (∑ m ν > 0.059) eV prior. These neutrino-mass constraints are substantially relaxed if the background dynamics are allowed to deviate from flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Candidate strongly lensed type Ia supernovae in the Zwicky Transient Facility archive

Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 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 (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) 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. The candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, 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 ± 3 and 34 ± 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia 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, such as the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

We present ∼115 new spectroscopically identified members of the GD-1 tidal stream observed with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a “cocoon,” which is a broad (FWHM ∼ 2:932 ∼ 460 pc) and kinematically hot (velocity dispersion, σ ∼ 5–8 km s −1 ) component that surrounds a narrower (FWHM ∼ 0:353 ∼ 55 pc) and colder [(σ = 3.09 ± 0.76 km s −1 )] thin stream component (based on a median per star velocity precision of 2.7 km s −1 ). The cocoon extends over at least a 30 ∘ segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: −2.54 ± 0.04 dex and −2.47 ± 0.06 dex, suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon comprises debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some of these. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman alpha forest

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 eff = 2.33 with 2% precision, H(z eff ) = (239.2 ± 4.8) (147.09 Mpc /r d ) km/s/Mpc. Similarly, we present a 2.4% measurement of the transverse comoving distance to the same redshift, DM(z eff ) = (5.84 ± 0.14) (r d /147.09 Mpc) Gpc. Together with other DESI BAO measurements at lower redshifts, these results are used in a companion paper to constrain cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The construction of large-scale structure catalogs for the Dark Energy Spectroscopic Instrument

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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Validating the galaxy and quasar catalog-level blinding scheme for the DESI 2024 analysis

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 series of sanity checks; the latter allows probing complexities in real data not captured in mocks. We find that the blinding strategy alters the data vector in a controlled way, and the BAO and RSD analysis choices are robust to blinding. 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

HOD-dependent systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

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 estimate the level of systematics induced in the DESI BAO analysis due the assumed Halo Occupation Distribution (HOD) model for the Emission Line Galaxy (ELG) tracer. We make use of mock galaxy catalogs constructed by fitting various HOD models to early DESI data, namely the One-Percent survey data. Our analysis includes typical HOD models for the ELG tracer used in the literature as well as extensions to the baseline models. Among the extensions, we consider various recipes for galactic conformity and assembly bias. We use 25 AbacusSummit simulations under the ΛCDM cosmology for each HOD model and perform independent analyses in Fourier space and in configuration space. To recover the BAO signal from our mocks we perform BAO reconstruction and apply the control variates technique to reduce sample variance noise. Our BAO analyses can recover the isotropic BAO parameter α iso within 0.1% and the Alcock Paczynski parameter α AP within 0.3%. Overall, we find that the 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, for the assumptions made.

79 ASTRONOMY AND ASTROPHYSICS↗