The cumulative emission resulting from hadronic cosmic-ray interactions in star-forming galaxies (SFGs) has been proposed as the dominant contribution to the astrophysical neutrino flux at TeV to PeV ...energies reported by IceCube. The same particle interactions also inevitably create γ-ray emission that could be detectable as a component of the extragalactic γ-ray background (EGB), which is now measured with the Fermi-LAT in the energy range from 0.1 to 820 GeV. New studies of the blazar flux distribution at γ-ray energies above 50 GeV place an upper bound on the residual non-blazar component of the EGB. We show that these results are in strong tension with models that consider SFGs as the dominant source of the diffuse neutrino backgrounds. A characteristic spectral index for parent cosmic rays in starburst galaxies of ΓSB 2.3 for dN / dE ∝ E − Γ SB is consistent with the observed scaling relation between γ-ray and IR luminosity for SFGs, the bounds from the non-blazar EGB, and the observed γ-ray spectra of individual starbursts, but underpredicts the IceCube data by approximately an order of magnitude.
Numerical simulations of structure formation have made remarkable progress in 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 subhalo population based on the results of the recent hydrodynamical simulation suite of Milky Way–sized galaxies F. Marinacci, R. Pakmor, and V. Springel, Mon. Not. R. Astron. Soc. 437, 1750 (2014).. 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 subhalos. 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 subhalos among the unassociated sources is very small. A null number of detectable subhalos in the Fermi-LAT 3FGL catalog would imply upper limits on the dark matter annihilation cross section into bb¯ of 2×10−26(5×10−25) cm3/s with MDM=50(1000) GeV. We find less than one extended subhalo in the Fermi-LAT 3FGL catalog. As a matter of fact, the differences in the spatial and mass distribution of subhalos between hydrodynamic and dark matter–only runs do not have significant impact on the detectability of dark subhalos in gamma rays.
A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has recently been published by the Fermi-LAT collaboration, reporting that the emission ...broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). The weight of the evidence is steadily accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare these observations with what is, perhaps, the simplest model for an MSP population, which is solely obtained by rescaling of the MSP luminosity function that is determined in the local MW disk via the respective stellar mass of the systems. Remarkably, we find that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to only contribute about a quarter of the detected emission, although a stronger contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSPs is not far below the present upper bound. We also discuss additional refinements of this simple model. Using the correlation between globular cluster gamma-ray luminosity and stellar encounter rate, we gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level of 5%, it could affect the signal's morphology. We also comment on the limitations of our model and on future perspectives for improved diagnostics.
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 light of new cross section measurements by the NA61 experiment of p+p→p¯+X and the first ever measurement of p+He→p¯+X by the LHCb experiment, we update the parametrization of proton-proton and proton-nucleon cross sections. We find that the LHCb pHe 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, pHe and pC data we estimate the uncertainty on the Lorentz invariant cross section for p+He→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 up to ±20% and slightly increase below antiproton energies of 5 GeV. This uncertainty is dominated by the p+p→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 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.
ABSTRACT
The γ-ray emission from stars is induced by the interaction of cosmic rays with stellar atmospheres and photon fields. This emission is expected to come in two components: a stellar disc ...emission, where γ-rays are mainly produced in atmospheric showers generated by hadronic cosmic rays, and an extended halo emission, where the high density of soft photons in the surroundings of stars create a suitable environment for γ-ray production via inverse Compton (IC) scattering by cosmic ray electrons. Besides the Sun, no other disc or halo from single stars has ever been detected in γ-rays. However, by assuming a cosmic ray spectrum similar to that observed on Earth, the predicted γ-ray emission of superluminous stars, e.g. Betelgeuse and Rigel, could be high enough to be detected by the Fermi Large Area Telescope (LAT) after its first decade of operations. In this work, we use 12 yr of Fermi-LAT observations along with IC models to study nine superluminous nearby stars, both individually and via stacking analysis. Our results show no significant γ-ray emission, but allow us to restrict the stellar γ-ray fluxes to be on average <3.3 × 10−11 ph cm−2 s−1 at a 3σ confidence level, which translates to an average local density of electrons in the surroundings of our targets to be less than twice of that observed for the Solar system.
The excess of γ rays in the data measured by the Fermi Large Area Telescope from the galactic center region is one of the most intriguing mysteries in astroparticle physics. This "galactic center ...excess" (GCE) has been measured with respect to different interstellar emission models (IEMs), source catalogs, data selections, and techniques. Although several proposed interpretations have appeared in the literature, there are no firm conclusions as to its origin. The main difficulty in solving this puzzle lies in modeling a region of such complexity and thus, precisely measuring the characteristics of the GCE. In this paper, we use 11 years of Fermi-LAT data, state of the art IEMs, and the newest 4FGL source catalog to provide precise measurements of the energy spectrum, spatial morphology, position, and sphericity of the GCE. We find that the GCE has a spectrum that is peaked at a few GeV and is well fit with a log parabola. The normalization of the spectrum changes by roughly 60% when using different IEMs, data selections, and analysis techniques. The spatial distribution of the GCE is compatible with a dark matter (DM) template produced with a generalized Navarro-Frenk-White density profile with slope γ=1.2–1.3. No energy evolution is measured for the GCE morphology between 0.6–30 GeV at a level larger than 10% of the γ average value, which is 1.25. The analysis of the GCE modeled with a DM template divided into quadrants shows that the spectrum and spatial morphology of the GCE is similar in different regions around the galactic center. Finally, the GCE centroid is compatible with the galactic center, with a best-fit position between l=−0.3°, 0.0 ° , b = − 0.1 ° , 0.0 ° , and it is compatible with a spherical symmetric morphology. In particular, fitting the DM spatial profile with an ellipsoid gives a major-to-minor axis ratio (aligned along the galactic plane) between 0.8–1.2.