No Arabic abstract
We examine the cosmic-ray protons (CRp) accelerated at collisionless shocks in galaxy clusters using cosmological structure formation simulations. We find that in the intracluster medium (ICM) within the virial radius of simulated clusters, only $sim7$% of shock kinetic energy flux is dissipated by the shocks that are expected to accelerate CRp, that is, supercritical, quasi-parallel ($Q_parallel$) shocks with sonic Mach number $M_sge2.25$. The rest is dissipated at subcritical shocks and quasi-perpendicular shocks, both of which may not accelerate CRp. Adopting the diffusive shock acceleration (DSA) model recently presented in Ryu et al. (2019), we quantify the DSA of CRp in simulated clusters. The average fraction of the shock kinetic energy transferred to CRp via DSA is assessed at $sim(1-2)times10^{-4}$. We also examine the energization of CRp through reacceleration using a model based on the test-particle solution. Assuming that the ICM plasma passes through shocks three times on average through the history of the universe and that CRp are reaccelerated only at supercritical $Q_parallel$-shocks, the CRp spectrum flattens by $sim0.05-0.1$ in slope and the total amount of CRp energy increases by $sim40-80$% from reacceleration. We then estimate diffuse $gamma$-ray and neutrino emissions, resulting from inelastic collisions between CRp and thermal protons. The predicted $gamma$-ray emissions from simulated clusters lie mostly below the upper limits set by Fermi-LAT for observed clusters. The neutrino fluxes towards nearby clusters would be $lesssim10^{-4}$ of the IceCube flux at $E_{ u}=1$ PeV and $lesssim10^{-6}$ of the atmospheric neutrino flux in the energy range of $E_{ u}leq1$ TeV.
Gamma-ray bursts are short-lived, luminous explosions at cosmological distances, thought to originate from relativistic jets launched at the deaths of massive stars. They are among the prime candidates to produce the observed cosmic rays at the highest energies. Recent neutrino data have, however, started to constrain this possibility in the simplest models with only one emission zone. In the classical theory of gamma-ray bursts, it is expected that particles are accelerated at mildly relativistic shocks generated by the collisions of material ejected from a central engine. We consider neutrino and cosmic-ray emission from multiple emission regions since these internal collisions must occur at very different radii, from below the photosphere all the way out to the circumburst medium, as a consequence of the efficient dissipation of kinetic energy. We demonstrate that the different messengers originate from different collision radii, which means that multi-messenger observations open windows for revealing the evolving GRB outflows.
Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $gamma$-ray flux from individual clusters in our sample.
Blazars are potential candidates of cosmic-ray acceleration up to ultrahigh energies ($Egtrsim10^{18}$ eV). For an efficient cosmic-ray injection from blazars, $pgamma$ collisions with the extragalactic background light (EBL) and cosmic microwave background (CMB) can produce neutrino spectrum peaks near PeV and EeV energies, respectively. We analyze the contribution of these neutrinos to the diffuse background measured by the IceCube neutrino observatory. The fraction of neutrino luminosity originating from individual redshift ranges is calculated using the distribution of BL Lacs and FSRQs provided in the textit{Fermi}-LAT 4LAC catalog. Furthermore, we use a luminosity dependent density evolution to find the neutrino flux from unresolved blazars. The results obtained in our model indicate that as much as $approx10%$ of the flux upper bound at a few PeV energies can arise from cosmic-ray interactions on EBL. The same interactions will also produce secondary electrons and photons, initiating electromagnetic cascades. The resultant photon spectrum is limited by the isotropic diffuse $gamma$-ray flux measured between 100 MeV and 820 GeV. The latter, together with the observed cosmic-ray flux at $E>10^{16.5}$ eV, can constrain the baryonic loading factor depending on the maximum cosmic-ray acceleration energy.
The Fermi gamma-ray satellite has recently detected gamma-ray emissions from radio galaxy cores. From these samples, we first examine the correlation between the luminosities at 5 GHz, L_{5GHz}, and at 0.1-10 GeV, L_{gamma}, of these gamma-ray loud radio galaxies. We find that the correlation is significant with L_{gamma} propto L_{5GHz}^{1.16} based on a partial correlation analysis. Using this correlation and the radio luminosity function (RLF) of radio galaxies, we further explore the contribution of gamma-ray loud radio galaxies to the unresolved extragalactic gamma-ray background (EGRB). The gamma-ray luminosity function is obtained by normalizing the RLF to reproduce the source count distribution of the Fermi gamma-ray loud radio galaxies. We find that gamma-ray loud radio galaxies will explain ~25% of the unresolved Fermi EGRB flux above 100 MeV and will also make a significant contribution to the EGRB in the 1-30 MeV energy band. Since blazars explain 22% of the EGRB above 100 MeV, radio loud active galactic nuclei (AGNs) population explains ~47% of the unresolved EGRB. We further make an interpretation on the origin of the EGRB. The observed EGRB spectrum at 0.2-100 GeV does not show an absorption signature by the extragalactic background light. Thus, the dominant population of the origin of EGRB at very high energy (>30 GeV) might be nearby gamma-ray emitting sources or sources with very hard gamma-ray spectrum.
Supernova remnants (SNRs) are believed to accelerate particles up to high energies through the mechanism of diffusive shock acceleration (DSA). Except for direct plasma simulations, all modeling efforts must rely on a given form of the diffusion coefficient, a key parameter that embodies the interactions of energetic charged particles with the magnetic turbulence. The so-called Bohm limit is commonly employed. In this paper we revisit the question of acceleration at perpendicular shocks, by employing a realistic model of perpendicular diffusion. Our coefficient reduces to a power-law in momentum for low momenta (of index $alpha$), but becomes independent of the particle momentum at high momenta (reaching a constant value $kappa_{infty}$ above some characteristic momentum $p_{rm c}$). We first provide simple analytical expressions of the maximum momentum that can be reached at a given time with this coefficient. Then we perform time-dependent numerical simulations to investigate the shape of the particle distribution that can be obtained when the particle pressure back-reacts on the flow. We observe that, for a given index $alpha$ and injection level, the shock modifications are similar for different possible values of $p_{rm c}$, whereas the particle spectra differ markedly. Of particular interest, low values of $p_{rm c}$ tend to remove the concavity once thought to be typical of non-linear DSA, and result in steep spectra, as required by recent high-energy observations of Galactic SNRs.