We present a measurement of Baryon Acoustic Oscillations (BAO) in the cross-correlation of triply ionized carbon C IV absorption with the positions of quasars (QSO) and Emission Line Galaxies (ELG). We use quasars and ELGs from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our data sample consists of 2.5 million quasars, 3.1 million ELGs, and the C IV absorption is measured along the line of sight of 1.5 million high redshift quasars with z > 1.3. We measure the isotropic BAO signal at 4.2σ for the CIV×QSO cross-correlation. This translates into a 3.0
Environmental effects are a primary driver of elevated gas-phase metallicities in galaxies around massive clusters, but the underlying physical mechanisms for this satellite metallicity enhancement (SME) are still unclear. Using the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, we present the first measurement of the average SME as a function of projected cluster-centric distance. The resulting profile reveals three distinct regimes: a steep decline from the cluster center, a plateau near the cluster boundary, and an extended downturn across several cluster radii. Remarkably, the complex shape and amplitude of this observed SME profile are successfully reproduced in the EAGLE cosmological simulation. Drawing insights from EAGLE, we develop a novel satellite chemical evolution model to decompose the observed SME into physical contributions from suppressed star formation, stellar mass loss, and enriched gas inflow. Our analysis shows that continuous accretion of enriched intracluster medium dominates the SME plateau within the cluster virial radius, while mass loss and quenching jointly drive the rapid metallicity decline in the cluster core. Our method disentangles the impacts of three environmental processes on galactic chemical enrichment in the cosmic web, providing a powerful framework for understanding cluster galaxy evolution with current and future spectroscopic surveys.
The determination of distances using time-series photometry of Type Ia supernovae (SNe Ia) relies on a 5
In preparation for the first cosmological measurements from the full shape of the Lyman-$\alpha$ (Ly $\alpha$) forest from Dark Energy Spectroscopic Instrument (DESI), we must carefully model all relevant systematics that might bias our analysis. It has been shown that random quasar redshift errors produce a smoothing effect on the mean quasar continuum in the Ly $\alpha$ forest region. This, in turn, gives rise to spurious features in the Ly $\alpha$ autocorrelation and its cross-correlation with quasars. Using synthetic data sets based on the DESI survey, we confirm that the impact on baryon acoustic oscillation measurements is small, but that a bias is introduced to parameters which depend on the full shape of our correlations. We combine a model of this contamination in the cross-correlation with a new model we introduce here for the autocorrelation. These are parametrized by three parameters, which, when included in a joint fit to both correlation functions, successfully eliminate any impact of redshift errors on our full-shape constraints. We also present a strategy for removing this contamination from real data, by removing $\sim$0.3 per cent of correlating pairs.
Here, we present the results of our search through the early data release (EDR) of the Dark Energy Spectroscopic Instrument (DESI) for extreme coronal line emitters (ECLEs) - a rare classification of galaxies displaying strong, high-ionization iron coronal emission lines within their spectra. With the requirement of a strong X-ray continuum to generate the coronal emission, ECLEs have been linked to both active galactic nuclei (AGNs) and tidal disruption events (TDEs). We focus our search on identifying TDE-linked ECLEs. We identify three such objects within the EDR sample, highlighting DESI's effectiveness for discovering new nuclear transients, and determine a galaxy-normalized TDE-linked ECLE rate of R-G = 5(-3 )(+5)& times; 10(-6 )galaxy(-1 )yr(-1) at a median redshift of z = 0.2 - broadly consistent with previous work. Additionally, we identify more than 200 AGNs displaying coronal emission lines, which serve as the primary astrophysical contaminants in searches for TDE-related events. We also include an outline of the custom python code developed for this search.
We perform a frequentist analysis using the standard profile likelihood method for clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). While Bayesian inferences for Effective Field Theory models of galaxy clustering can be highly sensitive to the choice of priors for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare Bayesian and frequentist constraints for the parameter set {σ_8, H_0, Ω_m, w_0, w_a} when fitting to the full-shape of the power spectrum multipoles, the post-reconstruction Baryon Acoustic Oscillation (BAO) measurements, as well as external datasets from the CMB and type Ia supernovae measurements. Bayesian prior effects are very significant for the w_0w_aCDM model; while the 1 σ frequentist confidence intervals encompass the maximum a posteriori (MAP), the Bayesian credible intervals almost always exclude the maximum likelihood estimate (MLE) and the MAP - indicating strong prior volume projection effects - unless supernovae data are included. We observe limited prior effects for the ΛCDM model, due to the reduced number of parameters. When DESI full-shape and BAO data are jointly fit, we obtain the following 1σ frequentist confidence intervals for ΛCDM (w_0w_aCDM): σ_8 = 0.867^+0.048_-0.041 , H_0 = 68.91^+0.80_-0.79 km s^-1Mpc^-1 , Ω_m = 0.3038±0.0110 (σ_8 = 0.793^+0.069_-0.048 , H_0 = 64.9^+4.8_-2.8 km s^-1Mpc^-1 , Ω_m = 0.369^+0.029_-0.059 , w_0 = -0.24^+0.17_-0.64 , w_a = -2.5^+1.9_), corresponding to 0.7σ, 0.3σ, 0.7σ (1.9σ, 3.4σ, 5.6σ, 5.5σ, 5.6σ) shifts between the MLE relative to the Bayesian posterior mean for ΛCDM (w_0w_aCDM) respectively.
The effective redshift distribution n(z) of galaxies is a critical component in the study of weak gravitational lensing. Here, we introduce a new method for determining n(z) for weak lensing surveys based on high-quality redshifts and neural-network-based importance weights. Additionally, we present the first unified photometric redshift calibration of the three leading stage-III weak lensing surveys, the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS), with state-of-the-art spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We verify our method using a new, data-driven approach and obtain n(z) constraints with statistical uncertainties of the order of sigma z & strns;similar to 0.01 and smaller. Our analysis is largely independent of previous photometric redshift calibrations and, thus, provides an important cross-check in light of recent cosmological tensions. Overall, we find excellent agreement with previously published results on the DES Y3 and HSC Y1 data sets, while there are some differences on the mean redshift with respect to the previously published KiDS-1000 results. We attribute the latter to mismatches in photometric noise properties in the COSMOS field compared to the wider KiDS self-organizing map-gold catalog. At the same time, the new n(z) estimates for KiDS do not significantly change estimates of cosmic structure growth from cosmic shear. Finally, we discuss how our method can be applied to future weak lensing calibrations with DESI data.
We present joint measurements of the pre- and post-reconstruction power spectra, P_ pre and P_ post, together with their cross-power spectrum, P_ cross, for the Luminous Red Galaxies (LRGs) in the DESI Data Release 1 (DR1). We jointly analyse these observables with an emulator-based full-shape modeling framework, thereby, for the first time, we extract complementary nonlinear information from the galaxy density field before and after reconstruction in real survey data. Specifically, including P_ post and P_ cross in addition to P_ pre (hereafter P_ all) yields an improvement of approximately 18-27% in the σ_8 constraint in both ΛCDM and wCDM, depending on the redshift bin, relative to the P_ pre-only analysis with the cosmic microwave background distance priors (hereafter CMB). In wCDM, the joint CMB+P_ all analysis can tighten the constraints on w by approximately 5-15% across the two LRG redshift bins, compared to the CMB+P_ pre combination. Further incorporating the Type Ia supernova dataset and comparing the cosmological constraints in wCDM from each individual power-spectrum component with those from the full combination, we find that P_ all consistently provides the tightest constraints. From the joint CMB+P_ all+DES-Dovekie dataset, we obtain Ω_m = 0.314 ± 0.0048 and w = -0.988 ± 0.023 for the sample, and Ω_m = 0.318 ± 0.0046 and w = -0.988 ± 0.025 for . These results demonstrate that combining pre- and post-reconstruction power spectra with their cross-correlation enables DESI to harvest additional nonlinear information, leading to tighter constraints on cosmological parameters.
The broad line region of quasars has long been probed by reverberation mapping techniques that measure time lags between continuum and broad emission line variations. Stacked reverberation mapping has been proposed as a less observationally expensive alternative to traditional methods. This ensemble approach also reduces biases from small-number statistics. The Dark Energy Spectroscopic Instrument (DESI) is conducting the most extensive spectroscopic survey of quasars to date. We create mock light curves emulating expected DESI quasar observations at redshifts 1.48<z<5.2 and luminosities 44.68 ≤log L_1350 λ/ erg s^-1≤ 45.99 to test stacked reverberation mapping feasibility using sparse spectroscopic data paired with well-sampled photometric data. The pipeline, using the lag estimation code JAVELIN, successfully recovers the simulated C IV lags within one sigma of the true values using spectroscopic light curves composed of only a few spectral epochs (2-10) with irregular cadences. We investigate how observational factors, including C IV flux error magnitude, number of stacked quasars, and spectral epoch count, affect performance. This work motivates a pathway for future stacked reverberation mapping projects with large scale spectroscopic surveys of quasars having ≥ 2 spectroscopic observations. Our results suggest an economical alternative for constraining and extending the radius-luminosity relation to higher redshifts and luminosities. Subsequently, this relation can be employed more reliably in single-epoch black hole mass measurements and quasar cosmology in these distant regimes.
Ly alpha emitters (LAEs) are valuable high-redshift cosmological probes traditionally identified using specialized narrowband photometric surveys. In ground-based spectroscopy, it can be difficult to distinguish the sharp LAE peak from residual sky emission lines using automated methods, leading to misclassified redshifts. We present a Bayesian spectral component separation technique to automatically determine spectroscopic redshifts for LAEs while marginalizing over sky residuals. We use visually inspected spectra of LAEs obtained using the Dark Energy Spectroscopic Instrument (DESI) to create a data-driven prior and can determine redshift by jointly inferring sky residual, LAE, and residual components for each individual spectrum. We demonstrate this method on 881 spectroscopically observed z = 2-4 DESI LAE candidate spectra and determine their redshifts with >90% accuracy when validated against visually inspected redshifts. Using the Delta chi(2) value from our pipeline as a proxy for detection confidence, we then explore potential survey design choices and implications for targeting LAEs with medium-band photometry. This method allows for scalability and accuracy in determining redshifts from DESI spectra, and the results provide recommendations for LAE targeting in anticipation of future high-redshift spectroscopic surveys.
We compute the 3D cross-correlation between the absorption of the $z\sim 2.3$ Lyman-alpha forest measured by the COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) survey, and 1642 foreground galaxies with spectroscopic redshifts from several different surveys, including 3D-HST, CLAMATO, zCOSMOS-deep, MOSDEF, and VUDS. For each survey, we compare the measured cross-correlation with models incorporating the galaxy linear bias as well as observed redshift dispersion and systematic redshift offset. The derived redshift dispersion and offsets are generally consistent with those expected from, e.g. spectroscopic redshifts measured with ultraviolet absorption lines or near-infrared emission lines observed with specific instruments, but we find hints of 'fingers-of-god' caused by overdensities in the field. We combine our foreground galaxy sample, and split them into thin bins of robustly estimated stellar mass in order to study the stellar mass-halo mass relationship. For subsamples with median stellar masses of $\log _{10}(M_* / \mathrm{M}_\odot) = [9.28,9.74,10.22]$, we find galaxy biases of $b_\mathrm{ g}\approx [2.1, 3.2,3.8]$, respectively. A comparison with mock measurements from the Bolshoi-Planck N-body simulation yields corresponding halo masses of $\log _{10}(M_* / \mathrm{M}_\odot) \approx [10.5,11.7,12.1]$ for these stellar mass bins. At the low-mass end, our results suggest enhanced star formation histories in mild tension with predictions from previous angular correlation and abundance matching-based observations, and the IllustrisTNG simulation.
We present a detailed study of the 3D morphology of two extended associations of multiple protoclusters at z=3.1. These protocluster `complexes', designated COSMOS-z3.1-A and COSMOS-z3.1-C, are the most prominent overdensities of z=3.1 Lyα emitters (LAEs) identified in the COSMOS field by the One-hundred-deg^2 DECam Imaging in Narrowbands (ODIN) survey. These protocluster complexes have been followed up with extensive spectroscopy from Keck, Gemini, and DESI. Using a probabilistic method that combines photometrically selected and spectroscopically confirmed LAEs, we reconstruct the 3D structure of these complexes on scales of ≈50 cMpc. We validate our reconstruction method using the IllustrisTNG300-1 cosmological hydrodynamical simulation and show that it consistently outperforms approaches relying solely on spectroscopic data. The resulting 3D maps reveal that both complexes are irregular and elongated along a single axis, emphasizing the impact of sightline on our perception of structure morphology. The complexes consist of multiple density peaks, ten in COSMOS-z3.1-A and four in COSMOS-z3.1-C. The former is confirmed to be a proto-supercluster, similar to Hyperion at z=2.4 but observed at an even earlier epoch. Multiple `tails' connected to the cores of the density peaks are seen, likely representing cosmic filaments feeding into these extremely overdense regions. The 3D reconstructions further provide strong evidence that Lyα blobs preferentially reside in the outskirts of the highest density regions. Descendant mass estimates of the density peaks suggest that COSMOS-z3.1-A and COSMOS-z3.1-C will evolve to become ultra-massive structures by z=0, with total masses log(M/M_⊙) ≳ 15.3, exceeding that of Coma.
Stellar streams are excellent probes of the gravitational potential in which they evolve. In the Milky Way (MW), globular cluster (GC) streams are routinely used to infer properties about time-dependent perturbations of the underlying potential. This implies that streams with Galactocentric radii small enough to be perturbed by the MW bar should offer constraints on it, such as its pattern speed, which currently has a wide range of values reported in the literature and is important when studying stellar kinematics. The GC M92 has a small pericentre and should be affected by the bar. It has a diffuse stellar stream, but confirming stream members has previously been hindered by a lack of spectroscopic data. In this paper, we use Dark Energy Spectroscopic Instrument (DESI) observations together with photometric and astrometric data to obtain spectroscopic members of the M92 stream for the first time. We identify a clear spatial distribution and gradients in distance moduli, proper motions, and radial velocities that confirm the stream's existence. We compare the observed stream to mock streams generated in different barred potentials and estimate the MW bar's pattern speed Ω= 29.1^+0.7_-0.4 km s^-1 kpc^-1 and Ω̇= 0.7^+3.5_-2.3 km s^-1 kpc^-1 Gyr^-1. This is the first time a stellar stream is used to probabilistically infer these bar properties, and it opens up an exciting realm of inner Galactic potential characterisation using stellar streams.
The standard cosmological analysis with the Lycx forest relies on a continuum fitting procedure that suppresses information on large scales and distorts the three-dimensional correlation function on all scales. In this work, we present the first cosmological forecasts without continuum fitting distortion in the Lycx forest, focusing on the recovery of large-scale information. Using idealized synthetic data, we compare the constraining power of the full shape of the Lycx forest auto-correlation and its cross-correlation with quasars using the baseline continuum fitting analysis versus the true continuum. We find that knowledge of the true continuum enables a similar to 10% reduction in uncertainties on the Alcock-Paczy & nacute;ski (AP) parameter and the matter density, ohm m. We also explore the impact of large-scale information by extending the analysis up to separations of 240 h-1Mpc along and across the line of sight. The combination of these analysis choices can recover significant large-scale information, yielding up to a similar to 15% improvement in AP constraints. This improvement is analogous to extending the Lycx forest survey area by similar to 40%.
Quasars acting as strong gravitational lenses offer a rare opportunity to probe the redshift evolution of scaling relations between supermassive black holes and their host galaxies, particularly the M _BH – M _host relation. Using these powerful probes, the mass of the host galaxy can be precisely inferred from the Einstein radius θ _E . Using 812,118 quasars from DESI DR1 (0.03 ≤ z ≤ 1.8), we searched for quasars lensing higher-redshift galaxies by identifying background emission-line features in their spectra. To detect these rare systems, we trained a convolutional neural network (CNN) on mock lenses constructed from real DESI spectra of quasars and emission-line galaxies (ELGs), achieving a high classification performance (AUC = 0.99). We also trained a regression network to estimate the redshift of the background ELG. Applying this pipeline, we identified seven high-quality (Grade A) lens candidates, each exhibiting a strong [O II ] doublet at a higher redshift than the foreground quasar; four candidates additionally show H β , [O III ] λ 4959, and [O III ] λ 5007 emission. These results significantly expand the sample of quasar lens candidates beyond the 12 identified and 3 confirmed in previous work and demonstrate the potential for scalable, data-driven discovery of quasars as strong lenses in upcoming spectroscopic surveys.
We present the first measurement of local-type primordial non-Gaussianity from the cross-correlation between 1.2 million spectroscopically confirmed quasars from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI) and the Planck PR4 CMB lensing reconstructions. The analysis is performed in three tomographic redshift bins covering 0.8 < z < 3.5, covering a sky fraction of similar to 20%. We adopt a catalog-based pseudo-C-l estimator and apply linear imaging weights validated on noiseless mocks. Compared to previous analyses using photometric quasar samples, our results benefit from the high purity of the DESI spectroscopic sample, the reduced noise of PR4 lensing, and the absence of excess large-scale power in the spectroscopic quasar auto-correlation. Fitting simultaneously for the non-Gaussianity parameter f(NL) and the linear bias amplitude in each redshift bin, we obtain f(NL) = 2(-34)(+28) for a response parameter p = 1.6, and f(NL) = 6(-24)(+20) for p = 1.0. These results improve the constraints on f(NL) by similar to 35% compared to the previous analysis based on the Legacy Imaging Survey DR9. Additionally, we derive an optimal weighting scheme to maximize the constraining power. In this case, and assuming p = 1.6, we obtain f(NL) = 19(-31)(+25). Our results demonstrate the statistical power of DESI quasars for probing inflationary physics, and highlight the promise of future DESI data releases.
Abstract The One-hundred-deg 2 DECam Imaging in Narrowbands (ODIN) survey is conducting the widest-field deep narrowband (NB) imaging of the equatorial and southern skies. ODIN uses three custom-built NB filters that sample Ly α -emitting galaxies (LAEs) within thin cosmic slices centered at z = 2.4, 3.1, and 4.5. In this work, we utilize extensive DESI spectroscopy of ODIN-selected galaxies in the COSMOS and XMM-LSS fields to validate our LAE selection. Exposures of 2–4 hr with DESI yielded redshift confirmation of 3075 ODIN LAE candidates with NB magnitudes brighter than 26 mag. Restricting to objects that yield high-confidence redshifts, the confirmation rates are (93%, 96%, 92%) at z = (2.4, 3.1, 4.5). The primary contaminants consist of active galactic nuclei at the expected Ly α- redshift range and lower redshifts (C iv , C iii ]), with the remainder being star-forming galaxies ([O ii ] and [O iii ]). We find minimal contamination from [O ii ] emitters in our sample (≲1%), implying that our rest-frame equivalent width (REW) > 20 Å NB excess photometry requirement is sufficient to remove them.
Protoclusters represent sites of accelerated galaxy formation and extreme astrophysical activity characteristic of dense environments. Identifying massive protoclusters and mapping their spatial structures are therefore crucial first steps in understanding how the large-scale environment influences galaxy evolution. We combine wide-field Lyα imaging from the ODIN survey with extensive DESI and ancillary spectroscopy across the extended COSMOS and XMM-LSS fields (≈14 deg^2) to search for massive protoclusters. We confirm six systems at z≈ 2.4 and z≈ 3.1, reconstruct their three-dimensional structures, estimate descendant halo masses, and, for one structure at z≈ 3.12, demonstrate that overlapping narrowband filters (NB497 and N501) provide accurate redshift tomography for emission-line galaxies. One protocluster at z≈ 2.45 overlaps with one of the LATIS tomographic fields, enabling direct comparison between galaxy and H i overdensities traced by Lyα forest absorption. Another at z≈ 3.12 hosts a massive quiescent galaxy (M_∗≈ 1.2 × 10^11M_⊙), indicating early quenching in a dense environment. By comparing Lyα emission properties across environments, we find that protocluster galaxies exhibit higher median line fluxes and a deficit of faint emitters relative to the field. The effect is strongest when both 2D and 3D density information are combined, indicating that galaxies in the densest protocluster cores are most affected by environmental processes. This effect is stronger at z≈3.1 than at z≈2.4, suggesting possible redshift evolution.
We present a study of the evolution of type 1 quasars at 1.0<z<3.5, covering the peak epoch of quasar activity. The quasar evolution has been extensively explored by a variety of previous works and the derived quasar luminosity functions (QLFs) are not well consistent with each other, presumably due to the complexities introduced by different quasar selection techniques and associated completeness corrections. We use a new strategy to construct QLFs based on a library of all known quasars. We focus on a wide region of ∼1700 deg^2 and a deep field of ∼265 deg^2 that have rich spectroscopic data primarily from SDSS and DESI. We then apply traditional color cuts in the rest-frame UV/optical to select quasar candidates and use the quasar library to identify them. Our final sample consists of 62,426 quasars at 1.0<z<3.5, with a high completeness (∼96
On 2025 August 18, the LIGO-Virgo-KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after similar to 5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of z = 0.084840 +/- 0.000006. If S250818k has an astrophysical origin, this places the transient within 2 sigma of the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of log10I approximate to 3.9-4.2 . Our analysis of the host galaxy's spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass similar to 1010 M circle dot, broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI's extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.