Source author record

A. Rassat

A. Rassat appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2015arXiv

3D galaxy clustering with future wide-field surveys: Advantages of a spherical Fourier-Bessel analysis

Upcoming spectroscopic galaxy surveys are extremely promising to help in addressing the major challenges of cosmology, in particular in understanding the nature of the dark universe. The strength of these surveys comes from their unprecedented depth and width. Optimal extraction of their three-dimensional information is of utmost importance to best constrain the properties of the dark universe. Although there is theoretical motivation and novel tools to explore these surveys using the 3D spherical Fourier-Bessel (SFB) power spectrum of galaxy number counts $C_\ell(k,k^\prime)$, most survey optimisations and forecasts are based on the tomographic spherical harmonics power spectrum $C^{(ij)}_\ell$. We performed a new investigation of the information that can be extracted from the tomographic and 3D SFB techniques by comparing the forecast cosmological parameter constraints obtained from a Fisher analysis in the context of planned stage IV wide-field galaxy surveys. The comparison was made possible by careful and coherent treatment of non-linear scales in the two analyses. Nuisance parameters related to a scale- and redshift-dependent galaxy bias were also included for the first time in the computation of both the 3D SFB and tomographic power spectra. Tomographic and 3D SFB methods can recover similar constraints in the absence of systematics. However, constraints from the 3D SFB analysis are less sensitive to unavoidable systematics stemming from a redshift- and scale-dependent galaxy bias. Even for surveys that are optimised with tomography in mind, a 3D SFB analysis is more powerful. In addition, for survey optimisation, the figure of merit for the 3D SFB method increases more rapidly with redshift, especially at higher redshifts, suggesting that the 3D SFB method should be preferred for designing and analysing future wide-field spectroscopic surveys.

preprint2014arXiv

Joint Planck and WMAP CMB Map Reconstruction

We present a novel estimate of the cosmological microwave background (CMB) map by combining the two latest full-sky microwave surveys: WMAP nine-year and Planck PR1. The joint processing benefits from a recently introduced component separation method coined "local-generalized morphological component analysis'' (LGMCA) based on the sparse distribution of the foregrounds in the wavelet domain. The proposed estimation procedure takes advantage of the IRIS 100 micron as an extra observation on the galactic center for enhanced dust removal. We show that this new CMB map presents several interesting aspects: i) it is a full sky map without using any inpainting or interpolating method, ii) foreground contamination is very low, iii) the Galactic center is very clean, with especially low dust contamination as measured by the cross-correlation between the estimated CMB map and the IRIS 100 micron map, and iv) it is free of thermal SZ contamination.

preprint2014arXiv

Planck CMB Anomalies: Astrophysical and Cosmological Secondary Effects and the Curse of Masking

Large-scale anomalies have been reported in CMB data with both WMAP and Planck data. These could be due to foreground residuals and or systematic effects, though their confirmation with Planck data suggests they are not due to a problem in the WMAP or Planck pipelines. If these anomalies are in fact primordial, then understanding their origin is fundamental to either validate the standard model of cosmology or to explore new physics. We investigate three other possible issues: 1) the trade-off between minimising systematics due to foreground contamination (with a conservative mask) and minimising systematics due to masking, 2) astrophysical secondary effects (the kinetic Doppler quadrupole and kinetic Sunyaev-Zel'dovich effect), and 3) secondary cosmological signals (the integrated Sachs-Wolfe effect). We address the masking issue by considering new procedures that use both WMAP and Planck to produce higher quality full-sky maps using the sparsity methodology (LGMCA maps). We show the impact of masking is dominant over that of residual foregrounds, and the LGMCA full-sky maps can be used without further processing to study anomalies. We consider four official Planck PR1 and two LGMCA CMB maps. Analysis of the observed CMB maps shows that only the low quadrupole and quadrupole-octopole alignment seem significant, but that the planar octopole, Axis of Evil, mirror parity and cold spot are not significant in nearly all maps considered. After subtraction of astrophysical and cosmological secondary effects, only the low quadrupole may still be considered anomalous, meaning the significance of only one anomaly is affected by secondary effect subtraction out of six anomalies considered. In the spirit of reproducible research all reconstructed maps and codes will be made available for download here http://www.cosmostat.org/anomaliesCMB.html.

preprint2014arXiv

PRISM: Recovery of the primordial spectrum from Planck data

The primordial power spectrum describes the initial perturbations that seeded the large-scale structure we observe today. It provides an indirect probe of inflation or other structure-formation mechanisms. In this letter, we recover the primordial power spectrum from the Planck PR1 dataset, using our recently published algorithm PRISM. PRISM is a sparsity-based inversion method, that aims at recovering features in the primordial power spectrum from the empirical power spectrum of the cosmic microwave background (CMB). This ill-posed inverse problem is regularised using a sparsity prior on features in the primordial power spectrum in a wavelet dictionary. Although this non-parametric method does not assume a strong prior on the shape of the primordial power spectrum, it is able to recover both its general shape and localised features. As a results, this approach presents a reliable way of detecting deviations from the currently favoured scale-invariant spectrum. We applied PRISM to 100 simulated Planck data to investigate its performance on Planck-like data. We also tested the algorithm's ability to recover a small localised feature at $k \sim 0.125$ Mpc$^{-1}$, which caused a large dip at $\ell \sim 1800$ in the angular power spectrum. We then applied PRISM to the Planck PR1 power spectrum to recover the primordial power spectrum. We find no significant departures from the fiducial Planck PR1 near scale-invariant primordial power spectrum with $A_s=2.215\times10^{-9}$ and $n_s = 0.9624$.

preprint2014arXiv

Sparse point-source removal for full-sky CMB experiments: application to WMAP 9-year data

Missions such as WMAP or Planck measure full-sky fluctuations of the cosmic microwave background and foregrounds, among which bright compact source emissions cover a significant fraction of the sky. To accurately estimate the diffuse components, the point-source emissions need to be separated from the data, which requires a dedicated processing. We propose a new technique to estimate the flux of the brightest point sources using a morphological separation approach: point sources with known support and shape are separated from diffuse emissions that are assumed to be sparse in the spherical harmonic domain. This approach is compared on both WMAP simulations and data with the standard local chi2 minimization, modelling the background as a low-order polynomial. The proposed approach generally leads to 1) lower biases in flux recovery, 2) an improved root mean-square error of up to 35% and 3) more robustness to background fluctuations at the scale of the source. The WMAP 9-year point-source-subtracted maps are available online.

preprint2013arXiv

On Preferred Axes in WMAP Cosmic Microwave Background Data after Subtraction of the Integrated Sachs-Wolfe Effect

There is currently a debate over the existence of claimed statistical anomalies in the cosmic microwave background (CMB), recently confirmed in Planck data. Recent work has focussed on methods for measuring statistical significance, on masks and on secondary anisotropies as potential causes of the anomalies. We investigate simultaneously the method for accounting for masked regions and the foreground integrated Sachs-Wolfe (ISW) signal. We search for trends in different years of WMAP CMB data with different mask treatments. We reconstruct the ISW field due to the 2 Micron All-Sky Survey (2MASS) and the NRAO VLA Sky Survey (NVSS) up to l=5, and we focus on the Axis of Evil (AoE) statistic and even/odd mirror parity, both of which search for preferred axes in the Universe. We find that removing the ISW reduces the significance of these anomalies in WMAP data, though this does not exclude the possibility of exotic physics. In the spirit of reproducible research, all reconstructed maps and codes will be made available for download at http://www.cosmostat.org/anomaliesCMB.html.

preprint2013arXiv

Removal of two large scale Cosmic Microwave Background anomalies after subtraction of the Integrated Sachs Wolfe effect

Though debated, the existence of claimed large-scale anomalies in the CMB is not totally dismissed. In parallel to the debate over their statistical significance, recent work focussed on masks and secondary anisotropies as potential sources of these anomalies. In this work we investigate simultaneously the impact of the method used to account for masked regions and the impact of the ISW effect, which is the large-scale secondary anisotropy most likely to affect the CMB anomalies. In this sense, our work is an update of both Francis & Peacock 2010 and Kim et al. 2012. Our aim is to identify trends in CMB data with different mask treatments. We reconstruct the ISW signal due to 2MASS and NVSS galaxies. We account for missing data using the sparse inpainting technique of Abrial et al. 2008 and sparse inpainting of the CMB, LSS and ISW and find that it constitutes a bias-free reconstruction method suitable to study large-scale statistical isotropy and the ISW effect. We focus on three large-scale CMB anomalies: the low quadrupole, the quad/oct alignment, and the octopole planarity. After sparse inpainting, the low quadrupole becomes more anomalous, whilst the quad/oct alignment becomes less anomalous. The significance of the low quadrupole is unchanged after subtraction of the ISW effect, while the trend is that the quad/oct alignment has reduced significance, yet other hypotheses remain possible as well (e.g. exotic physics). Our results also suggest that both of these anomalies may be due to the quadrupole alone. The octopole planarity significance is also reduced after inpainting and after ISW subtraction, however, we do not find that it was very anomalous to start with. In the spirit of reproducible research, we make all codes and resulting products which constitute main results of this paper public here: http://www.cosmostat.org/anomaliesCMB.html.

preprint2013arXiv

Sparsity and the Bayesian Perspective

Sparsity has been recently introduced in cosmology for weak-lensing and CMB data analysis for different applications such as denoising, component separation or inpainting (i.e. filling the missing data or the mask). Although it gives very nice numerical results, CMB sparse inpainting has been severely criticized by top researchers in cosmology, based on arguments derived from a Bayesian perspective. Trying to understand their point of view, we realize that interpreting a regularization penalty term as a prior in a Bayesian framework can lead to erroneous conclusions. This paper is by no means against the Bayesian approach, which has proven to be very useful for many applications, but warns about a Bayesian-only interpretation in data analysis, which can be misleading in some cases.

preprint2013arXiv

WMAP 9-year CMB estimation using sparsity

Recovering the Cosmic Microwave Background (CMB) from WMAP data requires galactic foreground emissions to be accurately separated out. Most component separation techniques rely on second order statistics such as Internal Linear Combination (ILC) techniques. In this paper, we present a new WMAP 9-year CMB map, with 15 arcmin resolution, which is reconstructed using a recently introduced sparse component separation technique, coined Local Generalized Morphological Component Analysis (LGMCA). LGMCA emphasizes on the sparsity of the components to be retrieved in the wavelet domain. We show that although derived from a radically different separation criterion ({i.e. sparsity), the LGMCA-WMAP 9 map and its power spectrum are fully consistent with their more recent estimates from WMAP 9.

preprint2012arXiv

3D Spherical Analysis of Baryon Acoustic Oscillations

Baryon Acoustic Oscillations (BAOs) are oscillatory features in the galaxy power spectrum and are a standard rod to measure the cosmological expansion. These have been studied in Cartesian space (Fourier or real space) or in Spherical Harmonic (SH) space in thin shells. Future wide-field surveys will cover both wide and deep regions of the sky and thus require a simultaneous treatment of the spherical sky and of an extended radial coverage. The Spherical Fourier-Bessel (SFB) decomposition is a natural basis for the analysis of fields in this geometry and facilitates the combination of BAO surveys with other cosmological probes readily described in this basis. We present here a new way to analyse BAOs by studying the BAO wiggles from the SFB power spectrum. In SFB space, the power spectrum generally has both a radial (k) and tangential (l) dependence and so do the BAOs. In the deep survey limit and ignoring evolution, the SFB power spectrum becomes radial and reduces to the Cartesian Fourier power spectrum. In the limit of a thin shell, all the information is contained in the tangential modes described by the 2D SH power spectrum. We find that the radialisation of the SFB power spectrum is still a good approximation even when considering an evolving and biased galaxy field with a finite selection function. This effect can be observed by all-sky surveys with depths comparable to current surveys. We find that the BAOs radialise more rapidly than the full SFB power spectrum. Our results suggest the first peak of the BAOs in SFB space becomes radial out to l ~ 10 for all-sky surveys with the same depth as SDSS or 2dF, and out to l ~ 70 for an all-sky stage IV survey. Subsequent BAO peaks also become radial, but for shallow surveys these may be in the non-linear regime. For modes that have become radial, measurements at different l's are useful in practice to reduce measurement errors.

preprint2012arXiv

3DEX: a code for fast spherical Fourier-Bessel decomposition of 3D surveys

High-precision cosmology requires the analysis of large-scale surveys in 3D spherical coordinates, i.e. spherical Fourier-Bessel decomposition. Current methods are insufficient for future data-sets from wide-field cosmology surveys. The aim of this paper is to present a public code for fast spherical Fourier-Bessel decomposition that can be applied to cosmological data or 3D data in spherical coordinates in other scientific fields. We present an equivalent formulation of the spherical Fourier-Bessel decomposition that separates radial and tangential calculations. We propose the use of the existing pixelisation scheme HEALPix for a rapid calculation of the tangential modes. 3DEX (3D EXpansions) is a public code for fast spherical Fourier-Bessel decomposition of 3D all-sky surveys that takes advantage of HEALPix for the calculation of tangential modes. We perform tests on very large simulations and we compare the precision and computation time of our method with an optimised implementation of the spherical Fourier-Bessel original formulation. For surveys with millions of galaxies, computation time is reduced by a factor 4-12 depending on the desired scales and accuracy. The formulation is also suitable for pre-calculations and external storage of the spherical harmonics, which allows for additional speed improvements. The 3DEX code can accommodate data with masked regions of missing data. 3DEX can also be used in other disciplines, where 3D data are to be analysed in spherical coordinates. The code and documentation can be downloaded at http://ixkael.com/blog/3dex.

preprint2012arXiv

Low-l CMB Analysis and Inpainting

Reconstruction of the CMB in the Galactic plane is extremely difficult due to the dominant foreground emissions such as Dust, Free-Free or Synchrotron. For cosmological studies, the standard approach consists in masking this area where the reconstruction is not good enough. This leads to difficulties for the statistical analysis of the CMB map, especially at very large scales (to study for e.g., the low quadrupole, ISW, axis of evil, etc). We investigate in this paper how well some inpainting techniques can recover the low-$\ell$ spherical harmonic coefficients. We introduce three new inpainting techniques based on three different kinds of priors: sparsity, energy and isotropy, and we compare them. We show that two of them, sparsity and energy priors, can lead to extremely high quality reconstruction, within 1% of the cosmic variance for a mask with Fsky larger than 80%.

preprint2012arXiv

Spherical 3D Isotropic Wavelets

Future cosmological surveys will provide 3D large scale structure maps with large sky coverage, for which a 3D Spherical Fourier-Bessel (SFB) analysis in spherical coordinates is natural. Wavelets are particularly well-suited to the analysis and denoising of cosmological data, but a spherical 3D isotropic wavelet transform does not currently exist to analyse spherical 3D data. The aim of this paper is to present a new formalism for a spherical 3D isotropic wavelet, i.e. one based on the SFB decomposition of a 3D field and accompany the formalism with a public code to perform wavelet transforms. We describe a new 3D isotropic spherical wavelet decomposition based on the undecimated wavelet transform (UWT) described in Starck et al. 2006. We also present a new fast Discrete Spherical Fourier-Bessel Transform (DSFBT) based on both a discrete Bessel Transform and the HEALPIX angular pixelisation scheme. We test the 3D wavelet transform and as a toy-application, apply a denoising algorithm in wavelet space to the Virgo large box cosmological simulations and find we can successfully remove noise without much loss to the large scale structure. We have described a new spherical 3D isotropic wavelet transform, ideally suited to analyse and denoise future 3D spherical cosmological surveys, which uses a novel Discrete Spherical Fourier-Bessel Transform. We illustrate its potential use for denoising using a toy model. All the algorithms presented in this paper are available for download as a public code called MRS3D at http://jstarck.free.fr/mrs3d.html

preprint2011arXiv

Euclid Definition Study Report

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

preprint2011arXiv

Measuring the Integrated Sachs-Wolfe Effect

One of the main challenges of modern cosmology is to understand the nature of dark energy. The Integrated Sachs-Wolfe (ISW) effect is sensitive to dark energy and presents an independent signature of dark energy in the absence of modified gravity and curvature. The ISW effect occurs on large scales, where cosmic variance is high and where there are large amounts of missing data in the CMB and large scale structure maps due to Galactic confusion. Moreover, existing methods in the literature often make strong assumptions about the statistics of the underlying fields or estimators. Together these effects can severely limit signal extraction. We review literature on the ISW effect, comparing statistical subtleties between existing methods, and identifying several limitations. We propose a novel method to detect and measure the ISW signal. This method assumes only that the primordial CMB field is Gaussian. It is based on a sparse inpainting method to reconstruct missing data and uses a bootstrap technique to avoid assumptions about the statistics of the estimator. It is a complete method, using three complementary statistical methods. We apply our method to Euclid-like simulations and show we can expect a \sim 7σmodel-independent detection of the ISW signal with WMAP7-like data, even with missing data. Other tests return \sim 4.7σdetection levels for a Euclid-like survey, with levels independent from whether the galaxy field is normally or lognormally distributed. We apply our method to the 2 Micron All Sky Survey (2MASS) and WMAP7 CMB data and find detections in the 1.1 - 2.0σrange, as expected from simulations. As a by-product, we reconstruct the full-sky temperature ISW field due to 2MASS data. We have presented a novel technique, based on sparse inpainting and bootstrapping, which accurately detects and reconstructs the ISW effect.