No Arabic abstract
The Lyman-$alpha$ forest is a highly non-linear field with a lot of information available in the data beyond the power spectrum. The flux probability distribution function (PDF) has been used as a successful probe of small-scale physics. In this paper we argue that measuring coefficients of the Legendre polyonomial expansion of the PDF offers several advantages over measuring the binned values as is commonly done. In particular, $n$-th coefficient can be expressed as a linear combination of the first $n$ moments, allowing these coefficients to be measured in the presence of noise and allowing a clear route for marginalisation over mean flux. Moreover, in the presence of noise, our numerical work shows that a finite number of coefficients are well measured with a very sharp transition into noise dominance. This compresses the available information into a small number of well-measured quantities. We find that the amount of recoverable information is a very non-linear function of spectral noise that strongly favors fewer quasars measured at better signal to noise.
We use the probability distribution function (PDF) of the lya forest flux at z=2-3, measured from high-resolution UVES/VLT data, and hydrodynamical simulations to obtain constraints on cosmological parameters and the thermal state of the intergalactic medium (IGM) at z 2-3. The observed flux PDF at z=3 alone results in constraints on cosmological parameters in good agreement with those obtained from the WMAP data, albeit with about a factor two larger errors. The observed flux PDF is best fit with simulations with a matter fluctuation amplitude of sigma_8=0.8-0.85 pm 0.07 and an inverted IGM temperature-density relation (gamma ~ 0.5-0.75), consistent with our previous results obtained using a simpler analysis. These results appear to be robust to uncertainties in the quasar (QSO) continuum placement. We further discuss constraints obtained by a combined analysis of the high-resolution flux PDF and the power spectrum measured from the Sloan Digital Sky Survey (SDSS) lya forest data. The joint analysis confirms the suggestion of an inverted temperature-density relation, but prefers somewhat higher values (sigma_8 ~ 0.9) of the matter fluctuation amplitude than the WMAP data and the best fit to the flux PDF alone. The joint analysis of the flux PDF and power spectrum (as well as an analysis of the power spectrum data alone) prefers rather large values for the temperature of the IGM, perhaps suggesting that we have identified a not yet accounted for systematic error in the SDSS flux power spectrum data or that the standard model describing the thermal state of the IGM at z ~ 2-3 is incomplete.
The Ly$alpha$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $langle z rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $gamma$, and temperature at mean-density, $T_0$, where $T(Delta) = T_0 Delta^{gamma-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $beta_mathrm{pLLS} sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$alpha$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $gamma=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($gamma leq 1$) are disfavored at a significance of over 4$sigma$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.
We present ANNz2, a new implementation of the public software for photometric redshift (photo-z) estimation of Collister and Lahav (2004), which now includes generation of full probability distribution functions (PDFs). ANNz2 utilizes multiple machine learning methods, such as artificial neural networks and boosted decision/regression trees. The objective of the algorithm is to optimize the performance of the photo-z estimation, to properly derive the associated uncertainties, and to produce both single-value solutions and PDFs. In addition, estimators are made available, which mitigate possible problems of non-representative or incomplete spectroscopic training samples. ANNz2 has already been used as part of the first weak lensing analysis of the Dark Energy Survey, and is included in the experiments first public data release. Here we illustrate the functionality of the code using data from the tenth data release of the Sloan Digital Sky Survey and the Baryon Oscillation Spectroscopic Survey. The code is available for download at https://github.com/IftachSadeh/ANNZ .
We introduce the position-dependent probability distribution function (PDF) of the smoothed matter field as a cosmological observable. In comparison to the PDF itself, the spatial variation of the position-dependent PDF is simpler to model and has distinct dependence on cosmological parameters. We demonstrate that the position-dependent PDF is characterized by variations in the local mean density, and we compute the linear response of the PDF to the local density using separate universe N-body simulations. The linear response of the PDF to the local density field can be thought of as the linear bias of regions of the matter field selected based on density. We provide a model for the linear response, which accurately predicts our simulation measurements. We also validate our results and test the separate universe consistency relation for the local PDF using global universe simulations. We find excellent agreement between the two, and we demonstrate that the separate universe method gives a lower variance determination of the linear response.
We examine the dust geometry and Ly{alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{alpha} emitting and starbursting systems with Hubble Space Telescope Ly{alpha}, H{alpha}, and H{beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{alpha}/H{alpha} and H{alpha}/H{beta} maps with spatial resolutions as low as $sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{alpha} photons only scatter. We show that the dust that affects the escape of Ly{alpha} is more restricted to the galaxies central regions, while the larger Ly{alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{alpha} scatters in the halo, and find that this characteristic scattering distance correlates with the measured size of the Ly{alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{alpha} and global dust properties.