The global contribution of unresolved gamma-ray point sources to the extragalactic gamma-ray background has been recently measured down to gamma-ray fluxes lower than those reached with standard source detection techniques, and by employing the statistical properties of the observed gamma-ray counts. We investigate and exploit the complementarity of the information brought by the one-point statistics of photon counts (using more than 10 years of Fermi-LAT data) and by the recent measurement of the angular power spectrum of the unresolved gamma-ray background (based on 8 years of Fermi-LAT data). We determine, under the assumption that the source-count distribution of the brightest unresolved objects is dominated by blazars, their gamma-ray luminosity function and spectral energy distribution down to fluxes almost two orders of magnitude smaller than the threshold for detecting resolved sources. The different approaches provide consistent predictions for the gamma-ray luminosity function of blazars, and they show a significant complementarity.

The hypothesis that pulsar wind nebulae (PWNe) can significantly contribute to the excess of the positron ($e^+$) cosmic-ray flux has been consolidated after the observation of a $\gamma$-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe, and at GeV energies for Geminga at a much larger extension. The $\gamma$-ray halos around these PWNe are interpreted as due to electrons ($e^-$) and $e^+$ accelerated and escaped by their PWNe, and inverse Compton scattering low-energy photons of the interstellar radiation fields. The extension of these halos suggests that the diffusion around these PWNe is suppressed by two orders of magnitude with respect to the average in the Galaxy. We implement a two-zone diffusion model for the propagation of $e^+$ accelerated by the Galactic population of PWNe. We consider pulsars from the ATNF catalog and build up simulations of the PWN Galactic population. In both scenarios, we find that within a two-zone diffusion model, the total contribution from PWNe and secondary $e^+$ is at the level of AMS-02 data, for an efficiency of conversion of the pulsar spin down energy in $e^\pm$ of $\eta\sim0.1$. For the simulated PWNe, a $1\sigma$ uncertainty band is determined, which is of at least one order of magnitude from 10 GeV up to few TeV. The hint for a decreasing $e^+$ flux at TeV energies is found, even if it is strongly connected to the chosen value of the radius of the low diffusion bubble around each source.

The HAWC Collaboration has discovered a $\gamma$-ray emission extended about 2 degrees around the Geminga and Monogem pulsar wind nebulae (PWNe) at $\gamma$-ray energies $E_\gamma >5$ TeV. We analyze, for the first time, almost 10 years of $\gamma$-ray data obtained with the Fermi Large Area Telescope at $E_\gamma >$ 8 GeV in the direction of Geminga and Monogem. Since these two pulsars are close the Galactic plane we run our analysis with 10 different interstellar emission models (IEMs) to study the systematics due to the modeling of this component. We detect a $\gamma$-ray halo around Geminga with a significance in the range $7.8-11.8\sigma$ depending on the IEM considered. This measurement is compatible with $e^+$ and $e^-$ emitted by the PWN, which inverse-Compton scatter (ICS) with photon fields located within a distance of about 100 pc from the pulsar, where the diffusion coefficient is estimated to be around $1.1 \times 10^{27}$ cm$^2$/s at 100 GeV. We include in our analysis the proper motion of the Geminga pulsar which is relevant for $\gamma$ rays produced for ICS in the Fermi-LAT energy range. We find that an efficiency of about $1\%$ for the conversion of the spin-down energy of the pulsar into $e^+$ and $e^-$ is required to be consistent with $\gamma$-ray data from Fermi-LAT and HAWC. The inferred contribution of Geminga to the $e^+$ flux is at most $20\%$ at the highest energy AMS-02 data. Our results are compatible with the interpretation that the cumulative emission from Galactic pulsars explains the positron excess.

Recently, a few-degrees extended $\gamma$-ray halo in the direction of Geminga pulsar has been detected by HAWC, Milagro and Fermi-LAT. These observations can be interpreted with positrons ($e^+$) and electrons ($e^-$) accelerated by Geminga pulsar wind nebula (PWN), released in a Galactic environment with a low diffusion coefficient ($D_0$), and inverse Compton scattering (ICS) with the interstellar radiation fields. We inspect here how the morphology of the ICS $\gamma$-ray flux depends on the energy, the pulsar age and distance, and the strength and extension of the low-diffusion bubble. In particular we show that $\gamma$-ray experiments with a peak of sensitivity at TeV energies are the most promising ones to detect ICS halos. We perform a study of the sensitivity of HAWC, HESS and the future CTA experiment finding that, with pulsar efficiencies, for the conversion of spin-down luminosity into $e^{\pm}$, of the order of few %, the first two experiments should have already detected a few tens of ICS halos while the latter will increase the number of detections by a factor of 4. We then consider a sample of sources associated to PWNe and detected in the HESS Galactic plane survey and in the second HAWC catalog. We use the information available in these catalogs for the $\gamma$-ray spatial morphology and flux of these sources to inspect the value of $D_0$ around them and the $e^{\pm}$ injection spectrum. All sources are detected as extended with a $\gamma$-ray emission extended about $15-80$ pc. Assuming that most of the $e^{\pm}$ accelerated by these sources have been released in the interstellar medium, the diffusion coefficient is $2-30 \cdot 10^{26}$ cm$^2$/s at 1 TeV, i.e. two orders of magnitude smaller than the value considered to be the average in the Galaxy. These observations imply that Galactic PWNe have low-diffusion bubbles with a size of at least 80 pc.

The data on the inclusive flux of cosmic positrons and electrons ($e^++e^{-}$) have been recently collected from GeV to tens of TeV energies by several experiments with unprecedented precision. In addition, the Fermi-LAT Collaboration has provided a new energy spectrum for the upper bounds on the $e^++e^{-}$ dipole anisotropy. This observable can bring information on the emission from local Galactic sources, notably measured with high precision at radio frequencies. We develop a framework in which $e^-$ and $e^+$ measured at Earth from GeV up to tens of TeV energies have a composite origin. A dedicated analysis is deserved to Vela YZ and Cygnus Loop Supernova Remnants (SNRs), for which we consider two different models for the injection of $e^-$. We investigate the consistency of these models using the three physical observables: the radio flux from Vela YZ and Cygnus Loop at all the available frequencies, the $e^++e^-$ flux from five experiments from the GeV to tens of TeV energy, the $e^++e^-$ dipole anisotropy upper limits from 50 GeV to about 1 TeV. We find that the radio flux for these nearby SNRs strongly constraints the properties of the injection electron spectrum, partially compatible with the looser constraints derived from the $e^+ + e^-$ flux data. We also perform a multi-wavelength multi-messenger analysis by fitting simultaneously the radio flux on Vela YZ and Cygnus Loop and the $e^+ + e^-$ flux, and checking the outputs against the $e^+ + e^-$ dipole anisotropy data. Remarkably, we find a model which is compatible with all the $e^++e^-$ flux data, the radio data for Vela YZ and Cygnus Loop, and with the anisotropy upper bounds. We show the severe constraints imposed by the most recent data on the $e^+ + e^-$ dipole anisotropy.

The origin of cosmic-ray positrons detected with an energy above 10 GeV is one of the most intriguing mysteries in Astroparticle Physics. Different interpretations have been invoked to solve this puzzle such as pulsar wind nebulae (PWNe), supernovae remnants or interactions of dark matter particles. The HAWC and Milagro experiments have measured an extended emission of photons above 10 TeV in the direction of Geminga and Monogem PWNe which can be interpreted by photons emitted through inverse Compton scattering by positrons and electrons accelerated by these sources. These HAWC and Milagro detections are extremely important to estimate the PWN contribution to the positron excess. Positrons and electrons emitted by PWNe and injected in the interstellar medium can produce photons with an energy between radio and X-ray through synchrotron radiation with the Galactic magnetic field. In this white paper we show that Astrophysics Probe mission concepts, such as AMEGO (the All-sky Medium Energy Gamma-ray Observatory) and AdEPT (The Advanced Energetic Pair Telescope), are suitable to detect the synchrotron halos from Monogem and Geminga PWNe and we will report the list of the most promising PWNe that AMEGO could be able to detect.

Gamma-ray searches for new physics such as dark matter are often driven by investigating the composition of the extragalactic gamma-ray background (EGB). Classic approaches to EGB decomposition manifest in resolving individual point sources and dissecting the intensity spectrum of the remaining unresolved component. Furthermore, statistical methods have recently been proven to outperform the sensitivity of classic source detection algorithms in finding point-source populations in the unresolved flux regime. In this article, we employ the 1-point photon count statistics of eight years of Fermi-LAT data to resolve the population of extragalactic point sources and to decompose the diffuse isotropic background contribution for Galactic latitudes |b|>30 deg. We use three adjacent energy bins between 1 and 10 GeV. For the first time, we extend the analysis to incorporate a potential contribution from annihilating dark matter smoothly distributed in the Galaxy. We investigate the sensitivity reach of 1-point statistics for constraining the thermally-averaged self-annihilation cross section <\sigma v> of dark matter, using different template models for the Galactic foreground emission. Given the official Fermi-LAT interstellar emission model, we set upper bounds on the DM self-annihilation cross section <\sigma v> that are comparable with the constraints obtained by other indirect detection methods, in particular by the stacking analysis of several dwarf spheroidal galaxies.

The Fermi Large Area Telescope (LAT) Collaboration has recently released the Third Catalog of Hard Fermi-LAT Sources (3FHL), which contains 1556 sources detected above 10 GeV with seven years of Pass 8 data. We investigate the source count distribution of 3FHL sources at Galactic latitudes $|b|>20^{\circ}$, where the sources are mostly blazars. We use two complementary techniques: 1) a source-detection efficiency correction method and 2) an analysis of pixel photon count statistics with the 1-point probability distribution function (1pPDF). With the first method, using realistic Monte Carlo simulations of the $\gamma$-ray sky, we calculate the efficiency of the LAT to detect point sources. This enables us to find the intrinsic source count distribution at photon fluxes down to $7.5\times10^{-12}$ ph cm$^{-2}$s$^{-1}$. With this method we detect a flux break at $(3.5\pm0.4) \times 10^{-11} \mathrm{ph}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ with a significance of at least $5.4 \sigma$. The power-law indexes of the source count distribution above and below the break are $2.09\pm0.04$ and $1.07\pm0.27$, respectively. This result is confirmed with the 1pPDF method, which has a sensitivity reach of $\sim10^{-11}$ ph cm$^{-2}$s$^{-1}$. Integrating the derived source count distribution above the sensitivity of our analysis, we find that $(42\pm8)\%$ of the extragalactic $\gamma$-ray background originates from blazars.

The flux of positrons and electrons ($e^+ + e^-$) has been measured by the $Fermi$ Large Area Telescope (LAT) in the energy range between 7~GeV and 2~TeV. We discuss a number of interpretations of Pass 8 $Fermi$-LAT $e^+ + e^-$ spectrum, combining electron and positron emission from supernova remnants (SNRs) and pulsar wind nebulae (PWNe), or produced by the collision of cosmic rays with the interstellar medium. We find that the $Fermi$-LAT spectrum is compatible with the sum of electrons from a smooth SNR population, positrons from cataloged PWNe, and a secondary component. If we include in our analysis constraints from AMS-02 positron spectrum, we obtain a slightly worse fit to the $e^+ + e^-$ $Fermi$-LAT spectrum, depending on the propagation model. As an additional scenario, we replace the smooth SNR component within 0.7 kpc with the { individual sources} found in Green's catalog of Galactic SNRs. We find that separate consideration of far and near sources helps to reproduce the $e^+ + e^-$ $Fermi$-LAT spectrum. However, we show that the fit degrades when the radio constraints on the positron emission from Vela SNR (which is the main contributor at high energies) are taken into account. We find that a break in the power-law injection spectrum at about 100 GeV can also reproduce the measured $e^+ + e^-$ spectrum} and, among the cosmic-ray propagation models that we consider, no reasonable break of the power-law dependence of the diffusion coefficient can modify the electron flux enough to reproduce the observed shape.

High energy cosmic ray electrons and positrons probe the local properties of our Galaxy. In fact, electromagnetic energy losses limit the typical propagation scale of GeV-TeV electrons and positrons to a few kpc. In the diffusion model, nearby and dominant sources may produce an observable dipole anisotropy in the cosmic ray fluxes. We present a detailed study on the role of anisotropies from nearby sources in the interpretation of the observed GeV-TeV cosmic ray electron and positrons fluxes. We compute predictions for the anisotropies from known astrophysical sources as supernova remnants and pulsar wind nebulae of the ATNF catalog. Our results are compared with current anisotropy upper limits from the Fermi- LAT, AMS-02 and PAMELA experiments.

The cosmic electrons and positrons have been measured with unprecedented statistics up to several hundreds GeV, thus permitting to explore the role that close single sources can have in shaping the flux at different energies. In the present analysis, we consider electrons and positrons in cosmic rays to be produced by spallations of hadron fluxes with the interstellar medium, by a smooth Supernova Remnant (SNR) population, by all the ATNF catalog pulsars, and by few discrete, local SNRs. We test several source models on the $e^++e^-$ and $e^+$ AMS-02 flux data. For the configurations compatible with the data, we compute the dipole anisotropy in $e^++e^-$, $e^+$, $e^+/e^-$ from single sources. Our study includes a dedicated analysis to the Vela SNR. We show that Fermi-LAT present data on dipole anisotropy of $e^++e^-$ start to explore some of the models for the Vela SNR selected by AMS-02 flux data. We also investigate how the observed anisotropy could result from a combination of local sources. Our analysis shows that the search of anisotropy in the lepton fluxes up to TeV energies can be an interesting tool for the inspection of properties of close SNRs, complementary to the high precision flux data.