Unified reconstruction of the Lyman-alpha power spectrum with Hamiltonian Monte Carlo
Abstract
The complex geometry of the Lyα forest data has motivated the use of various two-point statistics as alternatives to the three-dimensional power spectrum (P3D), which carries cosmological information in Fourier space. On large scales, the three-dimensional correlation function (3D) has provided robust measurements of the baryon acoustic oscillation (BAO) scale at 150~Mpc. On smaller scales, the one-dimensional power spectrum, P1D(k\|), has been the primary tool for extracting information. At the same time, the cross-spectrum, P×(θ, k\|), has been introduced to incorporate angular information without the complications caused by survey window functions. We propose an analytical forward-modeling framework to reconstruct P3D from all these observables, based on the mathematical relation between them and P3D. We demonstrate the performance of our method using a hypothetical mock data vector representative of future Dark Energy Spectroscopic Instrument (DESI) measurements. We show that the monopole of P3D can be reconstructed in 25 k bins between 0.07~Mpc-1 and 1.8~Mpc-1, achieving an average precision of σP/P=13\% across the bins. Our method can serve as an intermediary for consistency checks, though it is not intended to replace direct P3D estimation.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.