Pith. Sign in

REVIEW

Improving cosmological distance measurements at high redshift using Lyman-Alpha forest BAO reconstruction

Cosmology and Nongalactic Astrophysics · Submitted 27 Oct 2026 · 123 references

Pith's verdict: A genuine, cleanly executed extension of Zel'dovich BAO reconstruction to the Lyman-alpha forest, with honest upper-limit framing on the real-data gains; the open question is how much continuum fitting erodes the radial improvement.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 1 major objection · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.37664 (version 1)
Pith IDRVNLUZYO, Pith's permanent record of this paper and its review
CitationBoryana Hadzhiyska, Roger de Belsunce. "Improving cosmological distance measurements at high redshift using Lyman-Alpha forest BAO reconstruction." arXiv:2610.37664 (2026). https://pith.science/paper/RVNLUZYO
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Applying Zel'dovich reconstruction, the standard tool of galaxy BAO analyses, to the Lyman-alpha forest halves the radial damping of the BAO peak at z≈2.5 and tightens the radial and transverse distance constraints by roughly 25% and…

The reading

DESI's highest-redshift BAO measurement comes from the Lyman-alpha forest, absorption features in quasar spectra that trace matter at z≈2.3–4. Non-linear structure formation damps and shifts the acoustic peak, and current analyses capture the damping only through a Gaussian broadening in the template. This paper argues that the Zel'dovich reconstruction routinely applied to galaxy surveys can be extended to the forest, and that doing so halves the radial BAO damping scale and erases most of the non-linear shift that pipelines currently must model away with nuisance parameters. The demonstration runs on idealized mocks whose continuum is perfectly known and whose line-of-sight modes are intact, and the authors state plainly that their gains are upper limits, with continuum fitting flagged as the likeliest effect to erode them. If the gains survive realistic data treatment, DESI and DESI-II gain a materially sharper standard ruler in the matter-dominated era, which is where distance measurements best discriminate between competing explanations of the reported dark-energy tension.

What carries the argument

The central object is the reconstructed displacement field. The paper solves the linear redshift-space equation $\nabla\cdot\Psi + \beta_F\, \nabla\cdot[(\Psi\cdot\hat z)\hat z] = -S\star\delta_F/b_F$, where $b_F<0$ is the forest's linear density bias, $\beta_F = f b_\eta/b_F$ its redshift-space distortion parameter set by the velocity-gradient bias and growth rate $f$, and $S$ a Gaussian smoothing with $R=10\,h^{-1}\,$Mpc. Because $\delta_F/b_F$ carries the Kaiser factor $(1+\beta_F\mu^2)$, that factor cancels and the solution reduces exactly to the real-space Zel'dovich displacement of the smoothed linear density. The flux pixels, together with an unmoved reference copy that encodes the window function, are shifted by the redshift-space displacement $\Psi_s = \Psi + f(\Psi\cdot\hat z)\hat z$ following the RecSym convention, and the field is repainted. This addresses the two obstacles unique to the forest: a negative density response with an independent velocity-gradient response, and a near-continuous field sampled densely along but sparsely across the line of sight. A pair-count spectral estimator adapted to the forest measures the power spectrum on masked sightlines without grid-based aliasing, and linear control variates subtract the part of each measured spectrum predictable from the known initial conditions, exposing scale shifts at the $\pm0.05\%$ level.

Extended reading notes

Core claim

The paper's central claim is that the flux decrement of the Lyman-$\alpha$ forest can be fed through the Zel'dovich reconstruction equation, with the negative density bias $b_F$ and the independent velocity-gradient bias $b_\eta$ handled explicitly, to yield a displacement field that matches the linear one to sub-percent accuracy (cross-transfer function unity within 0.1% at all models) and that moving the flux pixels by this displacement halves the fitted radial damping scale, from 4.7 to 2.2 $h^{-1}\,$Mpc in the fiducial model. Using linear control variates to suppress sample variance by up to a factor of $\sqrt{50}$, the uncertainties on the radial and transverse scaling parameters $\alpha_\parallel$ and $\alpha_\perp$ shrink by approximately 25% and 15% in the fully sampled regime, the radial non-linear shift falls from $+0.3\%$ to $+0.10\pm0.07\%$, and $\alpha_\perp$ stays unbiased. Under a DESI-like sightline mask of 50 deg$^{-2}$, the gains persist at roughly 27% and 10%, whether the displacement is estimated from the forest itself or from an external halo tracer at the number density of DESI-II's Lyman-break galaxy and Lyman-$\alpha$ emitter targets, $\bar n \simeq 10^{-3}\,h^3\,$Mpc$^{-3}$. The figure of merit for the pair $(\alpha_\parallel,\alpha_\perp)$ improves by factors of 1.66, 1.65, and 1.76 for the fully sampled, DESI-like, and high-density configurations.

Load-bearing premise

The load-bearing premise is that mocks with perfectly known continuum, intact line-of-sight modes, no pixel noise, and no radial mask capture the information loss that matters for reconstruction, and the authors themselves state that continuum fitting — which projects out the long-wavelength modes along each sightline — will likely have the strongest impact and that their gains are therefore upper limits.

Editorial extensions

If this is right

  • The non-linear radial shift drops from +0.3% to +0.10±0.07%, meaning less of the BAO position must be absorbed by a large set of EFT nuisance parameters whose priors can dilute the constraint.
  • Because an external tracer at $\bar n \gtrsim 10^{-3}\,h^3\,$Mpc$^{-3}$ reproduces almost all of the self-reconstruction gain, the method does not need to depend on the observed flux and is partly shielded from continuum-fitting systematics.
  • Figure-of-merit improvements of 1.65–1.76 in the BAO scaling parameters are equivalent to enlarging the effective survey volume by roughly 65–76% for the distance measurement.
  • The same reconstruction also sharpens the Lyman-alpha forest–quasar cross-correlation, which carries much of the constraining power of current forest BAO analyses, with a radial gain near 1.37.
  • Gains hold across all four forest models spanning $\beta_F \approx 1.0$–2.4, so the result does not hinge on a fine-tuned redshift-space distortion parameter.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If a realistic continuum-fitting step projects out exactly the long-wavelength modes the reconstruction uses, the paper's Regime II result suggests the fallback of computing displacements from external tracers may rescue most of the gain in real DESI data, since the tracer displacement never touches the observed flux.
  • The paper's six boxes share a single cosmology, so the dependence of the optimal smoothing radius and bias parameters on cosmology is untested; a natural extension is to repeat the shifts with the control-variate procedure across several cosmologies to check that reconstruction does not imprint a cosmology-dependent scale.
  • The clean separation of a sharpened peak from a broadened one could let future forest BAO analyses drop their simulation-based bias priors entirely, testing the EFT nuisance-parameter bookkeeping against a fully independent route to the same acoustic scale.
Original abstract

The Dark Energy Spectroscopic Instrument (DESI) measures baryon acoustic oscillations (BAO) in the Lyman-$\alpha$ (Lyman-alpha) forest and provides a robust probe of the expansion history in the matter-dominated era. Yet, non-linear structure formation broadens and shifts the BAO peak. To undo these effects, we extend standard Zel'dovich reconstruction, routinely applied to galaxy surveys, to the Lyman-alpha forest. We validate our methodology on large-volume idealized Lyman-alpha forest mocks painted on top of the \textsc{AbacusSummit} $N$-body simulations, using linear control variates to suppress sample variance and thereby reduce the error bars by a factor of up to $\approx \sqrt{50}$. In fully sampled periodic boxes, reconstruction halves the radial BAO damping scale and reduces the uncertainties on the radial ($\apar$) and transverse ($\aperp$) scaling parameters by $\approx25\%$ and $\approx15\%$, respectively. It also reduces the non-linear shift from $\Delta\apar=+0.3\%$ to $+0.10\pm0.07\%$, while $\aperp$ remains unbiased. To make this approach readily applicable to observational data, we demonstrate the impact of a survey mask with a DESI-like density of $50\,\mathrm{deg}^{-2}$. For this mask, we compare two approaches to estimate the displacement field: either directly from the Lyman-alpha forest or using an external tracer with a number density of $\bar n \simeq 10^{-3}\,h^3\mathrm{Mpc}^{-3}$, such as Lyman-break galaxies and Lyman-$\alpha$ emitters---both key targets for DESI-II. Both approaches give comparable gains and reduce the uncertainties on $\apar$ and $\aperp$ by $\approx 27\%$ and $\approx 10\%$, respectively, with mild improvements for futuristic high-density surveys with $160\,\mathrm{deg}^{-2}$. Correspondingly, the figure of merit improves by factors of $1.66$ for fully sampled boxes, $1.65$ for a DESI-like, and $1.76$ for a high-density survey.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-28: 1 request, 8,166 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.50 · human participants 0.01 · controlled comparison 0.78 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.86 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The Dark Energy Spectroscopic Instrument (DESI) measures baryon acoustic oscillations (BAO) in the Lyman-$\alpha$ (Lyman-alpha) forest and provides a robust probe of the expansion history in the matter-dominated era. Yet, non-linear structure formation broadens and shifts the BAO peak. To undo these effects, we extend standard Zel'dovich reconstruction, routinely applied to galaxy surveys, to the Lyman-alpha forest. We validate our methodology on large-volume idealized Lyman-alpha forest mocks painted on top of the \textsc{AbacusSummit} $N$-body simulations, using linear control variates to suppress sample variance and thereby reduce the error bars by a factor of up to $\approx \sqrt{50}$. In fully sampled periodic boxes, reconstruction halves the radial BAO damping scale and reduces the uncertainties on the radial ($\apar$) and transverse ($\aperp$) scaling parameters by $\approx25\%$ and $\approx15\%$, respectively. It also reduces the non-linear shift from $\Delta\apar=+0.3\%$ to $+0.10\pm0.07\%$, while $\aperp$ remains unbiased. To make this approach readily applicable to observational data, we demonstrate the impact of a survey mask with a DESI-like density of $50\,\mathrm{deg}^{-2}$. For this mask, we compare two approaches to estimate the displacement field: either directly from the Lyman-alpha forest or using an external tracer with a number density of $\bar n \simeq 10^{-3}\,h^3\mathrm{Mpc}^{-3}$, such as Lyman-break galaxies and Lyman-$\alpha$ emitters---both key targets for DESI-II. Both approaches give comparable gains and reduce the uncertainties on $\apar$ and $\aperp$ by $\approx 27\%$ and $\approx 10\%$, respectively, with mild improvements for futuristic high-density surveys with $160\,\mathrm{deg}^{-2}$. Correspondingly, the figure of merit improves by factors of $1.66$ for fully sampled boxes, $1.65$ for a DESI-like, and $1.76$ for a high-density survey.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The Dark Energy Spectroscopic Instrument (DESI) measures baryon acoustic oscillations (BAO) in the Lyman-$\alpha$ (Lyman-alpha) forest and provides a robust probe of the expansion history in the matter-dominated era. background 0.99 0.21 0.05 asserted 1.00
s1 Yet, non-linear structure formation broadens and shifts the BAO peak. background 0.93 0.09 0.02 asserted 1.00
s2 To undo these effects, we extend standard Zel'dovich reconstruction, routinely applied to galaxy surveys, to the Lyman-alpha forest. contribution 0.67 0.87 0.06 asserted 0.98
s3 We validate our methodology on large-volume idealized Lyman-alpha forest mocks painted on top of the \textsc{AbacusSummit} $N$-body simulations, using linear control variates to suppress sample variance and thereby reduce the error bars by a factor of up to $\approx \sqrt{50}$. method 0.90 0.96 0.08 asserted 1.00
s4 In fully sampled periodic boxes, reconstruction halves the radial BAO damping scale and reduces the uncertainties on the radial ($\apar$) and transverse ($\aperp$) scaling parameters by $\approx25\%$ and $\approx15\%$, respectively. result 1.00 0.95 0.31 asserted 1.00
s5 It also reduces the non-linear shift from $\Delta\apar=+0.3\%$ to $+0.10\pm0.07\%$, while $\aperp$ remains unbiased. result 1.00 0.93 0.18 asserted 1.00
s6 To make this approach readily applicable to observational data, we demonstrate the impact of a survey mask with a DESI-like density of $50\,\mathrm{deg}^{-2}$. contribution 0.51 0.89 0.03 asserted 0.99
s7 For this mask, we compare two approaches to estimate the displacement field: either directly from the Lyman-alpha forest or using an external tracer with a number density of $\bar n \simeq 10^{-3}\,h^3\mathrm{Mpc}^{-3}$, such as Lyman-break galaxies and Lyman-$\alpha$ emitters---both key targets for DESI-II. method 0.91 0.80 0.04 asserted 0.87
s8 Both approaches give comparable gains and reduce the uncertainties on $\apar$ and $\aperp$ by $\approx 27\%$ and $\approx 10\%$, respectively, with mild improvements for futuristic high-density surveys with $160\,\mathrm{deg}^{-2}$. result 1.00 0.94 0.22 asserted 1.00
s9 Correspondingly, the figure of merit improves by factors of $1.66$ for fully sampled boxes, $1.65$ for a DESI-like, and $1.76$ for a high-density survey. result 0.98 0.93 0.33 asserted 1.00

What kind of paper the classifier reports this is: method or tool 0.93 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.