Engineering topics

BenZvi, S. (ORCID:0000000155374710)

Publications and source records attributed to BenZvi, S. (ORCID:0000000155374710).

At least 19 records

DOE OSTI2026

Extensive analysis of reconstruction algorithms for DESI 2024 baryon acoustic oscillations

Reconstruction of the baryon acoustic oscillation (BAO) signal has been a standard procedure in BAO analyses over the past decade and has helped to improve the BAO parameter precision by a factor of ∼2 on average. The Dark Energy Spectroscopic Instrument (DESI) BAO analysis for the first year (DR1) data uses the “standard” reconstruction framework, in which the displacement field is estimated from the observed density field by solving the linearized continuity equation in redshift space, and galaxy and random positions are shifted in order to partially remove non-linearities. There are several approaches to solving for the displacement field in real survey data, including the multigrid (MG), iterative Fast Fourier Transform (iFFT), and iterative Fast Fourier Transform particle (iFFTP) algorithms. In this work, we analyze these algorithms and compare them with various metrics including two-point statistics and the displacement itself using realistic DESI mocks. We focus on three representative DESI samples, the emission line galaxies (ELG), quasars (QSO), and the bright galaxy sample (BGS), which cover the extreme redshifts and number densities, and potential wide-angle effects. We conclude that the MG and iFFT algorithms agree within 0.4% in post-reconstruction power spectrum on BAO scales with the RecSym convention, which does not remove large-scale redshift space distortions (RSDs), in all three tracers. The RecSym convention appears to be less sensitive to displacement errors than the RecIso convention, which attempts to remove large-scale RSDs. However, iFFTP deviates from the first two; thus, we recommend against using iFFTP without further development. In addition, we provide the optimal settings for reconstruction for five years of DESI observation. The analyses presented in this work pave the way for DESI DR1 analysis as well as future BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

DESI DR1 Ly α 1D power spectrum: Validation of estimators

The Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) is the largest sample to date for small-scale Lyα forest cosmology, accessed through its one-dimensional power spectrum (P 1D ). The Lyα forest P 1D is extracted from quasar spectra that are highly inhomogeneous (both in wavelength and between quasars) in noise properties due to intrinsic properties of the quasar, atmospheric and astrophysical contamination, and also sensitive to low-level details of the spectral extraction pipeline. We employ two estimators in DR1 analysis to measure P 1D : the optimal estimator and the fast Fourier transform (FFT) estimator. To ensure robustness of our DR1 measurements, we validate these two power spectrum and covariance matrix estimation methodologies against the challenging aspects of the data. First, using a set of 20 synthetic 1D realizations of DR1, we derive the masking bias corrections needed for the FFT estimator and the continuum fitting bias needed for both estimators. We demonstrate that both estimators, including their covariances, are unbiased with these corrections using the Kolmogorov-Smirnov test. Second, we substantially extend our previous suite of CCD image simulations to include 675,000 quasars, allowing us to accurately quantify the pipeline's performance. This set of simulations reveals biases at the highest k values, corresponding to a resolution error of a few percent. We base the resolution systematics error budget of DR1 P 1D on these values, but do not derive corrections from them since the simulation fidelity is insufficient for precise corrections.

Lyman alpha forest↗

DOE OSTI2026

DESI Strong Lens Foundry. III. Keck Spectroscopy for Strong Lenses Discovered Using Residual Neural Networks

We present spectroscopic data of strong lenses and their source galaxies using the Keck Near-Infrared Echellette Spectrometer (NIRES) and the Dark Energy Spectroscopic Instrument (DESI), providing redshifts necessary for nearly all strong-lensing applications with these systems, especially the extraction of physical parameters from lensing modeling. These strong lenses were found in the DESI Legacy Imaging Surveys using residual neural networks and followed up by our Hubble Space Telescope program, with all systems displaying unambiguous lensed arcs. With NIRES, we target eight lensed sources at redshifts difficult to measure in the optical range and determine the source redshifts for six, between z s = 1.675 and 3.332. DESI observed one of the remaining source redshifts, as well as an additional source redshift within the six systems. The two systems with nondetections by NIRES were observed for a considerably shorter 600 s at high airmass. Combining NIRES infrared spectroscopy with optical spectroscopy from our DESI Strong Lensing Secondary Target Program, these results provide the complete lens and source redshifts for six systems, a resource for refining automated strong lens searches in future deep- and wide-field imaging surveys and addressing a range of questions in astrophysics and cosmology.

Agarwal, Shrihan [University of Chicago, IL (Unite↗

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

The DESI DR1 peculiar velocity survey: growth rate measurements from the maximum likelihood fields method

We present the constraint on the growth rate of structure from the combination of DESI DR1 BGS sample, Fundamental Plane, and Tully-Fisher peculiar velocity catalogues using the maximum likelihood fields method. The combined catalogue contains 415,523 galaxy redshifts and 76,616 peculiar velocity measurements. To handle the large amount of data in the DESI DR1 peculiar velocity catalogue, we significantly improve the computational efficiency by rewriting the algorithm with JAX. After removing outliers and Tully-Fisher galaxies that are affected by systematics, we find fσ 8 = 0.483 -0.043 +0.080 (stat) ± 0.018(sys), consistent within 1σ with the power spectrum and correlation function analysis using the same dataset. Combining all three measurements with appropriate correlations, the consensus measurement is fσ 8 (z eff = 0.07) = 0.450±0.055, consistent with Planck +ΛCDM cosmology (fσ 8 = 0.449±0.008). Combining with the high redshift growth rate of structure measurements from DESI ShapeFit, the constraint on the growth index is γ = 0.58±0.11, consistent with GR.

cosmic flows↗

DOE OSTI2026

H 0 without the sound horizon (or supernovae): A 2% measurement in DESI DR1

The sound horizon scale rs is a key source of information for measurements of H 0 from early-time data, and is therefore a common target of new physics proposed to solve the Hubble tension. We present a sub-2% measurement of the Hubble constant that is independent of this scale, using data from the first data release of the Dark Energy Spectroscopic Instrument (DESI DR1). Building on previous work, we remove dependency on the sound horizon size using a heuristic rescaling procedure at the power spectrum level. A key innovation is the inclusion of uncalibrated (agnostic to r s ) post-reconstruction BAO measurements from DESI DR1, as well as using the CMB acoustic scale θ* as a high-redshift anchor. Uncalibrated type-Ia supernovae are often included as an independent source of Ω m information; here we demonstrate the robustness of our results by additionally considering two supernova-independent alternative datasets. We find somewhat higher values of H 0 relative to our previous work: 69.2 +1.3 -1.4 , 70.3 +1.4 -1.2 , and 69.6 +1.3 -1.8 km s -1 Mpc -1 respectively when including measurements from i) Planck/ACT CMB lensing × unWISE galaxies, ii) the DES Year 3 6×2pt analysis, and iii) Planck/ACT CMB lensing + the DES Year 5 supernova analysis. These remarkably consistent constraints achieve better than 2% precision; they are among the most stringent sound horizon-independent measurements from LSS to date, and provide a powerful avenue for probing the origin of the Hubble tension.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

DESI EDR: Calibrating the Tully–Fisher Relationship with the DESI Peculiar Velocity Survey

We calibrate the Tully–Fisher relation (TFR) with data from the DESI Peculiar Velocity (PV) Survey taken during the Survey Validation (SV) period of the DESI galaxy redshift survey. Placing spectroscopic fibers on the centers and major axes of spatially extended spiral galaxies identified in the 2020 Siena Galaxy Atlas using the DESI Legacy Surveys, we measure the rotational velocities at 0.33R 26 for 1155 (1128 + 27 dwarf) spiral galaxies observed during SV. Using 39 spiral galaxies observed in the Coma cluster, we find a slope for the TFR of −8.32 ± 0.15 AB mag in the r band, with a scatter about the TFR of 1.12 ± 0.03 AB mag. We calibrate the zero-point of the TFR using galaxies with independent distances measured using type Ia supernovae (SNe Ia) via the cosmological distance ladder. From the SN Ia distances, we measure a zero-point of $-19.21^{+0.30}_{-0.31}$ AB mag in the r band. We produce a public catalog of the distances to these 1128 spiral galaxies observed during DESI SV as part of the DESI PV Survey with our calibrated TFR. This is, to our knowledge, the first catalog of TFR distances produced with velocities measured at a single point in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Evidence for Neutrino Emission from X-Ray Bright Seyfert Galaxies in the Southern Hemisphere Using Enhanced Starting Track Events with IceCube

IceCube recently reported the observation of TeV neutrinos from the nearby Seyfert galaxy NGC 1068, and the corresponding neutrino flux is significantly higher than the upper limit implied by observations of GeV–TeV gamma rays. This suggests that neutrinos are produced near the supermassive black hole, where the radiation density is high enough to obscure gamma rays. We use a set of muon neutrinos with interaction vertices inside the detector, which have good sensitivity to sources in the southern sky, from IceCube data recorded between 2011 and 2021. We then search for individual and collective neutrino signals from 14 Seyfert galaxies in the southern sky selected from the Swift Burst Alert Telescope AGN Spectroscopic Survey. Using the correlations between keV X-rays and TeV neutrinos predicted by disk–corona models, and assuming production characteristics similar to NGC 1068, a collective neutrino signal search reveals an excess of $6.{7}_{-3.2}^{+4.0}$ events, which is inconsistent with background expectations at the 3σ level of significance. In this Letter, we present new independent evidence that Seyfert galaxies contribute to the extragalactic flux of high-energy neutrinos.

Seyfert galaxies↗

DOE OSTI2026

Evidence for Neutrino Emission from X-Ray-bright Active Galactic Nuclei with IceCube

Recently, IceCube reported neutrino emission from the Seyfert galaxy NGC 1068. Using 13.1 yr of IceCube data, we present a follow-up search for neutrino sources in the northern sky. NGC 1068 remains the most significant neutrino source among 110 preselected gamma-ray emitters while also being spatially compatible with the most significant location in the northern sky. Its energy spectrum is characterized by an unbroken power-law with spectral index γ = 3.4 ± 0.2. Consistent with previous results, the observed neutrino flux exceeds its gamma-ray counterpart by at least 2 orders of magnitude. Motivated by this disparity and the high X-ray luminosity of the source, we selected 47 X-ray-bright Seyfert galaxies from the Swift/BAT spectroscopic survey that were not included in the list of gamma-ray emitters. When testing this collection for neutrino emission, we observe a 3.3σ excess from an ensemble of 11 sources, with NGC 1068 excluded from the sample. Our results strengthen the evidence that X-ray-bright cores of active galactic nuclei are neutrino emitters.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2026

Constraints on the Correlation of IceCube Neutrinos with a Tracer of Nearby Large-scale Structure

The IceCube Neutrino Observatory has observed extragalactic astrophysical neutrinos with an apparently isotropic distribution. Only a small fraction of the observed astrophysical neutrinos can be explained by known sources. Neutrino production is thought to occur in energetic environments that are ultimately powered by the gravitational collapse of dense regions of the large-scale mass distribution in the universe. Whatever their identity, neutrino sources likely trace this large-scale mass distribution. The clustering of neutrinos with a tracer of the large-scale structure may provide insight into the distribution of neutrino sources with respect to redshift and the identity of neutrino sources. We implement a two-point angular cross correlation of the Northern sky track events with an infrared galaxy catalog derived from the Wide-field Infrared Survey Explorer (WISE) and Two Micron All Sky Survey (2MASS) source catalogs, which trace the nearby large-scale structure. No statistically significant correlation is found between the neutrinos and this infrared galaxy catalog. We find that ≤54% of the diffuse muon neutrino flux can be attributed to sources correlated with the galaxy catalog with 90% confidence. Additionally, when assuming that the neutrino source comoving density evolves following a power law in redshift, dN s /dV ∝ (1 + z) k , we find that sources with negative evolution, in particular k < −1.75, are disfavored at the 90% confidence level.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2026

A Search for Millimeter-bright Blazars as Astrophysical Neutrino Sources

The powerful jets of blazars have been historically considered as likely sites of high-energy cosmic-ray acceleration. However, the particulars of the launched jet and the locations of leptonic and hadronic jet loading remain unclear. In the case when leptonic and hadronic particle injection occur jointly, a temporal correlation between synchrotron radiation and neutrino production is expected. We use a first catalog of millimeter wavelength (95–225 GHz) blazar light curves from the Atacama Cosmology Telescope for a time-dependent correlation with 12 yr of muon neutrino events from the IceCube South Pole Neutrino Observatory. Such millimeter emission traces activity of the bright jet base, which is often self-absorbed at lower frequencies and potentially gamma-ray opaque. We perform an analysis of the population, as well as analyses of individual, selected sources. We do not observe a significant signal from the stacked population. TXS 0506+056 is found as the most significant, individual source, though this detection is not globally significant in our analysis of selected active galactic nuclei. Our results suggest that the majority of millimeter-bright blazars are neutrino dim. In general, it is possible that many blazars have lighter, leptonic jets, or that only selected blazars provide exceptional conditions for neutrino production.

Abbasi, R. [Loyola University, Chicago, IL (United↗

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 OSTI2026

Time-integrated Southern-sky Neutrino Source Searches with 10 yr of IceCube Starting-track Events at Energies Down to 1 TeV

In the IceCube Neutrino Observatory, a signal of astrophysical neutrinos is obscured by backgrounds from atmospheric neutrinos and muons produced in cosmic-ray interactions. IceCube event selections used to isolate the astrophysical neutrino signal often focus on the morphology of the light patterns recorded by the detector. The analyses presented here use the new IceCube Enhanced Starting Track Event Selection (ESTES), which identifies events likely generated by muon–neutrino interactions within the detector geometry, focusing on neutrino energies of 1–500 TeV with a median angular resolution of 1.4°. Selecting for starting-track events filters out not only the atmospheric-muon background but also the atmospheric-neutrino background in the southern sky. This improves IceCube’s muon–neutrino sensitivity to southern-sky neutrino sources, especially for Galactic sources that are not expected to produce a substantial flux of neutrinos above 100 TeV. In this work, the ESTES sample was applied for the first time to search for astrophysical sources of neutrinos, including a search for diffuse neutrino emission from the Galactic plane. No significant excesses were identified from any of the analyses; however, constraining limits are set on the hadronic emission from TeV gamma-ray Galactic plane objects and models of the diffuse Galactic plane neutrino flux.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2025

All-sky Neutrino Point-source Search with IceCube Combined Track and Cascade Data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 2008 April 6 to 2022 May 23 to search the whole sky for neutrino sources, and using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing power and sensitivity in the northern sky, with all-sky cascades, with good energy resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube that targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5σ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest-significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068’s emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2025

A joint analysis of 3D clustering and galaxy × CMB-lensing cross-correlations with DESI DR1 galaxies

The spectroscopic data from DESI Data Release 1 (DR1) galaxies enables the analysis of 3D clustering by fitting galaxy power spectra and reconstructed correlation functions in redshift space. Given low measurements of the amplitude of structure from cosmic shear at z ∼ 1, redshift space distortions (RSD) + Baryon Acoustic Oscillation (BAO) signals from DESI galaxies combined with weak lensing can break degeneracies and provide a tight alternative constraint on the z ∼ 1 amplitude of structure. In this paper we perform joint analyses that combine full-shape + post-reconstruction information from the DESI DR1 BGS and LRG samples along with angular cross-correlations with Planck PR4 and ACT DR6 CMB lensing maps. We show that adding galaxy-lensing cross-correlations tightens clustering amplitude constraints, improving σ 8 uncertainties by 30% over RSD+BAO alone. We also include angular galaxy-galaxy and galaxy-lensing spectra using photometric samples from the DESI Legacy Survey to further improve constraints. Our headline results are σ 8 = 0.803 ± 0.017, Ω m = 0.3037 ± 0.0069, and S 8 = 0.808 ± 0.017. Given DESI's preference for higher σ 8 compared to lower values from BOSS, we perform a catalog-level comparison of LRG samples from both surveys. We test sensitivity to dark energy assumptions by relaxing our ΛCDM prior and allowing for evolving dark energy via the w 0 - w a parameterization. We find our S 8 constraints to be relatively unchanged despite a 3.5σ tension with the cosmological constant model when combining with the Union3 supernova likelihood. Finally we test general relativity (GR) by allowing the gravitational slip parameter (γ) to vary, and find γ = 1.17 ± 0.11 in mild (∼ 1.5σ) tension with the GR value of 1.0.

cosmological parameters from LSS↗

DOE OSTI2025

DESI Emission-line Galaxies: Clustering Dependence on Stellar Mass and [O II] Luminosity

We measure the projected two-point correlation functions of emission-line galaxies (ELGs) from the Dark Energy Spectroscopic Instrument One-Percent Survey and model their dependence on stellar mass and [O II] luminosity. We select ∼180,000 ELGs with redshifts of 0.8 < z < 1.6, and define 27 samples according to cuts in redshift and both galaxy properties. Following a framework that describes the conditional [O II] luminosity–stellar mass distribution as a function of halo mass, we simultaneously model the clustering measurements of all samples at fixed redshift. Based on the modeling result, most ELGs in our samples are classified as central galaxies, residing in halos of a narrow mass range with a typical median of ∼10 12.2−12.4 h −1 M ⊙ . We observe a weak dependence of clustering amplitude on stellar mass, which is reflected in the model constraints and is likely a consequence of the 0.5 dex measurement uncertainty in the stellar mass estimates. The model shows a trend between galaxy bias and [O II] luminosity at high redshift (1.2 < z < 1.6) that is otherwise absent at lower redshifts.

Cosmology↗

DOE OSTI2025

Missing components in ΛCDM from DESI Y1 baryonic acoustic oscillation measurements: Insights from redshift remapping

We explore transformations of the Friedman-Lemaître-Robertson-Walker (FLRW) metric and cosmological parameters that align with observational data while aiming to gain insights into potential extensions of standard cosmological models. We modified the FLRW metric by introducing a scaling factor, e 2Θ(a) –the cosmological scaling function (CSF), which alters the standard relationship between cosmological redshift and the cosmic scale factor without affecting angular measurements or cosmic microwave background (CMB) anisotropies. Using data from DESI Year 1, Pantheon+ supernovae, and the Planck CMB temperature power spectrum, we constrained both the CSF and cosmological parameters through a Markov chain Monte Carlo approach. Our results indicate that the CSF model fits observational data with a lower Hubble constant (although it is compatible with the value given by Planck 2018 within 1σ) and is predominantly dark matter dominated. Additionally, the CSF model produces temperature and lensing power spectra similar to those predicted by the standard model, though with lower values in the CSF model at large scales. We also checked that when fitting a CSF model without dark energy to the data, we obtain a more negative conformal function. This suggests that the CSF model may offer hints about missing elements and opens up a new avenue for exploring physical interpretations of cosmic acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Erratum: “Tripling the Census of Dwarf AGN Candidates Using DESI Early Data” (2025, ApJ, 982, 10)

We identified an error in the EmFit code, which was employed to estimate the emission-line fluxes and widths used in the published article. Specifically, the error spectrum was not corrected for Galactic reddening before being incorporated into the fitting process. Although this does not affect the measured fluxes or line widths, it led to underestimated uncertainties. Consequently, the signal-to-noise ratios (SNRs) for individual sources were systematically overestimated, resulting in a larger number of sources satisfying a given SNR threshold.

Pucha, Ragadeepika [University of Utah, Salt Lake ↗