No Arabic abstract
We present an in-depth and systematic variability study of a sample of 20 powerful blazars, including 12 BL Lacs and 8 flat spectrum radio quasars, applying various analysis tools such as flux distribution, symmetry analysis, and time series analysis on the decade-long Fermi/LAT observations. The results show that blazars with steeper gama-ray spectral indexes are found to be more variable; and the gama-ray flux distribution closely resembles lognormal probability distribution function. The statistical variability properties of the sources as studied by power spectral density analysis are consistent with emph{flicker noise} ($P( u)propto1/ u$) -- an indication of long-memory processes at work. Statistical analysis of the distribution of flux rise and decay rates in the light curves of the sources, aimed at distinguishing between particle acceleration and energy dissipation timescales, counter-intuitively suggests that both kinds of rates follow a similar distribution and the derived mean variability timescales are in the order of a few weeks. The corresponding emission region size is used to constrain location of gama-ray production sites in the sources to be a few parsecs. Additionally, using Lomb-Scargle periodogram and weighted wavelet z-transform methods and extensive Monte Carlo simulations, we detected year timescale quasi-periodic oscillations in the sources S5 0716+714, Mrk 421, ON +325, PKS 1424-418 and PKS 2155-304. The detection significance was computed taking proper account of the red-noise and other artifacts inherent in the observations. We explain the results in the light of current blazar models with relativistic shocks propagating down the jet viewed close to the line of sight.
We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.
Optical observations of a sample of 12 $gamma$-ray bright blazars from four optical data archives, AAVSO, SMARTS, Catalina, and Steward Observatory, are compiled to create densely sampled light curves spanning more than a decade. As a part of the bla
We present the time variability properties of a sample of six blazars, AO 0235+164, 3C 273, 3C 279, PKS 1510-089, PKS 2155-304, and 3C 454.3, at optical-IR as well as gamma-ray energies. These observations were carried out as a part of the Yale/SMARTS program during 2008-2010 that has followed the variations in emission of the bright Fermi-LAT-monitored blazars in the southern sky with closely-spaced observations at BVRJK bands. We find the optical/IR time variability properties of these blazars to be remarkably similar to those at the gamma-ray energies. The power spectral density (PSD) functions of the R-band variability of all six blazars are fit well by simple power-law functions with negative slope such that there is higher amplitude variability on longer timescales. No clear break is identified in the PSD of any of the sources. The average slope of the PSD of R-band variability of these blazars is similar to what was found by the Fermi team for the gamma-ray variability of a larger sample of bright blazars. This is consistent with leptonic models where the optical-IR and gamma-ray emission is generated by the same population of electrons through synchrotron and inverse-Compton processes, respectively. The prominent flares present in the optical-IR as well as the gamma-ray light curves of these blazars are predominantly symmetric, i.e., have similar rise and decay timescales, indicating that the long-term variability is dominated by the crossing time of radiation or a disturbance through the emission region rather than by the acceleration or energy-loss timescales of the radiating electrons. In the blazar 3C 454.3, which has the highest-quality light curves, the location of a large gamma-ray outburst during 2009 December is consistent with being in the jet at ~18 pc from the central engine. This poses strong constraints on the models of high energy emission in the jets of blazars.
We study the expected variability patterns of blazars within the two-zone acceleration model putting special emphasis on flare shapes and spectral lags. We solve semi-analytically the kinetic equations which describe the particle evolution in the acceleration and radiation zone. We then perturb the solutions by introducing Lorentzian variations in its key parameters and examine the flaring behavior of the system. We apply the above to the X-ray observations of blazar 1ES 1218+304 which exhibited a hard lag behavior during a flaring episode and discuss possibilities of producing it within the context of our model. The steady-state radio to X-rays emission of 1ES 1218+304 can be reproduced with parameters which lie well within the ones generally accepted from blazar modeling. Additionally, we find that the best way to explain its flaring behavior is by varying the rate of particles injected in the acceleration zone.
We present a new method for identifying blazar candidates by examining the locus, i.e. the region occupied by the Fermi gamma-ray blazars in the three-dimensional color space defined by the WISE infrared colors. This method is a refinement of our previous approach that made use of the two-dimensional projection of the distribution of WISE gamma-ray emitting blazars (the Strip) in the three WISE color-color planes (Massaro et al. 2012a). In this paper, we define the three-dimensional locus by means of a Principal Component (PCs) analysis of the colors distribution of a large sample of blazars composed by all the ROMA-BZCAT sources with counterparts in the WISE All-Sky Catalog and associated to gamma-ray source in the second Fermi LAT catalog (the WISE Fermi Blazars sample, WFB). Our new procedure yields a total completeness of c~81% and total efficiency of e~97%. We also obtain local estimates of the efficiency and completeness as functions of the WISE colors and galactic coordinates of the candidate blazars. The catalog of all WISE candidate blazars associated to the WFB sample is also presented, complemented by archival multi-frequency information for the alternative associations. Finally, we apply the new association procedure to all gamma-ray blazars in the 2FGL and provide a catalog containing all the gamma-ray candidates blazars selected according to our procedure.