Multiple methods for estimating the bispectrum of the cosmic microwave background with application to the MAXIMA data

We describe different methods for estimating the bispectrum of cosmic microwave background data. In particular, we construct a minimum-variance estimator for the flat-sky limit and compare results with previously studied frequentist methods. Application to the MAXIMA data set shows consistency with primordial Gaussianity. Weak quadratic non-Gaussianity is characterized by a tunable parameter   f _NL  , corresponding to non-Gaussianity at a level of  ∼10⁻⁵ f _NL  (the ratio of non-Gaussian to Gaussian terms), and we find limits of   f _NL= 1500 ± 950  for the minimum-variance estimator and   f _NL= 2700 ± 1650  for the usual frequentist estimator. These are the tightest limits on primordial non-Gaussianity, which include the full effects of the radiation transfer function.     

Multi-Conjugate Adaptive Optics with laser guide stars: performance in the infrared and visible

The assumption of the Gaussianity of primordial perturbations plays an important role in modern cosmology. The most direct test of this hypothesis consists of testing the Gaussianity of cosmic microwave background (CMB) maps. Counting the pixels with the temperatures in given ranges and thus estimating the one-point probability function of the field is the simplest of all the tests. Other usually more complex tests of Gaussianity generally use a great deal of the information already contained in the probability function. However, the most interesting outcome of such a test would be the signal of non-Gaussianity independent of the probability function. It is shown that the independent information has purely morphological character i.e. it depends on the geometry and topology of the level contours only. As an example we discuss in detail the quadratic model   v = u + α ( u ²-1)  ( u is a Gaussian field with   u¯ =0  and  〈 u ²〉=1  , α is a parameter) that may arise in slow-roll or two-field inflation models. We show that in the limit of small amplitude α the full information about the non-Gaussianity is contained in the probability function. If other tests are performed on this model they simply recycle the same information. A simple procedure allowing us to assess the sensitivity of any statistics to the morphological information is suggested. We provide an analytic estimate of the statistical limit for detecting the quadratic non-Gaussianity α _c as a function of the map size in the ideal situation when the scale of the field is resolved. This estimate is in a good agreement with the results of the Monte Carlo simulations of 256² and 1024² maps. The effect of resolution on the detection quadratic non-Gaussianity is also briefly discussed.     

Tests for primordial non-Gaussianity

We investigate the relative sensitivities of several tests for deviations from Gaussianity in the primordial distribution of density perturbations. We consider models for non-Gaussianity that mimic that which comes from inflation as well as that which comes from topological defects. The tests we consider involve the cosmic microwave background (CMB), large-scale structure, high-redshift galaxies, and the abundances and properties of clusters. We find that the CMB is superior at finding non-Gaussianity in the primordial gravitational potential (as inflation would produce), while observations of high-redshift galaxies are much better suited to find non-Gaussianity that resembles that expected from topological defects. We derive a simple expression that relates the abundance of high-redshift objects in non-Gaussian models to the primordial skewness.     

Limits on isocurvature perturbations from non-Gaussianity in WMAP temperature anisotropy

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of cosmic microwave background (CMB) temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable  σ,  S =σ²−〈σ²〉  , can give the same signal-to-noise ratio as   f _NL= 30  even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum α. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-year data. We do not find a significant signal of isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction  α < 0.070  (95 per cent CL) for a scale-invariant spectrum which is comparable to the limit obtained from the power spectrum.     

The effect of primordial non-Gaussianity on the topology of large-scale structure

We study the effect of primordial non-Gaussianity on the development of large-scale cosmic structure using high-resolution N -body simulations. In particular, we focus on the topological properties of the 'cosmic web', quantitatively characterized by the Minkowski functionals (MFs), for models with quadratic non-linearities with different values of the usual non-Gaussianity parameter f _NL. In the weakly non-linear regime (the amplitude of mass density fluctuations σ₀ < 0.1), we find that analytic formulae derived from perturbation theory agree with the numerical results within a few per cent of the amplitude of each MF when | f _NL| < 1000. In the non-linear regime, the detailed behaviour of the MFs as functions of threshold density deviates more strongly from the analytical curves, while the overall amplitude of the primordial non-Gaussian effect remains comparable to the perturbative prediction. When smaller-scale information is included, the influence of primordial non-Gaussianity becomes increasingly significant statistically due to decreasing sample variance. We find that the effect of the primordial non-Gaussianity with  | f _NL| = 50  is comparable to the sample variance of mass density fields with a volume of 0.125( h ⁻¹ Gpc)³ when they are smoothed by Gaussian filter at a scale of 5  h ⁻¹ Mpc. The detectability of this effect in actual galaxy surveys will strongly depend on residual uncertainties in cosmological parameters and galaxy biasing.     

Correlation of excursion sets for non-Gaussian cosmic microwave background temperature distributions

We present a method, based on the correlation function of excursion sets above a given threshold, to test the Gaussianity of the cosmic microwave background (CMB) temperature fluctuations in the sky. In particular, this method can be applied to discriminate between standard inflationary scenarios and those producing non-Gaussianity such as topological defects. We have obtained the normalized correlation of excursion sets, including different levels of noise, for two-point probability density functions constructed from the Gaussian, χ² _n and Laplace one-point probability density functions in two different ways. Considering subdegree angular scales, we find that this method can distinguish between different distributions even if the corresponding marginal probability density functions and/or the radiation power spectra are the same.     

Non-Gaussianity detections in the Bianchi VIIh corrected WMAP one-year data made with directional spherical wavelets

Many of the current anomalies reported in the WMAP text ( WMAP ) one-year data disappear after 'correcting' for the best-fitting embedded Bianchi type VII _h component, albeit assuming no dark energy component. We investigate the effect of this Bianchi correction on the detections of non-Gaussianity in the WMAP data that we previously made using directional spherical wavelets. We confirm that the deviations from Gaussianity in the kurtosis of spherical Mexican hat wavelet coefficients are eliminated once the data are corrected for the Bianchi component, as previously discovered by Jaffe et al. This is due to the reduction of the cold spot at Galactic coordinates  ( l , b ) = (209°, −57°)  , which Cruz et al. claimed to be the sole source of non-Gaussianity introduced in the kurtosis. Our previous detections of non-Gaussianity observed in the skewness of spherical wavelet coefficients are not reduced by the Bianchi correction. Indeed, the most significant detection of non-Gaussianity made with the spherical real Morlet wavelet at a significant level of 98.4 per cent remains (using a very conservative method to estimate the significance). Furthermore, we perform preliminary tests to determine if foregrounds or systematics are the source of this non-Gaussian signal, concluding that it is unlikely that these factors are responsible. We make our code to simulate Bianchi-induced temperature fluctuations publicly available.     

Wilkinson Microwave Anisotropy Probe 5-yr constraints on fnl with wavelets

We present a Gaussianity analysis of the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-yr cosmic microwave background (CMB) temperature anisotropy data maps. We use several third-order estimators based on the spherical Mexican hat wavelet. We impose constraints on the local non-linear coupling parameter f _nl using well-motivated non-Gaussian simulations. We analyse the WMAP maps at resolution of 6.9 arcmin for the Q , V , and W frequency bands. We use the KQ 75 mask recommended by the WMAP team which masks out 28 per cent of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 150 arcmin. With these coefficients, we compute the third order estimators which are used to perform a  χ²  analysis. The  χ²  statistic is used to test the Gaussianity of the WMAP data as well as to constrain the f _nl parameter. Our results indicate that the WMAP data are compatible with the Gaussian simulations, and the f _nl parameter is constrained to  −8 < f _nl < +111  at 95 per cent confidence level (CL) for the combined   V + W   map. This value has been corrected for the presence of undetected point sources, which add a positive contribution of  Δ f _nl= 3 ± 5  in the   V + W   map. Our results are very similar to those obtained by the WMAP team using the bispectrum.     

Geometrical estimators as a test of Gaussianity in the cosmic microwave background

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background (CMB). In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.     

Non-Gaussian cosmic microwave background temperature fluctuations from peculiar velocities of clusters

We use numerical simulations of a (480 Mpc  h ⁻¹)³ volume to show that the distribution of peak heights in maps of the temperature fluctuations from the kinematic and thermal Sunyaev–Zeldovich (SZ) effects will be highly non-Gaussian, and very different from the peak-height distribution of a Gaussian random field. We then show that it is a good approximation to assume that each peak in either SZ effect is associated with one and only one dark matter halo. This allows us to use our knowledge of the properties of haloes to estimate the peak-height distributions. At fixed optical depth, the distribution of peak heights resulting from the kinematic effect is Gaussian, with a width that is approximately proportional to the optical depth; the non-Gaussianity comes from summing over a range of optical depths. The optical depth is an increasing function of halo mass and the distribution of halo speeds is Gaussian, with a dispersion that is approximately independent of halo mass. This means that observations of the kinematic effect can be used to put constraints on how the abundance of massive clusters evolves, and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, on the other hand, comes primarily from the fact that, on average, the effect is larger in more massive haloes, and the distribution of halo masses is highly non-Gaussian. We also show that because haloes of the same mass may have a range of density and velocity dispersion profiles, the relation between halo mass and the amplitude of the thermal effect is not deterministic, but has some scatter.     

Cross-correlation between WMAP and 2MASS: non-Gaussianity induced by the SZ effect

We study the non-Gaussianity induced by the Sunyaev–Zel'dovich (SZ) effect in cosmic microwave background (CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy clusters, the CMB maps should have similar non-Gaussian features to the galaxy and cluster fields. Using the WMAP data and 2MASS galaxy catalogue, we show that the non-Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can be seen with the fourth-order cross-correlation between the wavelet variables of the WMAP maps and 2MASS clusters. The intensity of the fourth-order non-Gaussian features is found to be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that this non-Gaussianity can not be seen by the high-order autocorrelation of the WMAP . This is because the SZ signals in the autocorrelations of the WMAP data generally are weaker than the WMAP –2MASS cross-correlations by a factor f ², which is the ratio between the powers of the SZ-effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of high-order autocorrelations of CMB maps to cross-correlations of the CMB maps and galaxy field would be effective to constrain the powers of the SZ effect on various scales.     

Wavelet analysis and the detection of non-Gaussianity in the cosmic microwave background

We investigate the use of wavelet transforms in detecting and characterizing non-Gaussian structure in maps of the cosmic microwave background (CMB). We apply the method to simulated maps of the KaiserStebbins effect resulting from cosmic strings, on to which Gaussian signals of varying amplitudes are superposed. We find that the method significantly outperforms standard techniques based on measuring the moments of the pixel temperature distribution. We also compare the results with those obtained using techniques based on Minkowski functionals, and we again find the wavelet method to be superior. In particular, using the wavelet technique, we find that it is possible to detect non-Gaussianity even in the presence of a superposed Gaussian signal with 3 times the rms amplitude of the original cosmic string map. We also find that the wavelet technique is useful in characterizing the angular scales at which the non-Gaussian signal occurs.     

Cosmic microwave background bispectrum and slow-roll inflation

Recent tentative findings of non-Gaussian structure in the COBE -DMR data set have triggered renewed attention on candidate models from which such intrinsic signature could arise. In the framework of slow-roll inflation with built-in non-linearities in the inflaton field evolution, we present expressions for both the cosmic microwave background (CMB) skewness and the full angular bispectrum _ℓ1ℓ₂ℓ₃ in terms of the slow-roll parameters. We use an estimator for the angular bispectrum recently proposed in the literature and calculate its variance for an arbitrary ℓ _i multipole combination. We stress that a real detection of non-Gaussianity in the CMB would imply that an important component of the anisotropies arises from processes other than primordial quantum fluctuations. We further investigate the behaviour of the signal-to-(theoretical) noise ratio and demonstrate for generic inflationary models that it decreases in the limited range of small ℓs considered for increasing multipole ℓ, while the opposite applies for the standard _ℓs.     

The effects of primordial non-Gaussianity on the cosmological reionization

We investigate the effects of non-Gaussianity in the primordial density field on the reionization history. We rely on a semi-analytic method to describe the processes acting on the intergalactic medium (IGM), relating the distribution of the ionizing sources to that of dark matter haloes. Extending previous work in the literature, we consider models in which the primordial non-Gaussianity is measured by the dimensionless non-linearity parameter f _NL, using the constraints recently obtained from cosmic microwave background data. We predict the ionized fraction and the optical depth at different cosmological epochs assuming two different kinds of non-Gaussianity characterized by a scale-independent and a scale-dependent f _NL and comparing the results to those for the standard Gaussian scenario. We find that a positive f _NL enhances the formation of high-mass haloes at early epochs when reionization begins, and, as a consequence, the IGM ionized fraction can grow by a factor of up to 5 with respect to the corresponding Gaussian model. The increase of the filling factor has a small impact on the reionization optical depth and is of the order of ∼10 per cent if a scale-dependent non-Gaussianity is assumed. Our predictions for non-Gaussian models are in agreement with the latest Wilkinson Microwave Anisotropy Probe results within the error bars, but a higher precision is required to constrain the scale dependence of non-Gaussianity.     

Constraints on primordial non-Gaussianity from a needlet analysis of the WMAP-5 data

We look for a non-Gaussian signal in the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-year temperature anisotropy maps by performing a needlet-based data analysis. We use the foreground-reduced maps obtained by the WMAP team through the optimal combination of the W, V and Q channels, and perform realistic non-Gaussian simulations in order to constrain the non-linear coupling parameter f _NL. We apply a third-order estimator of the needlet coefficients skewness and compute the  χ²  statistics of its distribution. We obtain  −80 < f _NL < 120  at 95 per cent confidence level, which is consistent with a Gaussian distribution and comparable to previous constraints on the non-linear coupling. We then develop an estimator of f _NL based on the same simulations and we find consistent constraints on primordial non-Gaussianity. We finally compute the three-point correlation function in needlet space: the constraints on f _NL improve to  −50 < f _NL < 110  at 95 per cent confidence level.     

Observation number correlation in WMAP data

A remarkable similarity between the large-scale non-Gaussian pattern of cosmic microwave background (CMB) temperatures obtained by the Wilkinson Microwave Anisotropy Probe ( WMAP ) mission and the distribution features of observation numbers is noted. Motivated by such a similarity, in this work we check the WMAP data for the correlation between pixel temperature t and observation number N . Systematic effects of imbalance in the differential observations and significant t – N correlations in magnitude, distribution of non-Gaussianity and north–south asymmetry are found. Our results indicate that, for precision cosmology studies based on WMAP observations, the observation effect on released WMAP temperature maps requires further careful study.     

Scalar statistics on the sphere: application to the cosmic microwave background

A method to compute several scalar quantities of cosmic microwave background (CMB) maps on the sphere is presented. We consider here four type of scalars: the Hessian matrix scalars, the distortion scalars, the gradient-related scalars and the curvature scalars. Such quantities are obtained directly from the spherical harmonic coefficients   a _{ℓ m}   of the map. We also study the probability density function of these quantities for the case of a homogeneous and isotropic Gaussian field, which are functions of the power spectrum of the initial field. From these scalars it is possible to construct a new set of scalars which are independent of the power spectrum of the field. We test our results using simulations and find good agreement between the theoretical probability density functions and those obtained from simulations. Therefore, these quantities are proposed to investigate the presence of non-Gaussian features in CMB maps. Finally, we show how to compute the scalars in the presence of anisotropic noise and realistic masks.     

Primordial non-Gaussianity: local curvature method and statistical significance of constraints on fNL from WMAP data

We test the consistency of estimates of the non-linear coupling constant f _NL using non-Gaussian cosmic microwave background (CMB) maps generated by the method described in the work of Liguori, Matarrese & Moscardini. This procedure to obtain non-Gaussian maps differs significantly from the method used in previous works on the estimation of f _NL. Nevertheless, using spherical wavelets, we find results in very good agreement with Mukherjee & Wang, 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⁺²⁷⁰₋₂₆₀  at the 2σ level (Bayesian) and   f _NL=−10⁺³¹⁰₋₂₇₀  (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 the work of Cabella et al., the estimates may be combined to reduce the error bars. In this way, we obtain   f _NL=−5 ± 85  and   f _NL=−5 ± 175  at the 1σ and 2σ level respectively using the frequentist approach.