Blazar Jets Research Articles

VERITAS OBSERVATIONS OF RELATIVISTIC JETS

VERITAS is an array of four imaging atmospheric Cherenkov telescopes located in southern Arizona, and is now starting its fifth year of observations of the northern Sky at energies > 100 GeV. Observations of astrophysical objects in the TeV band are sensitive probes of highly energetic processes occurring in these sources. The majority of the active galaxies detected at TeV energies are blazars, sources where we view the jet nearly along its axis. TeV and multiwavelength observations of blazars help us to better understand the mechanisms of ultrarelativistic jet production by supermassive black holes, and constrain models of particle acceleration in blazar jets. Galactic sources at TeV energies include supernova remnants, pulsar wind nebulae, and binary systems, and TeV emission is a key diagnostic of highly energetic particles in these objects. This report presents recent results from VERITAS on Galactic and extragalactic sources.

THE HIGH-ENERGY PROPERTIES OF 3C 454.3: A FOUR-YEAR STUDY

The blazar 3C 454.3 has become the most active and brightest γ-ray source of the sky, earning the nickname of Crazy Diamond. The short-term variability in the γ-ray energy band and the extremely high peak fluxes reached during intense flaring episodes make 3C 454.3 one of the best targets to investigate the blazar jet properties. We will review four years of observational properties of this remarkable source, discussing both short- and long-term multi-wavelength campaigns, with particular emphasis on the recent flaring episode which occurred on 2010 November 20, when 3C 454.3 reached on a daily time-scale a gamma-ray flux (E > 100 MeV) higher than 6 × 10-5 photons cm-2 s-1, about six times the flux of the brightest γ-ray steady source, the Vela Pulsar.

High-frequency predictions for number counts and spectral properties of extragalactic radio sources. New evidence of a break at mm wavelengths in spectra of bright blazar sources

We present models to predict high frequency counts of extragalactic radio sources using physically grounded recipes to describe the complex spectral behaviour of blazars, that dominate the mm-wave counts at bright flux densities. We show that simple power-law spectra are ruled out by high-frequency (nu>100 GHz) data. These data also strongly constrain models featuring the spectral breaks predicted by classical physical models for the synchrotron emission produced in jets of blazars (Blandford & Konigl 1979; Konigl 1981). A model dealing with blazars as a single population is, at best, only marginally consistent with data coming from current surveys at high radio frequencies. Our most successful model assumes different distributions of break frequencies, nu_M, for BL Lacs and Flat-Spectrum Radio Quasars (FSRQs). The former objects have substantially higher values of nu_M, implying that the synchrotron emission comes from more compact regions; therefore, a substantial increase of the BL Lac fraction at high radio frequencies and at bright flux densities is predicted. Remarkably, our best model is able to give a very good fit to all the observed data on number counts and on distributions of spectral indices of extragalactic radio sources at frequencies above 5 and up to 220 GHz. Predictions for the forthcoming sub-mm blazar counts from Planck, at the highest HFI frequencies, and from Herschel surveys are also presented.

Rapid optical variability of TeV blazars

In this first systematic attempt to characterise the intranight optical variability (INOV) of TeV detected blazars, we have monitored a well defined set of 9 TeV blazars on total 26 nights during 2004-2010. In this R (or V)-band monitoring programme only one blazar was monitored per night for a minimum duration of 4 hours. Using the CCD, an INOV detection threshold of ~ 1-2 % was achieved in the densely sampled DLCs. We have further expanded the sample by including another 13 TeV blazars from literature. This enlarged sample of 22 TeV blazars, monitored on a total of 116 nights (including 55 nights newly reported here), has enabled us to arrive at the first estimate of the INOV duty cycle of TeV detected blazars. Applying the C-test, the INOV DC is found to be 59 %, which decreases to 47 % if only INOV fractional amplitudes above 3 % are considered. These observations also permit, for the first time, a comparison of the INOV characteristics of the two major subclasses of TeV detected BL Lacs, namely LBLs and HBLs, for which we find the INOV DCs to be ~ 63 % and ~ 38 %, respectively. This demonstrates that the INOV differential between LBLs and HBLs persists even when only their TeV detected subsets are considered. Despite dense sampling, the intranight light curves of the 22 TeV blazars have not revealed even a single feature on time scale substantially shorter than 1 hour, even though the inner jets of TeV blazars are believed to have exceptionally large bulk Lorentz factors (and correspondingly stronger time compression). An intriguing feature, clearly detected in the light curve of the HBL J1555+1111, is a 4 per cent `dip' on a 1 hour timescale. This unique feature could have arisen from absorption in a dusty gas cloud, occulting a superluminally moving optical knot in the parsec scale jet of this relatively luminous BL Lacs object.

Photopolarimetric Monitoring of Blazars in the Optical and Near-Infrared Bands with the Kanata Telescope. I. Correlations between Flux, Color, and Polarization

Abstract We report on the correlation between the flux, color, and polarization variations on time scales of days–months in blazars, and discuss their universal aspects. We performed monitoring of 42 blazars in the optical and near-infrared bands from 2008 to 2010 using TRISPEC attached to the ``Kanata'' 1.5-m telescope. We found that 28 blazars exhibited ``bluer-when-brighter'' trends in their whole or a part of time-series data sets. This corresponds to 88% of objects that were observed for $\gt\ $10 days. Thus, our observation unambiguously confirmed that the ``bluer-when-brighter'' trend is common in the emission from blazar jets. This trend was apparently generated by a variation component with a constant and relatively blue color and an underlying red component. Prominent short-term flares on time scales of days–weeks tended to exhibit a spectral hysteresis; their rising phases were bluer than their decay phases around the flare maxima. In contrast to the strong flux–color correlation, the correlation of the flux and polarization degree was relatively weak; only 10 objects showed significant positive correlations. Rotations of polarization were detected only in three objects: PKS 1510$-$089, 3C 454.3, and PKS 1749$+$096, and possibly in S5 0716$+$714. We also investigated the dependence of the degree of variability on the luminosity and the synchrotron peak frequency, $\nu _{\rm peak}$. As a result, we found that lower luminosity and higher $\nu _{\rm peak}$ objects had smaller variations in their amplitudes both in the flux, color, and polarization degree. Our observation suggests the presence of several distinct emitting sources, which have different variation time-scales, colors, and polarizations. We propose that the energy injection by, for example, internal shocks in relativistic shells is a major factor for blazar variations on time scales of both days and months.

Models for High-Energy Radiation from Blazars

We discuss on the modelling of blazar jets as emitters of multiwavelength radiation with the implementation of a lepto-hadronic treatment. Assuming that injection of non-thermal electrons and protons can take place at the base of the jet, the stationary particle distributions can be found using an inhomogeneous one-dimensional transport equation with cooling and convection. The goal of this approach is to replace the widely used one-zone purely leptonic approximation by a more realistic model. We argue that the rapid variability observed in emission from blazars can be obtained as a result of interaction of the jet with obstacles, i.e., molecular clouds and stars. Long term variability is likely related to changes in the injection and physical conditions in the acceleration region.

Study of the variability of blazars gamma-ray emission

The γ-ray emission of blazar jets shows a pronounced variability and this feature provides limits to the size and to the speed of the emitting region. We study the γ-ray variability of bright blazars using data from the first 18 months of activity of the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. From the daily light-curves of the blazars characterized by a remarkable activity, we firstly determine the minimum variability time-scale, giving an upper limit for the size of the emitting region of the sources, assumed to be spheroidal blobs in relativistic motion. These regions must be smaller than ∼10 −3 parsec. Another interesting time-scale is the duration of the outbursts. We conclude that they cannot correspond to radiation produced by a single blob moving relativistically along the jet, but they are either the signature of emission from a standing shock extracting energy from a modulated jet, or the superposition of a number of flares occurring on a shorter time-scale. We also derive lower limits on the bulk Lorentz factor needed to make the emitting region transparent for gamma-rays interacting through photon–photon collisions.

EMISSION FROM HOT DUST IN THE INFRARED SPECTRA OF GAMMA-RAY BRIGHT BLAZARS

A possible source of $\gamma$-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four $\gamma$-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510$-$089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature $\sim1200$ K, plus a much weaker optically thin component at $\sim660$ K. The total luminosity of the thermal dust emission is $7.9\pm0.2 \times 10^{45}$ erg s$^{-1}$. If the dust lies in an equatorial torus, the density of IR photons from the torus is sufficient to explain the $\gamma$-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is $7 \pm 2 \times 10^{45}$ erg s$^{-1}$. However, in CTA102 the far-IR spectra are too noisy to detect the $10 \mu$m silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510-089, and ON231, are, $2.3\times10^{45}$, and $6.6\times10^{43}$ erg s$^{-1}$, respectively. These upper limits do not rule out the possibility of inverse Compton up-scattering of IR photons to $\gamma$-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust.

TIME-DEPENDENT RADIATION TRANSFER IN THE INTERNAL SHOCK MODEL SCENARIO FOR BLAZAR JETS

We describe the time-dependent radiation transfer in blazar jets, within the internal shock model. We assume that the central engine, which consists of a black hole and an accretion disk, spews out relativistic shells of plasma with different velocity, mass, and energy. We consider a single inelastic collision between a faster (inner) and a slower (outer) moving shell. We study the dynamics of the collision and evaluate the subsequent emission of radiation via the synchrotron and synchrotron self Compton (SSC) processes after the interaction between the two shells has begun. The collision results in the formation of a forward shock (FS) and a reverse shock (RS) that convert the ordered bulk kinetic energy of the shells into magnetic field energy and accelerate the particles, which then radiate. We assume a cylindrical geometry for the emission region of the jet. We treat the self-consistent radiative transfer by taking into account the inhomogeneity in the photon density throughout the region. In this paper, we focus on understanding the effects of varying relevant input parameters on the simulated spectral energy distribution (SED) and spectral variability patterns.

Incidence of Mg ii absorbers towards blazars and the GRB/QSO puzzle

In order to investigate the origin of the excess of strong MgII systems towards GRB afterglows as compared to QSO sightlines, we have measured the incidence of MgII absorbers towards a third class of objects: the Blazars. This class includes the BL Lac object population for which a tentative excess of MgII systems had already been reported. We observed with FORS1 at the ESO-VLT 42 Blazars with an emission redshift 0.8<z_em<1.9, to which we added the three high z northern objects belonging to the 1Jy BL Lac sample. We detect 32 MgII absorbers in the redshift range 0.35-1.45, leading to an excess in the incidence of MgII absorbers compared to that measured towards QSOs by a factor ~2, detected at 3 sigma. The amplitude of the effect is similar to that found along GRB sightlines. Our analysis provides a new piece of evidence that the observed incidence of MgII absorbers might depend on the type of background source. In front of Blazars, the excess is apparent for both 'strong' (w_ r(2796) > 1.0 A) and weaker (0.3 < w_r(2796) < 1.0 A) MgII systems. The dependence on velocity separation with respect to the background Blazars indicates, at the ~1.5 sigma level, a potential excess for beta = v/c ~0.1. We show that biases involving dust extinction or gravitational amplification are not likely to notably affect the incidence of MgII systems towards Blazars. Finally we discuss the physical conditions required for these absorbers to be gas entrained by the powerful Blazar jets. More realistic numerical modelling of jet-ambient gas interaction is required to reach any firm conclusions as well as repeat observations at high spectral resolution of strong MgII absorbers towards Blazars in both high and low states.

RESOLVING DOPPLER-FACTOR CRISIS IN ACTIVE GALACTIC NUCLEI: NON-STEADY MAGNETIZED OUTFLOWS

Magnetically-driven non-stationary acceleration of jets in AGNs results in the leading parts of the flow been accelerated to much higher Lorentz factors than in the case of steady state acceleration with the same parameters. The higher Doppler-boosted parts of the flow may dominate the high energy emission of blazar jets. We suggest that highly variable GeV and TeV emission in blazars is produced by the faster moving leading edges of highly magnetized non-stationary ejection blobs, while the radio data trace the slower-moving bulk flow. Model predictions compare favorably with the latest Fermi gamma-ray and MOJAVE radio VLBI results.

Variability studies in blazar jets with SF analysis: caveats and problems

AbstractBlazars are radio-loud active galactic nuclei (AGN) dominated by relativistic jets seen at small angles to the line-of-sight. They exhibit dramatic flux variations across the electromagnetic spectrum. The fastest variations are observed in the X-ray and γ-ray bands on time-scales of hours or even minutes. Currently, a substantial part of the blazar literature has been based on the study of these temporal variations through the use of structure function (SF) analysis, the results of which are believed to put great constrains on the jet-physics.The SF is often invoked in the framework of shot-noise models to determine the temporal properties of individual shots within the jet as well as their geometrical sizes. We argue, that for blazar variability studies, the SF-results are sometimes erroneously interpreted leading to misconceptions about the actual source properties. Based on extensive simulations we caution that spurious breaks will appear in the SFs of almost all light-curves, even though these light-curves may contain no intrinsic characteristic time-scale.Finally, it is also commonly thought that SFs are immune to the sampling problems, such as data-gaps, which affects the estimators of the underlying power spectra density function such as the periodogram. However, we show that SFs are also troubled by gaps which can induce artefacts.

The picture of relativistic jet from Fermi-LAT and multi-band observations of blazar 3C 279

AbstractStrong and variable radiation detected over all accessible energy bands in blazar arises from a relativistic, Doppler-boosted jet pointing close to our line of sight. Flat Spectrum Radio Quasar 3C 279 was one of the brightest γ-ray blazars in the sky at the time of the discovery with EGRET. Since the successful launch of the Fermi Gamma-ray Space telescope in 2008, we have organized extensive multi-band observational campaign of 3C 279 from radio to γ-ray bands, also including optical polarimetric observations. The uninterrupted monitoring in the γ-ray band by Fermi-LAT together with the multi-band data provide us with new insights of the relativistic jet of blazar. Here, we present the results of the first-year multi-band campaign of 3C 279 including the discovery of a γ-ray flare event associated with a dramatic change of the optical polarization - as well as a discovery of an “orphan” X-ray flare, unassociated with prominent outbursts in other bands.

Blazar jet physics in the age of Fermi

AbstractThe impact of the Fermi Gamma-ray Space Telescope on blazar research is reviewed. This includes a brief description of the Fermi Large Area Telescope, a summary of the various classes of extragalactic sources found in the First Large Area Telescope AGN Catalog, and more detailed discussion of the flat spectrum radio quasar 3C454.3 and the BL Lac object PKS 2155-304. Some theoretical studies related to ongoing blazar research with Fermi are mentioned, including implications of γ-ray observations of radio galaxies on blazar unification scenarios, variability in colliding shells, and whether blazars are sources of ultra-high energy cosmic rays.

Additional Acceleration of Protons and Energetic Neutrino Production in a Filamentary Jet of the Blazar Markarian 501

Blazars have been regarded as one of the most powerful sources of the highest energy cosmic rays and also their byproducts, neutrinos. Provided that a magnetized filamentary system is established in a blazar jet as well, we could apply the mechanism of multi-stage diffusive shock acceleration to a feasible TeV emitter, Mrk 501 to evaluate the achievable maximum energy of protons. Taking conceivable energy restriction into account systematically, it seems adequate to say that EeV-protons are produced at this site by our present model. We also estimate neutrino fluxes generated by these accelerated protons, and discuss the detectability based on an updated kilometer-scale telescope, such as IceCube.

TIMING SIGNATURES OF THE INTERNAL-SHOCK MODEL FOR BLAZARS

We investigate the spectral and timing signatures of the internal-shock model for blazars. For this purpose, we develop a semi-analytical model for the time-dependent radiative output from internal shocks arising from colliding relativistic shells in a blazar jet. The emission through synchrotron and synchrotron-self Compton (SSC) radiation as well as Comptonization of an isotropic external radiation field are taken into account. We evaluate the discrete correlation function (DCF) of the model light curves in order to evaluate features of photon-energy dependent time lags and the quality of the correlation, represented by the peak value of the DCF. The almost completely analytic nature of our approach allows us to study in detail the influence of various model parameters on the resulting spectral and timing features. This paper focuses on a range of parameters in which the gamma-ray production is dominated by Comptonization of external radiation, most likely appropriate for gamma-ray bright flat-spectrum radio quasars (FSRQs) or low-frequency peaked BL Lac objects (LBLs). In most cases relevant for FSRQs and LBLs, the variability of the optical emission is highly correlated with the X-ray and high-energy (HE: > 100 MeV) gamma-ray emission. Our baseline model predicts a lead of the optical variability with respect to the higher-energy bands by 1 - 2 hours and of the HE gamma-rays before the X-rays by about 1 hour. We show that variations of certain parameters may lead to changing signs of inter-band time lags, potentially explaining the lack of persistent trends of time lags in most blazars.

Variability and stability in blazar jets on time-scales of years: optical polarization monitoring of OJ 287 in 2005-2009

OJ 287 is a BL Lac object at redshift z= 0.306 that has shown double-peaked bursts at regular intervals of ∼12 yr during the last ∼40 yr. We analyse optical photopolarimetric monitoring data from 2005 to 2009, during which the latest double-peaked outburst occurred. The aim of this study is twofold: firstly, we aim to analyse variability patterns and statistical properties of the optical polarization light curve. We find a strong preferred position angle in optical polarization. The preferred position angle can be explained by separating the jet emission into two components: an optical polarization core and chaotic jet emission. The optical polarization core is stable on time-scales of years and can be explained as emission from an underlying quiescent jet component. The chaotic jet emission sometimes exhibits a circular movement in the Stokes plane. We find six such events, all on the time-scales of 10–20 d. We interpret these events as a shock front moving forwards and backwards in the jet, swiping through a helical magnetic field. Secondly, we use our data to assess different binary black hole models proposed to explain the regularly appearing double-peaked bursts in OJ 287. We compose a list of requirements a model has to fulfil to explain the mysterious behaviour observed in OJ 287. The list includes not only characteristics of the light curve but also other properties of OJ 287, such as the black hole mass and restrictions on accretion flow properties. We rate all existing models using this list and conclude that none of the models is able to explain all observations. We discuss possible new explanations and propose a new approach to understanding OJ 287. We suggest that both the double-peaked bursts and the evolution of the optical polarization position angle could be explained as a sign of resonant accretion of magnetic field lines, a ‘magnetic breathing’ of the disc.

EXTENDED RADIO EMISSION IN MOJAVE BLAZARS: CHALLENGES TO UNIFICATION

We present the results of a study on the 1.4 GHz kpc-scale radio emission in the complete flux density limited MOJAVE sample, comprising 135 radio-loud AGNs. While extended emission is detected in the majority of the sources, about 7% of the sources exhibit only radio core emission. Many BL Lacs exhibit extended radio power and kpc-scale morphology typical of powerful FRII jets, while a substantial number of quasars possess radio powers intermediate between FRIs and FRIIs. This poses challenges to the simple radio-loud unified scheme, which links BL Lacs to FRIs and quasars to FRIIs. We find a significant correlation between extended radio emission and pc-scale jet speeds: the more radio powerful sources possess faster jets. This indicates that the 1.4 GHz (or low frequency) radio emission is indeed related to jet kinetic power. Various properties such as extended radio power and apparent pc-scale jet speeds vary smoothly between different blazar subclasses, suggesting that, at least in terms of radio jet properties, the distinction between quasars and BL Lac objects, at an emission-line equivalent width of 5 Angstrom is essentially an arbitrary one. Based on the assumption that the extended radio luminosity is affected by the kpc-scale environment, we define the ratio of extended radio power to absolute optical magnitude as a proxy for environmental effects. Trends with this parameter suggest that the pc-scale jet speeds and the pc-to-kpc jet misalignments are not affected by the large-scale environment, but are more likely to depend upon factors intrinsic to the AGN, or its local pc-scale environment. We suggest that some of the extremely misaligned MOJAVE blazar jets could be "hybrid" morphology sources, with an FRI jet on one side and an FRII jet on the other. (Abridged)

DISCOVERY OF VARIABILITY IN THE VERY HIGH ENERGY γ-RAY EMISSION OF 1ES 1218+304 WITH VERITAS

We present results from an intensive VERITAS monitoring campaign of the high-frequency peaked BL Lac object 1ES 1218+304 in 2008/2009. Although 1ES 1218+304 was detected previously by MAGIC and VERITAS at a persistent level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite common in many nearby blazars, the case of 1ES 1218+304 (redshift z = 0.182) is particularly interesting since it belongs to a group of blazars that exhibit unusually hard VHE spectra considering their redshifts. When correcting the measured spectra for absorption by the extragalactic background light, 1ES 1218+304 and a number of other blazars are found to have differential photon indices less than 1.5. The difficulty in modeling these hard spectral energy distributions in blazar jets has led to a range of theoretical gamma-ray emission scenarios, one of which is strongly constrained by these new VERITAS observations. We consider the implications of the observed light curve of 1ES 1218+304, which shows day scale flux variations, for shock acceleration scenarios in relativistic jets, and in particular for the viability of kiloparsec-scale jet emission scenarios.

Contribution from unresolved discrete sources to the extragalactic gamma-ray background (EGRB)

The origin of the extragalactic gamma-ray background (EGRB) is still an open question, even nearly forty years after its discovery. The emission could originate either from truly diffuse processes or from unresolved point sources. Although the majority of the 271 point sources detected by EGRET (Energetic Gamma Ray Experiment Telescope) are unidentified, of the identified sources, blazars are the dominant candidates. Therefore, unresolved blazars may be considered the main contributor to the EGRB, and many studies have been carried out to understand their distribution, evolution and contribution to the EGRB. Considering that γ-ray emission comes mostly from jets of blazars and that the jet emission decreases rapidly with increasing jet to line-of-sight angle, it is not surprising that EGRET was not able to detect many large inclination angle active galactic nuclei (AGNs). Though Fermi could only detect a few large inclination angle AGNs during the first three months of its survey, it is expected to detect many such sources in the near future. Since non-blazar AGNs are expected to have higher density as compared to blazars, these could also contribute significantly to the EGRB. In this paper, we discuss contributions from unresolved discrete sources including normal galaxies, starburst galaxies, blazars and off-axis AGNs to the EGRB.