Analytic compression of the effective field theory of the Lyman-alpha forest

The 1D flux power spectrum ($P_{\mathrm{1D}}$) of the Ly$α$ forest provides an exceptionally high-resolution probe of structure formation down to small scales ($k\approx1-10~\text{$h~$Mpc$^{-1}$}$). These scales carry the imprints of massive neutrinos, warm dark matter, and the running of the primordial power spectrum spectral index. The effective field theory (EFT) is a promising perturbative approach to systematically and efficiently describe the Ly$α$ forest, but it faces challenges in its application to $P_{\mathrm{1D}}$, as many EFT parameters become degenerate when projected along the line of sight. In addition, this projection generates new stochastic terms from the integration over small-scale modes. In this work, we address these issues by compressing the EFT model space using the Fisher matrix formalism and linearizing the resulting compression directions, enabling analytic template marginalization and significantly reducing the computational cost of likelihood evaluation. We use hydrodynamical simulations to obtain a baseline estimate of EFT parameters, and use the DESI DR1 $P_{\mathrm{1D}}$ measurements to derive compression directions. We then marginalize over deviations from the baseline using these compression directions and forecast the constraining power of our formalism. We find that even in conservative scenarios where each data redshift bin requires its own set of EFT parameters, the cosmological constraints saturate with the linear bias, two leading-order 1D stochastic terms, and three principal combinations of the remaining EFT templates. In this case, our forecasted precision of the amplitude ($Δ^2_p$) and the logarithmic slope ($n_p$) of the linear matter power spectrum at the pivot scale ($k_p=0.7~\text{Mpc}^{-1}$) is $10\%$ and $2.0\%$, respectively, which is similar to emulator-based analyses that include observational data systematics.