Future Gamma-ray Observations Research Articles

Gamma rays from dark matter spikes in EAGLE simulations

Intermediate Mass Black Holes (IMBHs) with a mass range between 100 M⊙ and 106 M⊙ are expected to be surrounded by high dark matter densities, so-called dark matter spikes.The high density of self-annihilating Weakly Interacting Massive Particles (WIMPs) in these spikes leads to copious gamma-ray production. Sufficiently nearby IMBHs could therefore appear as unidentified gamma-ray sources.However, the number of IMBHs and their distribution within our own Milky Way is currently unknown. In this work, we provide a mock catalogue of IMBHs and their dark matter spikes obtained from the EAGLE simulations, in which black holes with a mass of 105 M⊙/h are seeded into the centre of halos greater than 1010 M⊙/h to model black hole feedback influencing the formation of galaxies.The catalogue contains the coordinates and dark matter spike parameters for about 2500 IMBHs present in about 150 Milky Way-like galaxies. We expect about 15+9 -6 IMBHs within our own galaxy, mainly distributed in the Galactic Centre and the Galactic Plane. In the most optimistic scenario, we find that current and future gamma-ray observatories, such as Fermi-LAT, H.E.S.S. and CTAO, would be sensitive enough to probe the cross section of dark matter self-annihilation around IMBHs down to many orders of magnitude below the thermal relic cross section for dark matter particles with masses from GeV to TeV.We have made the IMBH mock catalogue and the source code for our analysis publicly available, providing the resources to study dark matter self-annihilation around IMBHs with current and upcoming gamma-ray observatories.

Multiwavelength Modeling for the Shallow Decay Phase of Gamma-Ray Burst Afterglows

We simulate the emission in the shallow decay phase of gamma-ray burst afterglows using a time-dependent code. We test four models: the energy injection model, evolving the injection efficiency of nonthermal electrons, evolving the amplification of the magnetic field, and the wind model with a relatively low bulk Lorentz factor. All of the four models can reproduce the typical X-ray afterglow lightcurve. The spectral shape depends on not only the parameter values at the time corresponding to the observer time but also the past evolution of the parameters. The model differences appear in the evolution of the broadband spectrum, especially in the inverse Compton component. Future gamma-ray observations with imaging atmospheric Cherenkov telescopes such as the Cherenkov Telescope Array will reveal the mechanism of the shallow decay phase.

Observe Gamma-Rays and Neutrinos Associated with Ultra-High Energy Cosmic Rays

IceCube measures a diffuse neutrino flux comparable to the Waxman-Bahcall bound, which suggests the possibility that the ultra-high energy cosmic rays (UHECRs) have a common origin with diffuse high energy neutrinos. We propose high energy gamma-ray and/or neutrino observations toward the arrival directions of UHECRs to search for the sources and test this possibility. We calculate the detection probability of gamma-ray/neutrino sources, and find that the average probability per UHECR of >10 EeV is ∼10% if the sensitivity of the gamma-ray or neutrino telescope is ∼10−12 erg cm−2 s−1 and the source number density is ∼10−5 Mpc−3. Future gamma-ray and neutrino observations toward UHECRs, e.g., by LHAASO-WCDA, CTA, IceCube/Gen2, are encouraged to constrain the density of UHECR sources or even identify the sources of UHECRs.

The RTApipe framework for the gamma-ray real-time analysis software development

In the multi-messenger era, coordinating observations between astronomical facilities is mandatory to study transient phenomena (e.g. Gamma-ray bursts) and is achieved by sharing information with the scientific community through networks such as the Gamma-ray Coordinates Network. The facilities usually develop real-time scientific analysis pipelines to detect transient events, alert the astrophysical community, and speed up the reaction time of science alerts received from other observatories. We present in this work the RTApipe framework, designed to facilitate the development of real-time scientific analysis pipelines for present and future gamma-ray observatories. This framework provides pipeline architecture and automatisms, allowing the researchers to focus on the scientific aspects and integrate existing science tools developed with different technologies. The pipelines automatically execute all the configured analyses during the data acquisition. This framework can be interfaced with science alerts networks to perform follow-up analysis of transient events shared by other facilities. The analyses are performed in parallel and can be prioritised. The workload is highly scalable on a cluster of machines. The framework provides the required services using containerisation technology for easy deployment. We present the RTA pipelines developed for the AGILE space mission and the prototype of the SAG system for the ground-based future Cherenkov Telescope Array observatory confirming that the RTApipe framework can be used to successfully develop pipelines for the gamma-ray observatories, both space and ground-based.

Formula omitted]: A new variable for γ/h discrimination in large gamma-ray ground arrays

In this letter, a new strategy to enhance the discrimination of high energy gamma rays from the huge charged cosmic rays background in large cosmic rays ground arrays is presented. This strategy is based on the introduction of a new simple variable, Pγhα, which combines the probability of tagging muons and/or very energetic particles in each single array station. The discrimination power of this new variable, particularly important for and above multi-TeV energies, is illustrated for a few specific examples in the case of a hypothetical water Cherenkov detector cosmic ray array, both in the case of low and high particle stations occupancy. The results are very encouraging and hopefully will be demonstrated in the present and future gamma-ray Observatories.

Photohadronic modelling of the 2010 gamma-ray flare from Mrk 421

ABSTRACT Blazars are a subclass of active galactic nuclei (AGNs) that have a relativistic jet with a small viewing angle towards the observer. Recent results based on hadronic scenarios have motivated an ongoing discussion of how a blazar can produce high energy neutrinos during a flaring state and which scenario can successfully describe the observed gamma-ray behaviour. Markarian 421 is one of the closest and brightest objects in the extragalactic gamma-ray sky and showed flaring activity over a 14-days period in 2010 March. In this work, we describe the performed analysis of Fermi-LAT data from the source focused on the MeV range (100 MeV–1 GeV), and study the possibility of a contribution coming from the pγ interactions between protons and MeV SSC target photons to fit the very high energy (VHE) gamma-ray emission. The fit results were compared with two leptonic models (one-zone and two-zone) using the Akaike Information Criteria (AIC) test, which evaluates goodness-of-fit alongside the simplicity of the model. In all cases, the photohadronic model was favoured as a better fit description in comparison to the one-zone leptonic model, and with respect to the two-zone model in the majority of cases. Our results show the potential of a photohadronic contribution to a lepto-hadronic origin of gamma-ray flux of blazars. Future gamma-ray observations above tens of TeV and below 100 MeV in energy will be crucial to test and discriminate between models.

A crucial test of the phantom closed cosmological model

ABSTRACT We suggest a crucial direct-observational test for measuring distinction between the standard ΛCDM model and recently proposed phantom dark energy positive curvature cosmological model. The test is based on the fundamental distance–flux–redshift relation for general Friedmann models. It does not depend on the CMBR data, on the large-scale structure growth models, and also on the value of the Hubble constant H0. Our crucial test can be performed by future gamma-ray burst observations with THESEUS space mission and by using gravitational-wave standard siren observations with modern advanced LIGO–Virgo and also forthcoming LISA detectors.

The low rate of Galactic pevatrons

Although supernova remnants remain the main suspects as sources of Galactic cosmic rays up to the knee, the supernova paradigm still has many loose ends. The weakest point in this construction is the possibility that individual supernova remnants can accelerate particles to the rigidity of the knee, ~ 106 GV. This scenario heavily relies upon the possibility to excite current driven non-resonant hybrid modes while the shock is still at the beginning of the Sedov phase. These modes can enhance the rate of particle scattering thereby leading to potentially very–high maximum energies. Here we calculate the spectrum of particles released into the interstellar medium from the remnants of different types of supernovae. We find that only the remnants of very powerful, rare core–collapse supernova explosions can accelerate light elements such as hydrogen and helium nuclei, to the knee rigidity, and that the local spectrum of cosmic rays directly constrains the rate of such events, if they are also source of PeV cosmic rays. On the other hand, for remnants of typical core–collapse supernova explosions, the Sedov phase is reached at late times, when the maximum energy is too low and the spectrum at very–high energies is very steep, being mostly produced during the ejecta dominated phase. For typical thermonuclear explosions, resulting in type Ia supernovae, we confirm previous findings that these objects can only produce cosmic rays up to ≲ 105 GeV. The implications for the overall cosmic ray spectrum observed at the Earth and for the detection of PeVatrons by future gamma–ray observatories are discussed.

Global fit of pseudo-Nambu-Goldstone Dark Matter

We perform a global fit within the pseudo-Nambu-Goldstone Dark Matter (DM) model emerging from an additional complex scalar singlet with a softly broken global U (1) symmetry. Leading to a momentum-suppressed DM-nucleon cross section at tree level, the model provides a natural explanation for the null results from direct detection experiments. Our global fit combines constraints from perturbative unitarity, DM relic abundance, Higgs invisible decay, electroweak precision observables and latest Higgs searches at colliders. The results are presented in both frequentist and Bayesian statisical frameworks. Furthermore, post-processing our samples, we include the likelihood from gamma-ray observations of Fermi -LAT dwarf spheroidal galaxies and compute the one-loop DM-nucleon cross section. We find two favoured regions characterised by their dominant annihilation channel: the Higgs funnel and annihilation into Higgs pairs. Both are compatible with current Fermi -LAT observations, and furthermore, can fit the slight excess observed in four dwarfs in a mass range between about 30–300 GeV. While the former region is hard to probe experimentally, the latter can partly be tested by current observations of cosmic-ray antiprotons as well as future gamma-ray observations.

Searching for dark matter in the Galactic halo with a wide field of view TeV gamma-ray observatory in the Southern Hemisphere

Despite mounting evidence that dark matter (DM) exists in the Universe, its fundamental nature remains unknown. We present sensitivity estimates to detect DM particles with a future very-high-energy (≳ TeV) wide field-of-view gamma-ray observatory in the Southern Hemisphere. This observatory would search for gamma rays from the annihilation or decay of DM particles in the Galactic halo. With a wide field of view, both the Galactic Center and a large fraction of the Galactic halo will be detectable with unprecedented sensitivity to DM in the mass range of ∼500 GeV to ∼2 PeV . These results, combined with those from other present and future gamma-ray observatories, will likely probe the thermal relic annihilation cross section of Weakly Interacting Massive Particles for all masses from ∼80 TeV down to the GeV range in most annihilation channels.

Exacerbating the Cosmological Constant Problem with Interacting Dark Energy Models

Future cosmological surveys will probe the expansion history of the Universe and constrain phenomenological models of dark energy. Such models do not address the fine-tuning problem of the vacuum energy, i.e., the cosmological constant problem (CCP), but can make it spectacularly worse. We show that this is the case for "interacting dark energy" models in which the masses of the dark matter states depend on the dark energy sector. If realized in nature, these models have far-reaching implications for proposed solutions to the CCP that require the number of vacua to exceed the fine-tuning of the vacuum energy density. We show that current estimates of the number of flux vacua in string theory, N_{vac}∼O(10^{272 000}), are far too small to realize certain simple models of interacting dark energy and solve the cosmological constant problem anthropically. These models admit distinctive observational signatures that can be targeted by future gamma-ray observatories, hence making it possible to observationally rule out the anthropic solution to the cosmological constant problem in theories with a finite number of vacua.

Constraining high-energy cosmic neutrino sources: Implications and prospects

We consider limits on the local ($z=0$) density ($n_0$) of extragalactic neutrino sources set by the nondetection of steady high-energy neutrino sources producing $\gtrsim50$ TeV muon multiplets in the present IceCube data, taking into account the redshift evolution, luminosity function and neutrino spectrum of the sources. We show that the lower limit depends moderately on source spectra and strongly on redshift evolution. We find $n_0\gtrsim{10}^{-8}-{10}^{-7}~{\rm Mpc}^{-3}$ for standard candle sources evolving rapidly, $n_s\propto{(1+z)}^3$, and $n_0\gtrsim{10}^{-6}-{10}^{-5}~{\rm Mpc}^{-3}$ for nonevolving sources. The corresponding upper limits on their neutrino luminosity are $L_{{\nu_\mu}}^{\rm eff}\lesssim10^{42}-10^{43}~{\rm erg}~{\rm s}^{-1}$ and $L_{{\nu_\mu}}^{\rm eff}\lesssim10^{41}-10^{42}~{\rm erg}~{\rm s}^{-1}$, respectively. Applying these results to a wide range of classes of potential sources, we show that powerful blazar jets associated with active galactic nuclei are unlikely to be the dominant sources. For almost all other steady candidate source classes (including starbursts, radio galaxies, and galaxy clusters and groups), an order of magnitude increase in the detector sensitivity at $\sim0.1-1$ PeV will enable a detection (as point sources) of the few brightest objects. Such an increase, which may be provided by next-generation detectors like IceCube-Gen2 and an upgraded KM3NET, can improve the limit on $n_0$ by more than two orders of magnitude. Future gamma-ray observations (by Fermi, HAWC and CTA) will play a key role in confirming the association of the neutrinos with their sources.

In quest of non-thermal signatures in early-type stars

A reduced fraction of luminous, early-type stars in binaries has provided some of the most interesting sources in modern high-energy astrophysics. A fingerprint of the capability of these systems to accelerate particles up to TeV energies is the associated detection of non-thermal, synchrotron emission often in the radio domain. Here we aim to identify new early-type, luminous stars where energetic, non-thermal processes are at work to enable future comparative studies based on an extended sample. Moreover, these objects also appear as very interesting targets for future gamma-ray observatories such as the Cherenkov Telescope Array. We have designed a methodology to search for new examples of these interesting sources in order to enlarge the extremely reduced population currently known. Our search procedure is described in this paper, together with a practical application using public databases and catalogues currently available (Luminous Stars in the Northern Milky Way, NRAO VLA Sky Survey, and Westerbork Northern Sky Survey). Optical and radio interferometric follow-up observations of selected candidate stars were conducted to better assess their properties. Although no new discoveries of this kind have been achieved yet, the observational data analyzed in this work does provide some interesting side results, such as characterizing the spectroscopic and photometric properties of two early-type stars originally selected as candidates to be explored, namely TYC4051-1277-1 and TYC3594-2269-1. Unexpectedly, this last one was also found to be a new DAO-type white dwarf star instead of a non-degenerate star. In addition, both stars seem to display different optical periods based on our photometric monitoring.

Singlet Majorana fermion dark matter: a comprehensive analysis in effective field theory

We explore a singlet Majorana fermion dark matter candidate using an effective field theory (EFT) framework, respecting the relations imposed by the standard model $SU(3)_C \times SU(2)_L \times U(1)_Y$ gauge invariance among different couplings. All operators of dimension-5 and dimension-6, forming a complete basis, are taken into account at the same time, keeping in view ultraviolet completions which can give rise to more than one operator at a time. If in addition CP-conservation is assumed, the remaining parameter space, where an EFT description is valid, is found to be rather restricted after imposing constraints from relic abundance, direct, indirect and collider searches. On including the CP-violating dimension-5 operator, $(\overline{\chi}i \gamma_5 \chi) (H^\dagger H)$, a significantly larger parameter space opens up. We use the profile likelihood method to map out the remaining landscape of such a DM scenario. The reach of future searches using ton-scale direct detection experiments, an $e^+ e^-$ collider like the proposed ILC and limits from future gamma-ray observations are also estimated.

Dark matter profiles and annihilation in dwarf spheroidal galaxies: prospectives for present and future γ-ray observatories - I. The classical dwarf spheroidal galaxies

Due to their large dynamical mass-to-light ratios, dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter (DM) in gamma-rays. We examine their detectability by present and future gamma-ray observatories. The key innovative features of our analysis are: (i) We take into account the angular size of the dSphs; while nearby objects have higher gamma ray flux, their larger angular extent can make them less attractive targets for background-dominated instruments. (ii) We derive DM profiles and the astrophysical J-factor (which parameterises the expected gamma-ray flux, independently of the choice of DM particle model) for the classical dSphs directly from photometric and kinematic data. We assume very little about the DM profile, modelling this as a smooth split-power law distribution, with and without sub-clumps. (iii) We use a Markov Chain Monte Carlo (MCMC) technique to marginalise over unknown parameters and determine the sensitivity of our derived J-factors to both model and measurement uncertainties. (iv) We use simulated DM profiles to demonstrate that our J-factor determinations recover the correct solution within our quoted uncertainties. Our key findings are: (i) Sub-clumps in the dSphs do not usefully boost the signal; (ii) The sensitivity of atmospheric Cherenkov telescopes to dSphs within 20 kpc with cored halos can be up to ~50 times worse than when estimated assuming them to be point-like. Even for the satellite-borne Fermi-LAT the sensitivity is significantly degraded on the relevant angular scales for long exposures, hence it is vital to consider the angular extent of the dSphs when selecting targets; (iii) No DM profile has been ruled out by current data, but using a prior on the inner dark matter cusp slope 0<=gamma<=1 provides J-factor estimates accurate to a factor of a few if an appropriate angular scale [abridged]

HIGH ENERGY GAMMA-RAY ABSORPTION AND CASCADE EMISSION IN NEARBY STARBURST GALAXIES

High energy gamma-ray emission from two nearby bright starburst galaxies, M82 and NGC 253, have recently been detected by Fermi, H.E.S.S., and VERITAS. Since starburst galaxies have a high star formation rate and plenty of dust in the central starburst region, infrared emissions are strong there. Gamma-ray photons are absorbed by the interstellar radiation field photons via electron and positron pair creation. The generated electron and positron pairs up scatter the interstellar photons to very high energy gamma-ray photons via cascade emission through inverse Compton scattering. In this paper, we evaluate the contribution of this cascade emission to the gamma-ray spectra of M82 and NGC 253. Although it would be difficult to see direct gamma- ray evidence of cosmic-rays with an energy > 10 TeV due to the gamma-ray attenuation, the resulting cascade emission would be indirect evidence. By including the cascade component, we find that the total flux above 1 TeV increases ~18% and ~45% compared with the absorbed flux assuming the maximum kinetic proton energy as 45.3 TeV and 512 TeV, respectively. Future gamma-ray observatories such as CTA would be able to see the indirect evidence of cosmic-ray with an energy > 10 TeV by comparing with theoretical emission models including this cascade effect.

Counterpart candidates to the unidentified Fermi source 0FGL J1848.6-0138

Aims. We aim to contribute to the identification of the counterpart for one of the bright sources of gamma-rays in the catalogue obtained and released by the Fermi collaboration.Methods. Our work is based on a extensive identification of sources from different wavelength catalogues and databases.Results. As a first result, we report the finding of a few counterpart candidates inside the 95% confidence error box of the Fermi LAT unidentified gamma-ray source 0FGL J1848.6-0138. The globular cluster GLIMPSE-C01 is remarkably distinctive being among the most peculiar objects consistent with both the position uncertainty in the gamma-ray source and a conceivable physical scenario for gamma-ray production. The Fermi-observed spectrum is compared with theoretical predictions in the literature and the association is found to be plausible but not yet certain because of its low X-ray to gamma-ray luminosity ratio. Other competing counterparts are also discussed. In particular, we pay special attention to a possible Pulsar Wind Nebula inside the Fermi error box, whose nature is yet to be confirmed. Conclusions. Both a globular cluster and an infrared source resembling a Pulsar Wind Nebula were found to be in positional agreement with 0FGL J1848.6-0138. In addition, other interesting objects in the field are also reported. Future gamma-ray observations will reduce the position uncertainty and we hope eventually confirm one of the counterpart candidates reported here. If GLIMPSE-C01 is confirmed together with the possible Fermi detection of the well known globular cluster 47 Tuc, then this would provide strong support to theoretical predictions that globular clusters are possible gamma-ray sources.

ON THE POSSIBLE ASSOCIATION OF ULTRA HIGH ENERGY COSMIC RAYS WITH NEARBY ACTIVE GALAXIES

(Abridged) Data collected by the Pierre Auger Observatory (Auger) provide evidence for anisotropy in the arrival directions of cosmic rays (CRs) with energies >57 EeV that suggests a correlation with the positions of AGN located within ~75 Mpc. A detailed study of the sample of AGN whose positions are located within 3.2 degrees of the CR events (and extending our analysis out to ~150 Mpc) shows that most of them are classified as Seyfert 2 and low-ionization nuclear emission-line region (LINER) galaxies whose properties do not differ substantially from other local AGN of the same types. Therefore, if the production of the highest energy CRs is persistent in nature, i.e., operates in a single object on long (>Myr) timescales, the claimed correlation between the CR events observed by Auger and local active galaxies should be considered as resulting from a chance coincidence. Additionally, most of the selected sources do not show significant jet activity, and hence, in most conservative scenarios, there are no reasons for expecting them to accelerate CRs up to the highest energies, ~10^20 eV. A future analysis has to take into account AGN morphology and may yield a correlation with a larger deflection angle and/or more distant sources. We further argue that the nearby radio galaxy NGC 5128 (Cen A) alone could be associated with at least 4 events due to its large radio extent, and PKS 1343-60 (Cen B), another nearby radio galaxy, can be associated with more than 1 event due to its proximity to the Galactic plane and, correspondingly, the stronger Galactic magnetic field the UHECRs encounter during propagation to the Earth. Future gamma-ray observations (by, e.g., Fermi Gamma-ray Space Telescope, and HESS) may provide additional clues to the nature of the accelerators of the UHECRs in the local Universe.

Gamma-ray probe of cosmic ray pressure in galaxy clusters and cosmological implications

Cosmic rays produced in cluster accretion and merger shocks provide pressure to the intracluster medium (ICM) and affect the mass estimates of galaxy clusters. Although direct evidence for cosmic-ray ions in the ICM is still lacking, they produce gamma-ray emission through the decay of neutral pions produced in their collisions with ICM nucleons. We investigate the capability of the Gamma-ray Large Area Space Telescope (GLAST) and imaging atmospheric Cerenkov telescopes (IACTs) for constraining the cosmic-ray pressure contribution to the ICM. We show that GLAST can be used to place stringent upper limits, a few per cent for individual nearby rich clusters, on the ratio of pressures of the cosmic rays and thermal gas. We further show that it is possible to place tight (<~10%) constraints for distant (z <~ 0.25) clusters in the case of hard spectrum, by stacking signals from samples of known clusters. The GLAST limits could be made more precise with the constraint on the cosmic-ray spectrum potentially provided by IACTs. Future gamma-ray observations of clusters can constrain the evolution of cosmic-ray energy density, which would have important implications for cosmological tests with upcoming X-ray and Sunyaev-Zel'dovich effect cluster surveys.

The Guaranteed Gamma-Ray Background

The diffuse extragalactic gamma-ray background (EGRB) above 100 MeV encodes unique information about high-energy processes in the universe. Numerous sources for the EGRB have been proposed, but the two systems which are certain to make some contribution are active galaxies (blazars) as well as normal galaxies. In this paper, we evaluate the contribution to the background from both sources. The active galaxy contribution arises from unresolved blazars. We compute this contribution using the Stecker-Salamon model. For normal galaxies, the emission is due to cosmic-ray interactions with diffuse gas. Our key assumption here is that the cosmic-ray flux in a galaxy is proportional to the supernova rate and thus the massive star formation rate, quantified observationally by the cosmic star formation rate (CSFR). In addition, the existence of stars today requires a considerably higher ISM mass in the past. Using the CSFR to compute both these effects, we find that normal galaxies are responsible for a significant portion (\sim 1/3) of the EGRB near 1 GeV, but make a smaller contribution at other energies. Finally, we present a minimal two-component model which includes contributions from both normal galaxies and blazars. We show that the spectrum of the diffuse radiation is a key constraint on this model: while neither the blazar spectra, nor the galactic spectra, are separately optimal fits to the observed spectrum, the combined emission provides an excellent fit. We close by noting key observational tests of this two-component model, which can be probed by future gamma-ray observatories such as GLAST.