Unified reconstruction of the Lyman-alpha power spectrum with Hamiltonian Monte Carlo
Phys. Rev. D 113, 043528 – Published 20 February, 2026
DOI: https://doi.org/10.1103/b391-g2rm
Abstract
The complex geometry of the Lyman-alpha () forest data has motivated the use of various two-point statistics as alternatives to the three-dimensional power spectrum (), which carries cosmological information in Fourier space. On large scales, the three-dimensional correlation function () has provided robust measurements of the baryon acoustic oscillation scale at 150 Mpc. On smaller scales, the one-dimensional power spectrum, , has been the primary tool for extracting information. At the same time, the cross spectrum, , has been introduced to incorporate angular information without the complications caused by survey window functions. We propose an analytical forward-modeling framework to reconstruct from all these observables, based on the mathematical relation between them and . We demonstrate the performance of our method using a hypothetical mock data vector representative of future Dark Energy Spectroscopic Instrument measurements. We show that the monopole of can be reconstructed in bins between and , achieving an average precision of across the bins. Our method can serve as an intermediary for consistency checks, though it is not intended to replace direct estimation.