Engineering topics

Weaver, B. A. [Natl. Solar Observ., Tucson]

Publications and source records attributed to Weaver, B. A. [Natl. Solar Observ., Tucson].

At least 19 records

DOE OSTI2026

A Unified Tracer Analysis of DESI DR2 Baryon Acoustic Oscillations

We improve upon previous efforts to optimally combine overlapping galaxy samples in the DESI baryon acoustic oscillation analysis. By weighting each galaxy by its linear bias, overlapping galaxies are combined into a single, unified catalog, naturally avoiding double counting of cosmic volume and including all auto- and cross- information at the catalog level. Improvements over the previous effort include the addition of QSO out to $z=1.6$ to account for all overlapping DR2 tracers and redshift-dependent bias treatment to improve reconstruction. We report distance measurements using this unified tracer, and find them to be highly consistent with the baseline DR2 BAO analysis. We also test for tracer-dependent systematics within the DESI data, and find no evidence of tracer-dependent systematics within $0.8

Sanders, N. [Ohio U., Athens] (ORCID:0009000800202↗

DOE OSTI2026

PAC in DESI. II. Galaxy-halo connection into the $10^{6}{\rm M}_{\odot}$ frontier

Understanding dwarf galaxy formation is crucial for testing dark matter models and reionization physics. However, constructing stellar-mass complete spectroscopic samples at low masses is increasingly difficult, and the potential existence of a local void complicates studies in an average environment. The Photometric object Around Cosmic webs (PAC) method, which combines deep photometric and spectroscopic data to measure the excess surface density $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ of photometric objects around spectroscopic tracers, offers a promising path forward. We model 349 $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ measurements from DESI Y1 BGS and DECaLS, reaching $M_*=10^{6.4}\,{\rm M}_{\odot}$, using a stellar mass-halo mass relation (SHMR)-based subhalo abundance matching framework applied to two high-resolution $N$-body simulations from the Jiutian suite. The resulting SHMR is constrained down to $M_{\rm h}\simeq10^{8.0}\,h^{-1}{\rm M}_{\odot}$, revealing a clear upturn at $\sim10^{10.0}\,h^{-1}{\rm M}_{\odot}$ toward lower masses, indicating rising star-formation efficiency (SFE) in small haloes. This feature persists under extensions of the model that allow mass-dependent scatter, reionization-induced suppression of the halo occupation fraction, galaxy assembly bias, and alternative cosmologies. Combining with the results from Paper I, we find that central red galaxies dominate the low-mass regime. Our results motivate a hypothesis in which SFE is significantly higher than previously thought prior to reionization, enabling relatively massive galaxies to form in small haloes. These systems are subsequently quenched by the UV background, producing the central red dwarf galaxies observed. Finally, we obtain $3σ$ and $5σ$ upper mass bounds of $10^{8.80}\,h^{-1}{\rm M}_{\odot}$ and $10^{10.24}\,h^{-1}{\rm M}_{\odot}$ on the smallest haloes required to exist.

Xu, Kun [Pennsylvania U.; Durham U., ICC; Tsung-Da↗

DOE OSTI2026

ODIN: Spectroscopic Validation of Lyα-emitting Galaxy Samples with DESI

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.

Pinarski, Ethan [Purdue U., West Lafayette] (ORCID↗

DOE OSTI2026

Validation of the DESI DR2 Ly$α$ forest full-shape analysis

We present the validation of the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) Lyman-$α$ (Ly$α$) forest full-shape analysis. This analysis combines three-dimensional Ly$α$ forest auto-correlations and cross-correlations with quasars to extract information from both the baryon acoustic oscillation (BAO) feature and the broadband clustering signal, with primary emphasis on the Alcock-Paczynski (AP) measurement. Compared to the DESI DR1 analysis, the DR2 validation uses substantially larger and more realistic mock datasets, including CoLoRe 2LPT and AbacusSummit Ly$α$ forest simulations. The modeling framework is also improved through analytic marginalization over small scales ($<10$$h^{-1}$Mpc) and the impact of ultraviolet background fluctuations. The validation program was completed prior to unblinding and defines quantitative requirements for the cosmological parameters of interest, which are evaluated using hundreds of mock realizations. We further test the analysis through independent fits to the auto- and cross-correlations, multiple catalog splits, and a broad suite of analysis and modeling variations applied to both mocks and blinded observational data. We find that the BAO and AP parameters satisfy all validation requirements and remain stable across all tests. In contrast, mock studies reveal a significant bias in the inferred growth-rate parameter $fσ_8$, leading us to exclude this measurement from the final analysis. The consistency across mocks, data splits, and robustness tests demonstrates that the DR2 Ly$α$ full-shape analysis provides a reliable and substantially improved broadband AP measurement over previous Ly$α$ forest studies.

Herbold, M. [Chicago U., KICP; Ohio State U.] (ORC↗

DOE OSTI2026

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗

DOE OSTI2026

Lyman-$α$ forest holography: 3D predictions from 1D measurements

Cosmological analyses of Lyman-$α$ forest clustering rely on either one-dimensional correlations along individual sightlines or three-dimensional correlations between different sightlines. Because these observables probe the matter distribution on very different scales, they have traditionally been analyzed independently. In this work, we bridge this gap using ForestFlow, an emulator trained on a suite of cosmological hydrodynamical simulations that provides a unified description of Lyman-$α$ forest clustering from linear to nonlinear scales. This framework enables us to determine the range of three-dimensional clustering models compatible with the DESI one-dimensional flux power spectrum ($P_{\rm 1D}$). The resulting predictions successfully reproduce the large-scale clustering measured by the DESI BAO analysis and provide physically motivated priors on nonlinear clustering that are used in a companion paper presenting the full-shape analysis of the DESI DR2 Lyman-$α$ forest. We validate our methodology using the large-volume, high-resolution hydrodynamical simulation ACCEL-2, demonstrating excellent agreement across the full range of scales considered. Finally, we combine constraints from the $P_{\rm 1D}$ and BAO analyses on the parameter combinations $b_δσ_8$ and $b_ηf σ_8$, finding that the two probes provide comparable constraining power while exhibiting complementary parameter degeneracies. Our results establish a direct connection between one- and three-dimensional Lyman-$α$ forest measurements through ForestFlow, an approach we term Lyman-$α$ holography by analogy with the reconstruction of higher-dimensional structure from lower-dimensional information.

Chaves-Montero, J. [Barcelona, IFAE] (ORCID:000000↗

DOE OSTI2026

Electron temperature relations and the direct N, O, Ne, S, and Ar abundances of 49 959 star-forming galaxies in DESI data release 2

We present the largest direct-method abundance catalogue of galaxies to date, containing measurements of 49 959 star-forming galaxies at z<0.96 from DESI (Dark Energy Spectroscopic Instrument) data release 2. By directly measuring electron temperatures across multiple ionization zones, we provide constraints on a number of electron temperature relations. Using the temperature measurements, we derive reliable abundances for N, O, Ne, S, and Ar, and measure the evolution of abundances and abundance ratios of as a function of metallicity and other galaxy properties. Our measurements include direct oxygen abundances for 49 507 galaxies, leading to the discovery of the two most metal-poor galaxies in the nearby Universe, with oxygen abundances of 12+log⁡(O/H)=6.77−0.03+0.03 dex (1.2 per cent Z⊙⁠) and 12+log⁡(O/H)=6.81−0.04+0.04 dex (1.3 per cent Z⊙⁠). We identify a rare outlier population of 24 galaxies with high-N/O ratios at low metallicity, reminiscent of galaxy abundances observed in the early Universe. We find the Ne/O ratio is constant at low metallicity but increases gradually at 12+log(O/H)>8.105±0.004 dex. We show that the S/O and Ar/O abundance ratios are strongly correlated, consistent with the expected additional Type Ia enrichment channel for S and Ar. In this work, we present an initial survey of the key properties of the sample, with this data set serving as a foundation for extensive future work on galaxy abundances at low redshift.

Scholte, D. [Edinburgh U., Inst. Astron.] (ORCID:0↗

DOE OSTI2026

Probing the matter-dominated expansion with multi-redshift Lyman-$α$ BAO from DESI DR2

We present a multi-redshift Baryon Acoustic Oscillations (BAO) analysis of the DESI Data Release 2 (DR2) Lyman-$α$ (Ly$α$) forest, splitting the forest auto-correlation and its cross-correlation with quasars into three redshift bins. We obtain BAO measurements at effective redshifts $z_{\rm eff} = 2.13$, $2.40$, and $2.81$ with $\sim2.0$--$2.5\%$ precision per bin in the radial and transverse directions, corresponding to $\sim1.1$--$1.2\%$ precision for the isotropic BAO measurement. Using the same data products and modeling framework as the DESI DR2 Ly$α$ BAO analysis, we validate the pipeline on $400$ synthetic datasets and find unbiased BAO recovery with well-calibrated uncertainties. The measurements show an increase in the isotropic dilation parameter $D_V/r_d$ from $30.26\pm0.39$ to $32.22\pm0.47$ and in the Alcock-Paczyński parameter $D_M/D_H$ from $3.96\pm0.15$ to $5.63^{+0.22}_{-0.24}$. The Hubble distance $D_H/r_d$ decreases from $9.40\pm0.20$ to $7.22\pm0.17$, providing a direct measurement of the expansion history consistent with $Λ$CDM and the expected matter-dominated scaling, with $H(z)\propto(1+z)^n$ giving $n=1.34\pm0.16$. The redshift split also provides a self-consistent measurement of clustering evolution: the Ly$α$ forest bias evolves as $(1+z)^γ$ with $γ_α=3.05\pm0.16$, the RSD parameter has a redshift evolution described by $γ_β=-0.97\pm0.26$, and the quasar bias evolves with $γ_Q=1.56\pm0.23$, consistent with independent quasar clustering measurements. Combining these three-bin BAO measurements with DESI DR2 galaxy and quasar BAO measurements yields cosmological constraints consistent with the single-bin Ly$α$ BAO analysis in flat $Λ$CDM and $w_0w_a$CDM and improves curvature constraints by $\sim12\%$ in $Λ$CDM$+Ω_\mathrm{K}$.

Herrera-Alcantar, Hiram K. [Paris, Inst. Astrophys↗

DOE OSTI2026

Alcock-Paczynski Blinding Scheme for the Ly-$α$ Forest Analysis

We present and validate a blinding method for the Lyman-$α$ (Ly$α$) forest analysis based on a modification of the Alcock-Paczynski test. In order to hide the background expansion history, the method employs a geometrical shift of each quasar (QSO) forest in wavelength space, once the quasar continuum has been fitted and the fluctuation field is extracted. The redshift positions for the QSO sample are also changed in a consistent manner. We show that the method remains effective when applied to real data, where contamination from metals and Lyman-$β$ is intrinsically mixed with the Lyman-$α$ forest. This limitation is primarily visible in the 1D correlation function, where other blinding strategies can mitigate the effect. To assess its effectiveness, the prescription is tested against a series of datasets of increasing complexity: from idealized low-noise mocks, to realistic DESI year one synthetic datasets, and finally to data from DESI first data release (DR1), using both the auto (Ly$α\times$Ly$α$) and cross (Ly$α\times$ QSO) correlations. We find that the method robustly shifts the BAO peak position from the 3D correlation functions to the expected value for cosmology changes of around 5% in the matter content, without altering the shape of the posteriors in the model parameters. In conclusion, this catalog-level blinding strategy is a viable method for cosmological inference with the Lyman-$α$ forest, particularly if a cross-analysis with other tracers using the same blinding strategy is pursued.

Perez-Sanchez, G. [Guanajuato U.] (ORCID:000900096↗

DOE OSTI2026

Weak Evolution of Cosmic Atomic Hydrogen over the Past 4.5 Billion Years

The cosmic star formation rate density (CSFRD) has declined sharply toward the present day, but the roles of the atomic and molecular gas reservoirs remain uncertain. We measure the cosmic HI density, $Ω_{\mathrm{HI}}$, over $0

Zhang, Chuan-Peng [NAOC, Beijing; Guizhou Norm. U.↗

DOE OSTI2026

DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs

High-redshift galaxy clustering provides a powerful probe of the growth of structure, testing models of dark matter, dark energy, and galaxy formation during the epoch when the Universe was rapidly evolving. Emission line galaxies (ELGs) and quasars (QSOs) are used as tracers of dark matter by the Dark Energy Spectroscopic Instrument (DESI) to probe this redshift regime. We present results from ELG and QSO mock catalogs created from the Uchuu N-body simulation and tuned to DESI Data Release 2 (DR2) clustering. Employing a modified subhalo abundance matching (SHAM) technique, we populate Uchuu halos and subhalos with QSOs between 0.8 < z < 2.1. For ELGs, we modify this method to select satellite galaxies with low velocities relative to their associated central halos, and populate a separate set of Uchuu halos and subhalos with ELGs between 0.8 < z < 1.6. In this paper, we reproduce the redshift evolution of number density and clustering statistics across the fitted range of scales. We also measure the large-scale clustering bias of both the data and mock samples. These results improve simulated lightcone construction from cosmological models and enhance our understanding of the galaxy-halo connection.

Vaisakh, R. [Southern Methodist U.] (ORCID:0009000↗

DOE OSTI2026

Clustering of high-redshift quasars with DESI DR2

We present clustering measurements for high-redshift quasars using data from the Dark Energy Spectroscopic Instrument Data Release 2. Our sample consists of quasars with $2.0 < z < 3.5$ in the luminosity range $M_{1450} \leq -19.94$ mag. We measure the mean quasar bias $b_Q(\bar{z} = 2.48) = 3.61 \pm 0.01$ for the full sample of $\sim 715,000$ quasars and quantify the redshift evolution of quasar bias by dividing the sample into four equal redshift bins. There is strong evolution of the quasar bias with redshift that is well fit by the function $b_Q(z) = a [(1 + z)^2 - 6.565] + b$ with $a=0.230 \pm 0.007$ and $b=2.394 \pm 0.035$, and this fit is also a good match to lower redshift measurements in the literature. This bias evolution is consistent with a characteristic halo mass of $\bar{M}_{\mathrm{h}} \sim 10^{12}\,\mathrm{M_\odot}$ that does not vary significantly with redshift. The inferred duty cycles for quasars in our sample are $f_{\mathrm{duty}} \sim 10^{-2}$, staying mostly constant over redshifts. We investigate the luminosity dependence of quasar clustering by dividing each of our four redshift bins into three luminosity bins. The size of our quasar sample permits the first statistically significant measurement of the luminosity dependence of quasar bias at these redshifts. We measure weak dependence of quasar bias on luminosity at fixed redshift, inconsistent with no dependence, but weaker than predicted by a model in which quasar luminosity is tightly correlated with halo mass. These clustering measurements provide a stringent test for models of active black hole light curves and the black hole-halo connection at high redshift.

Charles, M. [Ohio State U.; Chicago U., KICP] (ORC↗

DOE OSTI2026

Measuring local primordial non-Gaussianity from the clustering of DESI DR1 LRGs and QSOs

We report the first measurement of primordial non-Gaussianity (PNG), parameterized by $f_{\mathrm{NL}}$, in the configuration space two-point correlation function (2pcf). We employ simulation based modeling and a novel approach for the mitigation of imaging systematics. We apply this method to samples of luminous red galaxies (LRG) and quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) during the first year of its observations (DR1). The observed 68% CL interval on $f_{\mathrm{NL}}$ is $-3^{+22}_{-21}$ using LRGs, and $ 0^{+17}_{-16}$ using QSOs. The joint measurement yields $f_{\mathrm{NL}} = -3^{+12}_{-12}$ at $\ [68\%]$ CL. Our pipeline imposes a Gaussian prior on the value of $p$ (which defines the PNG bias via the Universality relation), with $p_{\rm LRG} = 1.0 \pm 0.1$ and $p_{\rm QSO} = 1.6\pm 0.1$. The observed constraining power of DESI tracers significantly exceeds that of previous large-scale structure (LSS) surveys, and encouragingly, approaches the sensitivity of CMB probes of PNG.

Brown, Z. [Kansas State U.]↗

DOE OSTI2026

The 3D clustering of Lyman Alpha Emitters measured with DESI

We present a clustering analysis of Lyman-$α$ emitters (LAEs) using spectroscopic observations from the Dark Energy Spectroscopic Instrument (DESI) of candidates selected from the Blanco/DECam Intermediate-Band Imaging Survey (IBIS). We measure the two-point correlation function and the power spectrum, including cross-correlations with DESI quasars. Using both analytical and halo occupation distribution (HOD) simulation-based modeling, we find a linear bias of $b \sim 2.31$--$2.62$ for LAEs over the redshift range $2.26 < z < 3.41$. The analytical modeling also provides constraints on the strength of radiative transfer effects, while the HOD analysis characterizes the LAE-halo connection across multiple models. Finally, we quantify the magnitude of non-perturbative clustering effects such as Fingers of God in the LAE population, providing essential input for the accurate modeling of LAE-based cosmological analyses in forthcoming high-redshift surveys such as DESI-II.

Ebina, H. [UC, Berkeley; LBL, Berkeley] (ORCID:000↗

DOE OSTI2026

DESI Data Release 2 ELGs: Property-dependent subsamples, imaging systematics, and clustering

Using emission-line galaxies (ELGs) from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we evaluate a property-dependent correction to imaging systematics. We derive systematic weights following the same linear regression method used for other DESI tracers, but do so separately on ELG subsamples to provide a physically-informed alternative to the fiducial, neural-network-based approach. In doing so, we show that the deeper imaging in the Dark Energy Survey (DES) footprint leads to a higher overall number density but a lack of targets with extreme $g-r$ and $r-z$ colors. ELGs in the DES region also show a distinct redshift distribution when subsampled by position in the $g-r$ vs. $r-z$ plane. To address these effects, we implement a separate treatment of the DES footprint within the DESI catalog production pipeline, which is generally well-motivated and, in some cases, imperative for accurate clustering measurements. With DES treated separately, we find that property-dependent systematic weights further mitigate spurious clustering signal in $\sim$10% of subsamples, while the fiducial scheme remains optimal for the full sample.

Hagen, T. [Utah U.]↗

DOE OSTI2026

Tomography of the gamma-ray sky from cross-correlation with DESI DR2 and unWISE galaxies

We study the origin of extragalactic gamma-ray emission observed by Fermi-LAT, using the cross-correlation of the gamma-ray sky with maps of large-scale structure provided by the DESI and unWISE surveys. Tomographic cross-correlation reveals the bias-weighted redshift distributions of gamma-ray sources. We first illustrate this method by cross-correlating detected gamma-ray point sources with large-scale structure. We find a significant cross-correlation and infer a point source redshift distribution broadly consistent with the distribution of identified optical counterparts previously reported in the literature, as well as a similar linear bias ($b \approx 2$) to massive galaxies that host bright active galactic nuclei. We then study the clustering of the Fermi unresolved gamma-ray background (UGRB), both in auto-correlation and in cross-correlation with large-scale structure. We detect the cross-correlation of the UGRB and LSS at $\sim 10σ$ in total, with highly significant detections from both DESI and unWISE. Our measurements suggest that the redshift distribution of the UGRB is broadly consistent with the redshift distribution of detected point sources. Additionally, we find a relatively weak amplitude for the cross-correlation with large-scale structure at z < 2, suggesting a significant fraction of the UGRB does not come from z < 2 large-scale structure. A natural candidate is contamination of from residual Galactic emission, and our best estimate of the contamination level derived from the UGRB auto-spectrum suggests that the mean bias of UGRB sources is indeed quite similar to the bias of detected Fermi point sources. However, we cannot exclude additional emission from gamma-ray sources at high redshift, z > 2, and we suggest that cross-correlation with tracers at z > 2, including CMB lensing, would be the ideal way to determine the fraction of z > 2 emission.

Krolewski, Alex [Waterloo U., Math. Dept.; Waterlo↗

DOE OSTI2026

Joint Modeling of GD-1 and C-19 as Old Streams

DESI observational data for the GD-1 and C-19 streams are compared to stream simulations in a common evolving multi-halo potential of a Milky Way-like galaxy based on a cosmological simulation. The goal is to find the best match of the stream velocity spread and the density power spectrum stream density to simulations having either CDM or WDM subhalos. The cocoon velocity width integrated over the length of the stream is independent of orbital blurring along the stream and the power spectrum integrates over the width of the stream, sidestepping the geometric details of the streams. Streams develop from star clusters inserted at $\simeq$1 Gyr after the Big Bang and evolved for 13 Gyr to their current orbital positions. Streams in a CDM subhalo population provide the best match to the velocity width, with streams younger than 10 Gyr ruled out as insufficiently hot. The progenitor star cluster masses, which determine the fraction of stars released at late times which comprise the stream core, are found to be $\simeq 8\times 10^4 M_\odot$ for GD-1 and $\simeq 4\times 10^4 M_\odot$ for C-19, although the mass depends on the star cluster half mass radius. Stream heating leads to stream lumpiness which is measurable for the relatively large and clean GD-1 dataset. The stream density power spectrum measured along the length of the DESI GD-1 sample is in good agreement with CDM simulations, with 1.7 to 1.9 times more power than WDM 7 keV and 5.5 keV simulations.

Carlberg, Raymond G. [Toronto U.] (ORCID:000000027↗

DOE OSTI2026

DESI DR2 Baryon Acoustic Oscillations from the Lyman Alpha Forest Multipoles

We present an alternative measurement of the Baryon Acoustic Oscillation (BAO) using the Legendre multipole representation of the Ly$α$ forest correlation functions from the second data release (DR2) of the Dark Energy Spectroscopic Instrument survey. Compressing the auto- and cross-correlation functions into Legendre multipoles yields a positive-definite covariance matrix without any smoothing -- unlike the baseline DR2 analysis -- thanks to a significantly reduced data vector size. We introduce the statistical corrections required to debias the finite-sample covariance matrix estimate and demonstrate that monopole and quadrupole terms for both auto- and cross-correlations can be used even when the correlation functions are distorted by continuum errors and contaminated by metals. This formalism has slightly diminished the constraining power of the BAO scale, while considerably weakening constraints on nuisance parameters. We measure the isotropic BAO scale with $0.93\%$ precision at $z_\mathrm{eff}=2.35$, the Hubble parameter $H(z_\mathrm{eff})=(239.5\pm3.4)~(147.09~\mathrm{Mpc}/r_d) ~\mathrm{km~s}^{-1}~\text{Mpc}^{-1}$, and the transverse comoving distance $D_M(z_\mathrm{eff})=(5.80 \pm 0.10)~(r_d/147.09~\mathrm{Mpc})$~Gpc for a given value of the sound horizon ($r_d$). Our BAO results are entirely consistent with the baseline DR2 analysis.