A multifrequency analysis of radio variability of blazars

We have carried out a multifrequency analysis of the radio variability of blazars, exploiting the data obtained during the extensive monitoring programs carried out at the University of Michigan Radio Astronomy Observatory (UMRAO, at 4.8, 8, and 14.5 GHz) and at the Metsahovi Radio Observatory (22 and 37 GHz). Two different techniques detect, in the Metsahovi light curves, evidences of periodicity at both frequencies for 5 sources (0224+671, 0945+408, 1226+023, 2200+420, and 2251+158). For the last three sources consistent periods are found also at the three UMRAO frequencies and the Scargle (1982) method yields an extremely low false-alarm probability. On the other hand, the 22 and 37 GHz periodicities of 0224+671 and 0945+408 (which were less extensively monitored at Metsahovi and for which we get a significant false-alarm probability) are not confirmed by the UMRAO database, where some indications of ill-defined periods about a factor of two longer are retrieved. We have also investigated the variability index, the structure function, and the distribution of intensity variations of the most extensively monitored sources. We find a statistically significant difference in the distribution of the variability index for BL Lac objects compared to flat-spectrum radio quasars (FSRQs), in the sense that the former objects are more variable. For both populations the variability index steadily increases with increasing frequency. The distribution of intensity variations also broadens with increasing frequency, and approaches a log-normal shape at the highest frequencies. We find that variability enhances by 20-30% the high frequency counts of extragalactic radio-sources at bright flux densities, such as those of the WMAP and Planck surveys.Comment: A&A accepted. 12 pages, 16 figure

Statistical analyses of long-term variability of AGN at high radio frequencies

We present a study of variability time scales in a large sample of Active Galactic Nuclei at several frequencies between 4.8 and 230 GHz. We investigate the differences of various AGN types and frequencies and correlate the measured time scales with physical parameters such as the luminosity and the Lorentz factor. Our sample consists of both high and low polarization quasars, BL Lacertae objects and radio galaxies. The basis of this work is the 22 GHz, 37 GHz and 87 GHz monitoring data from the Metsahovi Radio Observatory spanning over 25 years. In addition,we used higher 90 GHz and 230 GHz frequency data obtained with the SEST-telescope between 1987 and 2003. Further lower frequency data at 4.8 GHz, 8 GHz and 14.5 GHz from the University of Michigan monitoring programme have been used. We have applied three different statistical methods to study the time scales: The structure function, the discrete correlation function and the Lomb-Scargle periodogram. We discuss also the differences and relative merits of these three methods. Our study reveals that smaller flux density variations occur in these sources on short time scales of 1-2 years, but larger outbursts happen quite rarely, on the average only once in every 6 years. We do not find any significant differences in the time scales between the source classes. The time scales are also only weakly related to the luminosity suggesting that the shock formation is caused by jet instabilities rather than the central black hole.Comment: 19 pages, 12 figures, Accepted for publication in A&

30 years of multi-wavelength observations of 3C 273

We present a wide multi-wavelength database of most observations of the quasar 3C 273 obtained during the last 30 years. This database is the most complete set of observations available for an active galactic nucleus (AGN). It contains nearly 20'000 observations grouped together into 70 light curves covering 16 orders of magnitude in frequency from the radio to the gamma-ray domain. The database is constituted of many previously unpublished observations and of most publicly available data gathered in the literature and on the World Wide Web (WWW). It is complete to the best of our knowledge, except in the optical (UBV) domain where we chose not to add all observations from the literature. In addition to the photometric data, we present the spectra of 3C 273 obtained by the International Ultraviolet Explorer (IUE) satellite. In the X-ray domain, we used the spectral fit parameters from the literature to construct the light curves. Apart from describing the data, we show the most representative light curves and the average spectrum of 3C 273. The database is available on the WWW in a homogeneous and clear form and we wish to update it regularly by adding new observations.Comment: 12 pages, 6 figures, to be published in A&AS, data available at: http://obswww.unige.ch/3c273

Long-term Variability Properties and Periodicity Analysis for Blazars

In this paper, the compiled long-term optical and infrared measurements of some blazars are used to analyze the variation properties and the optical data are used to search for periodicity evidence in the lightcurve by means of the Jurkevich technique and the discrete correlation function (DCF) method. Following periods are found: 4.52-year for 3C 66A; 1.56 and 2.95 years for AO 0235+164; 14.4, 18.6 years for PKS 0735+178; 17.85 and 24.7 years for PKS 0754+100; 5.53 and 11.75 for OJ 287. 4.45, and 6.89 years for PKS 1215; 9 and 14.84 years for PKS 1219+285; 2.0, 13.5 and 22.5 for 3C273; 7.1 year for 3C279; 6.07 for PKS 1308+326; 3.0 and 16.5 years for PKS 1418+546; 2.0 and 9.35 years for PKS 1514-241; 18.18 for PKS 1807+698; 4.16 and 7.0 for 2155-304; 14 and 20 years for BL Lacertae. Some explanations have been discussed.Comment: 10 pages, 2 table, no figure, a proceeding paper for Pacific Rim Conference on Stellar Astrophysics, Aug. 1999, HongKong, Chin

Long-term optical variability of PKS 2155-304

Aims: The optical variability of the blazar PKS 2155-304 is investigated to characterise the red noise behaviour at largely different time scales from 20 days to O(>10 yrs). Methods: The long-term optical light curve of PKS 2155-304 is assembled from archival data as well as from so-far unpublished observations mostly carried out with the ROTSE-III and the ASAS robotic telescopes. A forward folding technique is used to determine the best-fit parameters for a model of a power law with a break in the power spectral density function (PSD). The best-fit parameters are estimated using a maximum-likelihood method with simulated light curves in conjunction with the Lomb Scargle Periodogram (LSP) and the first-order Structure Function (SF). In addition, a new approach based upon the so-called Multiple Fragments Variance Function (MFVF) is introduced and compared to the other methods. Simulated light curves have been used to confirm the reliability of these methods as well as to estimate the uncertainties of the best-fit parameters. Results: The light curve is consistent with the assumed broken power-law PSD. All three methods agree within the estimated uncertainties with the MFVF providing the most accurate results. The red-noise behaviour of the PSD in frequency f follows a power law with f^-{\beta}, {\beta}=1.8 +0.1/-0.2 and a break towards f^0 at frequencies lower than f_min=(2.7 +2.2/-1.6 yrs)^-1.Comment: 10 pages, 8 figures, the ROTSE-light curve can be downloaded from http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=J/A+A/531/A12

Helical motion and the origin of QPO in blazar-type sources

Recent observations and analysis of blazar sources provide strong evidence for (i) the presence of significant periodicities in their lightcurves and (ii) the occurrence of helical trajectories in their radio jets. In scenarios, where the periodicity is caused by differential Doppler boosting effects along a helical jet path, both of these facts may be naturally tied together. Here we discuss four possible driving mechanisms for the occurrence of helical trajectories: orbital motion in a binary system, Newtonian-driven jet precession, internal jet rotation and motion along a global helical magnetic field. We point out that for non-ballistic helical motion the observed period may appear strongly shortened due to classical travel time effects. Finally, the possible relevance of the above mentioned driving mechanisms is discussed for Mkn~501, OJ 287 and AO 0235+16.Comment: 6 pages, 1 figure; presented at the 5th Microquasar Workshop, Beijing, June 2004. Accepted for publication in the Chinese Journal of Astronomy and Astrophysic

Long-term radio variability of AGN: flare characteristics

We have studied the flare characteristics of 55 AGN at 8 different frequency bands between 4.8 and 230 GHz. Our extensive database enables us to study the various observational properties of flares in these sources and compare our results with theoretical models. We visually extracted 159 individual flares from the flux density curves and calculated different parameters, such as the peak flux density and duration, in all the frequency bands. The selection of flares is based on the 22 and 37 GHz data from Mets\"ahovi Radio Observatory and 90 and 230 GHz data from the SEST telescope. Additional lower frequency 4.8, 8, and 14.5 GHz data are from the University of Michigan Radio Observatory. We also calculated variability indices and compared them with earlier studies. The observations seem to adhere well to the shock model, but there is still large scatter in the data. Especially the time delays between different frequency bands are difficult to study due to the incomplete sampling of the higher frequencies. The average duration of the flares is 2.5 years at 22 and 37 GHz, which shows that long-term monitoring is essential for understanding the typical behaviour in these sources. It also seems that the energy release in a flare is independent of the duration of the flare.Comment: 11 pages, 9 figures, 2 tables, accepted for publication in A&

Blazar 3C 454.3 in Outburst and Quiescence During 2005-2007: Two Variable Synchrotron Emission Peaks

We monitored the flaring blazar 3C 454.3 during 2005 June-July with the Spitzer Infrared Spectrograph (IRS: 15 epochs), Infrared Array Camera (IRAC: 12 epochs) and Multiband Imaging Photometer (MIPS: 2 epochs). We also made Spitzer IRS, IRAC, and MIPS observations from 2006 December-2007 January when the source was in a low state, the latter simultaneous with a single Chandra X-ray observation. In addition, we present optical and sub-mm monitoring data. The 2005-2007 period saw 3 major outbursts. We present evidence that the radio-optical SED actually consists of two variable synchrotron peaks, the primary at IR and the secondary at sub-mm wavelengths. The lag between the optical and sub-mm outbursts may indicate that these two peaks arise from two distinct regions along the jet separated by a distance of 0.07-5 pc. The flux at 5-35 microns varied by a factor of 40 and the IR peak varied in frequency from <1E13 Hz to 4E13 Hz between the highest and lowest states in 2005 and 2006, respectively. Variability was well correlated across the mid-IR band, with no measurable lag. Flares that doubled in flux occurred on a time scale of 3 days. The IR SED peak moved to higher frequency as a flare brightened, then returned to lower frequency as it decayed. The fractional variability amplitude increased with frequency, which we attribute to decreasing synchrotron-self absorption optical depth. Mid-IR flares may signal the re-energization of a shock that runs into inhomogeneities along the pre-existing jet or in the external medium. The synchrotron peak frequencies during each major outburst may depend upon both the distance from the jet apex and the physical conditions in the shocks. Variation of the Doppler parameter along a curved or helical jet is another possibility. Frequency variability of the IR synchrotron peak may have important consequences for the interpretation of the blazar sequence, and the presence of a secondary peak may give insight into jet structure.Comment: 38 pages, 15 figures, submitted to ApJS, comments welcom

Optical Properties of faint FIRST Variable Radio Sources

A sample of 123 radio sources that exhibit significant variations at 1.4 GHz on a seven year base-line has been created using FIRST VLA B-array data from 1995 and 2002 on a strip at dec=0 near the south Galactic cap. This sample spans the range of radio flux densities from ~2 to 1000 mJy. It presents both in size and radio flux density range a unique starting point for variability studies of galaxies and quasars harboring lower luminosity Active Galactic Nuclei (AGN). We find, by comparing our variable and non-variable control samples to the Sloan Digital Sky Survey the following: 1) The quasar fraction of both the variable and non-variable samples declines as a function of declining radio flux density levels; 2) our variable sample contains a consistently higher fraction of quasars than the non-variable control sample, irrespective of radio flux; 3) the variable sources are almost twice as likely to be retrieved from the optical SDSS data than the non-variable ones; 4) based on relative numbers, we estimate that quasars are about five times more likely to harbor a variable radio source than are galaxies; and 5) there does not appear to be any significant optical color offset between the two samples, even though the suggestive trend for sources to be bluer when variable has been detected before and may be real. This leads us to conclude that both radio variability and radio flux density levels, in combination with accurate optical information, are important discriminators in the study of (radio) variability of galaxies. The latter start to dominate the source counts below about 20 mJy. In any case, variability appears to be an intrinsic property of radio sources, and is not limited to quasars. Radio variability at low flux density levels may offer a unique tool in AGN unification studies.Comment: 22 pages, 10 color figures. Accepted for publication in AJ, May issu

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 gamma-ray flares which are not accompanied by simultaneous X-ray flares (``orphan TeV flares'') is revisited.Comment: Invited Review at "The Multimessenger Approach to Gamma-Ray Sources", Barcelona, Spain, July 2006; submitted to Astrophysics and Space Science. 10 pages, including 6 eps figures. Uses Springer's ApSS macro