An Optimized Ly$\alpha$ Forest Inversion Tool Based on a Quantitative Comparison of Existing Reconstruction Methods

Abstract

We present a same-level comparison of the most prominent inversion methods for the reconstruction of the matter density field in the quasi-linear regime from the Ly$\alpha$ forest flux. Moreover, we present a pathway for refining the reconstruction in the framework of numerical optimization. We apply this approach to construct a novel hybrid method. The methods which are used so far for matter reconstructions are the Richardson-Lucy algorithm, an iterative Gauss-Newton method and a statistical approach assuming a one-to-one correspondence between matter and flux. We study these methods for high spectral resolutions such that thermal broadening becomes relevant. The inversion methods are compared on synthetic data (generated with the lognormal approach) with respect to their performance, accuracy, their stability against noise, and their robustness against systematic uncertainties. We conclude that the iterative Gauss-Newton method offers the most accurate reconstruction, in particular at small S/N, but has also the largest numerical complexity and requires the strongest assumptions. The other two algorithms are faster, comparably precise at small noise-levels, and, in the case of the statistical approach, more robust against inaccurate assumptions on the thermal history of the intergalactic medium (IGM). We use these results to refine the statistical approach using regularization. Our new approach has low numerical complexity and makes few assumptions about the history of the IGM, and is shown to be the most accurate reconstruction at small S/N, even if the thermal history of the IGM is not known. Our code will be made publicly available under https://github.com/hmuellergoe/reglyman.

arXiv comment21 pages, 13 figures, 3 tables, submitted to MNRAS, a python toolkit containing our implementations will be made publicly available under https://github.com/hmuellergoe/reglyman

References in corpus (71)

• Planck 2015 results. XIII. Cosmological parameters
• The Sloan Digital Sky Survey: Technical Summary
• The 2dF Galaxy Redshift Survey: Spectra and redshifts
• SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems
• The Cosmic Evolution Survey (COSMOS) -- Overview
• Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe
• How Filaments are Woven into the Cosmic Web
• Galaxy and Mass Assembly (GAMA): survey diagnostics and core data release
• Constraining Warm Dark Matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest
• 21cmFAST: A Fast, Semi-Numerical Simulation of the High-Redshift 21-cm Signal
• Equation of State of the Photoionized Intergalactic Medium
• A Survey of Star-Forming Galaxies in the z=1.4-2.5 `Redshift Desert': Overview
• Reionization of the Inhomogeneous Universe
• The Dark Energy Survey Data Release 1
• The Lyman-alpha Forest Power Spectrum from the Sloan Digital Sky Survey
• The VIMOS VLT Deep Survey final data release: a spectroscopic sample of 35016 galaxies and AGN out to z~6.7 selected with 17.5<=i_{AB}<=24.7
• P3M-SPH simulations of the Lyman-alpha Forest
• The Lyman Alpha Forest from Gravitational Collapse in the CDM+Lambda Model
• The Evolution of the Intergalactic Medium
• Evolution of Structure in the Intergalactic Medium and the Nature of the Ly-alpha Forest
• Neutrino masses and cosmology with Lyman-alpha forest power spectrum
• Inferring the dark matter power spectrum from the Lyman-alpha forest in high-resolution QSO absorption spectra
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• HeII Reionization and its Effect on the IGM
• A Direct Precision Measurement of the Intergalactic Lyman-alpha Opacity at 2<z<4.2
• The one-dimensional Ly-alpha forest power spectrum from BOSS
• Galaxy evolution in the metric of the Cosmic Web
• The Lyman-alpha forest opacity and the metagalactic hydrogen ionization rate at z~2-4
• A Refined Measurement of the Mean Transmitted Flux in the Ly-alpha Forest over 2 < z < 5 Using Composite Quasar Spectra
• The darkness that shaped the void: dark energy and cosmic voids
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• The Zeldovich approximation
• Experimental Constraints on Self-consistent Reionization Models
• First Data Release of the COSMOS Lyman-Alpha Mapping And Tomography Observations: 3D Lyman-$\alpha$ Forest Tomography at 2.05 < z < 2.55
• The UVES Spectral Quasar Absorption Database (SQUAD) Data Release 1: The first 10 million seconds
• The Temperature-Density Relation in the Intergalactic Medium at Redshift <z>=2.4
• Lyman-alpha Forest Tomography from Background Galaxies: The First Megaparsec-Resolution Large-Scale Structure Map at z>2
• Cosmological parameters sigma_8, the baryon density, and the UV background intensity from a calibrated measurement of H I Lyman-alpha absorption at z = 1.9
• The H I opacity of the intergalactic medium at redshifts 1.6 < z < 3.2
• Inversion of the Lyman-$\alpha$ forest: 3D investigation of the intergalactic medium
• Constraining the reionization history with QSO absorption spectra
• Shadow of a Colossus: A z=2.45 Galaxy Protocluster Detected in 3D Ly-a Forest Tomographic Mapping of the COSMOS Field
• A New Measurement of the Temperature Density Relation of the IGM From Voigt Profile Fitting
• A ``Skewed'' lognormal approximation to the probability distribution function of the large-scale density field
• Observational Requirements for Lyman-alpha Forest Tomographic Mapping of Large-Scale Structure at z ~ 2
• The Lyman-$\alpha$ forest as a diagnostic of the nature of the dark matter
• A first step towards a direct inversion of the Lyman forest in QSO spectra
• Proto-groups at 1.8<z<3 in the zCOSMOS-deep sample
• Semi analytic approach to understanding the distribution of neutral hydrogen in the universe: Comparison of simulations with observations
• An Emulator for the Lyman-alpha Forest
• Cosmic Web Reconstruction through Density Ridges: Catalogue
• Computational helioseismology in the frequency domain: acoustic waves in axisymmetric solar models with flows
• Protocluster Discovery in Tomographic Ly$\alpha$ Forest Flux Maps
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• Generating mock data sets for large-scale Lyman-α forest correlation measurements
• Sample variance and Lyman-alpha forest transmission statistics
• Modelling the IGM and the Lyalpha forest at high redshift from the dark matter distribution
• Nonparametric 3D map of the IGM using the Lyman-alpha forest
• A Detection of $z$~2.3 Cosmic Voids from 3D Lyman-$\alpha$ Forest Tomography in the COSMOS Field
• Inferring high redshift large-scale structure dynamics from the Lyman-alpha forest
• Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach
• On the observability of coupled dark energy with cosmic voids
• Acceleration techniques for regularized Newton methods applied to electromagnetic inverse medium scattering problems
• Large-scale 3D mapping of the intergalactic medium using the Lyman Alpha Forest
• Reconstructing the baryon acoustic oscillations using biased tracers
• The Effect of Large-Scale Power on Simulated Spectra of the Lya forest
• Simulating MOS science on the ELT: Ly$\alpha$ forest tomography
• Physical properties and small-scale structure of the Lyman-alpha forest: Inversion of the HE 1122-1628 UVES spectrum
• Cosmic density field reconstruction from Ly-alpha forest data
• Measuring the temperature and profiles of Lyman-$\alpha$ absorbers
• Trend Filtering -- II. Denoising Astronomical Signals with Varying Degrees of Smoothness

Cited by in corpus (3)

• Deep Forest: Neural Network reconstruction of the Lyman-alpha forest
• Searching for dilaton fields in the Ly$\alpha$ forest
• A Novel Estimator for the Equation of State of the IGM by Ly$\alpha$ Forest Tomography