Improving IGM temperature constraints using wavelet analysis on high-redshift quasars

The thermal state of the intergalactic medium contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young Universe. This thermal state is encoded in the small-scale structure of Lyman-α (Ly α) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines. At high redshifts, where the opacity is large, averaging mixes high signal-to-noise ratio transmission spikes with noisy absorption troughs. Wavelet amplitudes are an alternate statistic that maintains spatial information while quantifying fluctuations at the same spatial frequencies as the power spectrum, giving them the potential to more sensitively measure the small-scale structure. Previous Ly α forest studies using wavelet amplitude probability density functions (PDFs) used limited spatial frequencies and neglected strong correlations between PDF bins and across wavelets scales, resulting in suboptimal and unreliable parameter inference. Here we present a novel method for performing statistical inference using wavelet amplitude PDFs that spans the full range of spatial frequencies probed by the power spectrum and that fully accounts for these correlations. We applied this procedure to realistic mock data drawn from a simple thermal model parametrized by the temperature at mean density, T₀, and find that wavelets deliver 1σ constraints on T₀ that are on average 7 per cent more sensitive at z = 5 (12 per cent at z = 6) than those from the power spectrum. We consider the possibility of combing wavelet PDFs with the power, but find that this does not lead to improved sensitivity.