Formation of hard very high energy gamma-ray spectra of blazars due to internal photon–photon absorption

The energy spectra of TeV gamma-rays from blazars, after being corrected for intergalatic absorption in the extragalactic background light (EBL), appear unusually hard, a fact that poses challenges to the conventional models of particle acceleration in TeV blazars and/or to the EBL models. In this paper, we show that the internal absorption of gamma-rays caused by interactions with dense narrow-band radiation fields in the vicinity of compact gamma-ray production regions can lead to the formation of gamma-ray spectra of an almost arbitrary hardness. This allows significant relaxation of the current tight constraints on particle acceleration and radiation models, although at the expense of enhanced requirements to the available non-thermal energy budget. The latter, however, is not a critical issue, as long as it can be largely compensated by the Doppler boosting, assuming large (>10) Doppler factors of the relativistically moving gamma-ray production regions. The suggested scenario of formation of hard gamma-ray spectra predicts detectable synchrotron radiation of secondary electron–positron pairs which might require a revision of the current 'standard paradigm' of spectral energy distributions of gamma-ray blazars. If the primary gamma-rays are of hadronic origin related to pp or   p γ  interactions, the 'internal gamma-ray absorption' model predicts neutrino fluxes close to the detection threshold of the next generation high-energy neutrino detectors.     

Observational limits on inverse Compton processes in gamma-ray bursts

Inverse Compton (IC) scattering is one of two viable mechanisms that can produce prompt non-thermal soft gamma-ray emission in gamma-ray bursts. IC requires low-energy seed photons and a population of relativistic electrons that upscatter them. The same electrons will upscatter the gamma-ray photons to even higher energies in the TeV range. Using the current upper limits on the prompt optical emission, we show that under general conservative assumption the IC mechanism suffers from an 'energy crisis'. Namely, IC will overproduce a very high energy component that would carry much more energy than the observed prompt gamma-rays, or alternatively it will require a low-energy seed that is more energetic than the prompt gamma-rays. Our analysis is general, and it makes no assumptions on the specific mechanism that produces the relativistic electron population.     

TeV neutrinos and gamma-rays from pulsars

Recent studies suggest that pulsars could be strong sources of TeV muon neutrinos provided positive ions are accelerated by pulsar polar caps to PeV energies. In such a situation, muon neutrinos are produced through the Δ-resonance in interactions of pulsar-accelerated ions with its thermal radiation field. High-energy gamma-rays should also be produced simultaneously in pulsar environment as both charged and neutral pions are generated in the interactions of energetic hadrons with the ambient photon fields. Here, we estimate TeV gamma-ray flux at the Earth from a few nearby young pulsars. When compared with the observations, we find that proper consideration of the effect of polar cap geometry in flux calculation is important. Incorporating such an effect, we obtain the (revised) event rates at the Earth due to a few potential nearby pulsars. The results suggest that pulsars are unlikely to be detected by the upcoming neutrino telescopes. We also estimate TeV gamma-ray and neutrino fluxes from pulsar nebulae for the adopted model of particle acceleration.     

Modulation of the high energy γ-ray flux from PSR B1259-63/SS2883 due to the orbital variation of the maximum energy of accelerated electrons

The inverse Compton (IC) scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is generally believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. In such a system the acceleration takes place in the presence of a dense radiation field provided by a companion Be2-type star. Thus it is natural to expect an orbital phase dependence of the acceleration efficiency in the system. The HESS collaboration reported the tendency of reduction of TeV γ-rays around the periastron. In this paper we study a possible explanation of this effect by the “early” (sub-TeV) cutoffs in the energy spectrum of accelerated electrons due to the enhanced rate of Compton losses close to the periastron.     

TeV gamma-rays from accreting magnetars in massive binaries

This paper focuses on neutron stars (NS) of the magnetar type inside massive binary systems. We determine the conditions under which the matter from the stellar wind can penetrate the inner magnetosphere of the magnetar. At a certain distance from the NS surface, the magnetic pressure can balance the gravitational pressure of the accreting matter, creating a very turbulent, magnetized transition region. It is suggested that this region provides good conditions for the acceleration of electrons to relativistic energies. These electrons lose energy due to the synchrotron process and inverse Compton (IC) scattering of the radiation from the nearby massive stellar companion, producing high-energy radiation from X-rays up to ∼TeV γ-rays. The primary γ-rays can be further absorbed in the stellar radiation field, developing an IC  e^±  pair cascade. We calculate the synchrotron X-ray emission from primary electrons and secondary  e^±  pairs and the IC γ-ray emission from the cascade process. It is shown that quasi-simultaneous observations of the TeV γ-ray binary system LSI +61 303 in the X-ray and TeV γ-ray energy ranges can be explained with such an accreting magnetar model.     

河外TeV能段γ射线背景下限的估算

Fermi卫星对GeV能段的河外伽马射线背景(Extragalactic Gamma-ray Background, EGB)进行了较为精确的测量, 极大提高了对高能伽马射线背景的认识, 但是在TeV能段, 使用空间探测器进行观测非常困难, 只能依赖地面伽马射线探测器, 如成像大气切伦科夫望远镜. 目前, 对于TeV能段的河外伽马射线背景的认识还不完善. 使用有低活跃状态能谱的61个TeV源(包含2个星暴星系、6个射电星系以及53个耀变体)的累计流量给出河外TeV伽马射线背景的下限. 结果显示, 低能段(0.5--4.5TeV)流量由两个临近的耀变体Mrk 421和Mrk 501主导, 贡献了大约58%的累计背景流量; 而大于4.5TeV的能段, 由3个已观测到10TeV以上能段流量的极端耀变体H 1426+428、1ES 1959+650以及1ES 0229+200主导. 最后分别探究了星暴星系、射电星系以及耀变体对河外TeV伽马射线背景的贡献, 不同耀变体子类对河外TeV伽马射线背景的贡献以及不同红移区间TeV源对河外伽马射线背景的贡献.     

The cross-correlation between galaxies and groups: probing the galaxy distribution in and around dark matter haloes

Some massive binaries should contain energetic pulsars which inject relativistic leptons from their inner magnetospheres and/or pulsar wind regions. If the binary system is compact enough, then these leptons can initiate inverse Compton (IC) e^± pair cascades in the anisotropic radiation field of a massive star. γ-rays can be produced in the IC cascade during its development in a pulsar wind region and above a shock in a massive star wind region where the propagation of leptons is determined by the structure of a magnetic field around the massive star. For a binary system with specific parameters, we calculate phase-dependent spectra and fluxes of γ-rays escaping as a function of the inclination angle of the system and for different assumptions on injection conditions of the primary leptons (their initial spectra and location of the shock inside the binary). We conclude that the features of γ-ray emission from such massive binaries containing energetic pulsars should allow us to obtain important information on the acceleration of particles by the pulsars, and on interactions of a compact object with the massive star wind. Predicted γ-ray light curves and spectra in the GeV and TeV energy ranges from such binary systems within our Galaxy and Magellanic Clouds should be observed by future AGILE and GLAST satellites and low-threshold Cherenkov telescopes, such as MAGIC, HESS, VERITAS or CANGAROO III.     

Lightning black holes as unidentified TeV sources

Imaging Atmospheric Cherenkov Telescopes have revealed more than 100 TeV sources along the galactic plane, around 45% of them remain unidentified. However, radio observations revealed that dense molecular clumps are associated with 67% of 18 unidentified TeV sources. In this paper, we propose that an electron–positron magnetospheric accelerator emits detectable TeV gamma-rays when a rapidly rotating black hole enters a gaseous cloud. Since the general-relativistic effect plays an essential role in this magnetospheric lepton accelerator scenario, the emissions take place in the direct vicinity of the event horizon, resulting in a point-like gamma-ray image. We demonstrate that their gamma-ray spectra have two peaks around 0.1 GeV and 0.1 TeV and that the accelerators become most luminous when the mass accretion rate becomes about 0.01% of the Eddington accretion rate. We compare the results with alternative scenarios such as the cosmic-ray hadron scenario, which predicts an extended morphology of the gamma-ray image with a single power-law photon spectrum from GeV to 100 TeV.     

Spectral features in gamma-rays expected from millisecond pulsars

In the advent of next generation gamma-ray missions, we present general properties of spectral features of high-energy emission above 1 MeV expected for a class of millisecond, low magnetic field (∼10⁹ G) pulsars. We extend polar-cap model calculations of Rudak & Dyks by including inverse Compton scattering events in an ambient field of thermal X-ray photons and by allowing for two models of particle acceleration. In the range between 1 MeV and a few hundred GeV, the main spectral component is the result of curvature radiation of primary particles. The synchrotron component arising from secondary pairs becomes dominant only below 1 MeV. The slope of the curvature radiation spectrum in the energy range from 100 MeV to 10 GeV strongly depends on the model of longitudinal acceleration, whereas below ∼100 MeV all slopes converge to a unique value of 4/3 (in a ν ℱ _ν convention). The thermal soft X-ray photons, which come either from the polar cap or from the surface, are Compton upscattered to a very high energy domain and form a separate spectral component peaking at ∼1 TeV. We discuss the observability of millisecond pulsars by future high‐energy instruments and present two rankings relevant for GLAST and MAGIC. We point to the pulsar J0437−4715 as a promising candidate for observations.     

TeV γ-rays from old supernova remnants

We study the emission from an old supernova remnant (SNR) with an age of around 10⁵ yr and that from a giant molecular cloud (GMC) encountered by the SNR. When the SNR age is around 10⁵ yr, proton acceleration is efficient enough to emit TeV γ-rays both at the shock of the SNR and that in the GMC. The maximum energy of primarily accelerated electrons is so small that TeV γ-rays and X-rays are dominated by hadronic processes,  π⁰  -decay and synchrotron radiation from secondary electrons, respectively. However, if the SNR is older than several 10⁵ yr, there are few high-energy particles emitting TeV γ-rays because of the energy-loss effect and/or the wave-damping effect occurring at low-velocity isothermal shocks. For old SNRs or SNR–GMC interacting systems capable of generating TeV γ-ray emitting particles, we calculated the ratio of TeV γ-ray (1–10 TeV) to X-ray (2–10 keV) energy flux and found that it can be more than  ∼10²  . Such a source showing large flux ratio may be a possible origin of recently discovered unidentified TeV sources.     

Acceleration of particles in the vicinity of a massive black hole

Recent results of the gamma-ray Cherenkov astronomy definitely prove the existence of fast variability in the very high energy (V.H.E.) gamma-ray flux of some active galactic nuclei. The BL Lac PKS 2155-304 for instance showed variations down to a few minutes time scale. From standard light travel time argument, these variations put extremely strong constraints on the size of the TeV emitting zone, which has to be of the order of a few Schwarzschild radius, even for high values of the relativistic Doppler factor of the emitting jets. Such discovery is a challenge for particle acceleration scenarios, which have to imagine efficient acceleration processes at work in a very compact zone. Eventually, the immediate vicinity of the central black hole appears as the most conservative choice for the location of the TeV emission region of active galactic nuclei. In this paper, we propose a two-step mechanism for charged particle acceleration in the magnetosphere of a massive black hole surrounded by an accretion disk. Particles first gain energy by a stochastic process during the accretion phase. It is shown that effective proton acceleration up to energies 10¹⁷–10¹⁹ eV is possible in a low-luminosity magnetized accretion disk with 2D turbulent motion. The distribution function of energetic protons over energies is a power law function with typical index ≃−1. Here electrons are not very efficiently accelerated because of their drastic losses by synchrotron radiation. In a second time, part of the fast particles escape from the disk and are then entrained by the magnetic structure above the disk, in the rotating black hole magnetosphere. They thus gain additional energy by direct centrifugal mechanism, up to about 10²⁰ eV for the protons and to 10–100 TeV for the electrons when they cross the light cylinder surface. Such energetic particles can further radiate in the TeV spectral range observed by Cherenkov experiments as HESS, MAGIC and VERITAS. Energetic protons can produce γ-radiation in the energy band 1 GeV–100 TeV and above mainly by nuclei collisions with the disk matter, clouds, or ambient low energy photons. Energetic electrons can also reach the required spectral range by inverse Compton emission. However their acceleration is less efficient due to heavy radiation losses, and only gained by centrifugal process during the second phase of the whole mechanism we describe. Our present analysis would therefore favor hadronic scenarios for TeV emission of active galactic nuclei. It is tempting to relate long term variability over years of TeV active galactic nuclei to the first stochastic acceleration phase, which also provides the needed power law particle distributions, while short term variability over minutes is more likely due to perturbations of the second fast direct acceleration phase.     

Gamma-Ray Emission from Pulsar Outer Magnetospheres

We investigate a stationary pair production cascade in the outer magnetosphere of an isolated, spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate gamma-rays via curvature and inverse-Compton processes. Some of such gamma-rays collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the energy distribution of particles and gamma-rays at each point along the field lines. By solving the set of Maxwell and Boltzmann equations, we demonstrate that an external injection of charged particles at nearly Goldreich-Julian rate does not quench the gap but shifts its position and that the particle energy distribution cannot be described by a power-law. The injected particles are accelerated in the gap and escape from it with large Lorentz factors. We show that such escaping particles migrating outside of the gap contribute significantly to the gamma-ray luminosity for young pulsars and that the soft gamma-ray spectrum between 100 MeV and 3 GeV observed for the Vela pulsar can be explained by this component. We also discuss that the luminosity of the gamma-rays emitted by the escaping particles is naturally proportional to the square root of the spin-down luminosity.     

One-dimensional electric field structure of an outer gap accelerator – III. Location of the gap and the gamma-ray spectrum

We investigate a stationary particle acceleration zone in the outer magnetosphere of an obliquely rotating neutron star. The charge depletion as a result of global current causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate curvature gamma-rays, some of which collide with the X-rays to materialize as pairs in the gap. As a result of this pair-production cascade, the replenished charges partially screen the electric field, which is self-consistently solved together with the distribution of particles and gamma-rays. If no current is injected at either of the boundaries of the accelerator, the gap is located around the so-called null surface, where the local Goldreich–Julian charge density vanishes. However, we find that the gap position shifts outwards (or inwards) when particles are injected at the inner (or outer) boundary. We apply the theory to the seven pulsars whose X-ray fields are known from observations. We show that the gap should be located near to or outside of the null surface for the Vela pulsar and PSR B1951+32, so that their expected GeV spectrum may be consistent with observations. We then demonstrate that the intrinsically large TeV flux from the outer gap of PSR B0540–69 is absorbed by the magnetospheric infrared photons, causing it to be undetectable. We also point out that the electrodynamic structure and the resultant GeV emission properties of millisecond pulsars are similar to young pulsars.     

Possible effects of pair echoes on gamma-ray burst afterglow emission

High-energy emission from gamma-ray bursts (GRBs) is widely expected but had been sparsely observed until recently when the Fermi satellite was launched. If >TeV gamma-rays are produced in GRBs and can escape from the emission region, they are attenuated by the cosmic infrared background photons, leading to regeneration of ∼GeV–TeV secondary photons via inverse-Compton scattering. This secondary emission can last for a longer time than the duration of GRBs, and it is called a pair echo. We investigate how this pair echo emission affects spectra and light curves of high-energy afterglows, considering not only prompt emission but also afterglow as the primary emission. Detection of pair echoes is possible as long as the intergalactic magnetic field (IGMF) in voids is weak. We find (1) that the pair echo from the primary afterglow emission can affect the observed high-energy emission in the afterglow phase after the jet break and (2) that the pair echo from the primary prompt emission can also be relevant, but only when significant energy is emitted in the TeV range, typically     . Even non-detections of the pair echoes could place interesting constraints on the strength of IGMF. The more favourable targets to detect pair echoes may be the 'naked' GRBs without conventional afterglow emission, although energetic naked GRBs would be rare. If the IGMF is weak enough, it is predicted that the GeV emission extends to >30–300 s.     

Inverse Compton e± pair cascade model for the γ-ray production in massive binary LSI +61° 303

We apply an inverse Compton   e ^±  pair cascade model for γ-ray production in the massive binary system LSI +61° 303 assuming that electrons are accelerated already inside the inner part of the jet launched by the compact object. γ-ray spectra, affected by the cascade process, and lower energy spectra, from the synchrotron cooling of the highest energy electrons in the jet, are calculated as a function of the phase of this binary system. γ-ray spectra expected in such a model have different shape than those ones produced by electrons in the jet directly to observer. Moreover, the model predicts clear anticorrelation between γ-ray fluxes in the GeV (1–10 GeV) and TeV (>200 GeV) energy ranges with the peak of the TeV emission at the phase ∼0.5 (the peak half-width ranges between the phases ∼0.4–0.9, for the inclination of the binary system equal to 60°, and ∼0.4–0.1 for 30°). The fine features of TeV γ-ray emission (fluxes and spectral shapes) as a function of the phase of the binary system are consistent with recent observations reported by the MAGIC collaboration. Future simultaneous observations in the GeV energies (by the GLAST and AGILE telescopes) and in the TeV energies (by the MAGIC and VERITAS telescopes) should test other predictions of the considered model supporting or disproving the hypothesis of acceleration of electrons already in the inner part of the microquasar jets.     

Tail emission of prompt gamma-ray burst jets

Tail emission of the prompt gamma-ray burst (GRB) is discussed using a multiple emitting sub-shell (inhomogeneous jet, sub-jets or mini-jets) model, where the whole GRB jet consists of many emitting sub-shells. One may expect that such a jet with angular inhomogeneity should produce spiky tail emission. However, we found that the tail is not spiky but is decaying roughly monotonically. The global decay slope of the tail is not so much affected by the local angular inhomogeneity but affected by the global sub-shell energy distribution. The fact that steepening GRB tail breaks appeared in some events prefers the structured jets. If the angular size of the emitting sub-shell is around 0.01–0.02 rad, some bumps or fluctuations appear in the tail emission observed frequently in long GRBs. If the parameter differences of sub-shell properties are large, the tail has frequent changes of the temporal slope observed in a few bursts. Therefore, the multiple emitting sub-shell model has the advantage of explaining the small-scale structure in the observed rapid decay phase.     

Precessing jets interacting with interstellar material as the origin for the light curves of gamma-ray bursts

We present an internal shock model with external characteristics for explaining the complicated light curves of gamma-ray bursts. Shocks produce gamma-rays in the interaction between a precessing beam of relativistic particles and the interstellar medium. Each time the particle beam passes the same line of sight with the observer the interstellar medium is pushed outward. Subsequent interactions between the medium and the beam are delayed by the extra distance to be travelled for the particles before the shock can form. This results in a natural retardation and leads to an intrinsic asymmetry in the light curves produced for gamma-ray bursts. In addition, we account for the cooling of the electron–proton plasma in the shocked region, which gives rise to an exponential decay in the gamma-ray flux. The combination of these effects and the precessing jet of ultrarelativistic particles produces light curves that can be directly compared with observed gamma-ray burst light curves. We illustrate the model by fitting a number of observed gamma-ray bursts that are difficult to explain with only a precessing jet. We develop a genetic algorithm to fit several observed gamma-ray bursts with remarkable accuracy. We find that for different bursts the observed fluence, assuming isotropic emission, easily varies over four orders of magnitude from the energy generated intrinsically.     

Gamma-ray binaries

Recent observations have shown that some compact stellar binaries radiate the highest energy light in the universe. The challenge has been to determine the nature of the compact object and whether the very high energy gamma-rays are ultimately powered by pulsar winds or relativistic jets. Multiwavelength observations have shown that one of the three gamma-ray binaries known so far, PSR B1259−63, is a neutron star binary and that the very energetic gamma-rays from this source and from another gamma-ray binary, LS I +61 303, may be produced by the interaction of pulsar winds with the wind from the companion star. At this time it is an open question whether the third gamma-ray binary, LS 5039, is also powered by a pulsar wind or a microquasar jet, where relativistic particles in collimated jets would boost the energy of the wind from the stellar companion to TeV energies. I.F. Mirabel is on leave from CEA, France.     

Broad-band non-thermal emission from molecular clouds illuminated by cosmic rays from nearby supernova remnants

Molecular clouds are expected to emit non-thermal radiation due to cosmic ray interactions in the dense magnetized gas. Such emission is amplified if a cloud is located close to an accelerator of cosmic rays and if energetic particles can leave the accelerator site and diffusively reach the cloud. We consider here a situation in which a molecular cloud is located in the proximity of a supernova remnant which is efficiently accelerating cosmic rays and gradually releasing them in the interstellar medium. We calculate the multiwavelength spectrum from radio to gamma rays which is emerging from the cloud as the result of cosmic ray interactions. The total energy output is dominated by the gamma-ray emission, which can exceed the emission in other bands by an order of magnitude or more. This suggests that some of the unidentified TeV sources detected so far, with no obvious or very weak counterparts in other wavelengths, might be in fact associated with clouds illuminated by cosmic rays coming from a nearby source. Moreover, under certain conditions, the gamma-ray spectrum exhibits a concave shape, being steep at low energies and hard at high energies. This fact might have important implications for the studies of the spectral compatibility of GeV and TeV gamma-ray sources.     

Gamma-ray bursts, ultra-high-energy cosmic rays, and cosmic gamma-ray background

We argue that gamma-ray bursts (GRBs) may be the origin of the cosmic gamma-ray background radiation observed in the GeV range. It has theoretically been discussed that protons may carry a much larger amount of energy than electrons in GRBs, and this large energy can be radiated in the TeV range by synchrotron radiation of ultra-high-energy protons ( 10²⁰ eV). The possible detection of GRBs above 10 TeV suggested by the Tibet and HEGRA groups also supports this idea. If this is the case, most of TeV gamma-rays from GRBs are absorbed in intergalactic fields and eventually form GeV gamma-ray background, whose flux is in good agreement with the recent observation.