Showing 1–50 of 104 results for author: Donato, F

Geminga's pulsar halo: an X-ray view

Authors: Silvia Manconi, Jooyun Woo, Ruo-Yu Shang, Roman Krivonos, Claudia Tang, Mattia Di Mauro, Fiorenza Donato, Kaya Mori, Charles J. Hailey

Abstract: Geminga is the first pulsar around which a remarkable TeV gamma-ray halo extending over a few degrees was discovered by MILAGRO, HAWC and later by H.E.S.S., and by Fermi-LAT in the GeV band. More middle-aged pulsars have exhibited gamma-ray halos, and they are now recognized as an emerging class of Galactic gamma-ray sources. The emission appears in the late evolution stage of pulsars, and is most… ▽ More Geminga is the first pulsar around which a remarkable TeV gamma-ray halo extending over a few degrees was discovered by MILAGRO, HAWC and later by H.E.S.S., and by Fermi-LAT in the GeV band. More middle-aged pulsars have exhibited gamma-ray halos, and they are now recognized as an emerging class of Galactic gamma-ray sources. The emission appears in the late evolution stage of pulsars, and is most plausibly explained by inverse Compton scattering of CMB and interstellar photons by relativistic electrons and positrons escaping from the pulsar wind nebulae. These observations pose a number of theoretical challenges. Tackling these questions requires constraining the ambient magnetic field properties, which can be achieved through X-ray observations. If the gamma-ray halos originate from a distribution of highly energetic electrons, synchrotron losses in the ambient magnetic fields of the same particles are expected to produce a diffuse X-ray emission with a similar spatial extension. We present the most comprehensive X-ray study of the Geminga pulsar halo to date, utilising archival data from XMM-Newton and NuSTAR. Our X-ray analysis covers a broad bandwidth ($0.5\rm{-}79$ keV) and large field of view ($\sim 4^\circ$) for the first time. This is achieved by accurately measuring the background over the entire field of view, and taking into account both focused and stray-light X-ray photons with NuSTAR. We find no significant emission and set robust constraints on the X-ray halo flux. These are translated to stringent constraints on the ambient magnetic field strength and the diffusion coefficient by using a physical model considering particle injection, diffusion and cooling over the pulsar's lifetime, which is tuned by fitting multi-wavelength data. Our novel methodology for modelling and searching for synchrotron X-ray halos can be applied to other pulsar halo candidates. △ Less

Submitted 4 April, 2024; v1 submitted 16 March, 2024; originally announced March 2024.

Comments: 15 pages, 11 figures. Comments are welcomed. v2: Few comments added. Matches version submitted to A&A

Report number: LAPTH--012/24

The Galactic center excess at the highest energies: morphology and photon-count statistics

Authors: Silvia Manconi, Francesca Calore, Fiorenza Donato

Abstract: The nature of the GeV gamma-ray Galactic center excess (GCE) in the data of Fermi-Large Area Telescope (LAT) is still to be unveiled. We present a new analysis of the inner Galaxy Fermi-LAT data at energies above 10 GeV, based on an innovative method which combines the skyFACT adaptive template fitting with and the 1pPDF pixel-count statistics. We find a strong evidence for the GCE also at high en… ▽ More The nature of the GeV gamma-ray Galactic center excess (GCE) in the data of Fermi-Large Area Telescope (LAT) is still to be unveiled. We present a new analysis of the inner Galaxy Fermi-LAT data at energies above 10 GeV, based on an innovative method which combines the skyFACT adaptive template fitting with and the 1pPDF pixel-count statistics. We find a strong evidence for the GCE also at high energies, $σ> 5$ regardless of the GCE spatial template. Remarkably, our fits prefer the bulge morphological model over the dark matter one at high significance, and show no evidence for an additional dark matter template on top of the bulge component. Through the 1pPDF analysis, we find that the model best describing the gamma-ray data requires a smooth, diffuse GCE following a bulge morphology, together with sub-threshold point sources. The 1pPDF fit reconstructs a consistent population of faint point sources down at least to $10^{-12}$ ph cm$^{-2}$ s$^{-1}$. Between $10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and $10^{-11}$ ph cm$^{-2}$ s$^{-1}$ the 1pPDF measures a number of point sources significantly higher than the ones in the Fermi 4FGL catalog. The robustness of our results brings further support to the attempt of explaining, at least partially, the high-energy tail of the GCE in terms of a population of point sources, likely corresponding to millisecond pulsars. △ Less

Submitted 24 May, 2024; v1 submitted 7 February, 2024; originally announced February 2024.

Comments: 11 pages, 6 figures. Appendix adds 3 pages and 4 figures. v2: few comments added, results unchanged. Matches version accepted by PRD

Report number: LAPTH-008/24

A novel prediction for secondary positrons and electrons in the Galaxy

Authors: Mattia Di Mauro, Fiorenza Donato, Michael Korsmeier, Silvia Manconi, Luca Orusa

Abstract: The Galactic flux of cosmic-ray (CR) positrons in the GeV to TeV energy range is very likely due to different Galactic components. One of these is the inelastic scattering of CR nuclei with the atoms of the interstellar medium. The precise amount of this component determines the eventual contribution from other sources. We present here a new estimation of the secondary CR positron flux by incorpor… ▽ More The Galactic flux of cosmic-ray (CR) positrons in the GeV to TeV energy range is very likely due to different Galactic components. One of these is the inelastic scattering of CR nuclei with the atoms of the interstellar medium. The precise amount of this component determines the eventual contribution from other sources. We present here a new estimation of the secondary CR positron flux by incorporating the latest results for the production cross sections of $e^\pm$ from hadronic scatterings calibrated on collider data. All the reactions for CR nuclei up to silicon scattering on both hydrogen and helium are included. The propagation models are derived consistently by fits on primary and secondary CR nuclei data. Models with a small halo size ($L \leq 2$ kpc) are disfavored by the nuclei data although the current uncertainties on the beryllium nuclear cross sections may impact this result. The resulting positron flux shows a strong dependence on the Galactic halo size, increasing up to factor 1.5 moving $L$ from 8 to 2 kpc. Within the most reliable propagation models, the positron flux matches the data for energies below 1 GeV. We verify that secondary positrons contribute less than $70\%$ of the data above a few GeV, corroborating that an excess of positrons is already present at very low energies. At larger energies, our predictions are below the data with the discrepancy becoming more and more pronounced. Our results are provided together with uncertainties due to propagation and hadronic cross sections. The former uncertainties are below 5\% at fixed $L$, while the latter are about 7\% almost independently of the propagation scheme. In addition to the predictions of positrons, we provide new predictions also for the secondary CR electron flux. △ Less

Submitted 22 September, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: Matches published version. Few comments and clarifications added, results unchanged

Report number: LAPTH-014/23, TTK-23-07

Journal ref: Phys. Rev. D 108, 063024, 2023

New determination of the production cross section for $γ$ rays in the Galaxy

Authors: Luca Orusa, Mattia Di Mauro, Fiorenza Donato, Michael Korsmeier

Abstract: The flux of $γ$ 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 intera… ▽ More The flux of $γ$ 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 $γ$-ray production in hadronic collisions. We utilize the limited total cross section data for $π^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 $π^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 $η$, $K^+$, $K^-$, $K^0_S$, and $K^0_L$ mesons as well as $Λ$, $Σ$, and $Ξ$ baryons. We determine the total differential cross section $dσ(p+p\rightarrow γ+X)/dE_γ$ from 10 MeV to 100 TeV with an uncertainty of 10% below 10 GeV of $γ$-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). △ Less

Submitted 26 April, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 12 pages, 8 figures. This version includes also the fit to the LHCf data on $π^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

Journal ref: Phys. Rev. D (2023) 107, 083031

New determination of the production cross section for secondary positrons and electrons in the Galaxy

Authors: Luca Orusa, Mattia Di Mauro, Fiorenza Donato, Michael Korsmeier

Abstract: 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 a… ▽ More 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 $π^\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 $μ^\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σ/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. △ Less

Submitted 1 July, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 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

Journal ref: Phys.Rev.D 105 (2022) 12, 123021

Snowmass2021 Cosmic Frontier: The landscape of cosmic-ray and high-energy photon probes of particle dark matter

Authors: Tsuguo Aramaki, Mirko Boezio, James Buckley, Esra Bulbul, Philip von Doetinchem, Fiorenza Donato, J. Patrick Harding, Chris Karwin, Jason Kumar, Rebecca K. Leane, Shigeki Matsumoto, Julie McEnry, Tom Melia, Kerstin Perez, Stefano Profumo, Daniel Salazar-Gallegos, Andrew W. Strong, Brandon Roach, Miguel A. Sanchez-Conde, Tom Shutt, Atsushi Takada, Toru Tanimori, John Tomsick, Yu Watanabe, David A. Williams

Abstract: This white paper discusses the current landscape and prospects for experiments sensitive to particle dark matter processes producing photons and cosmic rays. Much of the gamma-ray sky remains unexplored on a level of sensitivity that would enable the discovery of a dark matter signal. Currently operating GeV-TeV observatories, such as Fermi-LAT, atmospheric Cherenkov telescopes, and water Cherenko… ▽ More This white paper discusses the current landscape and prospects for experiments sensitive to particle dark matter processes producing photons and cosmic rays. Much of the gamma-ray sky remains unexplored on a level of sensitivity that would enable the discovery of a dark matter signal. Currently operating GeV-TeV observatories, such as Fermi-LAT, atmospheric Cherenkov telescopes, and water Cherenkov detector arrays continue to target several promising dark matter-rich environments within and beyond the Galaxy. Soon, several new experiments will continue to explore, with increased sensitivity, especially extended targets in the sky. This paper reviews the several near-term and longer-term plans for gamma-ray observatories, from MeV energies up to hundreds of TeV. Similarly, the X-ray sky has been and continues to be monitored by decade-old observatories. Upcoming telescopes will further bolster searches and allow new discovery space for lines from, e.g., sterile neutrinos and axion-photon conversion. Furthermore, this overview discusses currently operating cosmic-ray probes and the landscape of future experiments that will clarify existing persistent anomalies in cosmic radiation and spearhead possible new discoveries. Finally, the article closes with a discussion of necessary cross section measurements that need to be conducted at colliders to reduce substantial uncertainties in interpreting photon and cosmic-ray measurements in space. △ Less

Submitted 24 June, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 61 pages, 21 figures, Snowmass2021 Cosmic Frontier White Paper

NAJADS: a self-contained framework for the direct determination of astrophysical J-factors

Authors: Anna Balaudo, Francesca Calore, Valentina De Romeri, Fiorenza Donato

Abstract: Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy-cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J-factor… ▽ More Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy-cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J-factor. In the DM phenomenology community, J-factors are typically computed through the semi-analytical modelling of the DM mass distribution, which is affected by large uncertainties. With the scope of addressing and possibly reducing these uncertainties, we present NAJADS, a self-contained framework to derive the DM J-factor directly from the raw simulations data. We show how this framework can be used to compute all-sky maps of the J-factor, automatically accounting for the complex 3D structure of the simulated halos and for the boosting of the signal due to the density fluctuations along the line of sight. After validating our code, we present a proof-of-concept application of NAJADS to a realistic halo from the IllustrisTNG suite, and exploit it to make a thorough comparison between our numerical approach and traditional semi-analytical methods. △ Less

Submitted 9 February, 2024; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: Preprint now matches published version in JCAP. 31 pages, 14 figures

Constraining positron emission from pulsar populations with AMS-02 data

Authors: Luca Orusa, Silvia Manconi, Fiorenza Donato, Mattia Di Mauro

Abstract: The cosmic-ray flux of positrons is measured with high precision by the space-borne particle spectrometer AMS-02. 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 $γ$-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe. In this work we unde… ▽ More The cosmic-ray flux of positrons is measured with high precision by the space-borne particle spectrometer AMS-02. 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 $γ$-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe. In this work we undertake massive simulations of Galactic pulsars populations, adopting different distributions for their position in the Galaxy, intrinsic physical properties, pair emission models, in order to overcome the incompleteness of the ATNF catalog. We fit the $e^+$ AMS-02 data together with a secondary component due to collisions of primary cosmic rays with the interstellar medium. We find that several mock galaxies have a pulsar population able to explain the observed $e^+$ flux, typically by few, bright sources. We determine the physical parameters of the pulsars dominating the $e^+$ flux, and assess the impact of different assumptions on radial distributions, spin-down properties, Galactic propagation scenarios and $e^+$ emission time. △ Less

Submitted 6 December, 2021; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 21 pages, 7 figures. Typos corrected, references added. Matches version published by JCAP

Report number: TTK-21-27

Journal ref: JCAP12(2021)014

Does the Geminga, Monogem and PSR J0622+3749 $γ$-ray halos imply slow diffusion around pulsars?

Authors: Sarah Recchia, Mattia Di Mauro, Felix A. Aharonian, Luca Orusa, Fiorenza Donato, Stefano Gabici, Silvia Manconi

Abstract: The HAWC Collaboration has reported the detection of an extended $γ$-ray emission around the Geminga and Monogem pulsars of a few degree extension. Very recently, the LHAASO Collaboration released also the data for an extended $γ$-ray emission around the pulsar PSR J0622+3749. This flux can be explained with electrons and positrons injected from these sources and their inverse Compton Scattering o… ▽ More The HAWC Collaboration has reported the detection of an extended $γ$-ray emission around the Geminga and Monogem pulsars of a few degree extension. Very recently, the LHAASO Collaboration released also the data for an extended $γ$-ray emission around the pulsar PSR J0622+3749. This flux can be explained with electrons and positrons injected from these sources and their inverse Compton Scattering on the interstellar radiation fields. So far the size of such $γ-$ray halos has been interpreted as the result of the diffusion coefficient around the sources being about two orders of magnitude smaller than the average in the Galaxy. However, this conclusion is driven by the assumption that particles propagate diffusively right away after the injection without taking into account the ballistic propagation. The propagation of cosmic-ray leptons in the proximity of the Geminga, Monogem and PSR J0622+3749 pulsars is examined here considering the transition from the quasi-ballistic, valid for the most recently injected particles, to the diffusive transport regime. For typical interstellar values of the diffusion coefficient, the quasi-ballistic regime dominates the lepton distribution up to distances of a few tens of parsec from the pulsar for particle energies above $\sim 10$ TeV. In this regime the resulting $γ-$ray source tends to be rather compact, despite particles travel a long distance. Indeed, for larger values of the diffusion coefficient, particles propagate ballistically up to larger distances with the result of a more point-like $γ-$ray source. When such transition is taken into account, a good fit to the HAWC and LHAASO $γ-$ray data around Geminga, Monogem and PSR J0622+3749 is obtained without the need to invoke a strong suppression of the diffusion coefficient. △ Less

Submitted 22 October, 2021; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 10 pages, 4 figures. This version includes other two sources (Monogem and PSR J0622+3749) with respect to the first version

Dissecting the inner Galaxy with $γ$-ray pixel count statistics

Authors: Francesca Calore, Fiorenza Donato, Silvia Manconi

Abstract: We combine adaptive template fitting and pixel count statistics in order to assess the nature of the Galactic center excess in Fermi-LAT data. We reconstruct the flux distribution of point sources well below the Fermi-LAT detection threshold, and measure their radial and longitudinal profiles in the inner Galaxy. We find that all point sources $and$ the bulge-correlated diffuse emission each contr… ▽ More We combine adaptive template fitting and pixel count statistics in order to assess the nature of the Galactic center excess in Fermi-LAT data. We reconstruct the flux distribution of point sources well below the Fermi-LAT detection threshold, and measure their radial and longitudinal profiles in the inner Galaxy. We find that all point sources $and$ the bulge-correlated diffuse emission each contributes $\mathcal{O}$(10\%) of the total inner Galaxy emission, and disclose a potential sub-threshold point-source contribution to the Galactic center excess. △ Less

Submitted 28 October, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 5+8 pages, 2+6 figures. Minor changes to match version published in PRL

Report number: LAPTH-008/21, TTK-21-06

Journal ref: Phys. Rev. Lett. 127, 161102 2021

Investigating $γ$-ray halos around three HAWC bright sources in Fermi-LAT data

Authors: Mattia Di Mauro, Silvia Manconi, Michela Negro, Fiorenza Donato

Abstract: Numerous extended sources around Galactic pulsars have shown significant $γ$-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 tha… ▽ More Numerous extended sources around Galactic pulsars have shown significant $γ$-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^\pm$ propagating beyond the SNR-PWN system into the interstellar medium. Following the discovery of few $γ$-ray sources by HAWC at energies E$>100$~TeV, we investigate the presence of an extended $γ$-ray emission in Fermi-LAT data around the three brightest sources detected by HAWC up to 100 TeV. We find an extended emission of $θ_{68} = 1.00^{+0.05}_{-0.07}$ deg around eHWC J1825-134 and $θ_{68} = 0.71\pm0.10$ deg eHWC J1907+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 $γ$-ray Fermi-LAT data provide invaluable insight into the very high-energy electron and positron parent populations. △ Less

Submitted 6 December, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 18 pages and 11 figures. Phys.Rev.D 104 (2021) 10, 103002

Novel interpretation of the latest AMS-02 cosmic-ray electron spectrum

Authors: Mattia Di Mauro, Fiorenza Donato, Silvia Manconi

Abstract: 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 spe… ▽ More 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σ$, 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. △ Less

Submitted 28 October, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 13 pages, 6 figures and supplemental material. Phys.Rev.D 104 (2021) 8, 083012

Cosmic-ray Antinuclei as Messengers of New Physics: Status and Outlook for the New Decade

Authors: P. von Doetinchem, K. Perez, T. Aramaki, S. Baker, S. Barwick, R. Bird, M. Boezio, S. E. Boggs, M. Cui, A. Datta, F. Donato, C. Evoli, L. Fabris, L. Fabbietti, E. Ferronato Bueno, N. Fornengo, H. Fuke, C. Gerrity, D. Gomez Coral, C. Hailey, D. Hooper, M. Kachelriess, M. Korsmeier, M. Kozai, R. Lea , et al. (37 additional authors not shown)

Abstract: The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essential… ▽ More The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years. △ Less

Submitted 18 August, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: 45 pages, 14 figures

Journal ref: JCAP08(2020)035

Contribution of pulsars to cosmic-ray positrons in light of recent observation of inverse-Compton halos

Authors: Silvia Manconi, Mattia Di Mauro, Fiorenza Donato

Abstract: 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 $γ$-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 $γ$-ray halos around these PWNe are interpreted as due to electrons (… ▽ More 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 $γ$-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 $γ$-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 $η\sim0.1$. For the simulated PWNe, a $1σ$ 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. △ Less

Submitted 30 January, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 10 pages, 7 figures. Figures 2, 3 and 5 updated. Results unchanged

Journal ref: Phys. Rev. D 102, 023015 (2020)

Testing gamma-ray models of blazars in the extragalactic sky

Authors: Silvia Manconi, Michael Korsmeier, Fiorenza Donato, Nicolao Fornengo, Marco Regis, Hannes Zechlin

Abstract: 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 statist… ▽ More 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. △ Less

Submitted 25 May, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: Matches the published version; extended discussion on detection efficiency, 4FGL fit improved, comments added. Main results and conclusions unchanged

Journal ref: Phys. Rev. D 101, 103026 (2020)

Evidences of low-diffusion bubbles around Galactic pulsars

Authors: Mattia Di Mauro, Silvia Manconi, Fiorenza Donato

Abstract: Recently, a few-degrees extended $γ$-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 r… ▽ More Recently, a few-degrees extended $γ$-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 $γ$-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 $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 6 December, 2021; v1 submitted 8 August, 2019; originally announced August 2019.

Comments: 21 pages, 11 figures. Phys.Rev.D 101 (2020) 10, 103035

Prospects for the detection of synchrotron halos around middle-age pulsars

Authors: Mattia Di Mauro, Silvia Manconi, Fiorenza Donato

Abstract: 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… ▽ More 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. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 7 pages and two figures. White paper submitted to Astro2020 (Astronomy and Astrophysics Decadal Survey)

Detection of a $γ$-ray halo around Geminga with the Fermi-LAT and implications for the positron flux

Authors: Mattia Di Mauro, Silvia Manconi, Fiorenza Donato

Abstract: The HAWC Collaboration has discovered a $γ$-ray emission extended about 2 degrees around the Geminga and Monogem pulsar wind nebulae (PWNe) at $γ$-ray energies $E_γ>5$ TeV. We analyze, for the first time, almost 10 years of $γ$-ray data obtained with the Fermi Large Area Telescope at $E_γ>$ 8 GeV in the direction of Geminga and Monogem. Since these two pulsars are close the Galactic plane we run o… ▽ More The HAWC Collaboration has discovered a $γ$-ray emission extended about 2 degrees around the Geminga and Monogem pulsar wind nebulae (PWNe) at $γ$-ray energies $E_γ>5$ TeV. We analyze, for the first time, almost 10 years of $γ$-ray data obtained with the Fermi Large Area Telescope at $E_γ>$ 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 $γ$-ray halo around Geminga with a significance in the range $7.8-11.8σ$ 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 $γ$ 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 $γ$-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. △ Less

Submitted 17 December, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 22 pages, 16 figures, This version matches the one published by PRD

Journal ref: Phys. Rev. D 100, 123015 (2019)

Multi-messenger constraints to the local emission of cosmic-ray electrons

Authors: Silvia Manconi, Mattia Di Mauro, Fiorenza Donato

Abstract: 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 Galact… ▽ More 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. △ Less

Submitted 12 April, 2019; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 20 pages, 9 figures. Accepted for publication in the JCAP. Changes in v3: discussion and results extended to include an evolutionary model for the injection of cosmic-ray electrons in SNRs

Journal ref: JCAP 1904 (2019) no.04, 024

Production cross sections of cosmic antiprotons in the light of new data from the NA61 and LHCb experiments

Authors: Michael Korsmeier, Fiorenza Donato, Mattia Di Mauro

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: 16 pages, 11 figures, matches published version

Report number: TTK-18-06

Journal ref: Phys. Rev. D 97, 103019 (2018)

Prospects to verify a possible dark matter hint in cosmic antiprotons with antideuterons and antihelium

Authors: Michael Korsmeier, Fiorenza Donato, Nicolao Fornengo

Abstract: Cosmic rays are an important tool to study dark matter annihilation in our Galaxy. Recently, a possible hint for dark matter annihilation was found in the antiproton spectrum measured by AMS-02, even though the result might be affected by theoretical uncertainties. A complementary way to test its dark matter interpretation would be the observation of low-energy antinuclei in cosmic rays. We determ… ▽ More Cosmic rays are an important tool to study dark matter annihilation in our Galaxy. Recently, a possible hint for dark matter annihilation was found in the antiproton spectrum measured by AMS-02, even though the result might be affected by theoretical uncertainties. A complementary way to test its dark matter interpretation would be the observation of low-energy antinuclei in cosmic rays. We determine the chances to observe antideuterons with GAPS and AMS-02, and the implications for the ongoing AMS-02 antihelium searches. We find that the corresponding antideuteron signal is within the GAPS and AMS-02 detection potential. If, more conservatively, the putative signal was considered as an upper limit on DM annihilation, our results would indicate the highest possible fluxes for antideuterons and antihelium compatible with current antiproton data. △ Less

Submitted 3 May, 2018; v1 submitted 22 November, 2017; originally announced November 2017.

Comments: 11 pages, 5 figures, matches published version

Report number: TTK-17-42

Journal ref: Phys. Rev. D 97, 103011 (2018)

Deriving the contribution of blazars to the Fermi-LAT Extragalactic $γ$-ray background at $E>10$ GeV with efficiency corrections and photon statistics

Authors: Mattia Di Mauro, Silvia Manconi, Hannes-S. Zechlin, Marco Ajello, Eric Charles, Fiorenza Donato

Abstract: 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-dete… ▽ More 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 $γ$-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 σ$. 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 $γ$-ray background originates from blazars. △ Less

Submitted 6 December, 2021; v1 submitted 8 November, 2017; originally announced November 2017.

Comments: 13 pages, 13 figures, 3 table. 2018 ApJ 856 106

Constraining Galactic dark matter with gamma-ray pixel counts statistics

Authors: Hannes-S. Zechlin, Silvia Manconi, Fiorenza Donato

Abstract: 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 sensi… ▽ More 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 <σ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 <σv> that are comparable with the constraints obtained by other indirect detection methods, in particular by the stacking analysis of several dwarf spheroidal galaxies. △ Less

Submitted 25 August, 2018; v1 submitted 4 October, 2017; originally announced October 2017.

Comments: 11 pages, 7 figures, 1 table; v2: major changes improving the selection of the ROI

Journal ref: Phys. Rev. D 98, 083022 (2018)

Prescriptions on antiproton cross section data for precise theoretical antiproton flux predictions

Authors: Fiorenza Donato, Michael Korsmeier, Mattia Di Mauro

Abstract: After the breakthrough from the satellite-borne PAMELA detector, the flux of cosmic-ray (CR) antiprotons has been provided with unprecedented accuracy by AMS-02 on the International Space Station. Its data spans an energy range from below 1 GeV up to 400 GeV and most of the data points contain errors below the amazing level of 5%. The bulk of the antiproton flux is expected to be produced by the s… ▽ More After the breakthrough from the satellite-borne PAMELA detector, the flux of cosmic-ray (CR) antiprotons has been provided with unprecedented accuracy by AMS-02 on the International Space Station. Its data spans an energy range from below 1 GeV up to 400 GeV and most of the data points contain errors below the amazing level of 5%. The bulk of the antiproton flux is expected to be produced by the scatterings of CR protons and helium off interstellar hydrogen and helium atoms at rest. The modeling of these interactions, which requires the relevant production cross sections, induces an uncertainty in the determination of the antiproton source term that can even exceed the uncertainties in the CR $\bar{p}$ data itself. The aim of the present analysis is to determine the uncertainty required for $p+p\rightarrow \bar{p} + X$ cross section measurements such that the induced uncertainties on the $\bar{p}$ flux are at the same level. Our results are discussed both in the center-of-mass reference frame, suitable for collider experiments, and in the laboratory frame, as occurring in the Galaxy. We find that cross section data should be collected with accuracy better that few percent with proton beams from 10 GeV to 6 TeV and a pseudorapidity $η$ ranging from 2 to almost 8 or, alternatively, with $p_T$ from 0.04 to 2 GeV and $x_R$ from 0.02 to 0.7. Similar considerations hold for the $p$He production channel. The present collection of data is far from these requirements. Nevertheless, they could, in principle, be reached by fixed target experiments with beam energies in the reach of CERN accelerators. △ Less

Submitted 4 June, 2018; v1 submitted 12 April, 2017; originally announced April 2017.

Comments: 15 pages, 13 figures, matches published version

Report number: TTK-17-08

Journal ref: Phys. Rev. D 96, 043007 (2017)

Theoretical interpretation of Pass 8 Fermi-LAT $e^+ + e^-$ data

Authors: Mattia Di Mauro, Silvia Manconi, Andrea Vittino, Fiorenza Donato, Nicolao Fornengo, Luca Baldini, Raffaella Bonino, Nicolò Di Lalla, Luca Latronico, Simone Maldera, Alberto Manfreda, Michela Negro, Melissa Pesce-Rollins, Carmelo Sgrò, Francesca Spada

Abstract: 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 inter… ▽ More 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. △ Less

Submitted 12 November, 2020; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: 15 pages, 12 figures. The Astrophysical Journal, Volume 845, Number 2

Report number: TUM-HEP-1078/17

Statistical Measurement of the Gamma-ray Source-count Distribution as a Function of Energy

Authors: H. -S. Zechlin, A. Cuoco, F. Donato, N. Fornengo, M. Regis

Abstract: Photon counts statistics have recently been proven to provide a sensitive observable for characterizing gamma-ray source populations and for measuring the composition of the gamma-ray sky. In this work, we generalize the use of the standard 1-point probability distribution function (1pPDF) to decompose the high-latitude gamma-ray emission observed with Fermi-LAT into: (i) point-source contribution… ▽ More Photon counts statistics have recently been proven to provide a sensitive observable for characterizing gamma-ray source populations and for measuring the composition of the gamma-ray sky. In this work, we generalize the use of the standard 1-point probability distribution function (1pPDF) to decompose the high-latitude gamma-ray emission observed with Fermi-LAT into: (i) point-source contributions, (ii) the Galactic foreground contribution, and (iii) a diffuse isotropic background contribution. We analyze gamma-ray data in five adjacent energy bands between 1 and 171 GeV. We measure the source-count distribution dN/dS as a function of energy, and demonstrate that our results extend current measurements from source catalogs to the regime of so far undetected sources. Our method improves the sensitivity for resolving point-source populations by about one order of magnitude in flux. The dN/dS distribution as a function of flux is found to be compatible with a broken power law. We derive upper limits on further possible breaks as well as the angular power of unresolved sources. We discuss the composition of the gamma-ray sky and capabilities of the 1pPDF method. △ Less

Submitted 21 February, 2017; originally announced February 2017.

Comments: 4 pages, 2 figures. Proceedings, 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma2016), Heidelberg, Germany, July 11-15, 2016

Journal ref: AIP Conference Proceedings 1792, 070020 (2017)

Anisotropies in GeV-TeV cosmic ray electrons and positrons

Authors: Silvia Manconi, Mattia Di Mauro, Fiorenza Donato

Abstract: 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 near… ▽ More 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. △ Less

Submitted 5 January, 2017; originally announced January 2017.

Comments: XXV ECRS 2016 Proceedings - eConf C16-09-04.3

Dipole anisotropy in cosmic electrons and positrons: inspection on local sources

Authors: Silvia Manconi, Mattia Di Mauro, Fiorenza Donato

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 26 pages, 14 figures. Comments are welcome

Journal ref: JCAP 1701 (2017) no.01, 006

Realistic estimation for the detectability of dark matter sub-halos with Fermi-LAT

Authors: Francesca Calore, Valentina De Romeri, Mattia Di Mauro, Fiorenza Donato, Federico Marinacci

Abstract: Numerical simulations of structure formation have recorded a remarkable progress in the recent years, in particular due to the inclusion of baryonic physics evolving with the dark matter component. We generate Monte Carlo realizations of the dark matter sub-halo population based on the results of the recent hydrodynamical simulation suite of Milky Way-sized galaxies. We then simulate the gamma-ray… ▽ More Numerical simulations of structure formation have recorded a remarkable progress in the recent years, in particular due to the inclusion of baryonic physics evolving with the dark matter component. We generate Monte Carlo realizations of the dark matter sub-halo population based on the results of the recent hydrodynamical simulation suite of Milky Way-sized galaxies. We then simulate the gamma-ray sky for both the setup of the 3FGL and 2FHL Fermi Large Area Telescope (LAT) catalogs, including the contribution from the annihilation of dark matter in the sub-halos. We find that the flux sensitivity threshold strongly depends on the particle dark matter mass, and more mildly also on its annihilation channel and the observation latitude. The results differ for the 3FGL and 2FHL catalogs, given their different energy thresholds. We also predict that the number of dark matter sub-halos among the unassociated sources is very small. A null number of detectable sub-halos in the Fermi-LAT 3FGL catalog would imply upper limits on the dark matter annihilation cross section into $b\bar{b}$ of $2 \cdot 10^{-26}$ ($5 \cdot 10^{-25}$) cm$^3$/s with $M_{\rm DM}$= 50 (1000) GeV. We find less than one extended sub-halo in the Fermi-LAT 3FGL catalog. As a matter of fact, the differences in the spatial and mass distribution of sub-halos between hydrodynamic and dark matter-only runs do not have significant impact on the gamma-ray dark matter phenomenology. △ Less

Submitted 10 November, 2016; originally announced November 2016.

Comments: 33 pages, 2 tables, 14 figures. Comments are welcome

Report number: LAPTH-064/16; IFT-UAM/CSIC-16-118; FTUAM-16-40

Journal ref: Phys. Rev. D 96, 063009 (2017)

Statistical Measurement of the Gamma-ray Source-count Distribution as a Function of Energy

Authors: Hannes-S. Zechlin, Alessandro Cuoco, Fiorenza Donato, Nicolao Fornengo, Marco Regis

Abstract: Statistical properties of photon count maps have recently been proven as a new tool to study the composition of the gamma-ray sky with high precision. We employ the 1-point probability distribution function of 6 years of Fermi-LAT data to measure the source-count distribution dN/dS and the diffuse components of the high-latitude gamma-ray sky as a function of energy. To that aim, we analyze the ga… ▽ More Statistical properties of photon count maps have recently been proven as a new tool to study the composition of the gamma-ray sky with high precision. We employ the 1-point probability distribution function of 6 years of Fermi-LAT data to measure the source-count distribution dN/dS and the diffuse components of the high-latitude gamma-ray sky as a function of energy. To that aim, we analyze the gamma-ray emission in five adjacent energy bands between 1 GeV and 171 GeV. It is demonstrated that the source-count distribution as a function of flux is compatible with a broken power law up to energies of ~50 GeV. The index below the break is between 1.95 and 2.0. For higher energies, a simple power-law fits the data, with an index of 2.2^{+0.7}_{-0.3} in the energy band between 50 GeV and 171 GeV. Upper limits on further possible breaks as well as the angular power of unresolved sources are derived. We find that point-source populations probed by this method can explain 83^{+7}_{-13}% (81^{+52}_{-19}%) of the extragalactic gamma-ray background between 1.04 GeV and 1.99 GeV (50 GeV and 171 GeV). The method has excellent capabilities for constraining the gamma-ray luminosity function and the spectra of unresolved blazars. △ Less

Submitted 21 July, 2016; v1 submitted 13 May, 2016; originally announced May 2016.

Comments: 6 pages, 2 figures, 2 tables; v2: minor changes, accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, Volume 826, Number 2, L31 (2016)

Unveiling the Gamma-ray Source Count Distribution Below the Fermi Detection Limit with Photon Statistics

Authors: Hannes-S. Zechlin, Alessandro Cuoco, Fiorenza Donato, Nicolao Fornengo, Andrea Vittino

Abstract: The source-count distribution as a function of their flux, dN/dS, is one of the main quantities characterizing gamma-ray source populations. We employ statistical properties of the Fermi Large Area Telescope (LAT) photon counts map to measure the composition of the extragalactic gamma-ray sky at high latitudes (|b|>30 deg) between 1 GeV and 10 GeV. We present a new method, generalizing the use of… ▽ More The source-count distribution as a function of their flux, dN/dS, is one of the main quantities characterizing gamma-ray source populations. We employ statistical properties of the Fermi Large Area Telescope (LAT) photon counts map to measure the composition of the extragalactic gamma-ray sky at high latitudes (|b|>30 deg) between 1 GeV and 10 GeV. We present a new method, generalizing the use of standard pixel-count statistics, to decompose the total observed gamma-ray emission into (a) point-source contributions, (b) the Galactic foreground contribution, and (c) a truly diffuse isotropic background contribution. Using the 6-year Fermi-LAT data set (P7REP), we show that the dN/dS distribution in the regime of so far undetected point sources can be consistently described with a power law of index between 1.9 and 2.0. We measure dN/dS down to an integral flux of ~2x10^{-11} cm^{-2}s^{-1}, improving beyond the 3FGL catalog detection limit by about one order of magnitude. The overall dN/dS distribution is consistent with a broken power law, with a break at 2.1^{+1.0}_{-1.3}x10^{-8} cm^{-2}s^{-1}. The power-law index n_1=3.1^{+0.7}_{-0.5} for bright sources above the break hardens to n_2=1.97+-0.03 for fainter sources below the break. A possible second break of the dN/dS distribution is constrained to be at fluxes below 6.4x10^{-11} cm^{-2}s^{-1} at 95% confidence level. The high-latitude gamma-ray sky between 1 GeV and 10 GeV is shown to be composed of ~25% point sources, ~69.3% diffuse Galactic foreground emission, and ~6% isotropic diffuse background. △ Less

Submitted 22 July, 2016; v1 submitted 22 December, 2015; originally announced December 2015.

Comments: 28 pages, 20 figures, 6 tables; v2: minor changes and additions, accepted for publication in ApJS

Report number: TUM-HEP 1034/15

Journal ref: The Astrophysical Journal Supplement Series, Volume 225, Number 2, 18 (2016)

Radio detection prospects for a bulge population of millisecond pulsars as suggested by Fermi LAT observations of the inner Galaxy

Authors: Francesca Calore, Mattia Di Mauro, Fiorenza Donato, Jason W. T. Hessels, Christoph Weniger

Abstract: Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in t… ▽ More Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in the inner 15 deg about the Galactic center (Fermi GeV excess). The excess can be interpreted as the collective emission of thousands of MSPs in the Galactic bulge, with a spherical distribution that strongly peaks towards the Galactic center. In order to fully establish the MSP interpretation, it is essential to find corroborating evidence in multi-wavelength searches, most notably through the detection of radio pulsation from individual bulge MSPs. Based on globular cluster observations and the gamma-ray emission from the inner Galaxy, we investigate the prospects for detecting MSPs in the Galactic bulge. While previous pulsar surveys failed to identify this population, we demonstrate that, in the upcoming years, new large-area surveys with focus on regions a few degrees north or south of the Galactic center should lead to the detection of dozens of bulge MSPs. Additionally, we show that, in the near future, deep targeted searches of unassociated Fermi sources should be able to detect the first few MSPs in the bulge. The prospects for these deep searches are enhanced by a tentative gamma-ray/radio correlation that we infer from high-latitude gamma-ray MSPs. Such detections would constitute the first clear discoveries of field MSPs in the Galactic bulge, with far-reaching implications for gamma-ray observations, the formation history of the central Milky Way and strategy optimization for future radio observations. △ Less

Submitted 12 September, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 24 pages, 17 figures, 5 tables. Minor clarifications. Matches version published in ApJ

Journal ref: Astrophys.J. 827 (2016) no.2, 143

Dark matter vs. astrophysics in the interpretation of AMS-02 electron and positron data

Authors: Mattia Di Mauro, Fiorenza Donato, Nicolao Fornengo, Andrea Vittino

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 13 July, 2016; v1 submitted 24 July, 2015; originally announced July 2015.

Comments: 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

Journal ref: JCAP 1605 (2016) 031

Status of cosmic-ray antideuteron searches

Authors: P. von Doetinchem, T. Aramaki, S. Boggs, S. Bufalino, L. Dal, F. Donato, N. Fornengo, H. Fuke, M. Grefe, C. Hailey, B. Hamilton, A. Ibarra, J. Mitchell, I. Mognet, R. A. Ong, R. Pereira, K. Perez, A. Putze, A. Raklev, P. Salati, M. Sasaki, G. Tarle, A. Urbano, A. Vittino, S. Wild , et al. (2 additional authors not shown)

Abstract: The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihil… ▽ More The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785). △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 9 pages, 4 figures, ICRC 2015 proceedings

Journal ref: PoS (ICRC 2015) 1218

Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuterons

Authors: T. Aramaki, S. Boggs, S. Bufalino, L. Dal, P. von Doetinchem, F. Donato, N. Fornengo, H. Fuke, M. Grefe, C. Hailey, B. Hamilton, A. Ibarra, J. Mitchell, I. Mognet, R. A. Ong, R. Pereira, K. Perez, A. Putze, A. Raklev, P. Salati, M. Sasaki, G. Tarle, A. Urbano, A. Vittino, S. Wild , et al. (2 additional authors not shown)

Abstract: Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons… ▽ More Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons benefit from very suppressed conventional backgrounds, offering a potential breakthrough in unexplored phase space for dark matter. This article is based on the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA in June 2014. It reviews broad classes of dark matter candidates that result in detectable cosmic-ray antideuteron fluxes, as well as the status and prospects of current experimental searches. The coalescence model of antideuteron production and the influence of antideuteron measurements at particle colliders are discussed. This is followed by a review of the modeling of antideuteron propagation through the magnetic fields, plasma currents, and molecular material of our Galaxy, the solar system, the Earth's geomagnetic field, and the atmosphere. Finally, the three ongoing or planned experiments that are sensitive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As cosmic-ray antideuteron detection is a rare event search, multiple experiments with orthogonal techniques and backgrounds are essential. Many theoretical and experimental groups have contributed to these studies over the last decade, this review aims to provide the first coherent discussion of the relevant dark matter theories that antideuterons probe, the challenges to predictions and interpretations of antideuteron signals, and the experimental efforts toward cosmic antideuteron detection. △ Less

Submitted 4 April, 2016; v1 submitted 28 May, 2015; originally announced May 2015.

Comments: 72 pages, 20 figures

The Connection Between the Positron Fraction Anomaly and the Spectral Features in Galactic Cosmic-Ray Hadrons

Authors: N. Tomassetti, F. Donato

Abstract: Recent data on Galactic cosmic-ray (CR) leptons and hadrons gave rise to two exciting problems: on the lepton side, the origin of the rise of the CR positron fraction e+/(e- + e+) at ~10 - 300 GeV of energy; on the hadron side, the nature of the spectral hardening observed in CR protons and nuclei at ~TeV energies. The lepton anomaly indicates the existence of a nearby e+/- source. It has been pro… ▽ More Recent data on Galactic cosmic-ray (CR) leptons and hadrons gave rise to two exciting problems: on the lepton side, the origin of the rise of the CR positron fraction e+/(e- + e+) at ~10 - 300 GeV of energy; on the hadron side, the nature of the spectral hardening observed in CR protons and nuclei at ~TeV energies. The lepton anomaly indicates the existence of a nearby e+/- source. It has been proposed that high-energy positrons can be produced inside nearby supernova remnants (SNRs) via interactions of CR hadrons with the ambient medium. A distinctive prediction of this mechanism is a high-energy rise of the boron-to-carbon ratio, which has not been observed. It also requires old SNRs at work (with ineffective magnetic field amplification and slow shock speed), that cannot account for the CR hadronic spectra observed up to the knee energies (~5 PeV). We propose a new picture where, in addition to such a nearby CR accelerator, the high-energy spectrum of CR hadrons is provided by the large-scale population of SNRs, on average younger, that can efficiently accelerate CRs up to the knee. Under this scenario, the spectral hardening of CR hadrons can be naturally interpreted as the transition between the two components. As we will show, our two-component model breaks the connection between the positron fraction and the boron-to-carbon ratio, which is now predicted to decrease with energy in accordance with the data. Forthcoming data from AMS will be crucial for testing this model. △ Less

Submitted 15 April, 2015; v1 submitted 21 February, 2015; originally announced February 2015.

Comments: 5 pages, 3 figures; matches the published version

Journal ref: Astrophysical Journal Letters 803 (2015) L15

Composition of the Fermi-LAT isotropic gamma-ray background intensity: Emission from extragalactic point sources and dark matter annihilations

Authors: Mattia Di Mauro, Fiorenza Donato

Abstract: A new estimation of the isotropic diffuse gamma-ray background (IGRB) observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) has been presented for 50 months of data, in the energy range 100 MeV-820 GeV and for different modelings of the Galactic foreground. We attempt here the interpretation of the Fermi-LAT IGRB data in terms of the gamma-ray unresolved e… ▽ More A new estimation of the isotropic diffuse gamma-ray background (IGRB) observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) has been presented for 50 months of data, in the energy range 100 MeV-820 GeV and for different modelings of the Galactic foreground. We attempt here the interpretation of the Fermi-LAT IGRB data in terms of the gamma-ray unresolved emission from different extragalactic populations. We find very good fits to the experimental IGRB, obtained with theoretical predictions for the emission from active galactic nuclei and star-forming galaxies. In addition, we probe a possible emission coming from the annihilation of weakly interacting dark matter (DM) particles in the halo of our Galaxy. We set stringent limits on its annihilation cross section into gamma rays, which are about the thermal relic value for a wide range of DM masses. We also identify regions in the DM mass and annihilation cross section parameter space which can significantly improve the fit to the IGRB data. Our analysis is conducted within the different IGRB data sets obtained from different models for the Galactic emission, which is shown to add a significant ambiguity on the IGRB interpretation. △ Less

Submitted 15 May, 2017; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: 17 pages, 9 figures, 5 tables. This version matches the published version, minor changes only

Journal ref: Phys.Rev. D91 (2015) 123001

Gamma rays from Galactic pulsars

Authors: Francesca Calore, Mattia Di Mauro, Fiorenza Donato

Abstract: Gamma rays from young pulsars and milli-second pulsars are expected to contribute to the diffuse gamma-ray emission measured by the {\it Fermi} Large Area Telescope (LAT) at high latitudes. We derive the contribution of the pulsars undetected counterpart by using information from radio to gamma rays and we show that they explain only a small fraction of the isotropic diffuse gamma-ray background. Gamma rays from young pulsars and milli-second pulsars are expected to contribute to the diffuse gamma-ray emission measured by the {\it Fermi} Large Area Telescope (LAT) at high latitudes. We derive the contribution of the pulsars undetected counterpart by using information from radio to gamma rays and we show that they explain only a small fraction of the isotropic diffuse gamma-ray background. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Comments: Proceeding of the Neutrino Oscillation Workshop, NOW 2014. 3 pages, 2 figures

A new evaluation of the antiproton production cross section for cosmic ray studies

Authors: Mattia di Mauro, Fiorenza Donato, Andreas Goudelis, Pasquale Dario Serpico

Abstract: 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 Collabor… ▽ More 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. △ Less

Submitted 26 July, 2018; v1 submitted 1 August, 2014; originally announced August 2014.

Comments: Two typos corrected (entry D4 in Table II and mishandled units in Eq. (10)). None of the results in the article are affected

Journal ref: Phys. Rev. D 90, 085017 (2014)

Fermi-LAT gamma-ray anisotropy and intensity explained by unresolved Radio-Loud Active Galactic Nuclei

Authors: Mattia Di Mauro, Alessandro Cuoco, Fiorenza Donato, Jennifer M. Siegal-Gaskins

Abstract: Radio-loud active galactic nuclei (AGN) are expected to contribute substantially to both the intensity and anisotropy of the isotropic gamma-ray background (IGRB). In turn, the measured properties of the IGRB can be used to constrain the characteristics of proposed contributing source classes. We consider individual subclasses of radio-loud AGN, including low-, intermediate-, and high-synchrotron-… ▽ More Radio-loud active galactic nuclei (AGN) are expected to contribute substantially to both the intensity and anisotropy of the isotropic gamma-ray background (IGRB). In turn, the measured properties of the IGRB can be used to constrain the characteristics of proposed contributing source classes. We consider individual subclasses of radio-loud AGN, including low-, intermediate-, and high-synchrotron-peaked BL Lacertae objects, flat-spectrum radio quasars, and misaligned AGN. Using updated models of the gamma-ray luminosity functions of these populations, we evaluate the energy-dependent contribution of each source class to the intensity and anisotropy of the IGRB. We find that collectively radio-loud AGN can account for the entirety of the IGRB intensity and anisotropy as measured by the Fermi Large Area Telescope (LAT). Misaligned AGN provide the bulk of the measured intensity but a negligible contribution to the anisotropy, while high-synchrotron-peaked BL Lacertae objects provide the dominant contribution to the anisotropy. In anticipation of upcoming measurements with the Fermi-LAT and the forthcoming Cherenkov Telescope Array, we predict the anisotropy in the broader energy range that will be accessible to future observations. △ Less

Submitted 2 December, 2014; v1 submitted 11 July, 2014; originally announced July 2014.

Comments: 27 pages, 29 figures. This version matches the published version, minor changes only

Journal ref: JCAP11(2014)021

Diffuse gamma-ray emission from galactic pulsars

Authors: F. Calore, M. Di Mauro, F. Donato

Abstract: Millisecond Pulsars are second most abundant source population discovered by the Fermi-LAT. They might contribute non-negligibly to the diffuse emission measured at high latitudes by Fermi-LAT, the IDGRB. Gamma-ray sources also contribute to the anisotropy of the IDGRB measured on small scales by Fermi-LAT. We aim to assess the contribution of the unresolved counterpart of the detected MSPs popula… ▽ More Millisecond Pulsars are second most abundant source population discovered by the Fermi-LAT. They might contribute non-negligibly to the diffuse emission measured at high latitudes by Fermi-LAT, the IDGRB. Gamma-ray sources also contribute to the anisotropy of the IDGRB measured on small scales by Fermi-LAT. We aim to assess the contribution of the unresolved counterpart of the detected MSPs population to the IDGRB and the maximal fraction of the measured anisotropy produced by this source class. We model the MSPs spatial distribution in the Galaxy and the gamma-ray emission parameters by considering radio and gamma-ray observational constraints. By simulating a large number of MSPs populations, we compute the average diffuse emission and the anisotropy 1-sigma upper limit. The emission from unresolved MSPs at 2 GeV, where the peak of the spectrum is located, is at most 0.9% of the measured IDGRB above 10 degrees in latitude. The 1-sigma upper limit on the angular power for unresolved MSP sources turns out to be about a factor of 60 smaller than Fermi-LAT measurements above 30 degrees. Our results indicate that this galactic source class represents a negligible contributor to the high-latitude gamma-ray sky and confirm that most of the intensity and geometrical properties of the measured diffuse emission are imputable to other extragalactic source classes. Nevertheless, given the MSP distribution, we expect them to contribute significantly to the gamma-ray diffuse emission at low latitudes. Since, along the galactic disk, the population of young Pulsars overcomes in number the one of MSPs, we compute the gamma-ray emission from the whole population of unresolved Pulsars in two low-latitude regions: the inner Galaxy and the galactic center. △ Less

Submitted 4 November, 2014; v1 submitted 10 June, 2014; originally announced June 2014.

Comments: 19 pages, 26 figures. It matches the published version, minor changes only

Journal ref: Astrophys.J. 796 (2014) 1, 14

Gamma-ray anisotropies from dark matter in the Milky Way: the role of the radial distribution

Authors: F. Calore, V. De Romeri, M. Di Mauro, F. Donato, J. Herpich, A. V. Maccio', L. Maccione

Abstract: The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simul… ▽ More The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations. △ Less

Submitted 18 June, 2014; v1 submitted 3 February, 2014; originally announced February 2014.

Comments: 7 pages, 8 figures. It matches the version published in MNRAS

Report number: LAPTH-009/14, PCCF RI 14-03

Journal ref: Monthly Notices of the Royal Astronomical Society 2014 442 (1): 1151-1156

Interpretation of AMS-02 electrons and positrons data

Authors: M. Di Mauro, F. Donato, N. Fornengo, R. Lineros, A. Vittino

Abstract: We perform a combined analysis of the recent AMS-02 data on electrons, positrons, electrons plus positrons and positron fraction, in a self-consistent framework where we realize a theoretical modeling of all the astrophysical components that can contribute to the observed fluxes in the whole energy range. The primary electron contribution is modeled through the sum of an average flux from distant… ▽ More We perform a combined analysis of the recent AMS-02 data on electrons, positrons, electrons plus positrons and positron fraction, in a self-consistent framework where we realize a theoretical modeling of all the astrophysical components that can contribute to the observed fluxes in the whole energy range. The primary electron contribution is modeled through the sum of an average flux from distant sources and the fluxes from the local supernova remnants in the Green catalog. The secondary electron and positron fluxes originate from interactions on the interstellar medium of primary cosmic rays, for which we derive a novel determination by using AMS-02 proton and helium data. Primary positrons and electrons from pulsar wind nebulae in the ATNF catalog are included and studied in terms of their most significant (while loosely known) properties and under different assumptions (average contribution from the whole catalog, single dominant pulsar, a few dominant pulsars). We obtain a remarkable agreement between our various modeling and the AMS-02 data for all types of analysis, demonstrating that the whole AMS-02 leptonic data admit a self-consistent interpretation in terms of astrophysical contributions. △ Less

Submitted 7 April, 2014; v1 submitted 3 February, 2014; originally announced February 2014.

Comments: 33 pages, 26 figures and 4 tables, v2: accepted for publication in JCAP, minor changes relative to v1

Report number: LAPTH--008/14, SACLAY--T14/010

Journal ref: JCAP04(2014)006

Diffuse $γ$-ray emission from unresolved BL Lac objects

Authors: M. Di Mauro, F. Donato, G. Lamanna, D. A. Sanchez, P. D. Serpico

Abstract: Blazars represent the most abundant class of high-energy extragalactic $γ$-ray sources. The subset of blazars known as BL Lac objects is on average closer to Earth and characterized by harder spectra at high energy than the whole sample. The fraction of BL Lacs that is too dim to be detected and resolved by current $γ$-ray telescopes is therefore expected to contribute to the high-energy isotropic… ▽ More Blazars represent the most abundant class of high-energy extragalactic $γ$-ray sources. The subset of blazars known as BL Lac objects is on average closer to Earth and characterized by harder spectra at high energy than the whole sample. The fraction of BL Lacs that is too dim to be detected and resolved by current $γ$-ray telescopes is therefore expected to contribute to the high-energy isotropic diffuse $γ$-ray background (IGRB). The IGRB has been recently measured over a wide energy range by the Large Area Telescope (LAT) on board the Gamma-ray Space Telescope ({\it Fermi}). We present a new prediction of the diffuse $γ$-ray flux due to the unresolved BL Lac blazar population. The model is built upon the spectral energy distribution and the luminosity function derived from the fraction of BL Lacs detected (and spectrally characterized) in the $γ$-ray energy range. We focus our attention on the ${\cal O}(100)$ GeV energy range, predicting the emission up to the TeV scale and taking into account the absorption on the extragalactic background light. In order to better shape the BL Lac spectral energy distribution, we combine the {\it Fermi}-LAT data with Imaging Atmospheric Cerenkov Telescopes measurements of the most energetic sources. Our analysis is carried on separately for low- and intermediate-synchrotron-peaked BL Lacs on one hand, and high-synchrotron-peaked BL Lacs on the other one: we find in fact statistically different features for the two. The diffuse emission from the sum of both BL Lac classes increases from about 10$\%$ of the measured IGRB at 100 MeV to $\sim$100$\%$ of the data level at 100 GeV. At energies greater than 100 GeV, our predictions naturally explain the IGRB data, accommodating their softening with increasing energy. Uncertainties are estimated to be within of a factor of two of the best-fit flux up to 500 GeV. △ Less

Submitted 24 May, 2016; v1 submitted 22 November, 2013; originally announced November 2013.

Comments: 14 pages, 6 figures. Version 4: Minor typos corrections, notably minus sign in the exponent of the square brackets of Eq.s A2, A6 and luminosity units in Table 2

Report number: LAPTH-067/13

Journal ref: 2014 ApJ 786 129

Diffuse $γ$-ray emission from misaligned active galactic nuclei

Authors: M. Di Mauro, F. Donato, F. Calore

Abstract: We calculate the diffuse $γ$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $γ$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-L… ▽ More We calculate the diffuse $γ$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $γ$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $γ$-ray luminosity function by means of the source count distribution of the MAGN detected by the {\it Fermi}-LAT. We finally estimate the diffuse $γ$-ray flux due to the whole MAGN population which ranges from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate also the room left to galactic DM at high latitudes ($>10^\circ$), by taking into account the results on the MAGN together with the other significant galactic and extragalactic $γ$-rays emitting sources. △ Less

Submitted 17 May, 2013; originally announced May 2013.

Comments: Proceedings of Rencontres de Moriond, La Thuile, March 9th - 16th, 2013

Diffuse $γ$-ray emission from misaligned active galactic nuclei

Authors: M. Di Mauro, F. Calore, F. Donato, M. Ajello, L. Latronico

Abstract: Active galactic nuclei (AGN) with jets seen at small viewing angles are the most luminous and abundant objects in the $γ$-ray sky. AGN with jets misaligned along the line-of-sight appear fainter in the sky, but are more numerous than the brighter blazars. We calculate the diffuse $γ$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\i… ▽ More Active galactic nuclei (AGN) with jets seen at small viewing angles are the most luminous and abundant objects in the $γ$-ray sky. AGN with jets misaligned along the line-of-sight appear fainter in the sky, but are more numerous than the brighter blazars. We calculate the diffuse $γ$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $γ$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $γ$-ray luminosity function by means of the source count distribution of the radio galaxies (RGs) detected by the {\it Fermi}-LAT. We finally calculate the diffuse $γ$-ray flux due to the whole MAGN population. Our results demonstrate that the MAGN can contribute from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate a theoretical uncertainty on the flux of almost an order of magnitude. △ Less

Submitted 18 November, 2013; v1 submitted 3 April, 2013; originally announced April 2013.

Comments: Accepted for publication in ApJ

Constraining dark matter annihilation with the isotropic $γ$-ray background: updated limits and future potential

Authors: Torsten Bringmann, Francesca Calore, Mattia Di Mauro, Fiorenza Donato

Abstract: The nature of the Isotropic $γ$-ray Background (IGRB) measured by the Large Area Telescope (LAT) on the Fermi $γ$-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 $γ$-ray emission… ▽ More The nature of the Isotropic $γ$-ray Background (IGRB) measured by the Large Area Telescope (LAT) on the Fermi $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 3 February, 2014; v1 submitted 13 March, 2013; originally announced March 2013.

Comments: 11 pages revtex4, 9 figures. Matches published version

Journal ref: Phys. Rev. D 89,023012 (2014)

Secondary Cosmic Ray Nuclei from Supernova Remnants and Constraints to the Propagation Parameters

Authors: N. Tomassetti, F. Donato

Abstract: The secondary-to-primary B/C ratio is widely used to study the cosmic ray (CR) propagation processes in the Galaxy. It is usually assumed that secondary nuclei such as Li-Be-B are entirely generated by collisions of heavier CR nuclei with the interstellar medium (ISM). We study the CR propagation under a scenario where secondary nuclei can also be produced or accelerated from galactic sources. We… ▽ More The secondary-to-primary B/C ratio is widely used to study the cosmic ray (CR) propagation processes in the Galaxy. It is usually assumed that secondary nuclei such as Li-Be-B are entirely generated by collisions of heavier CR nuclei with the interstellar medium (ISM). We study the CR propagation under a scenario where secondary nuclei can also be produced or accelerated from galactic sources. We consider the processes of hadronic interactions inside supernova remnants (SNRs) and re-acceleration of background CRs in strong shocks. Thus, we investigate their impact in the propagation parameter determination within present and future data. The spectra of Li-Be-B nuclei emitted from SNRs are harder than those due to CR collisions with the ISM. The secondary-to-primary ratios flatten significantly at ~TeV/n energies, both from spallation and re-acceleration in the sources. The two mechanisms are complementary to each other and depend on the properties of the local ISM around the expanding remnants. The secondary production in SNRs is significant for dense background media, n ~1 cm^-3, while the amount of re-accelerated CRs is relevant for SNRs expanding into rarefied media, n ~0.1 cm-3. Due to these effects, the the diffusion parameter 'delta' may be misunderstood by a factor of ~5-15%. Our estimations indicate that an experiment of the AMS-02 caliber can constrain the key propagation parameters while breaking the source-transport degeneracy, for a wide class of B/C-consistent models. Given the precision of the data expected from on-going experiments, the SNR production/acceleration of secondary nuclei should be considered, if any, to prevent a possible mis-determination of the CR transport parameters. △ Less

Submitted 22 July, 2012; v1 submitted 27 March, 2012; originally announced March 2012.

Comments: 13 pages, 9 figures; matches the published version

Journal ref: A&A 544, A16 (2012)

Conservative upper limits on WIMP annihilation cross section from Fermi-LAT $γ$-rays

Authors: F. Donato, F. Calore, V. De Romeri

Abstract: The spectrum of an isotropic extragalactic $γ$-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subtraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation… ▽ More The spectrum of an isotropic extragalactic $γ$-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subtraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain $conservative$ upper limits on their annihilation cross section \sigmav. Severe bounds on a possible Sommerfeld enhancement of the annihilation cross section are set as well. Finally, would {\sigmav} be inversely proportional to the WIMP velocity, very severe limits are derived for the velocity-independent part of the annihilation cross section. △ Less

Submitted 18 December, 2011; originally announced December 2011.

Comments: Proceedings of XII Taup Conference, Munich, September 2011

Astrophysics of Galactic charged cosmic rays

Authors: Antonella Castellina, Fiorenza Donato

Abstract: A review is given of the main properties of the charged component of galactic cosmic rays, particles detected at Earth with an energy spanning from tens of MeV up to about 10^19 eV. After a short introduction to the topic and an historical overview, the properties of cosmic rays are discussed with respect to different energy ranges. The origin and the propagation of nuclei in the Galaxy are dealt… ▽ More A review is given of the main properties of the charged component of galactic cosmic rays, particles detected at Earth with an energy spanning from tens of MeV up to about 10^19 eV. After a short introduction to the topic and an historical overview, the properties of cosmic rays are discussed with respect to different energy ranges. The origin and the propagation of nuclei in the Galaxy are dealt with from a theoretical point of view. The mechanisms leading to the acceleration of nuclei by supernova remnants and to their subsequent diffusion through the inhomogeneities of the galactic magnetic field are discussed and some clue is given on the predictions and observations of fluxes of antimatter, both from astrophysical sources and from dark matter annihilation in the galactic halo. The experimental techniques and instrumentations employed for the detection of cosmic rays at Earth are described. Direct methods are viable up to 10^14 eV, by means of experiments flown on balloons or satellites, while above that energy, due to their very low flux, cosmic rays can be studied only indirectly by exploiting the particle cascades they produce in the atmosphere. The possible physical interpretation of the peculiar features observed in the energy spectrum of galactic cosmic rays, and in particular the so-called "knee" at about 4x10^15 eV, are discussed. A section is devoted to the region between about 10^18 and 10^19 eV, which is believed to host the transition between galactic and extragalactic cosmic rays. The conclusion gives some perspectives on the cosmic ray astrophysics field. Thanks to a wealth of different experiments, this research area is living a very flourishing era. The activity is exciting both from the theoretical and the instrumental sides, and its interconnection with astronomy, astrophysics and particle physics experiences non-stop growth. △ Less

Submitted 13 October, 2011; originally announced October 2011.

Comments: Invited review, to appear in Vol.5 of "Planets, Stars and Stellar Systems" Editor-in-chief: Oswalt, T.D McLean, I.S.; Bond, H.E.; French, L.; Kalas, P.; Barstow, M.; Gilmore,G.F.; Keel, W. (Eds.) 1st Edition., 2012, 4760 p., Springer, ISBN 978-90-481-8817-8 Due: February 15, 2012