Journal Article

Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization

F. S. Kitaura and T. A. Enßlin

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 389, issue 2, pages 497-544
Published in print September 2008 | ISSN: 0035-8711
Published online September 2008 | e-ISSN: 1365-2966 | DOI:
Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization

Show Summary Details


We address the inverse problem of cosmic large-scale structure reconstruction from a Bayesian perspective. For a linear data model, a number of known and novel reconstruction schemes, which differ in terms of the underlying signal prior, data likelihood and numerical inverse extraregularization schemes are derived and classified. The Bayesian methodology presented in this paper tries to unify and extend the following methods: Wiener filtering, Tikhonov regularization, ridge regression, maximum entropy and inverse regularization techniques. The inverse techniques considered here are the asymptotic regularization, the Jacobi, Steepest Descent, Newton–Raphson, Landweber–Fridman and both linear and non-linear Krylov methods based on Fletcher–Reeves, Polak–Ribière and Hestenes–Stiefel conjugate gradients. The structures of the up-to-date highest performing algorithms are presented, based on an operator scheme, which permits one to exploit the power of fast Fourier transforms. Using such an implementation of the generalized Wiener filter in the novel argo software package, the different numerical schemes are benchmarked with one-, two- and three-dimensional problems including structured white and Poissonian noise, data windowing and blurring effects. A novel numerical Krylov scheme is shown to be superior in terms of performance and fidelity. These fast inverse methods ultimately will enable the application of sampling techniques to explore complex joint posterior distributions. We outline how the space of the dark matter density field, the peculiar velocity field and the power spectrum can jointly be investigated by a Gibbs-sampling process. Such a method can be applied for the redshift distortions correction of the observed galaxies and for time-reversal reconstructions of the initial density field.

Keywords: methods: data analysis; methods: numerical; methods: statistical; techniques: image processing; galaxies: distances and redshifts; large-scale structure of Universe

Journal Article.  24223 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.