TeV γ-ray fluxes from the long campaigns on Mrk 421 as constraints on the emission of TeV–PeV neutrinos and UHECRs

Long TeV γ-ray campaigns have been carried out to study the spectrum, variability and duty cycle of the BL Lac object Markarian 421. These campaigns have given some evidence of the presence of protons in the jet: (i) Its spectral energy distribution which shows two main peaks; one at low energies (∼1 keV) and the other at high energies (hundreds of GeV), has been described by using synchrotron proton blazar model. (ii) The study of the variability at GeV γ-rays and X-rays has indicated no significant correlation. (iii) TeV γ-ray detections without activity in X-rays, called “orphan flares” have been observed in this object.Recently, The Telescope Array Collaboration reported the arrival of 72 ultra-high-energy cosmic rays with some of them possibly related to the direction of Markarian 421. The IceCube Collaboration reported the detection of 37 extraterrestrial neutrinos in the TeV–PeV energy range collected during three consecutive years. In particular, no neutrino track events were associated with this source. In this paper, we consider the proton–photon interactions to correlate the TeV γ-ray fluxes reported by long campaigns with the neutrino and ultra-high-energy cosmic ray observations around this blazar. Considering the results reported by The IceCube and Telescope Array Collaborations, we found that only from ∼25% to 70% of TeV fluxes described with a power law function with exponential cutoff can come from the proton–photon interactions.     

Gamma rays from colliding winds of massive stars

Colliding winds of massive binaries have long been considered as potential sites of non-thermal high-energy photon production. This is motivated by the detection of non-thermal spectra in the radio band, as well as by correlation studies of yet unidentified EGRET γ-ray sources with source populations appearing in star formation regions. This work re-considers the basic radiative processes and its properties that lead to high energy photon production in long-period massive star systems. We show that Klein–Nishina effects as well as the anisotropic nature of the inverse Compton scattering, the dominating leptonic emission process, likely yield spectral and variability signatures in the γ-ray domain at or above the sensitivity of current or upcoming gamma ray instruments like GLAST-LAT. In addition to all relevant radiative losses, we include propagation (such as convection in the stellar wind) as well as photon absorption effects, which a priori can not be neglected. The calculations are applied to WR 140 and WR 147, and predictions for their detectability in the γ-ray regime are provided. Physically similar specimen of their kind like WR 146, WR 137, WR 138, WR 112 and WR 125 may be regarded as candidate sources at GeV energies for near-future γ-ray experiments. Finally, we discuss several aspects relevant for eventually identifying this source class as a γ-ray emitting population. Thereby we utilize our findings on the expected radiative behavior of typical colliding wind binaries in the γ-ray regime as well as its expected spatial distribution on the γ-ray sky.     

Discovery of GeV gamma-ray emission from the LMC B0443-6657 with the <Emphasis Type="Italic">Fermi</Emphasis> Large Area Telescope

We report the discovery of gamma-ray detection from the Large Magellanic Cloud (LMC) B0443-6657 using the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. LMC B0443-6657 is a flat-spectrum radio source, possibly associated with a supernova remnant in the Large Magellanic Cloud (LMC N4). Employing the LAT data of 8 years, our results show a significant excess (\(>9.4\sigma \)) of gamma rays in the range of 0.2–100 GeV above the gamma-ray background. A power-law function is found to adequately describe the 0.2–\(100\mbox{ GeV}\)\(\gamma \)-ray spectrum, which yields a photon flux of \(3.27\pm 0.53\ \text{photon}\,\mbox{cm}^{2}\,\mbox{s}^{-1}\) with a photon index of \(2.35\pm 0.11\), corresponding to an isotropic gamma-ray luminosity of \(5.3\times 10^{40}~\mbox{erg}\,\mbox{s}^{-1}\). The hadronic model predicts a low X-ray and TeV flux while the leptonic model predicts an observable flux in these two energy bands. The follow-up observations of the LMC B0443-6657 in X-ray or TeV band would distinguish the radiation models of gamma rays from this region.     

Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

Minimal spanning tree algorithm for γ-ray source detection in sparse photon images: cluster parameters and selection strategies

The minimal spanning tree (MST) algorithm is a graph-theoretical cluster-finding method. We previously applied it to γ-ray bidimensional images, showing that it is quite sensitive in finding faint sources. Possible sources are associated with the regions where the photon arrival directions clusterize. MST selects clusters starting from a particular “tree” connecting all the point of the image and performing a cut based on the angular distance between photons, with a number of events higher than a given threshold. In this paper, we show how a further filtering, based on some parameters linked to the cluster properties, can be applied to reduce spurious detections. We find that the most efficient parameter for this secondary selection is the magnitude M of a cluster, defined as the product of its number of events by its clustering degree. We test the sensitivity of the method by means of simulated and real Fermi-Large Area Telescope (LAT) fields. Our results show that $\sqrt{M}$ is strongly correlated with other statistical significance parameters, derived from a wavelet based algorithm and maximum likelihood (ML) analysis, and that it can be used as a good estimator of statistical significance of MST detections. We apply the method to a 2-year LAT image at energies higher than 3 GeV, and we show the presence of new clusters, likely associated with BL Lac objects.     

Simulation studies of MACE-I: Trigger rates and energy thresholds

The MACE (Major Atmospheric Cherenkov Experiment) is an upcoming Very High Energy (VHE) γ-ray telescope, based on imaging atmospheric Cherenkov technique, being installed at Hanle, a high altitude astronomical site in Ladakh, India. Here we present Monte Carlo simulation studies of trigger rates and threshold energies of MACE in the zenith angle range of 0°–60° for on-axis γ-ray coming from point source and various cosmic ray species. We have simulated the telescope’s response to γ-rays, proton, electron and alpha initiated atmospheric Extensive Air Showers (EAS) in the broad energy range of 5 GeV to 20 TeV. For γ-rays we consider power law and log parabolic spectra while other particles are simulated with their respective cosmic ray spectrum. Trigger rates and threshold energies are estimated for the trigger configuration of 4 Close Cluster Nearest Neighbour(CCNN) pixels as implemented in MACE hardware, in combination with single channel discriminator threshold ranging from 6–10 photo electrons (pe). We find that MACE can achieve the γ-ray trigger energy threshold of ∼ 17 GeV (4 CCNN, 9 pe) at 0° zenith angle for power law spectrum. The total trigger rate at 0° zenith is expected to be ∼650 Hz, with protons contributing ∼ 80% to it. For the zenith range of 0°-40° we find that the telescope can achieve γ-ray trigger threshold energies of ∼22 GeV at 20° zenith angle and ∼40 GeV at 40° zenith angle. Integral rates are also almost constant for this zenith angle range. At zenith angle of 60°, trigger energy threshold increases to ∼173 GeV and total integral rate falls down to ∼305 Hz.     

The cascading of cosmic-ray nuclei in various media

More than a thousand interactions of primary heavy nuclei of the cosmic radiation with charge ≥10 and energy >1 GeV nucleon^?1 in nuclear emulsion have been studied with emphasis on how the primary nucleus fragments. It has been determined that the cases of multiple successive fragmentations that have been observed do not occur more frequently than expected. The fragmentation ofZ>20 nuclei does depend on the target nucleus to some extent so it is important to try to separate the interactions in emulsion by theirN _h (number of evaporation prongs). The fragmentation of ₈ ¹⁶ O at 2.1 GeV nucleon^?1 measured at the Bevalac shows a similar dependence on target nucleus. By using data from these new interactions combined with published data we have simulated on a computer nuclear cascades in both emulsion and air. Results on these cascades are given for both primary silicon and primary iron nuclei. These results are used to discuss the fluctuations expected in extensive air showers produced by heavy primary nuclei.     

Cosmic rays from binary neutron stars

Within the more than 30 yr of cosmic ray astrophysics, neither their origin nor their precise mode of propagation have found undisputable explanations. Among the favoured boosters have been point sources, like supernovae and pulsars, as well as extended sources, like cosmic clouds and supernova remnants. Extended sources have been proposed by Fermi (1949), and pushed more recently by a number of investigators because of the huge available reservoirs, and because repetitive shock acceleration can generate power law spectra which are similar to the ones observed (Axfordet al., 1977; Bell, 1978; Blandford and Ostriker, 1978; Krymsky, 1977). Yet the shock acceleration model cannot easily be adjusted to achieve particle energies in excess of some critical energy, of order 10^4±1 GeV (Völket al., 1981). For this and several other reasons, the suggestion is revived that neutron stars are the dominant source of high-energy cosmic rays. To be more precise: the (relativistic) ionic component of the cosmic rays is argued to be injected by young binary neutron stars (?10⁵ yr) whose rotating magnetospheres act like grindstones in the wind of their companion (Kundt, 1976). The high-energy (?30 GeV) electron-positron component may be generated by young pulsars (?10⁵ yr) and by collision processes, and the electron component below 30 GeV predominantly by supernova remnants.     

Milagro observations of potential TeV emitters

This paper reports the results from three targeted searches of Milagro TeV sky maps: two extragalactic point source lists and one pulsar source list. The first extragalactic candidate list consists of 709 candidates selected from the Fermi-LAT 2FGL catalog. The second extragalactic candidate list contains 31 candidates selected from the TeVCat source catalog that have been detected by imaging atmospheric Cherenkov telescopes (IACTs). In both extragalactic candidate lists Mkn 421 was the only source detected by Milagro. This paper presents the Milagro TeV flux for Mkn 421 and flux limits for the brighter Fermi-LAT extragalactic sources and for all TeVCat candidates. The pulsar list extends a previously published Milagro targeted search for Galactic sources. With the 32 new gamma-ray pulsars identified in 2FGL, the number of pulsars that are studied by both Fermi-LAT and Milagro is increased to 52. In this sample, we find that the probability of Milagro detecting a TeV emission coincident with a pulsar increases with the GeV flux observed by the Fermi-LAT in the energy range from 0.1 GeV to 100 GeV.     

Bounding the time delay between high-energy neutrinos and gravitational-wave transients from gamma-ray bursts

We derive a conservative coincidence time window for joint searches of gravitational-wave (GW) transients and high-energy neutrinos (HENs, with energies ?100 GeV), emitted by gamma-ray bursts (GRBs). The last are among the most interesting astrophysical sources for coincident detections with current and near-future detectors. We take into account a broad range of emission mechanisms. We take the upper limit of GRB durations as the 95% quantile of the T₉₀’s of GRBs observed by BATSE, obtaining a GRB duration upper limit of ∼150 s. Using published results on high-energy (>100 MeV) photon light curves for 8 GRBs detected by Fermi LAT, we verify that most high-energy photons are expected to be observed within the first ∼150 s of the GRB. Taking into account the breakout-time of the relativistic jet produced by the central engine, we allow GW and HEN emission to begin up to 100 s before the onset of observable gamma photon production. Using published precursor time differences, we calculate a time upper bound for precursor activity, obtaining that 95% of precursors occur within ∼250 s prior to the onset of the GRB. Taking the above different processes into account, we arrive at a time window of t_HEN − t_GW ∈ [−500 s, +500 s]. Considering the above processes, an upper bound can also be determined for the expected time window of GW and/or HEN signals coincident with a detected GRB, t_GW − t_GRB ≈ t_HEN − t_GRB ∈ [−350 s, +150 s]. These upper bounds can be used to limit the coincidence time window in multimessenger searches, as well as aiding the interpretation of the times of arrival of measured signals.     

A local interstellar spectrum for galactic electrons

A heliopause spectrum at 122 AU from the Sun is presented for galactic electrons over an energy range from 1 MeV to 50 GeV that can be considered the lowest possible local interstellar spectrum (LIS). The focus of this work is on the spectral shape of the LIS below ∼1.0 GeV. The study is done by using a comprehensive numerical model for solar modulation in comparison with Voyager 1 observations at ∼112 AU from the Sun and PAMELA data at Earth. Below ∼1.0 GeV, this LIS exhibits a power law with E^{−(1.55 ± 0.05)}, where E is the kinetic energy of these electrons. However, reproducing the PAMELA electron spectrum averaged for 2009, requires a LIS with a different power law of the form E^{−(3.15 ± 0.05)} above ∼5 GeV. Combining the two power laws with a smooth transition from low to high energies yields a LIS over the full energy range that is relevant and applicable to the modulation of cosmic ray electrons in the heliosphere. The break occurs between ∼800 MeV and ∼2 GeV as a characteristic feature of this LIS. The power-law form below ∼1 GeV produces a challenge to the origin of these low energy galactic electrons. On the other hand, the results of this study can be used as a gauge for astrophysical modeling of the local interstellar spectrum for electrons.     

Joint Effects of Compton Backscattering and Low-Energy Cutoff on the Flattening of Solar Hard X-Ray Spectra at Lower Energies

Theoretical calculation has shown that the spectrum of the Compton backscattering component in solar hard X-ray flares has a peak around 30 keV for a primary power-law source. Thus the superposition of the Compton backscattering component could cause a photon spectrum received at the Earth to be flattened below the peak energy and steeper above the peak energy. On the other hand, because a thick-target bremsstrahlung photon with a given energy E only could be produced by a nonthermal electron with an energy larger than E, thus if a power-law electron spectrum is cutoff below an energy E _c, then the produced photon spectrum will become flattened below E _c. In this work we present a calculation of the joint effects of the Compton backscattering and the low-energy cutoff on the spectral characteristics of the received solar hard X-ray in the energy range 10–100 keV. The results show that the flattening caused purely by the Compton backscattering could be comparable with that by the low-energy cutoff for hard spectra. So, it is obvious that the joint effects of the low-energy cutoff and the Compton backscattering could result in the received photon spectra to be much more flattened at lower energies. On the other hand, compared to the primary photon spectrum, the received photon spectral index will increase about 0.15 due to the Compton backscattering at higher energy, which seems independent of the primary spectral index.     

High energy γ-ray emission from compact galactic sources in the context of observations with the next generation Cherenkov Telescope Arrays

The observational progress in the γ-ray astronomy in the last few years has led to the discovery of more than a thousand sources at GeV energies and more than a hundred sources at TeV energies. A few different classes of compact objects in the Galaxy have been established. They show many unexpected features at high energies the physics of which remains mainly unknown. At present it is clear that detailed investigation of these new phenomena can be performed only with the technical equipment which offer an order of magnitude better sensitivity, and a few times better energy, angular and time resolution in the broad energy range staring from a few tens of GeV up to a few hundreds TeV. Such facilities can be realized by the next generation of instruments such as the planned Cherenkov Telescope Array (CTA).The aim of this report is to summarize up to date observational results on the compact galactic sources in the GeV–TeV γ-ray energy range, discuss their theoretical implications, and indicate which hypothesis considered at present might be verified with the next generation of telescopes. We point out which of the observational features of the γ-ray sources are important to investigate with special care with the planned CTA in order to throw new light on physical processes involved. Their knowledge should finally allow us to answer the question on the origin of energetic particles in our Galaxy.     

Phenomenological explanation for the shape of the spectrum of cosmic rays with energies <Emphasis Type="Italic">E</Emphasis> > 10 GeV

Assuming that the energy gain by cosmic-ray (CR) particles is a stochastic process with stationary increments, we derive expressions for the shape of their energy spectrum up to energies E ～ 10¹⁸ eV. In the ultrarelativistic case under study, the energy is proportional to the momentum, whose time derivative is the force. According to the Fermi mechanism, a particle accelerates when it passes through a system of shock waves produced by supernova explosions. Since these random forces act on time scales much shorter than the particle lifetime, we assume them to be delta-correlated in time. In this case, due to the linear energy-momentum relationship, the mean square of the energy (increments) is proportional to the differential scale τ(E) ～ Eτ(≥E), where τ (≥E) is the cumulative time it takes for a particle to gain an energy ≥E. The probability of finding a particle with energy ≥E somewhere in the system is inversely proportional to the time it takes to gain the energy E. To estimate an upper limit for the space number density of CR particles, we use estimates of the CR volume energy density, a quantity known for our Galaxy. It is taken to be constant in the range 10 GeV ≤ E ≤ 3 × 10⁶ GeV, where the index of the energy spectrum was found to be ?8/3 ≈ ?2.67 against its empirical value of ?2.7. In the range 3 × 10⁶ GeV ≤ E < 10⁹ GeV, the upper limit for the volume energy density is estimated by using the results from the previous range to be ?28/9 ≈ ?3.11 against its empirical value of ?3.1. The numerical coefficients in the suggested shapes of the spectrum can be determined by comparison with observational data. Thus, the CR energy spectrumis the result of the random walks of ultrarelativistic particles in energy/momentum space caused by the Fermi mechanism.     

Crab pulsar spectrum: A non-extensive approach

Non extensive statistical physics has been applied to various problems in physics including astrophysics. In this paper we explore the possibility of using non-extensive approach to explain the recently observed pulsed γ-ray from Crab pulsar above 100 GeV observed by VERITAS γ-ray telescope.     

A classification of the available astrophysica data of particularHii regions

We have calculated the energy spectra of cosmic ray secondary antiprotons and positrons using the latest available data on inclusive reactions. Using the measured positron spectrum, we have found that the amount of matter traversed by the cosmic rays in the few GeV region to bem≈4.7±1.5 g cm^?2 of interstellar hydrogen. The computed antiproton to proton ratio is about 4×10^?4 for energies 5–10 GeV. This is sufficient to make observations of antiprotons feasible from balloon flights. We have also pointed out the type of information that can be obtained if accurate information of the spectra of these two components becomes available.     

Gamma rays from the Galactic Centre region: A review

During the last decades, increasingly precise astronomical observations of the Galactic Centre (GC) region at radio, infrared, and X-ray wavelengths laid the foundations to a detailed understanding of the high energy astroparticle physics of this most remarkable location in the Galaxy. Recently, observations of this region in high energy (HE, 10 MeV–100 GeV) and very high energy (VHE, > 100 GeV) γ-rays added important insights to the emerging picture of the Galactic nucleus as a most violent and active region where acceleration of particles to very high energies – possibly up to a PeV – and their transport can be studied in great detail. Moreover, the inner Galaxy is believed to host large concentrations of dark matter (DM), and is therefore one of the prime targets for the indirect search for γ-rays from annihilating or decaying dark matter particles. In this article, the current understanding of the γ-ray emission emanating from the GC is summarised and the results of recent DM searches in HE and VHE γ-rays are reviewed.     

Perspective of monochromatic gamma-ray line detection with the High Energy cosmic-Radiation Detection (HERD) facility onboard China’s space station

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China’s space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ( ∼1% for electrons and photons and 20% for nuclei) and a large geometry factor (>3 m² sr for electrons and diffuse photons and > [2]m² sr for nuclei). In this work we discuss the capability of HERD to detect monochromatic γ-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic γ-ray searches at energies between ∼ 10 to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.     

Thermal relic abundance and anisotropy due to modified gravity

Alternative cosmologies, based on extensions of General Relativity, predict modified thermal histories in the early universe during the pre Big Bang Nucleosynthesis (BBN) era. When the expansion rate is enhanced with respect to the standard case, thermal relics typically decouple with larger relic abundances. In this paper, we study the dynamical evolution of an f(R) model of gravity in a homogeneous and anisotropic background which is given by a Bianchi type-I model of the universe filled with dark matter, which is described by a perfect fluid with a barotropic equation of state. As an example of a consistent analysis of modified gravity, we apply the formalism to a simple background solution of R+βR ⁿ gravity. Our analysis shows that f(R) cosmology allows dark matter masses lesser than 100 GeV, in the regime ρ ^c ?ρ ^m . We finally discuss how these limits apply to some specific realizations of standard cosmologies: an f(R) gravity model, Einstein frame model.     

Resonant Compton scattering associated with pair creation

Studies of Compton scattering by relativistic electrons in a strong magnetic field have been restricted to either incident photon angles θ′ aligned along the magnetic field B or incident photon energies ω′ below the first pair creation threshold $\omega'_{PC}$ . When these restrictions are relaxed there is a resonance in Compton scattering associated with pair creation (PC), that is analogous to but independent of known resonances associated with gyromagnetic absorption (GA). As with the GA resonances, that may be labeled by the Landau quantum numbers of the relevant states, there is a sequence of PC resonances where the scattering cross section diverges. In this paper, the lowest divergence is studied for incident photon energies satisfying ω′²sin² θ′/(2eB)?1, assuming that the scattering electron is in its ground (Landau) state. This lowest resonance affects only parallel-polarized photons.