Related papers: The Information Content of Projected Galaxy Fields
Hamiltonian Monte Carlo (HMC) has been progressively incorporated within the statistician's toolbox as an alternative sampling method in settings when standard Metropolis-Hastings is inefficient. HMC generates a Markov chain on an augmented…
We derive and implement a full Bayesian large scale structure inference method aiming at precision recovery of the cosmological power spectrum from galaxy redshift surveys. Our approach improves over previous Bayesian methods by performing…
We present a novel approach for reconstructing the projected mass distribution of clusters of galaxies from sparse and noisy weak gravitational lensing shear data. The reconstructions are regularised using knowledge gained from numerical…
We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely…
We present a method to recover and study the projected gravitational tidal forces from a galaxy survey containing little or no redshift information. The method and the physical interpretation of the recovered tidal maps as a tracer of the…
We revisit the issue of non-parametric gravitational lens reconstruction and present a new method to obtain the cluster mass distribution using strong lensing data without using any prior information on the underlying mass. The method…
The combination of two- and three-point clustering statistics of galaxies and the underlying matter distribution has the potential to break degeneracies between cosmological parameters and nuisance parameters and can lead to significantly…
I present a new approach to recover the primordial density fluctuations and the cosmic web structure underlying a galaxy distribution. The method is based on sampling Gaussian fields which are compatible with a galaxy distribution and a…
This paper presents a new Expectation Propagation (EP) framework for image restoration using patch-based prior distributions. While Monte Carlo techniques are classically used to sample from intractable posterior distributions, they can…
We apply a mass reconstruction technique to simulated large-scale structure gravitational distortion maps, from 2.5' to 10 degree scales, for different cosmological scenarii. The projected mass is reconstructed using a non-parametric least…
We revise the Bayesian inference steps required to analyse the cosmological large-scale structure. Here we make special emphasis in the complications which arise due to the non-Gaussian character of the galaxy and matter distribution. In…
We investigate the effect of non-Gaussianity on the reconstruction of the initial mass field from the Ly$\alpha$ forest. We show that the transmitted flux of QSO absorption spectra are highly non-Gaussian in terms of the statistics, the…
We propose a method to estimate non-Gaussian error bars on the matter power spectrum from galaxy surveys in the presence of non-trivial survey selection functions. The estimators are often obtained from formalisms like FKP and PKL, which…
We present a method for fast optimal estimation of the temperature angular power spectrum from observations of the cosmic microwave background. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from the posterior…
In the $\Lambda$CDM model, cosmological observations from the late and recent universe reveal a puzzling $\sim 4.5\sigma$ tension in the current rate of universe expansion. In addition to the various scenarios suggested to resolve the…
Using a set of 28 high resolution, high signal to noise ratio (S/N) QSO Ly-alpha absorption spectra, we investigate the non-Gaussian features of the transmitted flux fluctuations, and their effect upon the power spectrum of this field. We…
Structure formation in our Universe creates non-Gaussian random fields that will soon be observed over almost the entire sky by the Euclid satellite, the Vera-Rubin observatory, and the Square Kilometre Array. An unsolved problem is how to…
The complex Gaussian distribution has been widely used as a fundamental spectral and noise model in signal processing and communication. However, its Gaussian structure often limits its ability to represent the diverse amplitude…
Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the…
A method for the recovery of the real space line-of-sight mass density field from Lyman absorption in QSO spectra is presented. The method makes use of a Lucy-type algorithm for the recovery of the HI density. The matter density is inferred…