Harvesting the Ly\alpha\ forest with convolutional neural networks

We develop a machine learning based algorithm using a convolutional neural network (CNN) to identify low HI column density Ly$\alpha$ absorption systems ($\log{N_{\mathrm{HI}}}/{\rm cm}^{-2}<17$) in the Ly$\alpha$ forest, and predict their physical properties, such as their HI column density ($\log{N}_{\mathrm{HI}}/{\rm cm}^{-2}$), redshift ($z_{\mathrm{HI}}$), and Doppler width ($b_{\mathrm{HI}}$). Our CNN models are trained using simulated spectra (S/N $\simeq10$), and we test their performance on high quality spectra of quasars at redshift $z\sim2.5-2.9$ observed with the High Resolution Echelle Spectrometer on the Keck I telescope. We find that $\sim78\%$ of the systems identified by our algorithm are listed in the manual Voigt profile fitting catalogue. We demonstrate that the performance of our CNN is stable and consistent for all simulated and observed spectra with S/N $\gtrsim10$. Our model can therefore be consistently used to analyse the enormous number of both low and high S/N data available with current and future facilities. Our CNN provides state-of-the-art predictions within the range $12.5\leq\log{N_{\mathrm{HI}}}/\mathrm{cm^{-2}}<15.5$ with a mean absolute error of $\Delta(\log{N}_{\mathrm{HI}}/{\rm cm}^{-2})=0.13$, $\Delta(z_{\mathrm{HI}})=2.7\times{10}^{-5}$, and $\Delta(b_{\mathrm{HI}})=4.1\ \mathrm{km\ s^{-1}}$. The CNN prediction costs $<3$ minutes per model per spectrum with a size of 120\,000 pixels using a laptop computer. We demonstrate that CNNs can significantly increase the efficiency of analysing Ly$\alpha$ forest spectra, and thereby greatly increase the statistics of Ly$\alpha$ absorbers.

Paper

References (12)

02GPyOpt: A Bayesian Optimization framework2020 · arXiv e-prints,
03Identifying Lyα systems with a CNN2019
04L85 APPENDIX A: THE IMPACT ON PREDICTING SPECTRA WITH DIFFERENT INITIAL SETUPS We test our pre-trained CNN model on 25 newly generated spectra with: (1) quasars at different redshifts2019 · MNRAS,
05MNRAS, 472, 1850 Garzilli A., Boyarsky A.2017 · Physics Letters B
06on Learning Representations,2015 · 3rd International Conference
07MNRAS, 395, 736 Bottrell C2009 · MNRAS
08VPFIT: Voigt profile fitting2009 · MNRAS,
09Overall, different FWHMs in this range show no impact to both the detection and physical property estimates. This is because the widths of the Lyα forest absorption lines (FWHM2000
10Kluwer Academic Publishers, USA1998
12Society of1994 · Photo-Optical

Similar papers

© 2026 NYSGPT2525 LLC