Characterization of the non-Gaussianity of radio and IR point sources at CMB frequencies

This study, using publicly available simulations, focuses on the characterization of the non-Gaussianity produced by radio point sources and by infrared (IR) sources in the frequency range of the cosmic microwave background from 30 to 350 GHz. We propose a simple prescription to infer the angular bispectrum from the power spectrum of point sources considering independent populations of sources, with or without clustering. We test the accuracy of our prediction using publicly available all-sky simulations of radio and IR sources and find very good agreement. We further characterize the configuration dependence and the frequency behaviour of the IR and radio bispectra. We show that the IR angular bispectrum peaks for squeezed triangles and that the clustering of IR sources enhances the bispectrum values by several orders of magnitude at scales l similar to 100. At 150 GHz the bispectrum of IR sources starts to dominate that of radio sources on large angular scales, and it dominates the whole multipole range at 350 GHz. Finally, we compute the bias on f(NL) induced by radio and IR sources. We show that the positive bias induced by radio sources is significantly reduced by masking the sources. We also show, for the first time, that the form of the IR bispectrum mimics a primordial 'local' bispectrum f NL. The IR sources produce a negative bias which becomes important for Planck-like resolution and at high frequencies (Delta f(NL) similar to -6 at 277 GHz and Delta f(NL) similar to -60-70 at 350 GHz). Most of the signal being due to the clustering of faint IR sources, the bias Delta f(NL)(IR) is not reduced by masking sources above a flux limit and may, in some cases, even be increased due to the reduction of the shot-noise term.