Field-level Lyman-α forest modeling in redshift space via augmented nonlocal Fluctuating Gunn-Peterson Approximation

¹ Instituto de Astrofísica de Canarias, Calle Via Láctea s/n, 38205 La Laguna, Tenerife, Spain
e-mail: fkitaura@iac.es
² Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
³ Department of Physics and Astronomy, Universitá degli Studi di Padova, Vicolo dell’Osservatorio 3, 35122 Padova, Italy
e-mail: francesco.sinigaglia@phd.unipd.it
⁴ INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy
⁵ Theoretical Astrophysics, Department of Earth and Space Science, Graduate School of Science, Osaka University, 1-1 Machikaneyama, Toyonaka, Osaka 560-0043, Japan
⁶ Kavli-IPMU (WPI), University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8583, Japan
⁷ Department of Physics & Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Pkwy, Las Vegas, NV 89154-4002, USA

Received: 17 May 2023
Accepted: 10 October 2023

Abstract

Context. Devising fast and accurate methods of predicting the Lyman-α forest at the field level, avoiding the computational burden of running large-volume cosmological hydrodynamic simulations, is of fundamental importance to quickly generate the massive set of simulations needed by the state-of-the-art galaxy and Lyα forest spectroscopic surveys.

Aims. We present an improved analytical model to predict the Lyα forest at the field level in redshift space from the dark matter field, expanding upon the widely used Fluctuating Gunn-Peterson Approximation (FGPA). Instead of assuming a unique universal relation over the whole considered cosmic volume, we introduce a dependence on the cosmic web environment (knots, filaments, sheets, and voids) in the model, thereby effectively accounting for nonlocal bias. Furthermore, we include a detailed treatment of velocity bias in the redshift space distortion modeling, allowing the velocity bias to be cosmic-web-dependent.

Methods. We first mapped the dark matter field from real to redshift space through a particle-based relation including velocity bias, depending on the cosmic web classification of the dark matter field in real space. We then formalized an appropriate functional form for our model, building upon the traditional FGPA and including a cutoff and a boosting factor mimicking a threshold and inverse-threshold bias effect, respectively, with model parameters depending on the cosmic web classification in redshift space. Eventually, we fit the coefficients of the model via an efficient Markov chain Monte Carlo scheme.

Results. We find evidence for a significant difference between the same model parameters in different environments, suggesting that for the investigated setup the simple standard FGPA is not able to adequately predict the Lyα forest in the different cosmic web regimes. We reproduce the summary statistics of the reference cosmological hydrodynamic simulation that we use for comparison, yielding an accurate mean transmitted flux, probability distribution function, 3D power spectrum, and bispectrum. In particular, we achieve maximum deviation and average deviation accuracy in the Lyα forest 3D power spectrum of ∼3% and ∼0.1% up to k ∼ 0.4 h Mpc⁻¹, and ∼5% and ∼1.8% up to k ∼ 1.4 h Mpc⁻¹.

Conclusions. Our new model outperforms previous analytical efforts to predict the Lyα forest at the field level in all the probed summary statistics, and has the potential to become instrumental in the generation of fast accurate mocks for covariance matrices estimation in the context of current and forthcoming Lyα forest surveys.