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&

Variability of the Spectral Energy Distribution of the Blazar S5 0716+714

The emission from blazars is known to be variable at all wavelengths. The flux variability is often accompanied by spectral changes. Spectral energy distribution (SED) changes must be associated with changes in the spectra of emitting electrons and/or the physical parameters of the jet. Meaningful modeling of blazar broadband spectra is required to understand the extreme conditions within the emission region. Not only is the broadband SED crucial, but also information about its variability is needed to understand how the highest states of emission occur and how they differ from the low states. This may help in discriminating between models. Here we present the results of our SED modeling of the blazar S5 0716+714 during various phases of its activity. The SEDs are classified into different bins depending on the optical brightness state of the source.Comment: 4 pages, 3 figures, contributed talk presented at the conference Multifrequency Variability of Blazars, Guangzhou, China, September 22-24, 2010. To appear in Journal of Astrophysics and Astronomy (JAA

Multi-epoch VLBA observations of 3C 66A

We present the results of six-epoch Very Long Baseline Array (VLBA) observations of 3C~66A. The high-resolution Very Long Baseline Interferometer (VLBI) maps obtained at multi-frequency (2.3, 8.4, and 22.2 GHz) simultaneously enabled us to identify the brightest compact component with the core. We find that the spectrum of the core can be reasonably fitted by the synchrotron self-absorption model. Our VLBA maps show that the jet of 3C~66A has two bendings at about 1.2 and 4 mas from the core. We also give possible identifications of our jet components with the components in previous VLBA observations by analysing their proper motions. We find consistent differences of the position from the core in one component between different frequencies at six epochs.Comment: 10 pages, 5 figures, received 30 January 2007, accepted 22 March 200

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

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&

Wavelet analysis of a large sample of AGN at high radio frequencies

We have studied the characteristic timescales of 80 AGNs at 22, 37 and 90 GHz examining the properties of the wavelet method and comparing them to traditional Fourier-based methods commonly used in astronomy. We used the continuous wavelet transform with the Morlet wavelet to study the characteristic timescales. We also gain information when the timescale is present in the flux curve and if it is persistent or not. Our results show that the sources are not periodic and changes in the timescales over a long time are common. The property of wavelets to be able to distinguish when the timescale has been present is superior to the Fourier-based methods. Therefore we consider it appropriate to use wavelets when the quasi-periodicities in AGNs are studied.Comment: 14 pages, 10 figures, A&A in pres

Frequency-dependent time-delays for strong outbursts in selected blazars from the Metsähovi and the University of Michigan Radio Astronomy Observatory monitoring databases – I

The combined data of the University of Michigan Radio Astronomy Observatory and Metsähovi Radio Observatory provide us with radio light curves for active galactic nuclei monitored by both observatories from 4.8 to 37 GHz covering time-intervals up to similar to 25 yr. We consider here such composite light curves for four gamma-ray blazars that have been nearly continuously monitored at both observatories: 0458-020, 0528+134, 1730-130 and 2230+114. We have decomposed the most prominent outbursts in the light curves of these four blazars into individual components using Gaussian model fitting, and estimated the epochs, amplitudes, and half-widths of these components as functions of frequency. We attempt to distinguish 'core outbursts', which show frequency-dependent time-delays and are associated with brightening of the core, from 'jet outbursts', which appear nearly synchronous at all frequencies and are accompanied by the emergence of new jet components and their subsequent evolution. The outbursts in 0528+134 and 2230+114 display fine structure and consist of individual sub-outbursts. Available 43-GHz Very Long Baseline Array images allow us to identify only one pure core outburst (in 2230+114) and one pure jet outburst (0458-020). Most of the outbursts analysed are mixed, in the sense that they display frequency-dependent time-delays (i.e. they are optically thick) and are associated with the eventual emergence of new jet components. The maxima of the jet and mixed outbursts probably correspond to epochs when newly ejected components become fully optically thin. These epochs are also marked by a significant increase in the angular velocities of the ejected components. There is evidence that the outbursts in 2230+114 repeat every 8.0 +/- 0.3 yr, with the positions of individual suboutbursts being preserved from one quasi-periodic eight-year cycle to another, even though their amplitudes vary by more than a factor of 2. Preliminary estimates of the total durations of possible activity cycles based on an analysis of total flux-density variations and all available very long baseline interferometry data are given for the remaining sources

A major radio outburst in III Zw 2 with an extremely inverted, millimeter-peaked spectrum

III Zw 2 is a spiral galaxy with an optical spectrum and faint extended radio structure typical of a Seyfert galaxy, but also with an extremely variable, blazar-like radio core. We have now discovered a new radio flare where the source has brightened more than twenty-fold within less than two years. A broad-band radio spectrum between 1.4 and 666 GHz shows a textbook-like synchrotron spectrum peaking at 43 GHz, with a self-absorbed synchrotron spectral index +2.5 at frequencies below 43 GHz and an optically thin spectral index -0.75 at frequencies above 43 GHz. The outburst spectrum can be well fitted by two homogenous, spherical components with equipartition sizes of 0.1 and 0.2 pc at 43 and 15 GHz, and with magnetic fields of 0.4 and 1 Gauss. VLBA observations at 43 GHz confirm this double structure and these sizes. Time scale arguments suggest that the emitting regions are shocks which are continuously accelerating particles. This could be explained by a frustrated jet scenario with very compact hotspots. Similar millimeter-peaked spectrum (MPS) sources could have escaped our attention because of their low flux density at typical survey frequencies and their strong variability.Comment: ApJ Letters, in press, (AAS)LaTeX, 3 figures, available at http://www2.mpifr-bonn.mpg.de/staff/hfalcke/publications.html#iiizw2 or in a few weeks at http://www.mpifr-bonn.mpg.de/staff/falcke/publications.html#iiizw

Frequency-dependent time delays for strong outbursts in selected blazars from the MetsÄhovi and UMRAO monitoring data bases – II

We analyse the radio light curves of the blazars 1308+326, 2223−052 and 2251+158 using University of Michigan Radio Observatory and MetsÄhovi Radio Observatory multifrequency monitoring data combined with high-resolution very long baseline interferometry (VLBI) observations in order to extract the properties of prominent outbursts. The outbursts are classified as ‘core’ and ‘jet’ events according to their behaviour at different frequencies and their associations with features appearing in the VLBI jet. We define the activity cycle for each blazar as the time interval between successive ‘core’ outbursts. The durations of the activity cycle derived in this way are ≥14 yr for 1308+326, ∼12 yr for 2223−052 and 12.4 ± 0.6 yr for 2251+158. We find an unusual frequency dependence for the time profiles for a major flare in 1308+326, which may provide evidence for acceleration or bending of the jet flow in the optically thick part of the emission region. Analysis of these activity cycles, combined with our earlier results, leads us to suggest that more luminous blazars possess shorter activity cycles, consistent with the accretion rates being higher in more powerful sources (relative to the Eddington rate).Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/72279/1/j.1365-2966.2007.12281.x.pd

A possible jet precession in the periodic quasar B0605-085

The quasar B0605-085 (OH 010) shows a hint for probable periodical variability in the radio total flux-density light curves. We study the possible periodicity of B0605-085 in the total flux-density, spectra and opacity changes in order to compare it with jet kinematics on parsec scales. We have analyzed archival total flux-density variability at ten frequencies (408 MHz, 4.8 GHz, 6.7 GHz, 8 GHz, 10.7 GHz, 14.5 GHz, 22 GHz, 37 GHz, 90 GHz, and 230 GHz) together with the archival high-resolution very long baseline interferometry data at 15 GHz from the MOJAVE monitoring campaign. Using the Fourier transform and discrete autocorrelation methods we have searched for periods in the total flux-density light curves. In addition, spectral evolution and changes of the opacity have been analyzed. We found a period in multi-frequency total flux-density light curves of 7.9+-0.5 yrs. Moreover, a quasi-stationary jet component C1 follows a prominent helical path on a similar time scale of 8 years. We have also found that the average instantaneous speeds of the jet components show a clear helical pattern along the jet with a characteristic scale of 3 mas. Taking into account average speeds of jet components, this scale corresponds to a time scale of about 7.7 years. Jet precession can explain the helical path of the quasi-stationary jet component C1 and the periodical modulation of the total flux-density light curves. We have fitted a precession model to the trajectory of the jet component C1, with a viewing angle phi=2.6+-2.2 degrees, aperture angle of the precession cone Omega=23.9+-1.9 degrees and fixed precession period (in the observers frame) P = 7.9 yrs.Comment: 14 pages, 16 figures, 5 tables, accepted for publication in A&