Search

Clusters of galaxies are the largest gravitationally-bound systems in the Universe. Their dynamics are dominated by dark matter (DM), which makes them among the best targets for indirect DM searches. We analyze 12 years of data collected by the Fermi Large Area Telescope (Fermi-LAT) in the direction of 49 clusters of galaxies selected for their proximity to the Earth and their high X-ray flux, which makes them the most promising targets. We first create physically motivated models for the DM density around each cluster considering different assumptions for the substructure distribution. Then we perform a combined search for a $\gamma$-ray signal in the {\it Fermi}-LAT data between 500 MeV and 1 TeV. We find a signal of $\gamma$ rays potentially associated with DM that is at a statistical significance of $2.5\sigma-3.0\sigma$ when considering a slope for the subhalo mass distribution $\alpha=1.9$ and minimum mass of $M_{\rm{min}}=10^{-6}$ $M_{\odot}$. The best-fit DM mass and annihilation cross-sections for a $b\bar{b}$ annihilation channel are $m_{\chi}=40-60$ GeV and $\langle \sigma v \rangle = (2-4) \times 10^{-25}$ cm$^3$/s. When we consider $\alpha=2.0$ and $M_{\rm{min}}=10^{-9}$ $M_{\odot}$, the best-fit of the cross section reduces to $\langle \sigma v \rangle = (4-10) \times 10^{-26}$ cm$^3$/s. For both DM substructure models there is a tension between the values of $\langle \sigma v \rangle$ that we find and the upper limits obtained with the non-detection of a $\gamma$-ray flux from Milky Way dwarf spheroidal galaxies. This signal is thus more likely associated with $\gamma$ rays produced in the intracluster region by cosmic rays colliding with gas and photon fields., Comment: 27 Pages, 13 Figures. Accepted for publication in the PRD journal

The flux of $\gamma$ rays is measured with unprecedented accuracy by the $\textit{Fermi}$ Large Area Telescope from 100 MeV to almost 1 TeV. In the future, the Cherenkov Telescope Array will have the capability to measure photons up to 100 TeV. To accurately interpret this data, precise predictions of the production processes, specifically the cross section for the production of photons from the interaction of cosmic-ray protons and helium with atoms of the ISM, are necessary. In this study, we determine new analytical functions describing the Lorentz-invariant cross section for $\gamma$-ray production in hadronic collisions. We utilize the limited total cross section data for $\pi^0$ production channels and supplement this information by drawing on our previous analyses of charged pion production to infer missing details. In this context, we highlight the need for new data on $\pi^0$ production. Our predictions include the cross sections for all production channels that contribute down to the 0.5% level of the final cross section, namely $\eta$, $K^+$, $K^-$, $K^0_S$, and $K^0_L$ mesons as well as $\Lambda$, $\Sigma$, and $\Xi$ baryons. We determine the total differential cross section $d\sigma(p+p\rightarrow \gamma+X)/dE_{\gamma}$ from 10 MeV to 100 TeV with an uncertainty of 10% below 10 GeV of $\gamma$-ray energies, increasing to 20% at the TeV energies. We provide numerical tables and a script for the community to access our energy-differential cross sections, which are provided for incident proton (nuclei) energies from 0.1 to $10^7$ GeV (GeV/n)., Comment: 12 pages, 8 figures. This version includes also the fit to the LHCf data on $\pi^0$ production. It matches version published by PRD. The updated tables of the energy differential cross sections of gamma rays can be found here: https://github.com/lucaorusa/gamma_cross_section

Satellite galaxies of the Milky Way with high mass-to-light ratios and little baryon content, i.e. dwarf spheroidal galaxies (dSphs), are among the most promising targets to detect or constrain the nature of dark matter (DM) through its final annihilation products into high-energy photons. Previously, the assumption that DM emission from dSphs is point-like has been used to set strong constraints on DM candidates using data from the Fermi Large Area Telescope (LAT). However, due to their high DM densities and proximity, dSphs actually have sufficient angular extension to be detected by the Fermi-LAT. Here, we assess, for the first time, the impact of accounting for angular extension in the search for gamma-ray DM signals towards known dSphs with Fermi-LAT. We show that, depending on the dSph under consideration, limits on the DM cross section can be weakened by up to a factor of 2-2.5, while the impact on the stacked, i.e. combined, limits is at most 1.5-1.8 depending on the annihilation channel. This result is of relevance when comparing dSphs limits to other multi-messenger DM constraints and for testing the DM interpretation of anomalous "excesses"., Comment: 16 pages, 8 figures and 2 tables. Paper accepted for Publication in the Physical Review D journal

The cosmic-ray fluxes of electrons and positrons ($e^{\pm}$) are measured with high precision by the space-borne particle spectrometer AMS-02. To infer a precise interpretation of the production processes for $e^{\pm}$ in our Galaxy, it is necessary to have an accurate description of the secondary component, produced by the interaction of cosmic-ray proton and helium with the interstellar medium atoms. We determine new analytical functions of the Lorentz invariant cross section for the production of $\pi^\pm$ and $K^\pm$ by fitting data from collider experiments. We also evaluate the invariant cross sections for several other channels, involving for example hyperon decays, contributing at the few \% level on the total cross section. For all these particles, the relevant 2 and 3 body decay channels are implemented, with the polarized $\mu^\pm$ decay computed with next-to-leading order corrections. The cross section for scattering of nuclei heavier than protons is modeled by fitting data on $p+C$ collisions. The total differential cross section $d\sigma/dT_{e^\pm}(p+p\rightarrow e^\pm+X)$ is predicted from 10 MeV up to 10 TeV of $e^\pm$ energy with an uncertainty of about 5-7\% in the energies relevant for AMS-02 positron flux, thus dramatically reducing the precision of the theoretical model with respect to the state of the art. Finally, we provide a prediction for the secondary Galactic $e^\pm$ source spectrum with an uncertainty of the same level. As a service for the scientific community, we provide numerical tables and a script to calculate energy-differential cross sections., Comment: 25 pages, 19 figures. Typos corrected, matches version published by PRD. The tables of the energy differential cross sections of electrons and positrons with a script to read them can be found here: https://github.com/lucaorusa/positron_electron_cross_section

Numerous extended sources around Galactic pulsars have shown significant $\ensuremath{\gamma}$-ray emission from GeV to TeV energies, revealing hundreds of TeV energy electrons scattering off of the underlying photon fields through inverse Compton scattering (ICS). HAWC TeV gamma-ray observations of few-degree extended emission around the pulsars Geminga and Monogem, and LAT GeV emission around Geminga, suggest that systems older than 100 kyr have multi-TeV ${e}^{\ifmmode\pm\else\textpm\fi{}}$ propagating beyond the SNR-PWN system into the interstellar medium. Following the discovery of few $\ensuremath{\gamma}$-ray sources by HAWC at energies $\mathrm{E}g100\text{ }\text{ }\mathrm{TeV}$, we investigate the presence of an extended $\ensuremath{\gamma}$-ray emission in Fermi-LAT data around the three brightest sources detected by HAWC up to 100 TeV. We find an extended emission of ${\ensuremath{\theta}}_{68}=1.0{0}_{\ensuremath{-}0.07}^{+0.05}\text{ }\text{ }\mathrm{deg}$ around eHWC J1825-134 and ${\ensuremath{\theta}}_{68}=0.71\ifmmode\pm\else\textpm\fi{}0.10\text{ }\text{ }\mathrm{deg}$ eHWC $\mathrm{J}1907+063$. The analysis with ICS templates on Fermi-LAT data point to diffusion coefficient values which are significantly lower than the average Galactic one. When studied along with HAWC data, the $\ensuremath{\gamma}$-ray Fermi-LAT data provide invaluable insight into the very high-energy electron and positron parent populations.

Dwarf spheroidal galaxies (dSphs) are ideal targets for Dark Matter (DM) indirect searches due to their high DM content and their negligible expected astrophysical background. In this presentation, we report on the combined analysis of the observations of 20 dSphs by Fermi-LAT, HAWC, H.E.S.S., MAGIC, and VERITAS collaborations in the search for DM, focusing on the Weakly Interacting Massive Particles (WIMPs) scenario. The combined analysis allows us to maximize the sensitivity by combining individual data sets from all five experiments for which the energy ranges of the search overlap. New constraints, spanning a range of DM masses from 5 GeV to 100 TeV, on the velocity-weighted cross section for DM self annihilation will be presented.

International audience; Extended gamma-ray halos around middle-aged pulsars, discovered at TeV energies from HAWC observations, are a new and potentially rich source class. The phenomenon is interpreted as the inverse-Compton scattering of ambient photons by electron/positron pairs accelerated by pulsars and their nebulae and confined in their vicinity. Physically, the dynamics of this pair confinement remains poorly understood and halos offer an opportunity to probe the neighborhood of pulsar/supernova remnant systems, e.g., the magnetic field structure and conditions for energetic particle transport. As a population, due to their large sizes and long lifetimes, halos can be expected to have a non-negligible contribution to the GeV-TeV emission from the Galaxy, in the form of currently unidentified sources and/or unresolved emission on large scales. In the GeV range, the first detections of halo candidates were achieved recently and the phenomenon essentially remains to be explored. We have performed a systematic search for halos around middle-aged pulsars using 12.8 years of Fermi-LAT data. We have set up an analysis suite to detect extended emission around a selection of ATNF pulsars likely to harbor halos at a detectable level. We present a list of promising halo candidates, together with dedicated studies in which we investigate the morphology and spectrum of a few selected targets. Combined with measurements at other wavelengths, this provides information on the physical processes underlying the formation of pulsar halos and contributes to a better assessment of halos as a population.

An excess of $\ensuremath{\gamma}$ rays in the data measured by the Fermi Large Area Telescope in the direction of the Galactic Center has been reported in several publications. This excess, labeled as the Galactic Center excess (GCE), is detected analyzing the data with different interstellar emission models, point source catalogs and analysis techniques. The characteristics of the GCE, recently measured with unprecedented precision, are all compatible with dark matter particles (DM) annihilating in the main halo of our Galaxy, even if other interpretations are still not excluded. We investigate the DM candidates that fit the observed GCE spectrum and spatial morphology. We assume a simple scenario with DM annihilating into a single channel but we inspect also more complicated models with two and three channels. We perform a search for a $\ensuremath{\gamma}$-ray flux from a list of 48 Milky Way dwarf spheroidal galaxies (dSphs) using state-of-the-art estimation of the DM density in these objects. Since we do not find any significant signal from the dSphs, we put upper limits on the annihilation cross section that result to be compatible with the DM candidate that fits the GCE. However, we find that the GCE DM signal is excluded at the 95% confidence level by the AMS-02 $\overline{p}$ flux data for all purely hadronic (semihadronic) channels unless the diffusive halo size $L$ is smaller than 1.7 kpc (2.6 kpc). Such a small diffusion halo is at the $2\ensuremath{\sigma}$ significance lower limit for the results inferred from fluxes of radioactive cosmic rays and is in some tension with results from analyses performed with radio and $\ensuremath{\gamma}$-ray data. Furthermore, AMS-02 ${e}^{+}$ data rule out the GCE DM interpretation with pure or partial annihilation into ${e}^{+}{e}^{\ensuremath{-}}$. The only DM candidate that fits the GCE spectrum and fulfills all constraints obtained with the combined dSphs analysis and the AMS-02 $\overline{p}$ and ${e}^{+}$ data annihilates purely (or very dominantly) into ${\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}$, has a mass of $\ensuremath{\sim}60\text{ }\text{ }\mathrm{GeV}$ and roughly a thermal cross section.

The $\gamma$-ray emission from stars is induced by the interaction of cosmic rays with stellar atmospheres and photon fields. This emission is expected to come in two components: a stellar disk emission, where $\gamma$-rays are mainly produced in atmospheric showers generated by hadronic cosmic rays, and an extended halo emission, where the high density of soft photons in the surroundings of stars create a suitable environment for $\gamma$-ray production via inverse Compton (IC) scattering by cosmic-ray electrons. Besides the Sun, no other disk or halo from single stars has ever been detected in $\gamma$-rays. However, by assuming a cosmic-ray spectrum similar to that observed on Earth, the predicted $\gamma$-ray emission of super-luminous stars, like e.g. Betelgeuse and Rigel, could be high enough to be detected by the Fermi Large Area Telescope (LAT) after its first decade of operations. In this work, we use 12 years of Fermi-LAT observations along with IC models to study 9 super-luminous nearby stars, both individually and via stacking analysis. Our results show no significant $\gamma$-ray emission, but allow us to restrict the stellar $\gamma$-ray fluxes to be on average $, Comment: 8 pages, 4 figures

The latest AMS-02 data on cosmic ray electrons show a break in the energy spectrum around 40 GeV, with a change in the slope of about 0.1. We perform a combined fit to the newest AMS-02 positron and electron flux data above 10 GeV using a semi-analytical diffusion model where sources includes production of pairs from pulsar wind nebulae (PWNe), electrons from supernova remnants (SNRs) and both species from spallation of hadronic cosmic rays with interstellar medium atoms. We demonstrate that within our setup the change of slope in the AMS-02 electron data is well explained by the interplay between the flux contributions from SNRs and from PWNe. In fact, the relative contribution to the data of these two populations changes by a factor of about 13 from 10 to 1000 GeV. The PWN contribution has a significance of at least $4\sigma$, depending on the model used for the propagation, interstellar radiation field and energy losses. We checked the stability of this result against low-energy effects by solving numerically the transport equation. as well as adding possible breaks in the injection spectrum of SNRs. The effect of the energy losses alone, when the inverse Compton scattering is properly computed within a fully numerical treatment of the Klein-Nishina cross section, cannot explain the break in the $e^-$ flux data, as recently proposed in the literature., Comment: 13 pages, 6 figures and supplemental material. Phys.Rev.D 104 (2021) 8, 083012

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., Comment: 10 pages, 7 figures. Figures 2, 3 and 5 updated. Results unchanged

Cold dark matter (DM) models for structure formation predict that DM subhalos are present in the Galaxy. In the standard paradigm of DM as weakly interacting massive particle, subhalos are expected to shine in gamma rays and to provide a signal detectable with current instruments, notably with the Large Area Telescope (LAT) aboard the Fermi~satellite. This is the main motivation behind searches for DM signals towards dwarf spheroidal galaxies and unidentified Fermi-LAT sources. A significant angular extension detected from unassociated sources located at relatively high latitudes is considered a "smoking gun" signature for identifying DM subhalos. In the present work, we systematically explore, by means of state-of-the-art models of cold DM halos in the Galaxy, the detectability of extended subhalos with Fermi-LAT. We simulate a DM signal exploring different assumptions of subhalos distribution in the Galaxy and DM profile, and reconstruct its flux through a realistic Fermi-LAT analysis pipeline. In the most optimistic case, we show that a detection of extended DM subhalos can be made for annihilation cross sections higher than $3 \times 10^{-26}$ cm$^3$/s (for a 100 GeV DM mass), still compatible with existing gamma-ray constraints, and that, in this case, the preference for extension of the source (vs point-like hypothesis) is significant. For fainter signals, instead, halos not only do not show significant extension, but they are not even detectable significantly as point-like sources., Comment: 12 pages, 11 figures. Comments as welcome!

The center of the Milky Way is offering one of the most striking mystery in Astroparticle Physics. An excess of gamma rays (GCE) has been measured by several groups in the data collected by the Fermi Large Area Telescope (LAT) towards the Galactic center region. The spectrum and spatial morphology of the GCE have been claimed by some groups to be compatible with a signal from the Galactic halo of dark matter (DM). Instead, other analyses have demonstrated that the GCE properties, e.g., its energy spectrum, highly depend on the choice of the Galactic interstellar emission (IEM) model source catalogs and analysis techniques. In this paper we investigate the sensitivity of Fermi-LAT to detect the characteristics of the GCE. In particular we simulate the GCE as given by DM and we verify that, with a perfect knowledge of the background components, its energy spectrum, position, spatial morphology and symmetry is properly measured. We also inspect two more realist cases for which there are imperfections in the IEM model. In the first we have an un-modeled gamma-ray source, constituted by the low-latitude component of the Fermi bubbles. In the second we simulate the data with one IEM template and analyze the data with an other. We verify that a mismodeling of the IEM introduces a systematics of about 10-15% in the GCE energy spectrum between 1-10 GeV and about 5% in the value of the slope for a NFW DM density profile, which is used to fit the GCE spatial morphology. Finally, we show how the GCE would be detected in case of alternative processes such as gamma-ray emission from a bulge population of pulsars or from electrons and positrons or protons injected from the Galactic center. We demonstrate that for each of these cases there is a distinctive smoking gun signature that would help to identify the real mechanism behind the origin of the GCE., Comment: Accepted for publication in PRD. 20 pages and 20 figures

The cosmic-ray flux of antiprotons is measured with high precision by the space-borne particle spectrometers AMS-02.Its interpretation requires a correct description of the dominant production process for antiprotons in our Galaxy, namely, the interaction of cosmic-ray proton and helium with the interstellar medium. In the light of new cross section measurements by the NA61 experiment of $p + p \rightarrow \bar{p} + X$ and the first ever measurement of $p + \mathrm{He} \rightarrow \bar{p} + X$ by the LHCb experiment, we update the parametrization of proton-proton and proton-nucleon cross sections.We find that the LHCb $p$He data constrain a shape for the cross section at high energies and show for the first time how well the rescaling from the $pp$ channel applies to a helium target. By using $pp$, $p$He and $p$C data we estimate the uncertainty on the Lorentz invariant cross section for $p + \mathrm{He} \rightarrow \bar{p} + X$. We use these new cross sections to compute the source term for all the production channels, considering also nuclei heavier than He both in cosmic rays and the interstellar medium. The uncertainties on the total source term are at the level of $\pm20$% and slightly increase below antiproton energies of 5 GeV. This uncertainty is dominated by the $p+p \rightarrow \bar{p} + X$ cross section, which translates into all channels since we derive them using the $pp$ cross sections. The cross sections to calculate the source spectra from all relevant cosmic-ray isotopes are provided in the Supplemental Material. We finally quantify the necessity of new data on antiproton production cross sections, and pin down the kinematic parameter space which should be covered by future data., Comment: 16 pages, 11 figures, matches published version

The Andromeda (M31) and Triangulum (M33) galaxies are the closest Local Group galaxies to the Milky Way, being only 785 and 870 kpc away. These two galaxies provide an independent view of high-energy processes that are often obscured in our own Galaxy, including possible signals of dark matter (DM) particle interactions. The Fermi Large Area Telescope (Fermi-LAT) preliminary eight year list of sources includes both M31, which is detected as extended with a size of about 0.4$^\circ$, and M33, which is detected as a point-like source. The spatial morphology of M31 $\gamma$-ray emission could trace a population of unresolved sources and energetic particles originating in sources not related to massive star formation. Alternatively, the $\gamma$-ray emission could also be an indication of annihilation or decay of DM particles. We investigate these two possibilities using almost 10 years of data from the Fermi LAT. An interpretation that involves only a DM $\gamma$-ray emission is in tension with the current limits from other searches, such as those targeting Milky Way dwarf spheroidal galaxies. When we include a template of astrophysical emission, tuned on $\gamma$-ray data or from observations of these galaxies in other wavelengths, we do not find any significant evidence for a DM contribution and we set limits for the annihilation cross section that probe the thermal cross section for DM masses up to a few tens of GeV in the $b\bar{b}$ and $\tau^+\tau^-$ channels. For models where the DM substructures have masses above $10^{-6}$ solar masses our limits probe the DM interpretation of the Fermi LAT Galactic center excess. We provide also the lower limit for the DM decay time assuming the same spatial models of the DM distribution in M31 and M33., Comment: 28 pages, 16 pages. Submitted to PRD. Comments are welcome!

One of the predictions of the $\Lambda$CDM cosmological framework is the hierarchical formation of structure, giving rise to dark matter (DM) halos and subhalos. When the latter are massive enough they retain gas (i.e., baryons) and become visible. This is the case of the dwarf satellite galaxies in the Milky Way (MW). Below a certain mass, halos may not accumulate significant amounts of baryons and remain completely dark. However, if DM particles are Weakly Interacting Massive Particles (WIMPs), we expect them to annihilate in subhalos, producing gamma rays which can be detected with the Fermi satellite. Using the three most recent point-source Fermi Large Area Telescope (LAT) catalogs (3FGL, 2FHL and 3FHL), we search for DM subhalo candidates among the unidentified sources, i.e., sources with no firm association to a known astrophysical object. We apply several selection criteria based on the expected properties of the DM-induced emission from subhalos, which allow us to significantly reduce the list of potential candidates. Then, by characterizing the minimum detection flux of the instrument and comparing our sample to predictions from the Via Lactea II (VL-II) N-body cosmological simulation, we place conservative and robust constraints on the $\langle\sigma v\rangle-m_{DM}$ parameter space. For annihilation via the $\tau^+\tau^-$ channel, we put an upper limit of $4\times 10^{-26}~(5\times 10^{-25})~cm^3~s^{-1}$ for a mass of 10 (100) GeV. A critical improvement over previous treatments is the repopulation we made to include low-mass subhalos below the VL-II mass resolution. With more advanced subhalo candidate filtering the sensitivity reach of our method can potentially improve these constraints by a factor 3 (2) for $\tau^+\tau^-$ ($b \bar{b}$) channel., Comment: 39 pages, 23 figures. JCAP accepted. Matches the accepted version. Supplementary material can be found as ancillary files

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 efficiencies of the order of a 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., Comment: 21 pages, 11 figures. Phys.Rev.D 101 (2020) 10, 103035

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., Comment: 22 pages, 16 figures, This version matches the one published by PRD

We introduce DRAGON2, the new version of the public software package designed to study cosmic-ray propagation in the Galaxy. Our aim is to illustrate the approach followed in the writing of the code and to present its most important features. We describe the properties of the numerical scheme that has been adopted in DRAGON2 to implement the different processes related to cosmic-ray transport and we investigate its correctness by comparing our numerical results with a set of analytical solutions. Starting from these validation tests, we study in detail the performances of the code by probing the different factors that influence its accuracy and its speed under a wide range of different conditions. Lastly, we investigate the new features introduced in DRAGON2 in the treatment of diffusion, energy losses and reacceleration and their impact on the predicted fluxes, in comparison also with the results given by the previous version of the code.

In recent years we witnessed several experiments measuring a large set of observables related to Cosmic-ray physics with an unprecedented level of precision. In order to be able to fully exploit this great amount of new data we must act to refine our theoretical predictions. This can be achieved by building more realistic models of Cosmic-ray Galactic transport. The DRAGON project has been pursued in order to model Cosmic-rays propagation under realistic conditions and to allow a comparison with a wide set of experimental data. Studies brought forth with DRAGON showed how a treatment of Cosmic rays energy losses as realistic as possible is pivotal. In DRAGON2, the new version of the code, a more accurate, second order scheme for Cosmic Rays energy losses is implemented. In addition, the new version of the code allows us to investigate the impact of different models for Interstellar Radiation Field or galactic magnetic field. We present comparison between the previous energy losses approach and the new one, as well as validation test by comparing our numerical results with a set of analytical solutions. We study in particular the interplay of diffusion, reacceleration, and energy losses in a realistic case, and their impact on leptonic spectrum.

Understanding the isotopic composition of cosmic rays (CRs) observed near Earth represents a milestone towards the identification of their origin. Local fluxes contain all the known stable and long-lived isotopes, reflecting the complex history of primaries and secondaries as they traverse the interstellar medium. For that reason, a numerical code which aims at describing the CR transport in the Galaxy must unavoidably rely on accurate modelling of the production of secondary particles. In this work we provide a detailed description of the nuclear cross sections and decay network as implemented in the forthcoming release of the galactic propagation code DRAGON2. We present the secondary production models implemented in the code and we apply the different prescriptions to compute quantities of interest to interpret local CR fluxes (e.g., nuclear fragmentation timescales, secondary and tertiary source terms). In particular, we develop a nuclear secondary production model aimed at accurately computing the light secondary fluxes (namely: Li, Be, B) above 1 GeV/n. This result is achieved by fitting existing empirical or semi-empirical formalisms to a large sample of measurements in the energy range 100 MeV/n to 100 GeV/n and by considering the contribution of the most relevant decaying isotopes up to iron. Concerning secondary antiparticles (positrons and antiprotons), we describe a collection of models taken from the literature, and provide a detailed quantitative comparison., 22 pages, 12 figures

A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has been recently published by the Fermi-LAT collaboration, reporting that the emission broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). At the same time, evidence is accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare the mentioned observations with what is, perhaps, the simplest model for an MSP population, solely obtained by rescaling of the MSP luminosity function determined in the local MW disk via the respective stellar mass of the systems. It is remarkable that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to contribute only about a quarter of the detected emission, although a dominant contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSP is not far below the present upper bound. We also discuss a few refinements of this simple model. In particular, we use the correlation between globular cluster gamma-ray luminosity and stellar encounter rate to gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level $\lesssim 5\%$, but it may be of some importance in explaining the signal's morphology in the inner region of the Galaxy. We also comment on some effects which may lead to violations of the simple scaling used, on alternative models, and on future perspectives for improved diagnostics., 12 pages, 5 figures, accepted for publication in ApJ

The Fermi-LAT Collaboration has recently reported a new measurement of the inclusive cosmic-ray positron and electron spectrum in the energy range between 7 GeV and 2 TeV, obtained with almost seven years of all sky data processed with the Pass 8 event reconstruction. Here we discuss several interpretations of these results within a model where electrons and positrons are emitted by supernova remnants (SNRs) and pulsar wind nebulae (PWNe) or produced by the collision of cosmic rays with the interstellar medium. In particular, we study the emission from PWNe in relation to the limits set by the AMS-02 measurement of the positron flux. We investigate in a quantitative way the interplay between far and local SNRs by focusing on the impact of the radio constraints on the emission of electrons from the Vela SNR (which is possibly the main contributor at high energies). Furthermore, we discuss the possibility to exploit the features of the Fermi-LAT data to deduce some distinctive properties of the SNR electron spectrum at the injection, in view of the uncertainties associated to the modelling of cosmic-ray propagation in the Galaxy.

Fermipy is an open-source python framework that facilitates analysis of data collected by the Fermi Large Area Telescope (LAT). Fermipy is built on the Fermi Science Tools, the publicly available software suite provided by NASA for the LAT mission. Fermipy provides a high-level interface for analyzing LAT data in a simple and reproducible way. The current feature set includes methods for extracting spectral energy distributions and lightcurves, generating test statistic maps, finding new source candidates, and fitting source position and extension. Fermipy leverages functionality from other scientific python packages including NumPy, SciPy, Matplotlib, and Astropy and is organized as a community-developed package following an open-source development model. We review the current functionality of Fermipy and plans for future development.

We report on preliminary results from the Fermi High-Latitude Extended Sources Catalog (FHES), a comprehensive search for spatially extended gamma-ray sources at high Galactic latitudes ($|b|>5^\circ$) based on data from the Fermi Large Area Telescope (LAT). While the majority of high-latitude LAT sources are extragalactic blazars that appear point-like within the LAT angular resolution, there are several physics scenarios that predict the existence of populations of spatially extended sources. If Dark Matter consists of Weakly Interacting Massive Particles, the annihilation or decay of these particles in subhalos of the Milky Way would appear as a population of unassociated gamma-ray sources with finite angular extent. Gamma-ray emission from blazars could also be extended (so-called pair halos) due to the deflection of electron-positron pairs in the intergalactic magnetic field (IGMF). The pairs are produced in the absorption of gamma rays in the intergalactic medium and subsequently up-scatter photons of background radiation fields to gamma-ray energies. Measurement of pair halos could provide constraints on the strength and coherence length scale of the IGMF. In a dedicated search, we find 21 extended sources and 16 sources not previously characterized as extended. Limits on the flux of the extended source components are used to derive constraints on the strength of the IGMF using spectral and spatial templates derived from Monte Carlo simulations of electromagnetic cascades. This allows us to constrain the IGMF to be stronger than $3\times10^{-16}\,$G for a coherence length $\lambda \gtrsim 10\,$kpc., Comment: Proc. 35th ICRC, Busan, South Korea, PoS(ICRC2017)647

Fermi-LAT unidentified sources (unIDs) have proven to be compelling targets for performing indirect dark matter (DM) searches. In a previous work, we found that among the 1235 unIDs in Fermi-LAT catalogs (3FGL, 2FHL and 3FHL) only 44 of those are DM subhalos candidates. We now implement a spectral analysis to test whether these remaining sources are compatible or not with DM origin. This analysis is executed using almost 10 years of Pass 8 Fermi-LAT data. None of the unIDs are found to significantly prefer DM-induced emission compared to other, more conventional, astrophysical sources. In order to discriminate between pulsar and DM sources, we developed a new method which is based on the source's spectral curvature, peak energy, and its detection significance. We also look for spatial extension, which may be a hint for a DM origin according to our N-body simulation studies of the subhalo population. In addition, we used Gaia DR2 data to search for a potential stellar counterpart to our best DM subhalo candidates and, although no firm associations could be found, one of them coincides with the Sagittarius stream. Finally, previous constraints on the DM annihilation cross section are updated with the new number of remaining DM subhalo candidates among unIDs. Our limits now rule out canonical thermal WIMPs up to masses of 10 GeV for $b\bar{b}$ and 20 GeV for $\tau^+\tau^-$ annihilation channels, in this way being as sensitive and complementary to those obtained from other targets and probes., Comment: 27 pages, 10 figures, JCAP accepted. Matches the accepted version

We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions. The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios. The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well. The new version facilitates the users to include their own physical models by means of a modular C++ structure., Comment: 62 pages, 27 figures, 6 tables. Matches published JCAP version

We present the first interpretation of the new isotropic diffuse gamma-ray background (IGRB), measured by the Fermi Large Area Telescope (LAT), based on a statistical analysis. We demonstrate that the gamma-ray emission from unresolved active galactic nuclei and star forming galaxies is consistent with the Fermi-LAT IGRB data within the uncertainties both on the choice of the Galactic diffuse emission model and on the gamma-ray emission mechanism of these sources. Furthermore, adding to the extragalactic sources the contribution from a smooth Galactic halo of annihilating weakly interacting dark matter (DM) particles, we are able to set stringent limits on the DM annihilation cross section. Finally, we demonstrate that the addition of DM can significantly improve the fit to IGRB data., 2014 Fermi Symposium proceedings - eConf C14102.1

We perform a detailed quantitative analysis of the recent AMS-02 electron and positron data. We investigate the interplay between the emission from primary astrophysical sources, namely Supernova Remnants and Pulsar Wind Nebulae, and the contribution from a dark matter annihilation or decay signal. Our aim is to assess the information that can be derived on dark matter properties when both dark matter and primary astrophysical sources are assumed to jointly contribute to the leptonic observables measured by the AMS-02 experiment. We investigate both the possibility to set robust constraints on the dark matter annihilation/decay rate and the possibility to look for dark matter signals within realistic models that take into account the full complexity of the astrophysical background. Our results show that AMS-02 data enable to probe efficiently vast regions of the dark matter parameter space and, in some cases, to set constraints on the dark matter annihilation/decay rate that are comparable or even stronger than the ones derived from other indirect detection channels., Comment: 27 pages, 15 figures. In v2, several improvements in the statistical analysis have been implemented and some clarifications have been added. Conclusions unchanged. v2 matches version published in JCAP

The cumulative emission resulting from hadronic cosmic-ray interactions in star-forming galaxies (SFGs) has been proposed as the dominant contribution to the astrophysical neutrino flux at TeV to PeV energies reported by IceCube. The same particle interactions also inevitably create $\gamma$-ray emission that could be detectable as a component of the extragalactic $\gamma$-ray background (EGB), which is now measured with the Fermi-LAT in the energy range from 0.1 to 820 GeV. New studies of the blazar flux distribution at $\gamma$-ray energies above 50 GeV place an upper bound on the residual non-blazar component of the EGB. We show that these results are in strong tension with models that consider SFGs as the dominant source of the diffuse neutrino backgrounds. A characteristic spectral index for parent cosmic rays in starburst galaxies of $\Gamma_{\rm SB} \simeq 2.3$ for $dN/dE \propto E^{-\Gamma_{\rm SB}}$ is consistent with the observed scaling relation between $\gamma$-ray and IR luminosity for SFGs, the bounds from the non-blazar EGB, and the observed $\gamma$-ray spectra of individual starbursts, but underpredicts the IceCube data by approximately an order of magnitude., Comment: 10 pages, 6 figures, 1 table, published in ApJ

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., 26 pages, 14 figures. Comments are welcome

The nature of the Isotropic $\gamma$-ray Background (IGRB) measured by the Large Area Telescope (LAT) on the Fermi $\gamma$-ray space Telescope ({\it Fermi}) remains partially unexplained. Non-negligible contributions may originate from extragalactic populations of unresolved sources such as blazars, star-forming galaxies or galactic milli-second pulsars. A recent prediction of the diffuse $\gamma$-ray emission from Active Galactic Nuclei (AGN) with a large viewing angle with respect to the line-of-sight (l.o.s.) has demonstrated that this faint but numerous population is also expected to contribute significantly to the total IGRB intensity. A more exotic contribution to the IGRB invokes the pair annihilation of dark matter (DM) weakly interacting massive particles (WIMPs) into $\gamma$-rays. In this work, we evaluate the room left for galactic DM at high latitudes ($>10^\circ$) by including photons from both prompt emission and inverse Compton scattering, emphasizing the impact of the newly discovered contribution from misaligned AGN (MAGN) for such an analysis. Summing up all significant galactic and extragalactic components of the IGRB, we find that an improved understanding of the associated astrophysical uncertainties is still mandatory to put stringent bounds on thermally produced DM. On the other hand, we also demonstrate that the IGRB has the potential to be one of the most competitive {\it future} ways to test the DM WIMP hypothesis, once the present uncertainties are even slightly reduced. In fact, if MAGN contribute even at 90% of the maximal level consistent with our current understanding, thermally produced WIMPs would be severely constrained as DM candidates for masses up to several TeV., Comment: 11 pages revtex4, 9 figures. Matches published version

Theoretical predictions for the cosmic antiproton spectrum currently fall short of the corresponding experimental level of accuracy. Among the main sources of uncertainty are the antiproton production cross sections in cosmic ray inelastic interactions. We analyse existing data on antiproton production in $pp$ scattering, including for the first time the measurements performed by the NA49 Collaboration. We compute the antiproton spectrum finding that in the energy range where data are available (antiproton energies of about 4-550 GeV) different approaches lead to almost equivalent results, with an uncertainty of 10-20\%. Extrapolations outside this region lead to different estimates, with the uncertainties reaching the 50\% level around $1$ TeV, degrading the diagnostic power of the antiproton channel at those energies. We also comment on the uncertainties in the antiproton production source term coming from nuclei heavier than protons and from neutrons produced in $pp$ scatterings, and point out the need for dedicated experimental campaigns for all processes involving antiproton production in collisions of light nuclei., Comment: Two typos corrected (entry D4 in Table II and mishandled units in Eq. (10)). None of the results in the article are affected

The isotropic diffuse $\gamma$-ray background (IGRB) has been detected by various experiments and recently the Fermi-LAT Collaboration has precisely measured its spectrum in a wide energy range. The origin of the IGRB is still unclear and we show in this paper the significative improvements that have been done, thanks to the new Fermi-LAT catalogs, to solve this mystery. We demonstrate that the $\gamma$-ray intensity and spectrum of the IGRB is fully consistent with the unresolved emission from extragalactic point sources, namely Active Galactic Nuclei and Star Forming Galaxies. We show also that the IGRB can be employed to derive sever constraints for the $\gamma$-ray emission from diffuse processes such as annihilation of Dark Matter (DM) particles. Our method is able to provide low bounds for the thermal annihilation cross section for a wide range of DM masses., Comment: 6 pages, 5 figures. Proceedings of the XIV International Conference on Topics in Astroparticle and Underground Physics (TAUP 2015), 7-11 September 2015, Torino, Italy