We present a general expression for a lognormal filter given an arbitrary non-linear galaxy bias. We derive this filter as the maximum a posteriori solution assuming a lognormal prior distribution for the matter field with a given mean field and modelling the observed galaxy distribution by a Poissonian process. We have performed a three-dimensional implementation of this filter with a very efficient Newton-Krylov inversion scheme. Furthermore, we have tested it with a dark matter N-body simulation assuming a unit galaxy bias relation and compared the results with previous density field estimators like the inverse weighting scheme and Wiener filtering. Our results show good agreement with the underlying dark matter field for overdensities even above δ∼ 1000 which exceeds by one order of magnitude the regime in which the lognormal is expected to be valid. The reason is that for our filter the lognormal assumption enters as a prior distribution function, but the maximum a posteriori solution is also conditioned on the data. We find that the lognormal filter is superior to the previous filtering schemes in terms of higher correlation coefficients and smaller Euclidean distances to the underlying matter field. We also show how it is able to recover the positive tail of the matter overdensity field distribution for a unit bias relation down to scales of about ≳2 Mpc h-1. © 2010 The Authors. Journal compilation © 2010 RAS.

Recovering the non-linear density field from the galaxy distribution with a Poisson-lognormal filter

METCALF, ROBERT BENTON
2010

Abstract

We present a general expression for a lognormal filter given an arbitrary non-linear galaxy bias. We derive this filter as the maximum a posteriori solution assuming a lognormal prior distribution for the matter field with a given mean field and modelling the observed galaxy distribution by a Poissonian process. We have performed a three-dimensional implementation of this filter with a very efficient Newton-Krylov inversion scheme. Furthermore, we have tested it with a dark matter N-body simulation assuming a unit galaxy bias relation and compared the results with previous density field estimators like the inverse weighting scheme and Wiener filtering. Our results show good agreement with the underlying dark matter field for overdensities even above δ∼ 1000 which exceeds by one order of magnitude the regime in which the lognormal is expected to be valid. The reason is that for our filter the lognormal assumption enters as a prior distribution function, but the maximum a posteriori solution is also conditioned on the data. We find that the lognormal filter is superior to the previous filtering schemes in terms of higher correlation coefficients and smaller Euclidean distances to the underlying matter field. We also show how it is able to recover the positive tail of the matter overdensity field distribution for a unit bias relation down to scales of about ≳2 Mpc h-1. © 2010 The Authors. Journal compilation © 2010 RAS.
2010
Kitaura, Francisco-Shu; Jasche, Jens; Metcalf, R. Benton
File in questo prodotto:
Eventuali allegati, non sono esposti

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/593972
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 65
  • ???jsp.display-item.citation.isi??? 64
social impact