The direct evaluation of manifestly optimal, cut-sky CMB power spectrum and bispectrum estimators is numerically very costly, due to the presence of inverse-covariance filtering operations. This justifies the investigation of alternative approaches. In this work, we mostly focus on an inpainting algorithm that was introduced in recent CMB analyses to cure cut-sky suboptimalities of bispectrum estimators. First, we show that inpainting can equally be applied to the problem of unbiased estimation of power spectra. We then compare the performance of a novel inpainted CMB temperature power spectrum estimator to the popular apodised pseudo-$C_l$ (PCL) method and demonstrate, both numerically and with analytic arguments, that inpainted power spectrum estimates significantly outperform PCL estimates. Finally, we study the case of cut-sky bispectrum estimators, comparing the performance of three different approaches: inpainting, apodisation and a novel low-l leaning scheme. Providing an analytic argument why the local shape is typically most affected we mainly focus on local type non-Gaussianity. Our results show that inpainting allows to achieve optimality also for bispectrum estimation, but interestingly also demonstrate that appropriate apodisation, in conjunction with low-l cleaning, can lead to comparable accuracy.

We undertake a thorough search for signatures of sharp oscillatory features in the WMAP9 power spectrum and bispectrum as well as in the Planck power spectrum. For the first time, we carry out searches in both the power spectrum and bispectrum simultaneously, employing well-defined look-elsewhere statistics to assess significances in a rigorous manner. Developing efficient methods to scan power spectrum likelihoods for oscillatory features, we present results for the phenomenological bare sine and cosine modulations, allowing validation against existing Planck Likelihood surveys, as well as templates that include the correct sharp feature scaling. In particular, we study degeneracies between feature and cosmological parameters. For frequencies beyond the scale set by the acoustic peaks, the dependencies are realised through uninteresting adjustments of the comoving distance to last scattering. Hence, it is sufficient to keep cosmological parameters fixed and employ fast Gaussian approximations to the likelihood as a function of the feature model amplitude. In cases where results can be compared to the literature, our method shows excellent agreement. We supplement results from the Planck Likelihood with an analysis based on the Planck SMICA component separation map that, working on the assumption that the component separation algorithm is reliable, allows for the inclusion of a larger sky fraction. In principle, this allows us to place the most stringent constraints to date on the amplitudes of feature models in the temperature power spectrum. Invoking the WMAP bispectrum, we perform a combined power spectrum and bispectrum survey. We use and slightly generalise statistics developed in previous work to reliably judge the significance of large feature model amplitude estimates. We conclude that our results are entirely consistent with a featureless realisation of a Gaussian CMB.

The simplest inflationary models present us with few observable parameters to discriminate between them. A detection of features in the spectra of primordial density perturbations could provide valuable insights and lead to stringent tests of models of the early Universe. So far, searches for oscillatory features have not produced statistically significant results. In this work we consider a combined search for features in the power spectrum and bispectrum. We show that possible dependencies between the estimates of feature model amplitudes based on the two- and three-point correlators are largely statistically independent under the assumption of the null hypothesis of a nearly Gaussian featureless CMB. Building on this conclusion we propose an optimal amplitude estimator for a combined search and study the look-elsewhere effect in feature model surveys. In particular we construct analytic models for the distribution of amplitude estimates that allow for a reliable assessment of the significance of potential findings. We also propose a well-behaved integrated statistic that is designed to detect evidence for models exhibiting features at multiple frequencies.

Strong foreground contamination in high resolution CMB data requires masking which introduces statistical anisotropies and renders a full maximum likelihood analysis numerically intractable. Standard analysis methods like the pseudo-C_l framework lead to information loss due to estimator suboptimalities. We set out and validate a methodology for numerically efficient estimators for a masked sky that recover nearly as much information as a full maximum likelihood procedure. In addition to the standard pseudo-C_l statistic, the approach introduces an augmented basis designed to account for the mode coupling due to the masking of the sky. We motivate the choice of this basis by describing the basic structure of the covariance matrix. We demonstrate that the augmented estimator can achieve near-optimal results in the presence of a WMAP-realistic mask.