In recent years, many gamma-ray sources have been identified, yet the unresolved component hosts valuable information on the faintest emission. In order to extract it, a cross-correlation with gravitational tracers of matter in the Universe has been shown to be a promising tool. We report here the first identification of a cross-correlation signal between gamma rays and the distribution of mass in the Universe probed by weak gravitational lensing. We use the Dark Energy Survey Y1 weak lensing catalogue and the Fermi Large Area Telescope 9-year gamma-ray data, obtaining a signal-to-noise ratio of 5.3. The signal is mostly localised at small angular scales and high gamma-ray energies, with a hint of correlation at extended separation. Blazar emission is likely the origin of the small-scale effect. We investigate implications of the large-scale component in terms of astrophysical sources and particle dark matter emission.

We report the identification of a positive cross-correlation signal between the unresolved $\gamma$-ray emission, measured by the \emph{Fermi} Large Area Telescope, and four different galaxy cluster catalogues. The selected catalogues peak at low-redshift and span different frequency bands, including infrared, optical and X-rays. The signal-to-noise ratio of the detected cross-correlation amounts to 3.5 in the most significant case. We investigate and comment about its possible origin, in terms of compact $\gamma$-ray emission from AGNs inside clusters or diffuse emission from the intracluster medium. The analysis has been performed by introducing an accurate estimation of the cross-correlation power-spectrum covariance matrix, built with mock realisations of the gamma and galaxy cluster maps. Different methods to produce the mock realizations starting from the data maps have been investigated and compared, identifying suitable techniques which can be generalized to other cross-correlation studies.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

With a decade of gamma-ray data from the Fermi-LAT telescope, we can now hope to answer how well we know the local Universe at gamma-ray frequencies. On the other hand, with gamma-ray data alone it is not possible to directly access the distance of the emission and to point out the origin of unresolved sources. This obstacle can be overcome by cross-correlating the gamma-ray data with catalogs of objects with well-determined redshifts and positions. In this work, we cross-correlate Fermi-LAT skymaps with the 2MPZ catalog to study the local $z<0.2$ gamma-ray Universe, where about ten percent of the total unresolved gamma-ray background is produced. We find the signal to be dominated by AGN emissions, while star forming galaxies provide a subdominant contribution. Possible hints for a particle DM signal are discussed.