Bayesian luminosity function estimation in multi-depth datasets with selection effects: A case study for $3<z<5$ Lyman $α$ emitters
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , , , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866914189825540096 |
|---|---|
| author | Tornotti, Davide Fossati, Matteo Fumagalli, Michele Gerosa, Davide Pizzuti, Lorenzo Battaia, Fabrizio Arrigoni |
| author_facet | Tornotti, Davide Fossati, Matteo Fumagalli, Michele Gerosa, Davide Pizzuti, Lorenzo Battaia, Fabrizio Arrigoni |
| contents | We present a hierarchical Bayesian framework designed to infer the luminosity function of any class of object by jointly modelling data from multiple surveys with varying depth, completeness, and sky coverage. Our method explicitly accounts for selection effects and measurement uncertainties (e.g. in luminosity) and can be generalized to any extensive quantity, such as mass. We validated the model using mock catalogues; from this we determined that deep data reaching $\gtrsim 1.5$ dex below a characteristic luminosity ($\tilde{L}^\star$) are essential to reducing biases at the faint end ($\lesssim 0.1$ dex) and that wide-area data help constrain the bright end. As a proof of concept, we considered a combined sample of 1176 Lyman $α$ emitters at redshift $3 < z < 5$ drawn from several MUSE surveys, ranging from ultra-deep ($\gtrsim 90$ hr) and narrow ($\lesssim 1$ arcmin$^2$) fields to shallow ($\lesssim 5$ hr) and wide ($\gtrsim 20$ arcmin$^2$) fields. With this complete sample, we constrain the luminosity function parameters $\log(Φ^\star/\mathrm{Mpc^{-3}}) = -2.86^{+0.15}_{-0.17}$, $\log(L^\star/\mathrm{erg\,s^{-1}}) = 42.72^{+0.10}_{-0.09}$, and $α= -1.81^{+0.09}_{-0.09}$, where the uncertainties represent the $90\%$ credible intervals. These values are in agreement with the results of studies based on gravitational lensing that reach $\log(L/\mathrm{erg\,s^{-1}}) \approx 41$, although differences in the faint-end slope underscore how systematic errors are starting to dominate. In contrast, wide-area surveys represent the natural extension needed to constrain the brightest Lyman $α$ emitters [$\log(L/\mathrm{erg\,s^{-1}}) \gtrsim 43$], where statistical uncertainties still dominate. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2506_10083 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Bayesian luminosity function estimation in multi-depth datasets with selection effects: A case study for $3<z<5$ Lyman $α$ emitters Tornotti, Davide Fossati, Matteo Fumagalli, Michele Gerosa, Davide Pizzuti, Lorenzo Battaia, Fabrizio Arrigoni Astrophysics of Galaxies Cosmology and Nongalactic Astrophysics We present a hierarchical Bayesian framework designed to infer the luminosity function of any class of object by jointly modelling data from multiple surveys with varying depth, completeness, and sky coverage. Our method explicitly accounts for selection effects and measurement uncertainties (e.g. in luminosity) and can be generalized to any extensive quantity, such as mass. We validated the model using mock catalogues; from this we determined that deep data reaching $\gtrsim 1.5$ dex below a characteristic luminosity ($\tilde{L}^\star$) are essential to reducing biases at the faint end ($\lesssim 0.1$ dex) and that wide-area data help constrain the bright end. As a proof of concept, we considered a combined sample of 1176 Lyman $α$ emitters at redshift $3 < z < 5$ drawn from several MUSE surveys, ranging from ultra-deep ($\gtrsim 90$ hr) and narrow ($\lesssim 1$ arcmin$^2$) fields to shallow ($\lesssim 5$ hr) and wide ($\gtrsim 20$ arcmin$^2$) fields. With this complete sample, we constrain the luminosity function parameters $\log(Φ^\star/\mathrm{Mpc^{-3}}) = -2.86^{+0.15}_{-0.17}$, $\log(L^\star/\mathrm{erg\,s^{-1}}) = 42.72^{+0.10}_{-0.09}$, and $α= -1.81^{+0.09}_{-0.09}$, where the uncertainties represent the $90\%$ credible intervals. These values are in agreement with the results of studies based on gravitational lensing that reach $\log(L/\mathrm{erg\,s^{-1}}) \approx 41$, although differences in the faint-end slope underscore how systematic errors are starting to dominate. In contrast, wide-area surveys represent the natural extension needed to constrain the brightest Lyman $α$ emitters [$\log(L/\mathrm{erg\,s^{-1}}) \gtrsim 43$], where statistical uncertainties still dominate. |
| title | Bayesian luminosity function estimation in multi-depth datasets with selection effects: A case study for $3<z<5$ Lyman $α$ emitters |
| topic | Astrophysics of Galaxies Cosmology and Nongalactic Astrophysics |
| url | https://arxiv.org/abs/2506.10083 |