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.

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.

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.

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]

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.

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\'{e} 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.

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.

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.

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 $\tau_{NL}$.

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.

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.

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.

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.

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.

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.

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.

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)

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.

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\sigma$ level in either the frequency maps or the internal templates. Of these, 1116 are detected either at $5\sigma$ directly in the frequency channels or at $5\sigma$ in the internal templates and $\geq3\sigma$ 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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$\sigma$ 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.

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.

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.

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.

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\sigma$. The bounds on $w$ depend rather strongly on the assumed value of $\cs2$.

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\sigma$. 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$.

In this paper we study the asymptotic behavior of the angular bispectrum of spherical random fields. Here, the asymptotic theory is developed in the framework of fixed-radius fields, which are observed with increasing resolution as the sample size grows. The results we present are then exploited in a set of procedures aimed at testing non-Gaussianity; for these statistics, we are able to show convergence to functionals of standard Brownian motion under the null hypothesis. Analytic results are also presented on the behavior of the tests in the presence of a broad class of non-Gaussian alternatives. The issue of testing for non-Gaussianity on spherical random fields has recently gained enormous empirical importance, especially in connection with the statistical analysis of cosmic microwave background radiation.

We propose a fast and efficient bispectrum statistic for Cosmic Microwave Background (CMB) temperature anisotropies to constrain the amplitude of the primordial non-Gaussian signal measured in terms of the non-linear coupling parameter f_NL. We show how the method can achieve a remarkable computational advantage by focussing on subsets of the multipole configurations, where the non-Gaussian signal is more concentrated. The detection power of the test, increases roughly linearly with the maximum multipole, as shown in the ideal case of an experiment without noise and gaps. The CPU-time scales as l_{max}^3 instead of l_{max}^5 for the full bispectrum which for Planck resolution l_{max} \sim 3000 means an improvement in speed of a factor 10^7 compared to the full bispectrum analysis with minor loss in precision. We find that the introduction of a galactic cut partially destroys the optimality of the configuration, which will then need to be dealt with in the future. We find for an ideal experiment with l_{max}=2000 that upper limits of f_{NL}<8 can be obtained at 1 sigma. For the case of the WMAP experiment, we would be able to put limits of |f_{NL}|<40 if no galactic cut were present. Using the real data with galactic cut, we obtain an estimate of -80<f_{NL}<80 and -160<f_{NL}<160 at 1 and 2 sigma respectively.

We discuss the derivation of the analytic properties of the cross-power spectrum estimator from multi-detector CMB anisotropy maps. The method is computationally convenient and it provides unbiased estimates under very broad assumptions. We also propose a new procedure for testing for the presence of residual bias due to inappropriate noise subtraction in pseudo-$C_{\ell}$ estimates. We derive the analytic behavior of this procedure under the null hypothesis, and use Monte Carlo simulations to investigate its efficiency properties, which appear very promising. For instance, for full sky maps with isotropic white noise, the test is able to identify an error of 1% on the noise amplitude estimate.

The angular bispectrum of spherical random fields has recently gained an enormous importance, especially in connection with statistical inference on cosmological data. In this paper, we provide expressions for its moments of arbitrary order and we use these results to establish a multivariate central limit theorem and higher order approximations. The results rely upon combinatorial methods from graph theory and a detailed investigation for the asymptotic behaviour of Clebsch-Gordan coefficients; the latter are widely used in representation theory and quantum theory of angular momentum.

We test the consistency of estimates of the non-linear coupling constant f_{NL} using non-Gaussian CMB maps generated by the method described in (Liguori, Matarrese and Moscardini 2003). This procedure to obtain non-Gaussian maps differs significantly from the method used in previous works on estimation of f_{NL}. Nevertheless, using spherical wavelets, we find results in very good agreement with (Mukherjee and Wang 2004), showing that the two ways of generating primordial non-Gaussian maps give equivalent results. Moreover, we introduce a new method for estimating the non-linear coupling constant from CMB observations by using the local curvature of the temperature fluctuation field. We present both Bayesian credible regions (assuming a flat prior) and proper (frequentist) confidence intervals on f_{NL}, and discuss the relation between the two approaches. The Bayesian approach tends to yield lower error bars than the frequentist approach, suggesting that a careful analysis of the different interpretations is needed. Using this method, we estimate f_{NL}=-10^{+270}_{-260} at the 2\sigma level (Bayesian) and f_{NL}=-10^{+310}_{-270} (frequentist). Moreover, we find that the wavelet and the local curvature approaches, which provide similar error bars, yield approximately uncorrelated estimates of f_{NL} and therefore, as advocated in (Cabella et al. 2004), the estimates may be combined to reduce the error bars. In this way, we obtain f_{NL}=-5\pm 85 and f_{NL}=-5\pm 175 at the 1\sigma and 2\sigma level respectively using the frequentist approach.

We use the local curvature to investigate the possible existence of non-Gaussianity/asymmetry in the WMAP data. Considering the full sky we find results which are consistent with the Gaussian assumption. However, strong non-Gaussian features emerge when considering the northern and southern galactic hemisphere separately, particularly on scales between 1 and 5 degrees. Quite interestingly, the maximum non-Gaussianity is found for hemispheres centered near the ecliptic poles, which might suggest the presence of some systematic effect. The direction of the asymmetry seems consistent with the findings by Eriksen et al. 2004.

We present a comparison between three approaches to test non-Gaussianity of cosmic microwave background data. The Minkowski functionals, the empirical process method and the skewness of wavelet coefficients are applied to maps generated from non-standard inflationary models and to Gaussian maps with point sources included. We discuss the different power of the pixel, harmonic and wavelet space methods on these simulated almost full-sky data (with Planck like noise). We also suggest a new procedure consisting of a combination of statistics in pixel, harmonic and wavelet space.

In (Hansen et al. 2002) we presented a new approach for measuring non-Gaussianity of the Cosmic Microwave Background (CMB) anisotropy pattern, based on the multivariate empirical distribution function of the spherical harmonics a_lm of a CMB map. The present paper builds upon the same ideas and proposes several improvements and extensions. More precisely, we exploit the additional information on the random phases of the a_lm to provide further tests based on the empirical distribution function. Also we take advantage of the effect of rotations in improving the power of our procedures. The suggested tests are implemented on physically motivated models of non-Gaussian fields; Monte-Carlo simulations suggest that this approach may be very promising in the analysis of non-Gaussianity generated by non-standard models of inflation. We address also some experimentally meaningful situations, such as the presence of instrumental noise and a galactic cut in the map.

We present a new, model-independent approach for measuring non-Gaussianity of the Cosmic Microwave Background (CMB) anisotropy pattern. Our approach is based on the empirical distribution function of the normalized spherical harmonic expansion coefficients a_lm of a nearly full-sky CMB map, like the ones expected from forthcoming satellite experiments. Using a set of Kolmogorov-Smirnov type tests, we check for Gaussianity and independency of the a_lm. We test the method on two non-Gaussian toy-models of the CMB, one generated in spherical harmonic space and one in pixel (real) space. We also provide some rigorous results, possibly of independent interest, on the exact distribution of the spherical harmonic coefficients normalized by an estimated angular power spectrum.

We present a new approach for statistical inference on noise properties of CMB anisotropy data. We consider a Maximum Likelihood parametric estimator to recover the full dependence structure of the noise process. We also consider a semiparametric procedure which is only sensitive to the low frequency behavior of the noise spectral density. Both approaches are statistically robust and computationally convenient in the case of long memory noise, even under nonstationary circumstances. We show that noise properties can be consistently derived by such procedures without resorting to currently used iterative noise-signal methods. More importantly, we show that optimal (GLS) CMB maps can be obtained from the observed timestream with the only knowledge of the noise memory parameter, the outcome of our estimators.