Showing 1–50 of 60 results for author: Marinucci, D

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Geometric Methods for Spherical Data, with Applications to Cosmology

Authors: Javier Carrón Duque, Domenico Marinucci

Abstract: This survey is devoted to recent developments in the statistical analysis of spherical data, with a view to applications in Cosmology. We will start from a brief discussion of Cosmological questions and motivations, arguing that most Cosmological observables are spherical random fields. Then, we will introduce some mathematical background on spherical random fields, including spectral representati… ▽ More This survey is devoted to recent developments in the statistical analysis of spherical data, with a view to applications in Cosmology. We will start from a brief discussion of Cosmological questions and motivations, arguing that most Cosmological observables are spherical random fields. Then, we will introduce some mathematical background on spherical random fields, including spectral representations and the construction of needlet and wavelet frames. We will then focus on some specific issues, including tools and algorithms for map reconstruction (\textit{i.e.}, separating the different physical components which contribute to the observed field), geometric tools for testing the assumptions of Gaussianity and isotropy, and multiple testing methods to detect contamination in the field due to point sources. Although these tools are introduced in the Cosmological context, they can be applied to other situations dealing with spherical data. Finally, we will discuss more recent and challenging issues such as the analysis of polarization data, which can be viewed as realizations of random fields taking values in spin fiber bundles. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 25 pages, 6 figures

Minkowski Functionals in $SO(3)$ for the spin-2 CMB polarisation field

Authors: Javier Carrón Duque, Alessandro Carones, Domenico Marinucci, Marina Migliaccio, Nicola Vittorio

Abstract: The study of the angular power spectrum of Cosmic Microwave Background (CMB) anisotropies, both in intensity and in polarisation, has led to the tightest constraints on cosmological parameters. However, this statistical quantity is not sensitive to any deviation from Gaussianity and statistical isotropy in the CMB data. Minkowski Functionals (MFs) have been adopted as one of the most powerful stat… ▽ More The study of the angular power spectrum of Cosmic Microwave Background (CMB) anisotropies, both in intensity and in polarisation, has led to the tightest constraints on cosmological parameters. However, this statistical quantity is not sensitive to any deviation from Gaussianity and statistical isotropy in the CMB data. Minkowski Functionals (MFs) have been adopted as one of the most powerful statistical tools to study such deviations, since they characterise the topology and geometry of the field of interest. In this paper, we extend the application of MFs to CMB polarisation data by introducing a new formalism, where we lift the spin $2$ polarisation field to a scalar function in a higher-dimensional manifold: the group of rotations of the sphere, $SO(3)$. Such a function is defined as $f = Q \cos(2ψ) - U \sin(2ψ)$. We analytically obtain the expected values for the MFs of $f$ in the case of Gaussian isotropic polarisation maps. Furthermore, we present a new pipeline which estimates these MFs from input HEALPix polarisation maps. We apply it to CMB simulations in order to validate the theoretical results and the methodology. The pipeline is to be included in the publicly available Python package $\texttt{Pynkowski}$ available at https://github.com/javicarron/pynkowski. △ Less

Submitted 7 June, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 28 pages, 2 figures

Multi-Clustering Needlet-ILC for CMB B-modes component separation

Authors: Alessandro Carones, Marina Migliaccio, Giuseppe Puglisi, Carlo Baccigalupi, Domenico Marinucci, Nicola Vittorio, Davide Poletti

Abstract: The Cosmic Microwave Background (CMB) primordial B-modes signal is predicted to be much lower than the polarized Galactic emission (foregrounds) in any region of the sky pointing to the need for sophisticated component separation methods. Among them, the blind Needlet-ILC (NILC) has great relevance given our current poor knowledge of the B-modes foregrounds. However the expected level of spatial v… ▽ More The Cosmic Microwave Background (CMB) primordial B-modes signal is predicted to be much lower than the polarized Galactic emission (foregrounds) in any region of the sky pointing to the need for sophisticated component separation methods. Among them, the blind Needlet-ILC (NILC) has great relevance given our current poor knowledge of the B-modes foregrounds. However the expected level of spatial variability of the foreground spectral properties complicates the NILC subtraction of the Galactic contamination. In order to reach the ambitious targets of future CMB experiments, we therefore propose a novel extension of the NILC approach, the Multi-Clustering NILC (MC-NILC), which performs NILC variance minimization on separate regions of the sky (clusters) properly chosen to have similar spectral properties of the B-modes foregrounds emission. Clusters are identified thresholding the ratio of B-modes maps at two separate frequencies which is used as tracer of the spatial distribution of the spectral indices of the Galactic emission in B modes. We consider ratios either of simulated foregrounds-only B modes (ideal case) or of cleaned templates of Galactic emission obtained from realistic simulations. In this work we present an application of MC-NILC to the future LiteBIRD satellite, which targets the observation of both reionization and recombination peaks of the primordial B-modes angular power spectrum with a total error on the tensor-to-scalar ratio $δr < 0.001$. We show that MC-NILC provides a CMB solution with residual foregrounds and noise contamination that is significantly reduced with respect to NILC and lower than the primordial signal targeted by LiteBIRD at all angular scales for the ideal case and at the reionization peak for a realistic ratio. Thus, MC-NILC will represent a powerful method to mitigate B-modes foregrounds for future CMB polarization experiments. △ Less

Submitted 25 January, 2024; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 19 pages, 17 figures, updated to match version accepted by MNRAS. Abstract abridged for arXiv

Journal ref: MNRAS, Volume 525, Issue 2, October 2023, pp.3117-3135

Minkowski Functionals of CMB polarisation intensity with Pynkowski: theory and application to Planck and future data

Authors: Alessandro Carones, Javier Carrón Duque, Domenico Marinucci, Marina Migliaccio, Nicola Vittorio

Abstract: The angular power spectrum of the Cosmic Microwave Background (CMB) anisotropies is a key tool to study the Universe. However, it is blind to the presence of non--Gaussianities and deviations from statistical isotropy, which instead can be detected with other statistics such as Minkowski Functionals (MFs). These tools have been applied to CMB temperature and $E$-mode anisotropies with no detection… ▽ More The angular power spectrum of the Cosmic Microwave Background (CMB) anisotropies is a key tool to study the Universe. However, it is blind to the presence of non--Gaussianities and deviations from statistical isotropy, which instead can be detected with other statistics such as Minkowski Functionals (MFs). These tools have been applied to CMB temperature and $E$-mode anisotropies with no detection of deviations from Gaussianity and isotropy. In this work, we extend the MFs formalism to the CMB polarisation intensity, $P^2=Q^2+U^2$. We use the Gaussian Kinematic Formula to derive the theoretical predictions of MFs for Gaussian isotropic fields. We develop a software that computes MFs on $P^2$ HEALPix maps and apply it to simulations to verify the robustness of both theory and methodology. We then estimate MFs of $P^2$ maps from Planck, both in pixel space and needlet domain, comparing them with realistic simulations which include CMB and instrumental noise residuals. We find no significant deviations from Gaussianity or isotropy in Planck CMB polarisation intensity. However, MFs could play an important role in the analysis of CMB polarisation measurements from upcoming experiments with improved sensitivity. Therefore we forecast the ability of MFs applied to $P^2$ maps to detect much fainter non-Gaussian anisotropic signals than with Planck data for two future complementary experiments: the LiteBIRD satellite and the ground-based Simons Observatory. We publicly release the software to compute MFs in arbitrary scalar HEALPix maps as a fully-documented Python package called $\texttt{Pynkowski}$ (https://github.com/javicarron/pynkowski). △ Less

Submitted 25 January, 2024; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 18 pages, 9 figures, updated to match version accepted by MNRAS. $\texttt{Pynkowski}$ can be found at https://github.com/javicarron/pynkowski

Journal ref: MNRAS, Volume 527, Issue 1, January 2024, pp.756-773

Analysis of NILC performance on B-modes data of sub-orbital experiments

Authors: Alessandro Carones, Marina Migliaccio, Domenico Marinucci, Nicola Vittorio

Abstract: The observation of primordial B-modes in the Cosmic Microwave Background (CMB) represents the main scientific goal of most of the future CMB experiments. This signal is predicted to be much lower than polarised Galactic emission (foregrounds) in any region of the sky pointing to the need for effective components separation methods, such as the Needlet-ILC (NILC). In this work, we explore the possi… ▽ More The observation of primordial B-modes in the Cosmic Microwave Background (CMB) represents the main scientific goal of most of the future CMB experiments. This signal is predicted to be much lower than polarised Galactic emission (foregrounds) in any region of the sky pointing to the need for effective components separation methods, such as the Needlet-ILC (NILC). In this work, we explore the possibility of employing NILC for B-mode maps reconstructed from partial-sky data of sub-orbital experiments, addressing the complications that such an application yields: E-B leakage, needlet filtering and beam convolution. We consider two complementary simulated datasets from future experiments: the balloon-borne SWIPE telescope of the Large Scale Polarization Explorer, which targets the observation of both reionisation and recombination peaks of the primordial B-mode angular power spectrum, and the ground-based Small Aperture Telescope of Simons Observatory, which is designed to observe only the recombination bump. We assess the performance of two alternative techniques to correct for the CMB E-B leakage: the recycling technique (Liu et al. 2019) and the ZB method (Zhao & Baskaran 2010). We find that they both reduce the E-B leakage residuals at a negligible level given the sensitivity of the considered experiments, except for the recycling method in the SWIPE patch at $\ell < 20$. Thus, we implement two extensions of the pipeline, the iterative B-decomposition and the diffusive inpainting, which enable us to recover the input CMB B-mode power for $\ell \geq 5$. We demonstrate that needlet filtering and beam convolution do not affect the B-mode reconstruction. Finally, with an appropriate masking strategy, we find that NILC foregrounds subtraction allows to achieve sensitivities for the tensor-to-scalar ratio compatible to the targets of the considered CMB experiments. △ Less

Submitted 25 January, 2024; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 23 pages, 17 figures, updated to match version accepted by A&A. Abstract abridged for Arxiv

Journal ref: A&A 677, A147 (2023)

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

Authors: Euclid Collaboration, S. Ilić, N. Aghanim, C. Baccigalupi, J. R. Bermejo-Climent, G. Fabbian, L. Legrand, D. Paoletti, M. Ballardini, M. Archidiacono, M. Douspis, F. Finelli, K. Ganga, C. Hernández-Monteagudo, M. Lattanzi, D. Marinucci, M. Migliaccio, C. Carbone, S. Casas, M. Martinelli, I. Tutusaus, P. Natoli, P. Ntelis, L. Pagano, L. Wenzl , et al. (185 additional authors not shown)

Abstract: The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological para… ▽ More The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data. △ Less

Submitted 10 September, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 20 pages, 8 figures, 5 tables, 1 appendix; updated to match version accepted by journal

Journal ref: A&A 657, A91 (2022)

A novel Cosmic Filament catalogue from SDSS data

Authors: Javier Carrón Duque, Marina Migliaccio, Domenico Marinucci, Nicola Vittorio

Abstract: In this work we present a new catalogue of Cosmic Filaments obtained from the latest Sloan Digital Sky Survey (SDSS) public data. In order to detect filaments, we implement a version of the Subspace-Constrained Mean-Shift algorithm, boosted by Machine Learning techniques. This allows us to detect cosmic filaments as one-dimensional maxima in the galaxy density distribution. Our filament catalogue… ▽ More In this work we present a new catalogue of Cosmic Filaments obtained from the latest Sloan Digital Sky Survey (SDSS) public data. In order to detect filaments, we implement a version of the Subspace-Constrained Mean-Shift algorithm, boosted by Machine Learning techniques. This allows us to detect cosmic filaments as one-dimensional maxima in the galaxy density distribution. Our filament catalogue uses the cosmological sample of SDSS, including Data Release 16, so it inherits its sky footprint (aside from small border effects) and redshift coverage. In particular, this means that, taking advantage of the quasar sample, our filament reconstruction covers redshifts up to $z=2.2$, making it one of the deepest filament reconstructions to our knowledge. We follow a tomographic approach and slice the galaxy data in 269 shells at different redshift. The reconstruction algorithm is applied to 2D spherical maps. The catalogue provides the position and uncertainty of each detection for each redshift slice. We assess the quality of the detections with several metrics, which show improvement with respect to previous public catalogues obtained with similar methods. We also detect a highly significant correlation between our filament catalogue and galaxy cluster catalogues built from microwave observations of the Planck Satellite and the Atacama Cosmology Telescope. △ Less

Submitted 27 April, 2022; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 23 pages, 20 figures, version accepted for publication in A&A

Journal ref: A&A 659, A166 (2022)

Needlet thresholding methods in component separation

Authors: F. Oppizzi, A. Renzi, M. Liguori, F. K. Hansen, D. Marinucci, C. Baccigalupi, D. Bertacca, D. Poletti

Abstract: Foreground components in the Cosmic Microwave Background (CMB) are sparse in a needlet representation, due to their specific morphological features (anisotropy, non-Gaussianity). This leads to the possibility of applying needlet thresholding procedures as a component separation tool. In this work, we develop algorithms based on different needlet-thresholding schemes and use them as extensions of e… ▽ More Foreground components in the Cosmic Microwave Background (CMB) are sparse in a needlet representation, due to their specific morphological features (anisotropy, non-Gaussianity). This leads to the possibility of applying needlet thresholding procedures as a component separation tool. In this work, we develop algorithms based on different needlet-thresholding schemes and use them as extensions of existing, well-known component separation techniques, namely ILC and template-fitting. We test soft- and hard-thresholding schemes, using different procedures to set the optimal threshold level. We find that thresholding can be useful as a denoising tool for internal templates in experiments with few frequency channels, in conditions of low signal-to-noise. We also compare our method with other denoising techniques, showing that thresholding achieves the best performance in terms of reconstruction accuracy and data compression while preserving the map resolution. The best results in our tests are in particular obtained when considering template-fitting in an LSPE like experiment, especially for B-mode spectra. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 23 pages, 8 figures

Journal ref: Journal of Cosmology and Astroparticle Physics, 3, 054, 29 pp. (2020)

Point Source Detection and False Discovery Rate Control on CMB Maps

Authors: Javier Carrón Duque, Alessandro Buzzelli, Yabebal Fantaye, Domenico Marinucci, Armin Schwartzman, Nicola Vittorio

Abstract: We discuss a new procedure to search for point sources in Cosmic Microwave background maps; in particular, we aim at controlling the so-called False Discovery Rate, which is defined as the expected value of false discoveries among pixels which are labelled as contaminated by point sources. We exploit a procedure called STEM, which is based on the following four steps: 1) needlet filtering of the o… ▽ More We discuss a new procedure to search for point sources in Cosmic Microwave background maps; in particular, we aim at controlling the so-called False Discovery Rate, which is defined as the expected value of false discoveries among pixels which are labelled as contaminated by point sources. We exploit a procedure called STEM, which is based on the following four steps: 1) needlet filtering of the observed CMB maps, to improve the signal to noise ratio; 2) selection of candidate peaks, i.e., the local maxima of filtered maps; 3) computation of \emph{p-}values for local maxima; 4) implementation of the multiple testing procedure, by means of the so-called Benjamini-Hochberg method. Our procedures are also implemented on the latest release of Planck CMB maps. △ Less

Submitted 18 February, 2019; originally announced February 2019.

MSC Class: Primary 62M40; Secondary 62M30; 62M15; 60G60; 42C40

Journal ref: Astronomy and Computing, 28, 100310 (2019)

Needlet estimation of cross-correlation between CMB lensing maps and LSS

Authors: Federico Bianchini, Alessandro Renzi, Domenico Marinucci

Abstract: In this paper we develop a novel needlet-based estimator to investigate the cross-correlation between cosmic microwave background (CMB) lensing maps and large-scale structure (LSS) data. We compare this estimator with its harmonic counterpart and, in particular, we analyze the bias effects of different forms of masking. In order to address this bias, we also implement a MASTER-like technique in th… ▽ More In this paper we develop a novel needlet-based estimator to investigate the cross-correlation between cosmic microwave background (CMB) lensing maps and large-scale structure (LSS) data. We compare this estimator with its harmonic counterpart and, in particular, we analyze the bias effects of different forms of masking. In order to address this bias, we also implement a MASTER-like technique in the needlet case. The resulting estimator turns out to have an extremely good signal-to-noise performance. Our analysis aims at expanding and optimizing the operating domains in CMB-LSS cross-correlation studies, similarly to CMB needlet data analysis. It is motivated especially by next generation experiments (such as Euclid) which will allow us to derive much tighter constraints on cosmological and astrophysical parameters through cross-correlation measurements between CMB and LSS. △ Less

Submitted 26 November, 2016; v1 submitted 18 July, 2016; originally announced July 2016.

Comments: 18 pages, 8 figures, accepted by JCAP, updated to match published version

Journal ref: Journal of Cosmology and Astroparticle Physics, Issue 11, article id. 050 (2016)

Planck 2015 results. XVI. Isotropy and statistics of the CMB

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, Y. Akrami, P. K. Aluri, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill , et al. (220 additional authors not shown)

Abstract: We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separat… ▽ More We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date. △ Less

Submitted 22 January, 2016; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: Paper XVI of the Planck 2015 release. This is the version accepted by A&A. An additional section discussing the sensitivity of various anomalies to sky coverage has been included

Journal ref: A&A 594, A16 (2016)

Planck 2015 results. XVII. Constraints on primordial non-Gaussianity

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill , et al. (217 additional authors not shown)

Abstract: The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result f… ▽ More The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged] △ Less

Submitted 19 July, 2016; v1 submitted 5 February, 2015; originally announced February 2015.

Comments: 68 pages, 32 figures, 31 tables. Subsection 6.3 on "Primordial curvature reconstruction" added; analysis in subsection 8.3 extended (paragraph "High frequency resonance model estimator" added); introduction to Section 10 and Appendix B added. This paper is one of a set associated with the 2015 data release from Planck. Matches version accepted for publication in Astronomy & Astrophysics

Planck 2015 results. I. Overview of products and scientific results

Authors: Planck Collaboration, R. Adam, P. A. R. Ade, N. Aghanim, Y. Akrami, M. I. R. Alves, M. Arnaud, F. Arroja, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, S. Basak, P. Battaglia, E. Battaner, R. Battye, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, B. Bertincourt , et al. (330 additional authors not shown)

Abstract: The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including… ▽ More The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity. △ Less

Submitted 9 August, 2015; v1 submitted 5 February, 2015; originally announced February 2015.

Comments: 40 pages, 24 figures

A Simple Proposal for Radial 3D Needlets

Authors: Claudio Durastanti, Yabebal T. Fantaye, Frode K. Hansen, Domenico Marinucci, Isaac Z. Pesenson

Abstract: We present here a simple construction of a wavelet system for the three-dimensional ball, which we label \emph{Radial 3D Needlets}. The construction envisages a data collection environment where an observer located at the centre of the ball is surrounded by concentric spheres with the same pixelization at different radial distances, for any given resolution. The system is then obtained by weightin… ▽ More We present here a simple construction of a wavelet system for the three-dimensional ball, which we label \emph{Radial 3D Needlets}. The construction envisages a data collection environment where an observer located at the centre of the ball is surrounded by concentric spheres with the same pixelization at different radial distances, for any given resolution. The system is then obtained by weighting the projector operator built on the corresponding set of eigenfunctions, and performing a discretization step which turns out to be computationally very convenient. The resulting wavelets can be shown to have very good localization properties in the real and harmonic domain; their implementation is computationally very convenient, and they allow for exact reconstruction as they form a tight frame systems. Our theoretical results are supported by an extensive numerical analysis. △ Less

Submitted 5 August, 2014; originally announced August 2014.

Comments: 12 pages, 4 figure panels, submitted to PRD

Journal ref: Phys. Rev. D 90, 103532 (2014)

The Needlet CMB Trispesctrum

Authors: Antonino Troja, Simona Donzelli, Davide Maino, Domenico Marinucci

Abstract: We propose a computationally feasible estimator for the needlet trispectrum, which develops earlier work on the bispectrum by Donzelli et al. (2012). Our proposal seems to enjoy a number of useful properties, in particular a) the construction exploits the localization properties of the needlet system, and hence it automatically handles masked regions; b) the procedure incorporates a quadratic corr… ▽ More We propose a computationally feasible estimator for the needlet trispectrum, which develops earlier work on the bispectrum by Donzelli et al. (2012). Our proposal seems to enjoy a number of useful properties, in particular a) the construction exploits the localization properties of the needlet system, and hence it automatically handles masked regions; b) the procedure incorporates a quadratic correction term to correct for the presence of instrumental noise and sky-cuts; c) it is possible to provide analytic results on its statistical properties, which can serve as a guidance for simulations. The needlet trispectrum we present here provides the natural building blocks for the efficient estimation of nonlinearity parameters on CMB data, and in particular for the third order constants $g_{NL}$ and $τ_{NL}$. △ Less

Submitted 2 July, 2014; originally announced July 2014.

Comments: 3 pages, 1 figure, proceedings of IAU 306 Symposium Statistical Challenges in 21st Century Cosmology

Cosmological Applications of the Gaussian Kinematic Formula

Authors: Yabebal Fantaye, Frode Hansen, Davide Maino, Domenico Marinucci

Abstract: The Gaussian Kinematic Formula (GKF, see Adler and Taylor (2007,2011)) is an extremely powerful tool allowing for explicit analytic predictions of expected values of Minkowski functionals under realistic experimental conditions for cosmological data collections. In this paper, we implement Minkowski functionals on multipoles and needlet components of CMB fields, thus allowing a better control of c… ▽ More The Gaussian Kinematic Formula (GKF, see Adler and Taylor (2007,2011)) is an extremely powerful tool allowing for explicit analytic predictions of expected values of Minkowski functionals under realistic experimental conditions for cosmological data collections. In this paper, we implement Minkowski functionals on multipoles and needlet components of CMB fields, thus allowing a better control of cosmic variance and extraction of information on both harmonic and real domains; we then exploit the GKF to provide their expected values on spherical maps, in the presence of arbitrary sky masks, and under nonGaussian circumstances. All our results are validated by numerical experiments, which show a perfect agreement between theoretical predictions and Monte Carlo simulations. △ Less

Submitted 20 June, 2014; originally announced June 2014.

Comments: 19 pages, 10 figure panels

Journal ref: Phys. Rev. D 91, 063501 (2015)

The Stochastic Properties of $\ell^1$-Regularized Spherical Gaussian Fields

Authors: Valentina Cammarota, Domenico Marinucci

Abstract: Convex regularization techniques are now widespread tools for solving inverse problems in a variety of different frameworks. In some cases, the functions to be reconstructed are naturally viewed as realizations from random processes; an important question is thus whether such regularization techniques preserve the properties of the underlying probability measures. We focus here on a case which has… ▽ More Convex regularization techniques are now widespread tools for solving inverse problems in a variety of different frameworks. In some cases, the functions to be reconstructed are naturally viewed as realizations from random processes; an important question is thus whether such regularization techniques preserve the properties of the underlying probability measures. We focus here on a case which has produced a very lively debate in the cosmological literature, namely Gaussian and isotropic spherical random fields, and we prove that Gaussianity and isotropy are not conserved in general under convex regularization over a Fourier dictionary, such as the orthonormal system of spherical harmonics. △ Less

Submitted 13 September, 2013; originally announced September 2013.

MSC Class: 60G60; 33C53; 43A90

Journal ref: Applied and Computational Harmonic Analysis, 38, 2, 262-283 (2015)

Sparse Inpainting and Isotropy

Authors: Stephen M. Feeney, Domenico Marinucci, Jason D. McEwen, Hiranya V. Peiris, Benjamin D. Wandelt, Valentina Cammarota

Abstract: Sparse inpainting techniques are gaining in popularity as a tool for cosmological data analysis, in particular for handling data which present masked regions and missing observations. We investigate here the relationship between sparse inpainting techniques using the spherical harmonic basis as a dictionary and the isotropy properties of cosmological maps, as for instance those arising from cosmic… ▽ More Sparse inpainting techniques are gaining in popularity as a tool for cosmological data analysis, in particular for handling data which present masked regions and missing observations. We investigate here the relationship between sparse inpainting techniques using the spherical harmonic basis as a dictionary and the isotropy properties of cosmological maps, as for instance those arising from cosmic microwave background (CMB) experiments. In particular, we investigate the possibility that inpainted maps may exhibit anisotropies in the behaviour of higher-order angular polyspectra. We provide analytic computations and simulations of inpainted maps for a Gaussian isotropic model of CMB data, suggesting that the resulting angular trispectrum may exhibit small but non-negligible deviations from isotropy. △ Less

Submitted 31 January, 2014; v1 submitted 2 August, 2013; originally announced August 2013.

Comments: 18 pages, 6 figures. v3: matches version published in JCAP; formatting changes and single typo correction only. Code available from http://zuserver2.star.ucl.ac.uk/~smf/code.html

Report number: NSF-KITP-13-146

Journal ref: JCAP01(2014)050

On the Limiting Behaviour of Needlets Polyspectra

Authors: Valentina Cammarota, Domenico Marinucci

Abstract: This paper provides quantitative Central Limit Theorems for nonlinear transforms of spherical random fields, in the high frequency limit. The sequences of fields that we consider are represented as smoothed averages of spherical Gaussian eigenfunctions and can be viewed as random coefficients from continuous wavelets/needlets; as such, they are of immediate interest for spherical data analysis. In… ▽ More This paper provides quantitative Central Limit Theorems for nonlinear transforms of spherical random fields, in the high frequency limit. The sequences of fields that we consider are represented as smoothed averages of spherical Gaussian eigenfunctions and can be viewed as random coefficients from continuous wavelets/needlets; as such, they are of immediate interest for spherical data analysis. In particular, we focus on so-called needlets polyspectra, which are popular tools for nonGaussianity analysis in the astrophysical community, and on the area of excursion sets. Our results are based on Stein-Maliavin approximations for nonlinear transforms of Gaussian fields, and on an explicit derivation on the high-frequency limit of their variances, which may have some independent interest. △ Less

Submitted 17 July, 2013; originally announced July 2013.

MSC Class: 60G60; 62M15; 42C15

Journal ref: Annales de l'Institut Henri Poincare -Probabilites et Statistiques, 51, 1159-1189 (2015)

Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond , et al. (215 additional authors not shown)

Abstract: The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNL^local= 2.7+/-5.8, fNL^equil= -42+… ▽ More The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNL^local= 2.7+/-5.8, fNL^equil= -42+/-75, and fNL^ortho= -25+\-39 (68% CL statistical). NG is detected in the data; using skew-C_l statistics we find a nonzero bispectrum from residual point sources, and the ISW-lensing bispectrum at a level expected in the LambdaCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C_l, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, 3-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the 4-point function in the local model tauNL < 2800 (95% CL). These constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure. △ Less

Submitted 3 December, 2013; v1 submitted 20 March, 2013; originally announced March 2013.

Comments: 59 pages; LaTeX; 26 figures. This paper is one of a set associated with the 2013 data release from Planck. Updated version

Journal ref: A&A 571, A24 (2014)

Planck 2013 results. XXIII. Isotropy and statistics of the CMB

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, R. Battye, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi, L. Bonavera , et al. (230 additional authors not shown)

Abstract: The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependenc… ▽ More The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking. △ Less

Submitted 27 January, 2014; v1 submitted 20 March, 2013; originally announced March 2013.

Comments: 49 pages, 43 figures, 33 tables. Accepted for publication in A&A

Planck 2013 results. XV. CMB power spectra and likelihood

Authors: Planck collaboration, P. A. R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoit, A. Benoit-Levy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill , et al. (235 additional authors not shown)

Abstract: We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the… ▽ More We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged) △ Less

Submitted 25 March, 2013; v1 submitted 20 March, 2013; originally announced March 2013.

Comments: This paper is one of a set associated with the 2013 release of data from the Planck mission. 62 pages, 64 figures, 12 tables

Journal ref: A&A 571, A15 (2014)

Planck 2013 results. I. Overview of products and scientific results

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. I. R. Alves, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, H. Aussel, C. Baccigalupi, A. J. Banday, R. B. Barreiro, R. Barrena, M. Bartelmann, J. G. Bartlett, N. Bartolo, S. Basak, E. Battaner, R. Battye, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli , et al. (376 additional authors not shown)

Abstract: The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and p… ▽ More The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak. △ Less

Submitted 5 June, 2014; v1 submitted 20 March, 2013; originally announced March 2013.

High-frequency asymptotics for Lipschitz-Killing curvatures of excursion sets on the sphere

Authors: Domenico Marinucci, Sreekar Vadlamani

Abstract: In this paper, we shall be concerned with geometric functionals and excursion probabilities for some nonlinear transforms evaluated on Fourier components of spherical random fields. In particular, we consider both random spherical harmonics and their smoothed averages, which can be viewed as random wavelet coefficients in the continuous case. For such fields, we consider smoothed polynomial transf… ▽ More In this paper, we shall be concerned with geometric functionals and excursion probabilities for some nonlinear transforms evaluated on Fourier components of spherical random fields. In particular, we consider both random spherical harmonics and their smoothed averages, which can be viewed as random wavelet coefficients in the continuous case. For such fields, we consider smoothed polynomial transforms; we focus on the geometry of their excursion sets, and we study their asymptotic behaviour, in the high-frequency sense. We focus on the analysis of Euler-Poincaré characteristics, which can be exploited to derive extremely accurate estimates for excursion probabilities. The present analysis is motivated by the investigation of asymmetries and anisotropies in cosmological data. The statistics we focus on are also suitable to deal with spherical random fields which can only be partially observed, the canonical example being provided by the masking effect of the Milky Way on Cosmic Microwave Background (CMB) radiation data. △ Less

Submitted 12 January, 2016; v1 submitted 11 March, 2013; originally announced March 2013.

Comments: Published at http://dx.doi.org/10.1214/15-AAP1097 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AAP-AAP1097

Journal ref: Annals of Applied Probability 2016, Vol. 26, No. 1, 462-506

On the linear term correction for needlets/wavelets non-Gaussianity estimators

Authors: Simona Donzelli, Frode K. Hansen, Michele Liguori, Domenico Marinucci, Sabino Matarrese

Abstract: We derive the linear correction term for needlet and wavelet estimators of the bispectrum and the non-linearity parameter fNL on cosmic microwave background radiation data. We show that on masked WMAP-like data with anisotropic noise, the error bars improve by 10-20% and almost reach the optimal error bars obtained with the KSW estimator (Komatsu et al 2005). In the limit of full-sky and isotropic… ▽ More We derive the linear correction term for needlet and wavelet estimators of the bispectrum and the non-linearity parameter fNL on cosmic microwave background radiation data. We show that on masked WMAP-like data with anisotropic noise, the error bars improve by 10-20% and almost reach the optimal error bars obtained with the KSW estimator (Komatsu et al 2005). In the limit of full-sky and isotropic noise, this term vanishes. We apply needlet and wavelet estimators to the WMAP 7-year data and obtain our best estimate fNL=37.5 \pm 21.8. △ Less

Submitted 7 February, 2012; originally announced February 2012.

Comments: 10 pages, submitted to ApJ

Journal ref: Astrophysical Journal 755 (2012) 19

Needlet-Whittle Estimates on the Unit Sphere

Authors: Claudio Durastanti, Xiaohong Lan, Domenico Marinucci

Abstract: We study the asymptotic behaviour of needlets-based approximate maximum likelihood estimators for the spectral parameters of Gaussian and isotropic spherical random fields. We prove consistency and asymptotic Gaussianity, in the high-frequency limit, thus generalizing earlier results by Durastanti et al. (2011) based upon standard Fourier analysis on the sphere. The asymptotic results are then ill… ▽ More We study the asymptotic behaviour of needlets-based approximate maximum likelihood estimators for the spectral parameters of Gaussian and isotropic spherical random fields. We prove consistency and asymptotic Gaussianity, in the high-frequency limit, thus generalizing earlier results by Durastanti et al. (2011) based upon standard Fourier analysis on the sphere. The asymptotic results are then illustrated by an extensive Monte Carlo study. △ Less

Submitted 1 March, 2013; v1 submitted 7 February, 2012; originally announced February 2012.

Comments: 48 pages, 2 figures

MSC Class: 62M15; 62M30; 60G60; 42C40

Journal ref: Electronic Journal of Statistics, 7, 597-646, (2013)

Detection of new point sources in WMAP 7 year data using internal templates and needlets

Authors: Sandro Scodeller, Frode K. Hansen, Domenico Marinucci

Abstract: We have developed a new needlet based method to detect point sources in cosmic microwave background (CMB) maps and have applied it to the WMAP 7 year data. We use both the individual frequency channels as well as internal templates, the difference between pairs of frequency channels, with the advantage that the CMB component is eliminated. Using the area of the sky outside the Kq85 galactic mask,… ▽ More We have developed a new needlet based method to detect point sources in cosmic microwave background (CMB) maps and have applied it to the WMAP 7 year data. We use both the individual frequency channels as well as internal templates, the difference between pairs of frequency channels, with the advantage that the CMB component is eliminated. Using the area of the sky outside the Kq85 galactic mask, we detect a total of 2102 point sources at the $5σ$ level in either the frequency maps or the internal templates. Of these, 1116 are detected either at $5σ$ directly in the frequency channels or at $5σ$ in the internal templates and $\geq3σ$ at the corresponding position in the frequency channels. Of the 1116 sources, 603 are detections which have not been reported so far in WMAP data. We have made a catalogue of these sources available with position and flux estimated in the WMAP channels where they are seen. In total, we identified 1029 of the 1116 sources with counterparts at 5GHz and 69 at other frequencies. △ Less

Submitted 10 July, 2012; v1 submitted 27 January, 2012; originally announced January 2012.

Comments: Accepted version, added integral number counts, 29 pages (+16 pages for the catalogue),10 figures

Journal ref: 2012, ApJ, 753, 27

Planck Early Results. V. The Low Frequency Instrument data processing

Authors: A. Zacchei, D. Maino, C. Baccigalupi, M. Bersanelli, A. Bonaldi, L. Bonavera, C. Burigana, R. C. Butler, F. Cuttaia, G. de Zotti, J. Dick, M. Frailis, S. Galeotta, J. González-Nuevo, K. M. Górski, A. Gregorio, E. Keihänen, R. Keskitalo, J. Knoche, H. Kurki-Suonio, C. R. Lawrence, S. Leach, J. P. Leahy, M. López-Caniego, N. Mandolesi , et al. (124 additional authors not shown)

Abstract: We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic… ▽ More We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices required to compute statistical uncertainties on LFI and Planck products are also produced. Main beams are estimated down to the approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane. △ Less

Submitted 7 December, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: A&A Volume 536, December 2011 Planck early results

Planck early results. III. First assessment of the Low Frequency Instrument in-flight performance

Authors: A. Mennella, M. Bersanelli, R. C. Butler, A. Curto, F. Cuttaia, R. J. Davis, J. Dick, M. Frailis, S. Galeotta, A. Gregorio, H. Kurki-Suonio, C. R. Lawrence, S. Leach, J. P. Leahy, S. Lowe, D. Maino, N. Mandolesi, M. Maris, E. Martínez-González, P. R. Meinhold, G. Morgante, D. Pearson, F. Perrotta, G. Polenta, T. Poutanen , et al. (136 additional authors not shown)

Abstract: The scientific performance of the Planck Low Frequency Instrument (LFI) after one year of in-orbit operation is presented. We describe the main optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the main systematic effects is presented. For each of the performance parameters, we outline the methods used t… ▽ More The scientific performance of the Planck Low Frequency Instrument (LFI) after one year of in-orbit operation is presented. We describe the main optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the main systematic effects is presented. For each of the performance parameters, we outline the methods used to obtain them from the flight data and provide a comparison with pre-launch ground assessments, which are essentially confirmed in flight. △ Less

Submitted 19 December, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Published version

Journal ref: A&A Vol 536, A3 (Dec 2011)

Adaptive Nonparametric Regression on Spin Fiber Bundles

Authors: Claudio Durastanti, Daryl Geller, Domenico Marinucci

Abstract: The construction of adaptive nonparametric procedures by means of wavelet thresholding techniques is now a classical topic in modern mathematical statistics. In this paper, we extend this framework to the analysis of nonparametric regression on sections of spin fiber bundles defined on the sphere. This can be viewed as a regression problem where the function to be estimated takes as its values alg… ▽ More The construction of adaptive nonparametric procedures by means of wavelet thresholding techniques is now a classical topic in modern mathematical statistics. In this paper, we extend this framework to the analysis of nonparametric regression on sections of spin fiber bundles defined on the sphere. This can be viewed as a regression problem where the function to be estimated takes as its values algebraic curves (for instance, ellipses) rather than scalars, as usual. The problem is motivated by many important astrophysical applications, concerning for instance the analysis of the weak gravitational lensing effect, i.e. the distortion effect of gravity on the images of distant galaxies. We propose a thresholding procedure based upon the (mixed) spin needlets construction recently advocated by Geller and Marinucci (2008,2010) and Geller et al. (2008,2009), and we investigate their rates of convergence and their adaptive properties over spin Besov balls. △ Less

Submitted 22 September, 2010; originally announced September 2010.

Comments: 40 pages

MSC Class: 62G08; 62G20; 42B35; 42C10; 42C40; 46E35 42C40; 46E35

Journal ref: Journal of Multivariate Analysis, 104, pp. 16-38, (2012)

Mixed Needlets

Authors: Daryl Geller, Domenico Marinucci

Abstract: The construction of needlet-type wavelets on sections of the spin line bundles over the sphere has been recently addressed in Geller and Marinucci (2008), and Geller et al. (2008,2009). Here we focus on an alternative proposal for needlets on this spin line bundle, in which needlet coefficients arise from the usual, rather than the spin, spherical harmonics, as in the previous constructions. We la… ▽ More The construction of needlet-type wavelets on sections of the spin line bundles over the sphere has been recently addressed in Geller and Marinucci (2008), and Geller et al. (2008,2009). Here we focus on an alternative proposal for needlets on this spin line bundle, in which needlet coefficients arise from the usual, rather than the spin, spherical harmonics, as in the previous constructions. We label this system mixed needlets and investigate in full their properties, including localization, the exact tight frame characterization, reconstruction formula, decomposition of functional spaces, and asymptotic uncorrelation in the stochastic case. We outline astrophysical applications. △ Less

Submitted 18 June, 2010; originally announced June 2010.

Comments: 26 pages

MSC Class: 42C40; 60G60; 33C55; 62M15; 83F05; 58J05

Journal ref: Journal of Mathematical Analysis and its Applications, Vol. 375, pp. 610-630 (2011)

Introducing Mexican needlets for CMB analysis: Issues for practical applications and comparison with standard needlets

Authors: S. Scodeller, O. Rudjord, F. K. Hansen, D. Marinucci, D. Geller, A. Mayeli

Abstract: Over the last few years, needlets have a emerged as a useful tool for the analysis of Cosmic Microwave Background (CMB) data. Our aim in this paper is first to introduce in the CMB literature a different form of needlets, known as Mexican needlets, first discussed in the mathematical literature by Geller and Mayeli (2009a,b). We then proceed with an extensive study of the properties of both standa… ▽ More Over the last few years, needlets have a emerged as a useful tool for the analysis of Cosmic Microwave Background (CMB) data. Our aim in this paper is first to introduce in the CMB literature a different form of needlets, known as Mexican needlets, first discussed in the mathematical literature by Geller and Mayeli (2009a,b). We then proceed with an extensive study of the properties of both standard and Mexican needlets; these properties depend on some parameters which can be tuned in order to optimize the performance for a given application. Our second aim in this paper is then to give practical advice on how to adjust these parameters in order to achieve the best properties for a given problem in CMB data analysis. In particular we investigate localization properties in real and harmonic spaces and propose a recipe on how to quantify the influence of galactic and point source masks on the needlet coefficients. We also show that for certain parameter values, the Mexican needlets provide a close approximation to the Spherical Mexican Hat Wavelets (whence their name), with some advantages concerning their numerical implementation and the derivation of their statistical properties. △ Less

Submitted 26 August, 2011; v1 submitted 30 April, 2010; originally announced April 2010.

Comments: 40 pages, 11 figures, published version, main modification: added section on more realistic galactic and point source masks

Journal ref: S. Scodeller et al. 2011 ApJ 733 121

A statistical procedure for the identification of positrons in the PAMELA experiment

Authors: O. Adriani, G. C. Barbarino, G. A. Bazilevskaya, R. Bellotti, M. Boezio, E. A. Bogomolov, L. Bonechi, M. Bongi, V. Bonvicini, S. Borisov, S. Bottai, A. Bruno, F. Cafagna, D. Campana, R. Carbone, P. Carlson, M. Casolino, G. Castellini, L. Consiglio, M. P. De Pascale, C. De Santis, N. De Simone, V. Di Felice, A. M. Galper, W. Gillard , et al. (36 additional authors not shown)

Abstract: The PAMELA satellite experiment has measured the cosmic-ray positron fraction between 1.5 GeV and 100 GeV. The need to reliably discriminate between the positron signal and proton background has required the development of an ad hoc analysis procedure. In this paper, a method for positron identification is described and its stability and capability to yield a correct background estimate is shown… ▽ More The PAMELA satellite experiment has measured the cosmic-ray positron fraction between 1.5 GeV and 100 GeV. The need to reliably discriminate between the positron signal and proton background has required the development of an ad hoc analysis procedure. In this paper, a method for positron identification is described and its stability and capability to yield a correct background estimate is shown. The analysis includes new experimental data, the application of three different fitting techniques for the background sample and an estimate of systematic uncertainties due to possible inaccuracies in the background selection. The new experimental results confirm both solar modulation effects on cosmic-rays with low rigidities and an anomalous positron abundance above 10 GeV. △ Less

Submitted 20 January, 2010; originally announced January 2010.

Comments: 27 pages, 11 figures, 3 tables

Journal ref: Astropart.Phys.34:1-11,2010

Ergodicity and Gaussianity for Spherical Random Fields

Authors: Domenico Marinucci, Giovanni Peccati

Abstract: We investigate the relationship between ergodicity and asymptotic Gaussianity of isotropic spherical random fields, in the high-resolution (or high-frequency) limit. In particular, our results suggest that under a wide variety of circumstances the two conditions are equivalent, i.e. the sample angular power spectrum may converge to the population value if and only if the underlying field is asym… ▽ More We investigate the relationship between ergodicity and asymptotic Gaussianity of isotropic spherical random fields, in the high-resolution (or high-frequency) limit. In particular, our results suggest that under a wide variety of circumstances the two conditions are equivalent, i.e. the sample angular power spectrum may converge to the population value if and only if the underlying field is asymptotically Gaussian, in the high frequency sense. These findings may shed some light on the role of Cosmic Variance in Cosmic Microwave Background (CMB) radiation data analysis. △ Less

Submitted 12 November, 2009; originally announced November 2009.

Comments: 25 pages; PACS : 02.50-r, 98.70-Vc, 98.80-k

MSC Class: 60G60; 60G15; 60K40; 62M15; 42C10

Journal ref: Journal of Mathematical Physics, n.4, 043301, (2010)

Spin Needlets Spectral Estimation

Authors: Daryl Geller, Xiaohong Lan, Domenico Marinucci

Abstract: We consider the statistical analysis of random sections of a spin fibre bundle over the sphere. These may be thought of as random fields that at each point p in $S^2$ take as a value a curve (e.g. an ellipse) living in the tangent plane at that point $T_{p}S^2$, rather than a number as in ordinary situations. The analysis of such fields is strongly motivated by applications, for instance polariz… ▽ More We consider the statistical analysis of random sections of a spin fibre bundle over the sphere. These may be thought of as random fields that at each point p in $S^2$ take as a value a curve (e.g. an ellipse) living in the tangent plane at that point $T_{p}S^2$, rather than a number as in ordinary situations. The analysis of such fields is strongly motivated by applications, for instance polarization experiments in Cosmology. To investigate such fields, spin needlets were recently introduced by Geller and Marinucci (2008) and Geller et al. (2008). We consider the use of spin needlets for spin angular power spectrum estimation, in the presence of noise and missing observations, and we provide Central Limit Theorem results, in the high frequency sense; we discuss also tests for bias and asymmetries with an asymptotic justification. △ Less

Submitted 20 July, 2009; originally announced July 2009.

MSC Class: 60G60; 62M15; 42C40; 33C55; 58J05

Journal ref: Electronic Journal of Statistics, Vol. 3, (2009), pp.1497-1530

Directional Variations of the Non-Gaussianity Parameter f_NL

Authors: Øystein Rudjord, Frode K. Hansen, Xiaohong Lan, Michele Liguori, Domenico Marinucci, Sabino Matarrese

Abstract: We investigate local variations of the primordial non-Gaussianity parameter f_NL in the WMAP data, looking for possible influence of foreground contamination in the full-sky estimate of f_NL. We first improve the needlet bispectrum estimate in (Rudjord et al. 2009) on the full-sky to f_NL= 73 +/- 31 using the KQ75 mask on the co-added V+W channel. We find no particular values of f_NL estimates clo… ▽ More We investigate local variations of the primordial non-Gaussianity parameter f_NL in the WMAP data, looking for possible influence of foreground contamination in the full-sky estimate of f_NL. We first improve the needlet bispectrum estimate in (Rudjord et al. 2009) on the full-sky to f_NL= 73 +/- 31 using the KQ75 mask on the co-added V+W channel. We find no particular values of f_NL estimates close to the galactic plane and conclude that foregrounds are unlikely to affect the estimate of f_NL in the V and W bands even for the smaller KQ85 mask. In the Q band however, we find unexpectedly high values of f_NL in local estimates close to the galactic mask, as well as significant discrepancies between Q band estimates and V/W band estimates. We therefore conclude that the Q band is too contaminated to be used for non-Gaussianity studies even with the larger KQ75 mask. We further noted that the local f_NL estimates on the V+W channel are positive on all equatorial bands from the north to the south pole. The probability for this to happen in a universe with f_NL = 0 is less than one percent. △ Less

Submitted 28 May, 2010; v1 submitted 17 June, 2009; originally announced June 2009.

Comments: 6 pages, 2 figures, included section on non-Gaussian maps, references added

Journal ref: Astrophys.J.708:1321-1325,2010

An Estimate of the Primordial Non-Gaussianity Parameter f_NL Using the Needlet Bispectrum from WMAP

Authors: Oystein Rudjord, Frode K. Hansen, Xiaohong Lan, Michele Liguori, Domenico Marinucci, Sabino Matarrese

Abstract: We use the full bispectrum of spherical needlets applied to the WMAP data of the cosmic microwave background as an estimator for the primordial non-Gaussianity parameter f_NL. We use needlet scales up to l_max=1000 and the KQ75 galactic cut and find f_NL=84 +/- 40 corrected for point source bias. We also introduce a set of consistency tests to validate our results against the possible influence… ▽ More We use the full bispectrum of spherical needlets applied to the WMAP data of the cosmic microwave background as an estimator for the primordial non-Gaussianity parameter f_NL. We use needlet scales up to l_max=1000 and the KQ75 galactic cut and find f_NL=84 +/- 40 corrected for point source bias. We also introduce a set of consistency tests to validate our results against the possible influence of foreground residuals or systematic errors. In particular, fluctuations in the value of f_NL obtained from different frequency channels, different masks and different multipoles are tested against simulated maps. All variations in f_NL estimates are found statistically consistent with simulations. △ Less

Submitted 16 June, 2009; v1 submitted 20 January, 2009; originally announced January 2009.

Comments: 9 pages, 2 figures

Journal ref: Astrophys.J.701:369-376,2009

Spin Wavelets on the Sphere

Authors: Daryl Geller, Domenico Marinucci

Abstract: In recent years, a rapidly growing literature has focussed on the construction of wavelet systems to analyze functions defined on the sphere. Our purpose in this paper is to generalize these constructions to situations where sections of line bundles, rather than ordinary scalar-valued functions, are considered. In particular, we propose {\em needlet-type spin wavelets} as an extension of the nee… ▽ More In recent years, a rapidly growing literature has focussed on the construction of wavelet systems to analyze functions defined on the sphere. Our purpose in this paper is to generalize these constructions to situations where sections of line bundles, rather than ordinary scalar-valued functions, are considered. In particular, we propose {\em needlet-type spin wavelets} as an extension of the needlet approach recently introduced by Narcowich, Petrushev and Ward, and then considered for more general manifolds by Geller and Mayeli. We discuss localization properties in the real and harmonic domains, and investigate stochastic properties for the analysis of spin random fields. Our results are strongly motivated by cosmological applications, in particular in connection to the analysis of Cosmic Microwave Background polarization data. △ Less

Submitted 15 December, 2008; v1 submitted 18 November, 2008; originally announced November 2008.

Comments: 37 pages

MSC Class: 42C40; 60G60; 33C55; 14C21; 83F05; 58J05

Journal ref: Journal of Fourier Analysis and Applications, online first (2010)

Spin Needlets for Cosmic Microwave Background Polarization Data Analysis

Authors: Daryl Geller, Frode K. Hansen, Domenico Marinucci, Gerard Kerkyacharian, Dominique Picard

Abstract: Scalar wavelets have been used extensively in the analysis of Cosmic Microwave Background (CMB) temperature maps. Spin needlets are a new form of (spin) wavelets which were introduced in the mathematical literature by Geller and Marinucci (2008) as a tool for the analysis of spin random fields. Here we adopt the spin needlet approach for the analysis of CMB polarization measurements. The outcome… ▽ More Scalar wavelets have been used extensively in the analysis of Cosmic Microwave Background (CMB) temperature maps. Spin needlets are a new form of (spin) wavelets which were introduced in the mathematical literature by Geller and Marinucci (2008) as a tool for the analysis of spin random fields. Here we adopt the spin needlet approach for the analysis of CMB polarization measurements. The outcome of experiments measuring the polarization of the CMB are maps of the Stokes Q and U parameters which are spin 2 quantities. Here we discuss how to transform these spin 2 maps into spin 2 needlet coefficients and outline briefly how these coefficients can be used in the analysis of CMB polarization data. We review the most important properties of spin needlets, such as localization in pixel and harmonic space and asymptotic uncorrelation. We discuss several statistical applications, including the relation of angular power spectra to the needlet coefficients, testing for non-Gaussianity on polarization data, and reconstruction of the E and B scalar maps. △ Less

Submitted 10 December, 2008; v1 submitted 18 November, 2008; originally announced November 2008.

Comments: Accepted for publication in Phys. Rev. D

Journal ref: Phys.Rev.D78:123533,2008

Needlet Detection of Features in WMAP CMB Sky and the Impact on Anisotropies and Hemispherical Asymmetries

Authors: Davide Pietrobon, Alexandre Amblard, Amedeo Balbi, Paolo Cabella, Asantha Cooray, Domenico Marinucci

Abstract: We apply spherical needlets to the Wilkinson Microwave Anisotropy Probe 5-year cosmic microwave background (CMB) dataset, to search for imprints of non-isotropic features in the CMB sky. We use the needlets localization properties to resolve peculiar features in the CMB sky and to study how these features contribute to the anisotropy power spectrum of the CMB. In addition to the now well-known "… ▽ More We apply spherical needlets to the Wilkinson Microwave Anisotropy Probe 5-year cosmic microwave background (CMB) dataset, to search for imprints of non-isotropic features in the CMB sky. We use the needlets localization properties to resolve peculiar features in the CMB sky and to study how these features contribute to the anisotropy power spectrum of the CMB. In addition to the now well-known "cold spot" of the CMB map in the southern hemisphere, we also find two hot spots at greater than 99% confidence level, again in the southern hemisphere and closer to the Galactic plane. While the cold spot contributes to the anisotropy power spectrum in the multipoles between l=6 to l=33, the hot spots are found to be dominating the anisotropy power in the range between l=6 and l=18. Masking both the cold and the two hot spots results in a reduction by about 15% in the amplitude of the angular power spectrum of CMB around l=10. The resulting changes to the cosmological parameters when the power spectrum is estimated masking these features (in addition to the WMAP team's KQ85 mask) are within the 1$σ$ errors published with the WMAP mask only. We also study the asymmetry between the angular power spectra evaluated on the northern and southern hemispheres. When the features detected by needlets are masked, we find that the difference in the power, measured in terms of the anisotropy variance between l=4 and l=18, is reduced by a factor 2. We make available a mask related to needlet features for more detailed studies on asymmetries in the CMB anisotropy sky. △ Less

Submitted 10 November, 2008; v1 submitted 29 August, 2008; originally announced September 2008.

Comments: 12 pages, 12 figures; masks available at http://www.fisica.uniroma2.it/~cosmo/masks

Journal ref: Phys.Rev.D78:103504,2008

Adaptive density estimation for directional data using needlets

Authors: P. Baldi, G. Kerkyacharian, D. Marinucci, D. Picard

Abstract: This paper is concerned with density estimation of directional data on the sphere. We introduce a procedure based on thresholding on a new type of spherical wavelets called {\it needlets}. We establish a minimax result and prove its optimality. We are motivated by astrophysical applications, in particular in connection with the analysis of ultra high energy cosmic rays. This paper is concerned with density estimation of directional data on the sphere. We introduce a procedure based on thresholding on a new type of spherical wavelets called {\it needlets}. We establish a minimax result and prove its optimality. We are motivated by astrophysical applications, in particular in connection with the analysis of ultra high energy cosmic rays. △ Less

Submitted 31 July, 2008; originally announced July 2008.

Comments: 30 pages, 9 figures

MSC Class: 62G07; 62G20; 65T60

Journal ref: Annals of Statistics, Vol. 37, (2009), pp. 3362-3395

Statistical challenges in the analysis of Cosmic Microwave Background radiation

Authors: Paolo Cabella, Domenico Marinucci

Abstract: An enormous amount of observations on Cosmic Microwave Background radiation has been collected in the last decade, and much more data are expected in the near future from planned or operating satellite missions. These datasets are a goldmine of information for Cosmology and Theoretical Physics; their efficient exploitation posits several intriguing challenges from the statistical point of view.… ▽ More An enormous amount of observations on Cosmic Microwave Background radiation has been collected in the last decade, and much more data are expected in the near future from planned or operating satellite missions. These datasets are a goldmine of information for Cosmology and Theoretical Physics; their efficient exploitation posits several intriguing challenges from the statistical point of view. In this paper we review a number of open problems in CMB data analysis and we present applications to observations from the WMAP mission. △ Less

Submitted 15 May, 2009; v1 submitted 11 July, 2008; originally announced July 2008.

Comments: Published in at http://dx.doi.org/10.1214/08-AOAS190 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS190

Journal ref: Annals of Applied Statistics 2009, Vol. 3, No. 1, 61-95

On The Dependence Structure of Wavelet Coefficients for Spherical Random Fields

Authors: Xiaohong Lan, Domenico Marinucci

Abstract: We consider the correlation structure of the random coefficients for a wide class of wavelet systems on the sphere (Mexican needlets) which were recently introduced in the literature by Geller and Mayeli (2007). We provide necessary and sufficient conditions for these coefficients to be asymptotic uncorrelated in the real and in the frequency domain. Here, the asymptotic theory is developed in t… ▽ More We consider the correlation structure of the random coefficients for a wide class of wavelet systems on the sphere (Mexican needlets) which were recently introduced in the literature by Geller and Mayeli (2007). We provide necessary and sufficient conditions for these coefficients to be asymptotic uncorrelated in the real and in the frequency domain. Here, the asymptotic theory is developed in the high resolution sense. Statistical applications are also discussed, in particular with reference to the analysis of cosmological data. △ Less

Submitted 24 April, 2009; v1 submitted 27 May, 2008; originally announced May 2008.

Comments: Revised version for Stochastic Processes and their Applications

MSC Class: 60G60 (Primary); 62M40; 42C40; 42C10 (Secondary)

Journal ref: Stochastic Processes and their Applications, Vol. 119 (2009), pp. 3749-3766

The needlets bispectrum

Authors: Xiaohong Lan, Domenico Marinucci

Abstract: The purpose of this paper is to join two different threads of the recent literature on random fields on the sphere, namely the statistical analysis of higher order angular power spectra on one hand, and the construction of second-generation wavelets on the sphere on the other. To this aim, we introduce the needlets bispectrum and we derive a number of convergence results. Here, the limit theory… ▽ More The purpose of this paper is to join two different threads of the recent literature on random fields on the sphere, namely the statistical analysis of higher order angular power spectra on one hand, and the construction of second-generation wavelets on the sphere on the other. To this aim, we introduce the needlets bispectrum and we derive a number of convergence results. Here, the limit theory is developed in the high resolution sense. The leading motivation of these results is the need for statistical procedures for searching non-Gaussianity in Cosmic Microwave Background radiation. △ Less

Submitted 21 May, 2008; v1 submitted 27 February, 2008; originally announced February 2008.

Comments: Published in at http://dx.doi.org/10.1214/08-EJS197 the Electronic Journal of Statistics (http://www.i-journals.org/ejs/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-EJS-EJS_2008_197 MSC Class: 62G20 (Primary) 62M15; 60B15; 60G60 (Secondary)

Journal ref: Electronic Journal of Statistics 2008, Vol. 2, 332-367

Spherical Needlets for CMB Data Analysis

Authors: D. Marinucci, D. Pietrobon, A. Balbi, P. Baldi, P. Cabella, G. Kerkyacharian, P. Natoli, D. Picard, N. Vittorio

Abstract: We discuss Spherical Needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlets coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of t… ▽ More We discuss Spherical Needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlets coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of these properties, we believe needlets may turn out to be especially useful in the analysis of Cosmic Microwave Background (CMB) data on the incomplete sky, as well as of other cosmological observations. As a final advantage, we stress that the implementation of needlets is computationally very convenient and may rely completely on standard data analysis packages such as HEALPix. △ Less

Submitted 5 July, 2007; originally announced July 2007.

Comments: 7 pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 383, Issue 2, pp. 539-545, January 2008

Subsampling needlet coefficients on the sphere

Authors: P. Baldi, G. Kerkyacharian, D. Marinucci, D. Picard

Abstract: In a recent paper, we analyzed the properties of a new kind of spherical wavelets (called needlets) for statistical inference procedures on spherical random fields; the investigation was mainly motivated by applications to cosmological data. In the present work, we exploit the asymptotic uncorrelation of random needlet coefficients at fixed angular distances to construct subsampling statistics e… ▽ More In a recent paper, we analyzed the properties of a new kind of spherical wavelets (called needlets) for statistical inference procedures on spherical random fields; the investigation was mainly motivated by applications to cosmological data. In the present work, we exploit the asymptotic uncorrelation of random needlet coefficients at fixed angular distances to construct subsampling statistics evaluated on Voronoi cells on the sphere. We illustrate how such statistics can be used for isotropy tests and for bootstrap estimation of nuisance parameters, even when a single realization of the spherical random field is observed. The asymptotic theory is developed in detail in the high resolution sense. △ Less

Submitted 12 June, 2009; v1 submitted 28 June, 2007; originally announced June 2007.

Comments: Published in at http://dx.doi.org/10.3150/08-BEJ164 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

Report number: IMS-BEJ-BEJ164

Journal ref: Bernoulli 2009, Vol. 15, No. 2, 438-463

Dark Energy Constraints from Needlets Analysis of Wmap3 and NVSS Data

Authors: Davide Pietrobon, Amedeo Balbi, Domenico Marinucci

Abstract: We cross-correlate the new 3 year Wilkinson Microwave Anistropy Probe (WMAP3) cosmic microwave background (CMB) data with the NRAO VLA Sky Survey (NVSS) radio galaxy data, and find further evidence of late integrated Sachs-Wolfe (ISW) effect taking place at late times in cosmic history. Our detection makes use of a novel statistical method based on a new construction of spherical wavelets, calle… ▽ More We cross-correlate the new 3 year Wilkinson Microwave Anistropy Probe (WMAP3) cosmic microwave background (CMB) data with the NRAO VLA Sky Survey (NVSS) radio galaxy data, and find further evidence of late integrated Sachs-Wolfe (ISW) effect taking place at late times in cosmic history. Our detection makes use of a novel statistical method based on a new construction of spherical wavelets, called needlets. The null hypothesis (no ISW) is excluded at more than 99.7% confidence. When we compare the measured cross-correlation with the theoretical predictions of standard, flat cosmological models with a generalized dark energy component parameterized by its density, $\omde$, equation of state $w$ and speed of sound $\cs2$, we find $0.3\leq\omde\leq0.8$ at 95% c.l., independently of $\cs2$ and $w$. If dark energy is assumed to be a cosmological constant ($w=-1$), the bound on density shrinks to $0.41\leq\omde\leq 0.79$. Models without dark energy are excluded at more than $4σ$. The bounds on $w$ depend rather strongly on the assumed value of $\cs2$. △ Less

Submitted 27 November, 2006; originally announced November 2006.

Comments: 11th Marcel Grossmann Proceeding

Asymptotics for spherical needlets

Authors: P. Baldi, G. Kerkyacharian, D. Marinucci, D. Picard

Abstract: We investigate invariant random fields on the sphere using a new type of spherical wavelets, called needlets. These are compactly supported in frequency and enjoy excellent localization properties in real space, with quasi-exponentially decaying tails. We show that, for random fields on the sphere, the needlet coefficients are asymptotically uncorrelated for any fixed angular distance. This prop… ▽ More We investigate invariant random fields on the sphere using a new type of spherical wavelets, called needlets. These are compactly supported in frequency and enjoy excellent localization properties in real space, with quasi-exponentially decaying tails. We show that, for random fields on the sphere, the needlet coefficients are asymptotically uncorrelated for any fixed angular distance. This property is used to derive CLT and functional CLT convergence results for polynomial functionals of the needlet coefficients: here the asymptotic theory is considered in the high-frequency sense. Our proposals emerge from strong empirical motivations, especially in connection with the analysis of cosmological data sets. △ Less

Submitted 20 April, 2009; v1 submitted 23 June, 2006; originally announced June 2006.

Comments: Published in at http://dx.doi.org/10.1214/08-AOS601 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOS-AOS601 MSC Class: 62M40; 60F05; 60F17 (Primary) 62G20 (Secondary)

Journal ref: Annals of Statistics 2009, Vol. 37, No. 3, 1150-1171

Integrated Sachs-Wolfe effect from the cross correlation of WMAP3 year and the NRAO VLA sky survey data: New results and constraints on dark energy

Authors: Davide Pietrobon, Amedeo Balbi, Domenico Marinucci

Abstract: We cross-correlate the new 3 year Wilkinson Microwave Anistropy Probe (WMAP) cosmic microwave background (CMB) data with the NRAO VLA Sky Survey (NVSS) radio galaxy data, and find further evidence of late integrated Sachs-Wolfe (ISW) effect taking place at late times in cosmic history. Our detection makes use of a novel statistical method \cite{Baldi et al. 2006a, Baldi et al. 2006b} based on a… ▽ More We cross-correlate the new 3 year Wilkinson Microwave Anistropy Probe (WMAP) cosmic microwave background (CMB) data with the NRAO VLA Sky Survey (NVSS) radio galaxy data, and find further evidence of late integrated Sachs-Wolfe (ISW) effect taking place at late times in cosmic history. Our detection makes use of a novel statistical method \cite{Baldi et al. 2006a, Baldi et al. 2006b} based on a new construction of spherical wavelets, called needlets. The null hypothesis (no ISW) is excluded at more than 99.7% confidence. When we compare the measured cross-correlation with the theoretical predictions of standard, flat cosmological models with a generalized dark energy component parameterized by its density, $\omde$, equation of state $w$ and speed of sound $\cs2$, we find $0.3\leq\omde\leq0.8$ at 95% c.l., independently of $\cs2$ and $w$. If dark energy is assumed to be a cosmological constant ($w=-1$), the bound on density shrinks to $0.41\leq\omde\leq 0.79$. Models without dark energy are excluded at more than $4σ$. The bounds on $w$ depend rather strongly on the assumed value of $\cs2$. We find that models with more negative equation of state (such as phantom models) are a worse fit to the data in the case $\cs2=1$ than in the case $\cs2=0$. △ Less

Submitted 15 September, 2006; v1 submitted 20 June, 2006; originally announced June 2006.

Comments: 8 pages, 7 figures. Matches version published on Physical Review D

Journal ref: Phys.Rev. D74 (2006) 043524