Correlations in the three-dimensional Lyman-alpha forest contaminated by high column density absorbers

TL;DR

This paper investigates how high column density hydrogen absorbers contaminate 3D Lyman-alpha forest correlations, proposing a simple physical model to correct for residual effects in cosmological analyses.

Contribution

It introduces a linear Voigt profile model for residual contamination from high column density absorbers, improving upon existing approximations for analyzing the 3D Lyman-alpha flux power spectrum.

Findings

Residual contamination is scale-dependent and significant at small scales.

The linear Voigt model accurately fits simulation data for residuals up to |k|<1 h/Mpc.

High column density systems are mostly removed in surveys, reducing their impact on small scales.

Abstract

Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (HI; having column densities $N(\mathrm{HI}) > 1.6\times10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a function of the column density of the HCD absorbers and redshift by measuring 3D flux power spectra in cosmological hydrodynamical simulations from the Illustris project. Survey pipelines exclude regions containing the largest damping wings. We find that, even after this procedure, there is a scale-dependent correction to the 3D Lyman-alpha forest flux power spectrum from residual contamination. We model this residual using a simple physical model of the HCD absorbers as linearly biased tracers of the matter density distribution, convolved with their Voigt profiles and integrated over the column density distribution function. We recommend the use of this model over existing models used in data analysis, which approximate the damping wings as top-hats and so miss shape information in the extended wings. The simple 'linear Voigt model' is statistically consistent with our simulation results for a mock residual contamination up to small scales ($|k| < 1\,h\,\mathrm{Mpc}^{-1}$). It does not account for the effect of the highest column density absorbers on the smallest scales (e.g., $|k| > 0.4\,h\,\mathrm{Mpc}^{-1}$ for small damped Lyman-alpha absorbers; HCD absorbers with $N(\mathrm{HI}) \sim 10^{21}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$). However, these systems are in any case preferentially removed from survey data. Our model is appropriate for an accurate analysis of the baryon acoustic oscillations feature. It is additionally essential for reconstructing the full shape of the 3D flux power spectrum.