No Arabic abstract
The Tibet AS$gamma$ experiment provided the first measurement of the total diffuse gamma-ray emission from the Galactic disk in the sub-PeV energy range. Based on the analysis of TeV sources included in the HGPS catalog, we predict the expected contribution of unresolved sources in the two angular windows of the Galactic plane observed by Tibet AS$gamma$. We show that the sum of this additional diffuse component due to unresolved sources and the truly diffuse emission, due to cosmic ray interaction with the interstellar medium, well saturates the Tibet data, without the need to introduce a progressive hardening of the cosmic-ray spectrum toward the Galactic center.
Tibet-AS$gamma$ collaboration has recently reported a measurement of diffuse $gamma$-ray flux from the outer Galactic disk in the energy range reaching PeV. We complement this measurement with the Fermi/LAT measurement of the diffuse flux from the same sky region and study the pion decay model of the combined Fermi/LAT+Tibet-AS$gamma$ spectrum. We find that within such a model the average cosmic ray spectrum in the outer Galactic disk has the same characteristic features as the local cosmic ray spectrum. In particular, it experiences a hardening at several hundred GV rigidity and a knee feature in the PV rigidity range. The slope of the average cosmic ray spectrum above the break is close to the locally observed slope of the helium spectrum $gammasimeq 2.5$, but is harder than the slope of the local proton spectrum in the same rigidity range. Although the combination of Fermi/LAT and Tibet-AS$gamma$ data points to the presence of the knee in the average cosmic ray spectrum, the quality of the data is not yet sufficient for the study of knee shape and cosmic ray composition.
Blazars represent the most abundant class of high-energy extragalactic $gamma$-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 $gamma$-ray telescopes is therefore expected to contribute to the high-energy isotropic diffuse $gamma$-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 $gamma$-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 $gamma$-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.
Tibet AS-gamma collaboration has recently reported detection of gamma-rays with energies up to Peta-electronvolt from parts of the Galactic plane. We notice that the analysis of gamma-ray flux by the Tibet AS-gamma experiment also implies an upper bound on the diffuse gamma-ray flux from high Galactic latitudes(|b|>20 degrees) in the energy range between 100 TeV and 1 PeV. This bound is up to an order-of-magnitude stronger than previously derived bounds from GRAPES3, KASCADE and CASA-MIA experiments. We discuss the new TibetAS-gamma limit on high Galactic latitude gamma-ray flux in the context of possible mechanisms of multi-messenger (gamma-ray and neutrino) emission from nearby cosmic ray sources, dark matter decays and large scale cosmic ray halo of the Milky Way.
The Tibet ASgamma experiment just reported their measurement of sub-PeV diffuse gamma ray emission from the Galactic disk, with the highest energy up to 957 TeV. These gamma-rays are most likely the hadronic origin by cosmic ray interaction with interstellar gas in the Galaxy. This measurement provides direct evidence to the hypothesis that the Galactic cosmic rays can be accelerated beyond PeV energies. In this work, we try to explain the sub-PeV diffuse gamma-ray spectrum within cosmic rays diffusive propagation model. We find there is a tension between the sub-PeV diffuse gamma rays and the local cosmic ray spectrum. To describe the sub-PeV diffuse gamma-ray flux, it generally requires larger local cosmic-ray flux than measurement in the knee region. We further calculate the PeV neutrino flux from the cosmic ray propagation model. Even all of these sub-PeV diffuse gamma rays originate from the propagation, the Galactic neutrinos only account for less than ~15% of observed flux, most of which are still from extragalactic sources.
Very recently, diffuse gamma rays with $0.1,{rm PeV}<E_gamma <1,rm PeV$ have been discovered from the Galactic disk by the Tibet air shower array and muon detector array (Tibet AS+MD array). While the measured sub-PeV flux may be compatible with the hadronic origin in the conventional Galactic cosmic ray propagation model, we find that it is in possible tension with the non-detection of Galactic neutrino emissions by the IceCube neutrino telescope. We further find that the presence of an extra cosmic ray component of relatively hard spectrum, which is probably related to the Cygnus Cocoon region and other PeV cosmic-ray sources in the Galactic disk, would alleviate the tension. This scenario implies the existence of an extreme accelerator of either protons or electrons beyond PeV in the Cygnus region, and predicts the continuation of the gamma-ray spectrum of Cygnus Cocoon up to 1 PeV with a possible hardening beyond $sim 30-100,$TeV.