Modeling the emission processes in blazars

Blazars are the most violent steady/recurrent sources of high-energy gamma-ray emission in the known Universe. They are prominent emitters of electromagnetic radiation throughout the entire electromagnetic spectrum. The observable radiation most likely originates in a relativistic jet oriented at a small angle with respect to the line of sight. This review starts out with a general overview of the phenomenology of blazars, including results from a recent multiwavelength observing campaign on 3C279. Subsequently, issues of modeling broadband spectra will be discussed. Spectral information alone is not sufficient to distinguish between competing models and to constrain essential parameters, in particular related to the primary particle acceleration and radiation mechanisms in the jet. Short-term spectral variability information may help to break such model degeneracies, which will require snap-shot spectral information on intraday time scales, which may soon be achievable for many blazars even in the gamma-ray regime with the upcoming GLAST mission and current advances in Atmospheric Cherenkov Telescope technology. In addition to pure leptonic and hadronic models of gamma-ray emission from blazars, leptonic/hadronic hybrid models are reviewed, and the recently developed hadronic synchrotron mirror model for TeV γ-ray flares which are not accompanied by simultaneous X-ray flares (“orphan TeV flares”) is revisited. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Gamma-ray emission from pulsars

We have attempted to devise a scheme by which it may be possible to identify pulsars which are likely to be γ-ray pulsars. We apply this test to a representative population of pulsars and identify the likely candidates for γ emission. We also discuss some individual cases including the Crab and Vela pulsars.     

Gamma-ray emission from individual classical novae

Classical novae are important producers of radioactive nuclei, such as ⁷Be, ¹³N, ¹⁸F, ²²Na and ²⁶Al. The disintegration of these nuclei produces positrons (except for ⁷Be) that through annihilation with electrons produce photons of energies 511 keV and below. Furthermore, ⁷Be and ²²Na decay producing photons with energies of 478 and 1275 keV, respectively, well in the γ-ray domain. Therefore, novae are potential sources of γ-ray emission. We have developed two codes in order to analyse carefully the γ-ray emission of individual classical novae: a hydrodynamical one, which follows both the accretion and the explosion stages, and a Monte Carlo one, able to treat both production and transfer of γ-ray photons. Both codes have been coupled in order to simulate realistic explosions. The properties of γ-ray spectra and γ-ray light curves (for the continuum and for the lines at 511, 478 and 1275 keV) have been analysed, with a special emphasis on the difference between carbon–oxygen and oxygen–neon novae. Predictions of detectability of individual novae by the future SPI spectrometer on board the INTEGRAL satellite are made. Concerning ²⁶Al, its decay produces photons of 1809 keV but this occurs on a time-scale much longer than the typical time interval between nova outbursts in the Galaxy, making it undetectable in individual novae. The accumulated emission of ²⁶Al from many Galactic novae has not been modelled in this paper.     

Gamma-ray luminosity function of blazars and the cosmic gamma-ray background: evidence for the luminosity-dependent density evolution

We present a comprehensive study of the gamma-ray luminosity function (GLF) of blazars and their contribution to the extragalactic diffuse gamma-ray background (EGRB). Radio and gamma-ray luminosity correlation is introduced with a modest dispersion, consistent with observations, to take into account the radio detectability, which is important for blazar identification. Previous studies considered only pure luminosity evolution (PLE) or pure density evolution, but here we introduce the luminosity-dependent density evolution (LDDE) model, which is favored on the basis of the evolution of the X-ray luminosity function (XLF) of AGNs. The model parameters are constrained by likelihood analyses of the observed redshift and gamma-ray flux distributions of the EGRET blazars. Interestingly, we find that the LDDE model gives a better fit to the observed distributions than the PLE model, indicating that the LDDE model is also appropriate for gamma-ray blazars and that the jet activity is universally correlated with the accretion history of AGNs. We then find that only 25–50% of the EGRB can be explained by unresolved blazars with the best-fit LDDE parameters. Unresolved blazars can account for all the EGRB only with a steeper index of the faint-end slope of the GLF, which is marginally consistent with the EGRET data but inconsistent with XLF data. Therefore, unresolved AGNs cannot be the dominant source of the EGRB, unless there is a new population of gamma-ray emitting AGNs that evolves differently from the XLF of AGNs. Predictions for the GLAST mission are made, and we find that the best-fit LDDE model predicts about 3000 blazars in the entire sky, which is considerably fewer (by a factor of more than 3) than a previous estimate.     

Gamma-ray emission expected from Kepler’s SNR

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler’s SNR and, in particular, to predict the γ-eay spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E _sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler’s SN is a type Ia SN we determine for any given explosion energy E _sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected γ-eay flux. Exploring the expected distance range we find that for a typical explosion energy E _sn=10⁵¹ erg the expected energy flux of TeV γ-rays varies from 2×10⁻¹¹ to 10⁻¹³ erg/(cm² s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the γ-eay emission is dominated by π ⁰-decay γ-rays due to nuclear CRs. Therefore Kepler’s SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of γ-rays would mean that the actual source distance is larger than 7 kpc.     

Multifrequency emission analysis of TeV blazars H 2356-309 and 1ES 1218+304

The multiband nonthermal emissions in radio,X-ray,and very high-energy (VHE)γ-ray bands from two distant blazars,H 2356-309 and 1ES 1218+304,have been detected,and,especially from recent observations with the Suzaku,MAGIC and VERITAS telescopes,clearly reveal nonthermal power-law spectra. We study the broadband nonthermal spectra of the two sources by using a combination of a onezone homogeneous synchrotron self-Compton (SSC) model and an inhomogeneous conical jet model,where the new external background light (EBL) model is taken into account. The results show that (1) the nonthermal emissions of the two blazars,ranging from X-rays to VHE γ-rays,are from the homogeneous zone whereas the emissions in the radio bands can be explained as the radiation from the inhomogeneous conical jet; (2) a strict lower-limit EBL model can be used to explain their observed spectra well.     

Gamma-ray and synchrotron emission from neutralino annihilation in the Large Magellanic Cloud

We calculate the expected flux of γ-ray and radio emission from the LMC due to neutralino annihilation. Using rotation curve data to probe the density profile and assuming a minimum disk, we describe the dark matter halo of the LMC using models predicted by N-body simulations. We consider a range of density profiles including the NFW profile, a modified NFW profile proposed by Hayashi et al. (2003) to account for the effects of tidal stripping, and an isothermal sphere with a core. We find that the γ-ray flux expected from these models may be detectable by GLAST for a significant part of the neutralino parameter space. The prospects for existing and upcoming Atmospheric Cherenkov Telescopes (ACTs) are less optimistic, as unrealistically long exposures are required for detection. However, the effects of adiabatic compression due to the baryonic component may improve the chances for detection by ACTs. The maximum flux we predict is well below EGRET's measurements and thus EGRET does not constrain the parameter space. The expected synchrotron emission generally lies below the observed radio emission from the LMC in the frequency range of 19.7–8550 MHz. As long as σv<2×10⁻²⁶ cm³ s⁻¹ for a neutralino mass of 50 GeV, the observed radio emission is not primarily due to neutralinos and is consistent with the assumption that the main source is cosmic rays. We find that the predicted fluxes, obtained by integrating over the entire LMC, are not very strongly dependent on the inner slope of the halo profile, varying by less than an order of magnitude for the range of profiles we considered.     

Gamma-ray emission from pulsars: strength of the acceleration field in the outer gap

In this study we present and re-analyse the historical, 1889–1998, light curve (LC) of the eclipsing symbiotic binary AR Pav. For the first time, we show that the timing of mid-points of eclipses observed during a quiescent phase obeys a quadratic ephemeris, with an initial orbital period P ₀=605.18 d and a rate of period change     .
We determined a distance to the system of 5.8±1.5 kpc, the mass ratio of the giant to the hot star, M _g M _h=0.4±0.1, the mass of the giant, M _g=1.8+1/−0.5 M_⊙ and its radius, R _g=167±15 R_⊙.
During quiescence, the LC has characteristic features similar to those observed in cataclysmic variables (CVs). It can be well reproduced by a model of a large accretion disc surrounding the hot star. However, it is probable that the geometry of the transferred material in the Roche lobe of the accretor in AR Pav is different from that of CVs.
During active phases the shape of the LC changes considerably. A complex wave-like variation developed as a function of the orbital phase with an amplitude of ∼1 mag. It is interpreted in terms of a collisionally heated emission region located on the giant surface and arising from the hot star eruption.     

The evolutionary sequence of Fermi blazars

Using γ-ray data (α _γ , F _γ ) detected by Fermi Large Area Telescope (LAT) and black hole mass which has been compiled from literatures for 116 Fermi blazars, we calculated intrinsic γ-ray luminosity, intrinsic bolometric luminosity, intrinsic Eddington ratio and studied the relationships between all above parameters and redshift, between α _γ and L _γ . Furthermore, we obtained the histograms of key parameters. Our results are the following: (1) The main reason for the evolutionary sequence of three subclasses (HBLs, LBLs, FSRQs) may be Eddington ratio rather than black hole mass; (2) FSRQs occupy in the earlier, high-luminosity, high Eddington ratio, violent phase of the galactic evolution sequence, while BL Lac objects occur in the low luminosity, low Eddington ratio, late phase of the galactic evolution sequence; (3) These results imply that the evolutionary track of Fermi blazars is FSRQs ? LBLs ? HBLs.     

Optical variability of gamma-ray loud blazars

Blazars display flux variability on diverse timescales ranging from minutes to months. In our blazar monitoring project carried out on the 1.2meter telescope at Gurushikhar, Mount Abu, India, we selected someγ-ray loud blazars to study the variability over both the short and the long time scales. In this paper we have reported results based on photometric monitoring on 8 nights during the observing run in the first half of the year 2000 for 5 blazars: S5 0716+714, OJ 287, S4 0954+658, Mrk421 and PKS 2155–304 in B and R passbands. Microvariability and rapid variability has been observed in these sources suggesting that variability is a common property in these objects. This revised version was published online in July 2006 with corrections to the Cover Date.     

Gamma-ray telescopes

The last half-century has seen dramatic developments in γ?ray telescopes, from their initial conception and development through to their blossoming into full maturity as a potent research tool in astronomy. Gamma-ray telescopes are leading research in diverse areas such as γ?ray bursts, blazars, Galactic transients, and the Galactic distribution of ²⁶Al.     

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.