We report the measurement of the angular power spectrum of cross-correlation between the unresolved component of the Fermi-LAT gamma-ray sky-maps and the CMB lensing potential map reconstructed by the Planck satellite. The matter distribution in the
Universe determines the bending of light coming from the last scattering surface. At the same time, the matter density drives the growth history of astrophysical objects, including their capability at generating non-thermal phenomena, which in turn give rise to gamma-ray emissions. The Planck lensing map provides information on the integrated distribution of matter, while the integrated history of gamma-ray emitters is imprinted in the Fermi-LAT sky maps. We report here the first evidence of their correlation. We find that the multipole dependence of the cross-correlation measurement is in agreement with current models of the gamma-ray luminosity function for AGN and star forming galaxies, with a statistical evidence of 3.0$sigma$. Moreover, its amplitude can in general be matched only assuming that these extra-galactic emitters are also the bulk contribution of the measured isotopic gamma-ray background (IGRB) intensity. This leaves little room for a big contribution from galactic sources to the IGRB measured by Fermi-LAT, pointing toward a direct evidence of the extragalactic origin of the IGRB.
Based on realistic simulations, we propose an hybrid method to reconstruct the lensing potential power spectrum, directly on PLANCK-like CMB frequency maps. It implies using a large galactic mask and dealing with a strong inhomogeneous noise. For l <
100, we show that a full-sky inpainting method, already described in a previous work, still allows a minimal variance reconstruction, with a bias that must be accounted for by a Monte-Carlo method, but that does not couple to the deflection field. For l>100 we develop a method based on tiling the cut-sky with local 10x10 degrees overlapping tangent planes (referred to in the following as patches). It requires to solve various issues concerning their size/position, non-periodic boundaries and irregularly sampled data after the sphere-to-plane projection. We show how the leading noise term of the quadratic lensing estimator applied onto an apodized patch can still be taken directly from the data. To not loose spatial accuracy, we developed a tool that allows the fast determination of the complex Fourier series coefficients from a bi-dimensional irregularly sampled dataset, without performing an interpolation. We show that the multi-patch approach allows the lensing power spectrum reconstruction with a very small bias, thanks to avoiding the galactic mask and lowering the noise inhomogeneities, while still having almost a minimal variance. The data quality can be assessed at each stage and simple bi-dimensional spectra build, which allows the control of local systematic errors.
