Journal Article

Multiscale inference of matter fields and baryon acoustic oscillations from the Lyα forest

Francisco-Shu Kitaura, Simona Gallerani and Andrea Ferrara

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 420, issue 1, pages 61-74
Published in print February 2012 | ISSN: 0035-8711
Published online January 2012 | e-ISSN: 1365-2966 | DOI:
Multiscale inference of matter fields and baryon acoustic oscillations from the Lyα forest

More Like This

Show all results sharing this subject:

  • Astronomy and Astrophysics


Show Summary Details


We present a novel Bayesian method for the joint reconstruction of cosmological matter density fields, peculiar velocities and power spectra in the quasi-non-linear regime. We study its applicability to the Lyα forest based on multiple quasar absorption spectra. Our approach to this problem includes a multiscale, non-linear, two-step scheme since the statistics describing the matter distribution depends on scale, being strongly non-Gaussian on small scales (<0.1 h−1 Mpc) and closely lognormal on scales ≳10 h−1 Mpc. The first step consists on performing one-dimensional highly resolved matter density reconstructions along the line of sight towards z∼ 2–3 quasars based on an arbitrary non-Gaussian univariate model for matter statistics. The second step consists on Gibbs sampling based on conditional probability distribution functions. The matter density field is sampled in real space with Hamiltonian sampling using the Poisson/Gamma lognormal model, while redshift distortions are corrected with linear Lagrangian perturbation theory. The power spectrum of the lognormal transformed variable which is Gaussian distributed (and thus close to the linear regime) can consistently be sampled with the inverse Gamma distribution function. We test our method through numerical N-body simulations with a computational volume large enough (>1 h−3 Gpc3) to show that the linear power spectra are nicely recovered over scales larger than ≳20 h−1 Mpc, i.e. the relevant range where features imprinted by the baryon acoustics oscillations appear.

Keywords: methods: data analysis; methods: statistical; quasars: absorption lines; large-scale structure of Universe

Journal Article.  10532 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.