Browsing by Autor "J. Aguilar"

Now showing 1 - 20 of 20

Atacama Cosmology Telescope DR6 and DESI: Structure growth measurements from the cross-correlation of DESI legacy imaging galaxies and CMB lensing from ACT DR6 and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>P</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:math> PR4
(American Physical Society, 2025) Frank J. Qu; Qianjun Hang; Gerrit S. Farren; Boris Bolliet; J. Aguilar; S. P. Ahlen; Shadab Alam; D. Brooks; Yan-Chuan Cai; E. Calabrese
We measure the growth of cosmic density fluctuations on large scales and across the redshift range <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mn>0.3</a:mn><a:mo><</a:mo><a:mi>z</a:mi><a:mo><</a:mo><a:mn>0.8</a:mn></a:math> through galaxy clustering and the cross-correlation of the ACT data release 6 cosmic microwave background (CMB) lensing map and galaxies from the Dark Energy Spectroscopic Instrument Legacy Survey, using three galaxy samples spanning the redshifts of <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mn>0.3</c:mn><c:mo>≲</c:mo><c:mi>z</c:mi><c:mo>≲</c:mo><c:mn>0.45</c:mn></c:math>, <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mn>0.45</e:mn><e:mo>≲</e:mo><e:mi>z</e:mi><e:mo>≲</e:mo><e:mn>0.6</e:mn></e:math>, <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mn>0.6</g:mn><g:mo>≲</g:mo><g:mi>z</g:mi><g:mo>≲</g:mo><g:mn>0.8</g:mn></g:math>. We adopt a scale cut where nonlinear effects are negligible, so that the cosmological constraints are derived from the linear regime. We determine the amplitude of matter fluctuations over all three redshift bins using Atacama Cosmology Telescope (ACT) data alone to be <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:msub><i:mi>S</i:mi><i:mn>8</i:mn></i:msub><i:mo>≡</i:mo><i:msub><i:mi>σ</i:mi><i:mn>8</i:mn></i:msub><i:mo stretchy="false">(</i:mo><i:msub><i:mi mathvariant="normal">Ω</i:mi><i:mi>m</i:mi></i:msub><i:mo>/</i:mo><i:mn>0.3</i:mn><i:msup><i:mo stretchy="false">)</i:mo><i:mn>0.5</i:mn></i:msup><i:mo>=</i:mo><i:mn>0.772</i:mn><i:mo>±</i:mo><i:mn>0.040</i:mn></i:math> in a joint analysis combining the three redshift bins and ACT lensing alone. Using a combination of ACT and data we obtain <n:math xmlns:n="http://www.w3.org/1998/Math/MathML" display="inline"><n:msub><n:mi>S</n:mi><n:mn>8</n:mn></n:msub><n:mo>=</n:mo><n:mn>0.765</n:mn><n:mo>±</n:mo><n:mn>0.032</n:mn></n:math>. The lowest redshift bin used is the least constraining and exhibits a <p:math xmlns:p="http://www.w3.org/1998/Math/MathML" display="inline"><p:mo>∼</p:mo><p:mn>2</p:mn><p:mi>σ</p:mi></p:math> tension with the other redshift bins; thus we also report constraints excluding the first redshift bin, giving <r:math xmlns:r="http://www.w3.org/1998/Math/MathML" display="inline"><r:msub><r:mi>S</r:mi><r:mn>8</r:mn></r:msub><r:mo>=</r:mo><r:mn>0.785</r:mn><r:mo>±</r:mo><r:mn>0.033</r:mn></r:math> for the combination of ACT and . This result is in excellent agreement at the <t:math xmlns:t="http://www.w3.org/1998/Math/MathML" display="inline"><t:mn>0.3</t:mn><t:mi>σ</t:mi></t:math> level with measurements from galaxy lensing, but is <v:math xmlns:v="http://www.w3.org/1998/Math/MathML" display="inline"><v:mn>1.8</v:mn><v:mi>σ</v:mi></v:math> lower than predictions based on primary CMB data. Understanding whether this hint of discrepancy in the growth of structure at low redshifts arises from a fluctuation, from systematics in data, or from new physics is a high priority for forthcoming CMB lensing and galaxy cross-correlation analyses.
Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data
(American Physical Society, 2025) Bernardita Ried Guachalla; Emmanuel Schaan; Boryana Hadzhiyska; Simone Ferraro; J. Aguilar; S. P. Ahlen; Nicholas Battaglia; D. Bianchi; Richard Bond; David J. Brooks
The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel'dovich effect (kSZ). We precisely measure this effect ($S/N\ensuremath{\approx}10$) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over $\ensuremath{\ge}4,000\text{ }\text{ }{\mathrm{deg}}^{2}$. We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in $g$, $r$, and $z$ bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed $(\mathrm{PTE})=0.005$). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris $z=0.5$, $\mathrm{PTE}=0.37$) over weaker-feedback cases (IllustrisTNG $z=0.8$, $\mathrm{PTE}=0.045$), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast $S/N\ensuremath{\sim}50$ for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. These measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.
Clustering of DESI galaxies split by thermal Sunyaev-Zeldovich effect
(2025) M. Rashkovetskyi; Daniel J. Eisenstein; J. Aguilar; S. P. Ahlen; Abhijeet Anand; D. Bianchi; David Brooks; F. J. Castander; T. Claybaugh; Andrei Cuceu
The thermal Sunyaev-Zeldovich (tSZ) effect is associated with galaxy clusters - extremely large and dense structures tracing the dark matter with a higher bias than isolated galaxies. We propose to use the tSZ data to separate galaxies from redshift surveys into distinct subpopulations corresponding to different densities and biases independently of the redshift survey systematics. Leveraging the information from different environments, as in density-split and density-marked clustering, is known to tighten the constraints on cosmological parameters, like <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>Ω</mml:mi> <mml:mi>m</mml:mi> </mml:msub> </mml:math> , <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>σ</mml:mi> <mml:mn>8</mml:mn> </mml:msub> </mml:math> and neutrino mass. We use data from the Dark Energy Spectroscopic Instrument (DESI) and the Atacama Cosmology Telescope (ACT) in their region of overlap to demonstrate informative tSZ splitting of Luminous Red Galaxies (LRGs). We discover a significant increase in the large-scale clustering of DESI LRGs corresponding to detections starting from 1-2 sigma in the ACT DR6 + Planck tSZ Compton- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>y</mml:mi> </mml:math> map, below the cluster candidate threshold (4 sigma). We also find that such galaxies have higher line-of-sight coordinate (and velocity) dispersions and a higher number of close neighbors than both the full sample and near-zero tSZ regions. We produce simple simulations of tSZ maps that are intrinsically consistent with galaxy catalogs and do not include systematic effects, and find a similar pattern of large-scale clustering enhancement with tSZ effect significance. Moreover, we observe that this relative bias pattern remains largely unchanged with variations in the galaxy-halo connection model in our simulations. This is promising for future cosmological inference from tSZ-split clustering with semi-analytical models. Thus, we demonstrate that valuable cosmological information is present in the lower signal-to-noise regions of the thermal Sunyaev-Zeldovich map, extending far beyond the individual cluster candidates.
CMB lensing and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>Ly</mml:mi><mml:mi>α</mml:mi></mml:math> forest cross bispectrum from DESI’s first-year quasar sample
(American Physical Society, 2024) Naim Göksel Karaçaylı; Paul Martini; David H. Weinberg; Simone Ferraro; Roger de Belsunce; J. Aguilar; S. P. Ahlen; E. Armengaud; D. Brooks; T. Claybaugh
The squeezed cross-bispectrum ${B}_{\ensuremath{\kappa},\mathrm{Ly}\ensuremath{\alpha}}$ between the gravitational lensing in the cosmic microwave background and the 1D $\mathrm{Ly}\ensuremath{\alpha}$ forest power spectrum can constrain bias parameters and break degeneracies between ${\ensuremath{\sigma}}_{8}$ and other cosmological parameters. We detect ${B}_{\ensuremath{\kappa},\mathrm{Ly}\ensuremath{\alpha}}$ with $4.8\ensuremath{\sigma}$ significance at an effective redshift ${z}_{\mathrm{eff}}=2.4$ using Planck PR3 lensing map and over 280,000 quasar spectra from the Dark Energy Spectroscopic Instrument's first-year data. We test our measurement against metal contamination and foregrounds such as Galactic extinction and clusters of galaxies by deprojecting the thermal Sunyaev-Zeldovich effect. We compare our results to a tree-level perturbation theory calculation and find reasonable agreement between the model and measurement.
Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions
(Oxford University Press, 2025) Zachery Brown; R. Demina; A. G. Adame; S. Àvila; E. Chaussidon; Sihan Yuan; Violeta González-Pérez; J García-Bellido; B. Levi; J. Aguilar
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 $f_{\mathrm{NL}}$, 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 $f_{\mathrm{NL}}$ 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 haloes. 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 $f_{\mathrm{NL}}$. In LSS surveys, the effect of imaging systematics on two-point statistics is often degenerate with the PNG signal. Our proposed model employs three-point statistics primarily to break this degeneracy. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our method. Our forecast indicates the ability to constrain $f_{\mathrm{NL}}$ to a precision of $\sigma _{f_{\mathrm{NL}}} \approx 22$ with one year of DESI survey data, as well as the ability to constrain the imaging systematic weights in situ.
Construction of the damped <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>Ly</mml:mi> <mml:mi>α</mml:mi> </mml:mrow> </mml:math> absorber catalog for DESI DR2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mi>Ly</mml:mi> <mml:mi>α</mml:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:math> BAO
(American Physical Society, 2025) A. Brodzeller; Molly Wolfson; D. Martínez Santos; Ming-Feng Ho; T. Tan; Matthew M. Pieri; Andrei Cuceu; M. Karim; J. Aguilar; S. P. Ahlen
We present the Damped <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mrow> <a:mi>Ly</a:mi> <a:mi>α</a:mi> </a:mrow> </a:math> Toolkit for automated detection and characterization of damped <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"> <c:mrow> <c:mi>Ly</c:mi> <c:mi>α</c:mi> </c:mrow> </c:math> absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and column density. With an optimized quality cut on detection significance ( <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"> <e:mi mathvariant="normal">Δ</e:mi> <e:msubsup> <e:mi>χ</e:mi> <e:mi>r</e:mi> <e:mn>2</e:mn> </e:msubsup> <e:mo>></e:mo> <e:mn>0.03</e:mn> </e:math> ), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with <h:math xmlns:h="http://www.w3.org/1998/Math/MathML" display="inline"> <h:mi mathvariant="normal">S</h:mi> <h:mo>/</h:mo> <h:mrow> <h:mi mathvariant="normal">N</h:mi> <h:mo stretchy="false">></h:mo> <h:mn>2</h:mn> </h:mrow> </h:math> that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"> <m:mi mathvariant="normal">S</m:mi> <m:mo>/</m:mo> <m:mrow> <m:mi mathvariant="normal">N</m:mi> <m:mo stretchy="false">></m:mo> <m:mn>2</m:mn> </m:mrow> </m:math> spectra with predicted <r:math xmlns:r="http://www.w3.org/1998/Math/MathML" display="inline"> <r:mrow> <r:msub> <r:mrow> <r:mi>log</r:mi> </r:mrow> <r:mn>10</r:mn> </r:msub> <r:mo stretchy="false">(</r:mo> <r:msub> <r:mrow> <r:mi>N</r:mi> </r:mrow> <r:mrow> <r:mi mathvariant="monospace">HI</r:mi> </r:mrow> </r:msub> <r:mo stretchy="false">)</r:mo> <r:mo>></r:mo> <r:mn>20.3</r:mn> </r:mrow> </r:math> and detection significance <w:math xmlns:w="http://www.w3.org/1998/Math/MathML" display="inline"> <w:mi mathvariant="normal">Δ</w:mi> <w:msubsup> <w:mi>χ</w:mi> <w:mi>r</w:mi> <w:mn>2</w:mn> </w:msubsup> <w:mo>></w:mo> <w:mn>0.03</w:mn> </w:math> . We compare the Damped <z:math xmlns:z="http://www.w3.org/1998/Math/MathML" display="inline"> <z:mrow> <z:mi>Ly</z:mi> <z:mi>α</z:mi> </z:mrow> </z:math> Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the <bb:math xmlns:bb="http://www.w3.org/1998/Math/MathML" display="inline"> <bb:mrow> <bb:mi>Ly</bb:mi> <bb:mi>α</bb:mi> </bb:mrow> </bb:math> forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with <db:math xmlns:db="http://www.w3.org/1998/Math/MathML" display="inline"> <db:mrow> <db:msub> <db:mrow> <db:mi>log</db:mi> </db:mrow> <db:mn>10</db:mn> </db:msub> <db:mo stretchy="false">(</db:mo> <db:msub> <db:mrow> <db:mi>N</db:mi> </db:mrow> <db:mrow> <db:mi mathvariant="monospace">HI</db:mi> </db:mrow> </db:msub> <db:mo stretchy="false">)</db:mo> <db:mo>></db:mo> <db:mn>20.3</db:mn> </db:mrow> </db:math> from quasar spectra with <ib:math xmlns:ib="http://www.w3.org/1998/Math/MathML" display="inline"> <ib:mi mathvariant="normal">S</ib:mi> <ib:mo>/</ib:mo> <ib:mrow> <ib:mi mathvariant="normal">N</ib:mi> <ib:mo stretchy="false">></ib:mo> <ib:mn>2</ib:mn> </ib:mrow> </ib:math> . We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.
DESI DR2 results. I. Baryon acoustic oscillations from the Lyman alpha forest
(American Physical Society, 2025) M. Abdul Karim; J. Aguilar; S. Ahlen; C. Allende Prieto; O. Alves; A. Anand; U. Andrade; E. Armengaud; A. Aviles; S. Bailey
We present the baryon acoustic oscillation (BAO) measurements with the Lyman-<a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>α</a:mi></a:math> (<c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mi>Ly</c:mi><c:mi>α</c:mi></c:math>) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the autocorrelation of the <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mi>Ly</e:mi><e:mi>α</e:mi></e:math> forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of 2 larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mi>Ly</g:mi><g:mi>α</g:mi></g:math> absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:msub><i:mi>z</i:mi><i:mi>eff</i:mi></i:msub><i:mo>=</i:mo><i:mn>2.33</i:mn></i:math>. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"><k:mi>Ly</k:mi><k:mi>α</k:mi></k:math> BAO analysis for the first time. We measure the ratios <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"><m:msub><m:mi>D</m:mi><m:mi>H</m:mi></m:msub><m:mo stretchy="false">(</m:mo><m:msub><m:mi>z</m:mi><m:mi>eff</m:mi></m:msub><m:mo stretchy="false">)</m:mo><m:mo>/</m:mo><m:msub><m:mi>r</m:mi><m:mi>d</m:mi></m:msub><m:mo>=</m:mo><m:mn>8.632</m:mn><m:mo>±</m:mo><m:mn>0.098</m:mn><m:mo>±</m:mo><m:mn>0.026</m:mn></m:math> and <q:math xmlns:q="http://www.w3.org/1998/Math/MathML" display="inline"><q:msub><q:mi>D</q:mi><q:mi>M</q:mi></q:msub><q:mo stretchy="false">(</q:mo><q:msub><q:mi>z</q:mi><q:mi>eff</q:mi></q:msub><q:mo stretchy="false">)</q:mo><q:mo>/</q:mo><q:msub><q:mi>r</q:mi><q:mi>d</q:mi></q:msub><q:mo>=</q:mo><q:mn>38.99</q:mn><q:mo>±</q:mo><q:mn>0.52</q:mn><q:mo>±</q:mo><q:mn>0.12</q:mn></q:math>, where <u:math xmlns:u="http://www.w3.org/1998/Math/MathML" display="inline"><u:msub><u:mi>D</u:mi><u:mi>H</u:mi></u:msub><u:mo>=</u:mo><u:mi>c</u:mi><u:mo>/</u:mo><u:mi>H</u:mi><u:mo stretchy="false">(</u:mo><u:mi>z</u:mi><u:mo stretchy="false">)</u:mo></u:math> is the Hubble distance, <y:math xmlns:y="http://www.w3.org/1998/Math/MathML" display="inline"><y:msub><y:mi>D</y:mi><y:mi>M</y:mi></y:msub></y:math> is the transverse comoving distance, <ab:math xmlns:ab="http://www.w3.org/1998/Math/MathML" display="inline"><ab:msub><ab:mi>r</ab:mi><ab:mi>d</ab:mi></ab:msub></ab:math> is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.
DESI DR2 results. II. Measurements of baryon acoustic oscillations and cosmological constraints
(American Physical Society, 2025) M. Abdul Karim; J. Aguilar; S. P. Ahlen; Shadab Alam; L. Allen; Carlos Allende Prieto; O. Alves; Abhijeet Anand; U. Andrade; E. Armengaud
We present baryon acoustic oscillation (BAO) measurements from more than 14 million galaxies and quasars drawn from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2), based on three years of operation. For cosmology inference, these galaxy measurements are combined with DESI Lyman- <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi>α</a:mi> </a:math> forest BAO results presented in a companion paper (M. Abdul-Karim , companion paper, .). The DR2 BAO results are consistent with DESI DR1 and the Sloan Digital Sky Survey, and their distance-redshift relationship matches those from recent compilations of supernovae (SNe) over the same redshift range. The results are well described by a flat <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"> <c:mi mathvariant="normal">Λ</c:mi> </c:math> cold dark matter ( <f:math xmlns:f="http://www.w3.org/1998/Math/MathML" display="inline"> <f:mi mathvariant="normal">Λ</f:mi> <f:mi>CDM</f:mi> </f:math> ) model, but the parameters preferred by BAO are in mild, <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"> <i:mn>2.3</i:mn> <i:mi>σ</i:mi> </i:math> tension with those determined from the cosmic microwave background (CMB), although the DESI results are consistent with the acoustic angular scale <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"> <k:msub> <k:mi>θ</k:mi> <k:mo>*</k:mo> </k:msub> </k:math> that is well measured by Planck. This tension is alleviated by dark energy with a time-evolving equation of state parametrized by <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"> <m:msub> <m:mi>w</m:mi> <m:mn>0</m:mn> </m:msub> </m:math> and <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"> <o:msub> <o:mi>w</o:mi> <o:mi>a</o:mi> </o:msub> </o:math> , which provides a better fit to the data, with a favored solution in the quadrant with <q:math xmlns:q="http://www.w3.org/1998/Math/MathML" display="inline"> <q:msub> <q:mi>w</q:mi> <q:mn>0</q:mn> </q:msub> <q:mo>></q:mo> <q:mo>−</q:mo> <q:mn>1</q:mn> </q:math> and <s:math xmlns:s="http://www.w3.org/1998/Math/MathML" display="inline"> <s:msub> <s:mi>w</s:mi> <s:mi>a</s:mi> </s:msub> <s:mo><</s:mo> <s:mn>0</s:mn> </s:math> . This solution is preferred over <u:math xmlns:u="http://www.w3.org/1998/Math/MathML" display="inline"> <u:mi mathvariant="normal">Λ</u:mi> <u:mi>CDM</u:mi> </u:math> at <x:math xmlns:x="http://www.w3.org/1998/Math/MathML" display="inline"> <x:mn>3.1</x:mn> <x:mi>σ</x:mi> </x:math> for the combination of DESI BAO and CMB data. When also including SNe, the preference for a dynamical dark energy model over <z:math xmlns:z="http://www.w3.org/1998/Math/MathML" display="inline"> <z:mi mathvariant="normal">Λ</z:mi> <z:mi>CDM</z:mi> </z:math> ranges from <cb:math xmlns:cb="http://www.w3.org/1998/Math/MathML" display="inline"> <cb:mn>2.8</cb:mn> <cb:mo>−</cb:mo> <cb:mn>4.2</cb:mn> <cb:mi>σ</cb:mi> </cb:math> depending on which SNe sample is used. We present evidence from other data combinations which also favor the same behavior at high significance. From the combination of DESI and CMB we derive 95% upper limits on the sum of neutrino masses, finding <eb:math xmlns:eb="http://www.w3.org/1998/Math/MathML" display="inline"> <eb:mo>∑</eb:mo> <eb:msub> <eb:mi>m</eb:mi> <eb:mi>ν</eb:mi> </eb:msub> <eb:mo><</eb:mo> <eb:mn>0.064</eb:mn> <eb:mtext> </eb:mtext> <eb:mtext> </eb:mtext> <eb:mi>eV</eb:mi> </eb:math> assuming <gb:math xmlns:gb="http://www.w3.org/1998/Math/MathML" display="inline"> <gb:mi mathvariant="normal">Λ</gb:mi> <gb:mi>CDM</gb:mi> </gb:math> and <jb:math xmlns:jb="http://www.w3.org/1998/Math/MathML" display="inline"> <jb:mo>∑</jb:mo> <jb:msub> <jb:mi>m</jb:mi> <jb:mi>ν</jb:mi> </jb:msub> <jb:mo><</jb:mo> <jb:mn>0.16</jb:mn> <jb:mtext> </jb:mtext> <jb:mtext> </jb:mtext> <jb:mi>eV</jb:mi> </jb:math> in the <lb:math xmlns:lb="http://www.w3.org/1998/Math/MathML" display="inline"> <lb:msub> <lb:mi>w</lb:mi> <lb:mn>0</lb:mn> </lb:msub> <lb:msub> <lb:mi>w</lb:mi> <lb:mi>a</lb:mi> </lb:msub> </lb:math> model. Unless there is an unknown systematic error associated with one or more datasets, it is clear that <nb:math xmlns:nb="http://www.w3.org/1998/Math/MathML" display="inline"> <nb:mi mathvariant="normal">Λ</nb:mi> <nb:mi>CDM</nb:mi> </nb:math> is being challenged by the combination of DESI BAO with other measurements and that dynamical dark energy offers a possible solution.
Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey
(Oxford University Press, 2024) C Lamman; Daniel J. Eisenstein; J. E. Forero-Romero; J. Aguilar; S. P. Ahlen; S. Bailey; D. Bianchi; David J. Brooks; T. Claybaugh; Axel de la Macorra
ABSTRACT We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A bacusS ummitN-body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.
Detection of the pairwise kinematic Sunyaev-Zel’dovich effect and pairwise velocity with DESI DR1 galaxies and ACT DR6 and <i>Planck</i> CMB data
(American Physical Society, 2026) Yulin Gong; Patricio A. Gallardo; Rachel Bean; Jenna Moore; E. M. Vavagiakis; Nicholas Battaglia; Boryana Hadzhiyska; Yun-Hsin Hsu; J. Aguilar; S. Ahlen
Galaxy-multiplet clustering from DESI DR2
(Oxford University Press, 2025) Hanyue Wang; Daniel J. Eisenstein; J. Aguilar; S. P. Ahlen; D. Bianchi; David H. Brooks; T. Claybaugh; Axel de la Macorra; Arjun Dey; Biprateep Dey
ABSTRACT We present an efficient estimator for higher order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter haloes than individual galaxies. When comparing the observed clustering statistics with the mock catalogues generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point autocorrelation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher order statistic on galaxy–halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy haloes in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy–dark matter connection.
Identifying Anomalous DESI Galaxy Spectra with a Variational Autoencoder
(Oxford University Press, 2026) C. Nicolaou; Rowina S Nathan; O Lahav; A. Palmese; A. Saintonge; J. Aguilar; S. P. Ahlen; C. Allende Prieto; S. Bailey; S. BenZvi
ABSTRACT The tens of millions of spectra being captured by the Dark Energy Spectroscopic Instrument (DESI) provide tremendous discovery potential. In this work we show how Machine Learning, in particular Variational Autoencoders (VAE), can detect anomalies in a sample of approximately 200 000 DESI spectra comprising galaxies, quasars and stars. We demonstrate that the VAE can compress the dimensionality of a spectrum by a factor of 100, while still retaining enough information to accurately reconstruct spectral features. We detect anomalous spectra as those with high reconstruction error and those which are isolated in the VAE latent representation. The anomalies identified fall into two categories: spectra with artefacts and spectra with unique physical features. Awareness of the former could improve the DESI spectroscopic pipeline; whilst the latter could help us discover new and unusual objects. To further curate the list of outliers identified, we use the Astronomaly package which employs Active Learning to provide personalized outlier recommendations for visual inspection. In this work we also explore the VAE latent space, finding that different object classes and subclasses are separated despite being unlabelled. We inject controlled synthetic anomalies and analyse their locations in the latent space to illustrate how the VAE responds to atypical spectral features; and we demonstrate the interpretability of this latent space by identifying tracks within it that correspond to various spectral characteristics. In upcoming work we hope to apply the methods presented here to search for both systematics and astrophysically interesting objects in much larger datasets of DESI spectra.
Measurements of the thermal Sunyaev-Zel’dovich effect with ACT and DESI luminous red galaxies
(American Physical Society, 2025) R. Henry Liu; Simone Ferraro; Emmanuel Schaan; Rongpu Zhou; J. Aguilar; S. P. Ahlen; Nicholas Battaglia; Davide Bianchi; David Brooks; T. Claybaugh
Cosmic Microwave Background (CMB) photons scatter off the free-electron gas in galaxies and clusters, allowing us to use the CMB as a backlight to probe the gas in and around low-redshift galaxies. The thermal Sunyaev-Zel'dovich effect, sourced by hot electrons in high-density environments, measures the thermal pressure of the target objects, shedding light on halo thermodynamics and galaxy formation, and providing a path toward understanding the baryon distribution around cosmic structures. We use a combination of high-resolution CMB maps from the Atacama Cosmology Telescope and photometric luminous red galaxy catalogs from the Dark Energy Spectroscopic Instrument to measure the thermal Sunyaev-Zel'dovich signal in four redshift bins from $z=0.4$ to $z=1.2$, with a combined detection significance of $19\ensuremath{\sigma}$ when stacking on the fiducial CMB Compton-$y$ map. We discuss possible sources of contamination, finding that residual dust emission associated with the target galaxies is important and limits current analyses. We discuss several mitigation strategies and quantify the residual modeling uncertainty. This work complements closely related measurements of the kinematic Sunyaev-Zel'dovich and weak lensing of the same galaxies.
Nearby stellar substructures in the Galactic halo from DESI Milky Way Survey Year 1 Data Release
(Oxford University Press, 2025) Bokyoung Kim; S. E. Koposov; Ting S. Li; Sophia Lilleengen; Andrew P. Cooper; Andreia Carrillo; Monica Valluri; A. H. Riley; J. Han; J. Aguilar
ABSTRACT We report five nearby ($d_{\mathrm{helio}} &lt; 5$ kpc) stellar substructures in the Galactic halo from a subset of 138 661 stars in the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey Year 1 Data Release. With an unsupervised clustering algorithm, HDBSCAN*, these substructures are independently identified in Integrals of Motion ($E_{\rm tot}$, $L_{\rm z}$, $\log {J_r}$, $\log {J_z}$) space and Galactocentric cylindrical velocity space ($V_{R}$, $V_{\phi }$, $V_{z}$). We associate all identified clusters with known nearby substructures (Helmi streams, M18-Cand10/MMH-1, Sequoia, Antaeus, and ED-2) previously reported in various studies. With metallicities precisely measured by DESI, we confirm that the Helmi streams, M18-Cand10, and ED-2 are chemically distinct from local halo stars. We have characterized the chemodynamic properties of each dynamic group, including their metallicity dispersions, to associate them with their progenitor types (globular cluster or dwarf galaxy). Our approach for searching substructures with HDBSCAN* reliably detects real substructures in the Galactic halo, suggesting that applying the same method can lead to the discovery of new substructures in future DESI data. With more stars from future DESI data releases and improved astrometry from the upcoming Gaia Data Release 4, we will have a more detailed blueprint of the Galactic halo, offering a significant improvement in our understanding of the formation and evolutionary history of the Milky Way Galaxy.
PAC in DESI. I. Galaxy stellar mass function into the 106 M⊙ frontier
(Oxford University Press, 2025) Kun Xu; Yipeng Jing; Shaun Cole; Carlos S. Frenk; Sownak Bose; Willem Elbers; Hesheng Wang; Yirong Wang; Shannon Moore; J. Aguilar
ABSTRACT The Photometric objects Around Cosmic webs (PAC) method integrates cosmological photometric and spectroscopic surveys, offering valuable insights into galaxy formation. PAC measures the excess surface density of photometric objects, $\bar{n}_2w_{{\rm {p}}}$, with specific physical properties around spectroscopic tracers. In this study, we improve the PAC method to make it more rigorous and eliminate the need for redshift bins. We apply the enhanced PAC method to the DESI Y1 BGS Bright spectroscopic sample and the deep Dark Energy Camera Legacy Survey (DECaLS) photometric sample, obtaining $\bar{n}_2w_{{\mathrm {p}}}$ measurements across the complete stellar mass range, from $10^{5.3}$ to $10^{11.5}\,{\rm M}_{\odot }$ for blue galaxies, and from $10^{6.3}$ to $10^{11.9}\,{\rm M}_{\odot }$ for red galaxies. We combine $\bar{n}_2w_{{\rm {p}}}$ with $w_{{\rm {p}}}$ measurements from the BGS sample, which is not necessarily complete in stellar mass. Assuming that galaxy bias is primarily determined by stellar mass and colour, we derive the galaxy stellar mass functions (GSMFs) down to $10^{5.3}\,{\rm M}_{\odot }$ for blue galaxies and $10^{6.3}\,{\rm M}_{\odot }$ for red galaxies, while also setting lower limits for smaller masses. The blue and red GSMFs are well described by single and double Schechter functions, respectively, with low-mass end slopes of $\alpha _{\rm {blue}}=-1.54^{+0.02}_{-0.02}$ and $\alpha _{\rm {red}}=-2.50^{+0.08}_{-0.08}$, resulting in the dominance of red galaxies below $10^{7.6}\,{\rm M}_{\odot }$. Stage-IV cosmological photometric surveys, capable of reaching 2–3 mag deeper than DECaLS, present an opportunity to explore the entire galaxy population in the local universe with PAC. This advancement allows us to address critical questions regarding the nature of dark matter, the physics of reionization, and the formation of dwarf galaxies.
Spectroscopic search for optical emission lines from dark matter decay
(American Physical Society, 2024) Hanyue Wang; Daniel J. Eisenstein; J. Aguilar; S. P. Ahlen; S. Bailey; David J. Brooks; T. Claybaugh; Axel de la Macorra; Peter Doel; J. E. Forero-Romero
We search for narrow-line optical emission from dark matter decay by stacking dark-sky spectra from the Dark Energy Spectroscopic Instrument (DESI) at the redshift of nearby galaxies from DESI's Bright Galaxy and Luminous Red Galaxy samples. Our search uses regions separated by 5 to 20 arcsec from the centers of the galaxies, corresponding to an impact parameter of approximately 50 kpc. No unidentified spectral line shows up in the search, and we place a line flux limit of ${10}^{\ensuremath{-}19}\text{ }\text{ }\mathrm{ergs}/\mathrm{s}/{\mathrm{cm}}^{2}/{\mathrm{arcsec}}^{2}$ on emissions in the wavelength range of $2000--9000\stackrel{\ensuremath{\circ}}{\mathrm{A}}$. This places the tightest constraints yet on the two-photon decay of dark matter in the mass range of 5 to 12 eV, with a particle lifetime exceeding $3\ifmmode\times\else\texttimes\fi{}{10}^{25}\text{ }\text{ }\mathrm{s}$. This detection limit also implies that the line surface brightness contributed from all dark matter along the line of sight is at least 2 orders of magnitude lower than the measured extragalactic background light (EBL), ruling out the possibility that narrow optical-line emission from dark matter decay is a major source of the EBL.
Study of the connected four-point correlation function of galaxies from the DESI Data Release 1 luminous red galaxy sample
(American Physical Society, 2025) J. Hou; R. N. Cahn; J. Aguilar; S. Ahlen; D. Bianchi; D. Brooks; T. Claybaugh; P. Doel; S. Ferraro; J. E. Forero-Romero
We present a measurement of the non-Gaussian four-point correlation function (4PCF) from the DESI DR1 luminous red galaxy (LRG) sample. For the gravitationally induced parity-even 4PCF, we detect a signal with a significance of <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mrow> <a:mn>14.7</a:mn> <a:mi>σ</a:mi> </a:mrow> </a:math> using our fiducial setup. We assess the robustness of this detection through a series of validation tests, including auto and cross-correlation analyses, sky partitioning across multiple patch combinations, and variations in radial scale cuts. Due to the low completeness of the sample, we find that differences in fiber assignment implementation schemes can significantly impact estimation of the covariance and introduce biases in the data vector. After correcting for these effects, all tests yield consistent results. This is one of the first measurements of the connected 4PCF on the DESI LRG sample; the good agreement between the simulation and the data implies that the amplitude of the density fluctuation inferred from the connected 4PCF is consistent with the Planck <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"> <c:mi mathvariant="normal">Λ</c:mi> <c:mi>CDM</c:mi> </c:math> cosmology. The methodology and diagnostic framework established in this work provide a foundation for interpreting parity-odd 4PCF.
The atomic gas sequence and mass–metallicity relation from dwarfs to massive galaxies
(Oxford University Press, 2024) D. Scholte; A. Saintonge; John Moustakas; Barbara Catinella; Hu Zou; Biprateep Dey; J. Aguilar; S. P. Ahlen; Abhijeet Anand; Robert Blum
ABSTRACT Galaxy scaling relations provide insights into the processes that drive galaxy evolution. The extension of these scaling relations into the dwarf galaxy regime is of particular interest. This is because dwarf galaxies represent a crucial stage in galaxy evolution, and understanding them could also shed light on their role in reionizing the early Universe. There is currently no consensus on the processes that dominate the evolution of dwarfs. In this work, we constrain the atomic gas sequence (stellar mass versus atomic gas fraction) and mass–metallicity relation (stellar mass versus gas-phase metallicity) from dwarf ($10^{6.5} \, {\rm M}_{\odot }$) to massive ($10^{11.5} \, {\rm M}_{\odot }$) galaxies in the local Universe. The combined optical and 21-cm spectroscopic observations of the Dark Energy Spectroscopic Instrument and Arecibo Legacy Fast ALFA surveys allow us to constrain both scaling relations simultaneously. We find a slope change of the atomic gas sequence at a stellar mass of ${\sim} 10^{9} \, \textrm{M}_{\odot }$. We also find that the shape and scatter of the atomic gas sequence and mass–metallicity relation are strongly linked for both dwarfs and more massive galaxies. Consequently, the low-mass slope change of the atomic gas sequence is imprinted onto the mass–metallicity relation of dwarf galaxies. The mass scale of the measured slope change is consistent with a predicted escape velocity threshold below which low-mass galaxies experience significant supernova-driven gas loss, as well as with a reduction in cold gas accretion onto more massive galaxies.
The rate of extreme coronal line emitting galaxies in the Sloan Digital Sky Survey and their relation to tidal disruption events
(Oxford University Press, 2024) J. A. Callow; Or Graur; P. Clark; A. Palmese; J. Aguilar; S. P. Ahlen; S. BenZvi; David H. Brooks; T. Claybaugh; Axel de la Macorra
ABSTRACT High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalized rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included.
Validation of the DESI DR2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>Ly</mml:mi> <mml:mi>α</mml:mi> </mml:math> BAO analysis using synthetic datasets
(American Physical Society, 2025) L. Casas; H. K. Herrera-Alcantar; J. Chaves-Montero; A. Cuceu; A. Font-Ribera; M. Lokken; M. Abdul-Karim; C. Ramírez-Pérez; D. Alonso; J. Aguilar
International audience