M. F. Aller

MOJAVE: MONITORING OF JETS IN ACTIVE GALACTIC NUCLEI WITH VLBA EXPERIMENTS. V. MULTI-EPOCH VLBA IMAGES

We present images from a long-term program (MOJAVE: Monitoring of Jets in active galactic nuclei (AGNs) with VLBA Experiments) to survey the structure and evolution of parsec-scale jet phenomena associated with bright radio-loud active galaxies in the northern sky. The observations consist of 2424 15 GHz Very Long Baseline Array (VLBA) images of a complete flux-density-limited sample of 135 AGNs above declination –20°, spanning the period 1994 August to 2007 September. These data were acquired as part of the MOJAVE and 2 cm Survey programs, and from the VLBA archive. The sample-selection criteria are based on multi-epoch parsec-scale (VLBA) flux density, and heavily favor highly variable and compact blazars. The sample includes nearly all the most prominent blazars in the northern sky, and is well suited for statistical analysis and comparison with studies at other wavelengths. Our multi-epoch and stacked-epoch images show 94% of the sample to have apparent one-sided jet morphologies, most likely due to the effects of relativistic beaming. Of the remaining sources, five have two-sided parsec-scale jets, and three are effectively unresolved by the VLBA at 15 GHz, with essentially all of the flux density contained within a few tenths of a milliarcsecond.

The University of Michigan radio astronomy data base. I. Structure function analysis and the relation between BL Lacertae objects and quasi-stellar objects

A structure function analysis on a well-observed group of variable extragalactic sources, monitored in total and polarized flux in the centimeter waveband since 1965, has been performed. Total flux structure functions of BL Lac objects and QSOs exhibit very similar power-law slopes, suggesting that the same processes are responsible for the variability in both types of source. In view of the recent success of the shocked-jet model for radio source variability, it is argued that the current work provides circumstantial evidence for the ubiquity of such shocked flows

Simultaneous Planck, Swift, and Fermi observations of X-ray and γ-ray selected blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and -ray bands, and we compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi Large Area Telescope (LAT), whereas 30 to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the -ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, withh i 0 up to about 70 GHz, above which it steepens toh i 0:65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency ( S ) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples withh S i = 10 13:1 0:1 Hz, while the mean inverse-Compton peak frequency,h IC i, ranges from 10 21 to 10 22 Hz. The distributions of S and of IC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with -ray selected blazars peaking at 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection e ect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the -ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and S predicted by the blazar sequence.

The Relation Between AGN Gamma-Ray Emission and Parsec-Scale Radio Jets

We have compared the radio emission from a sample of parsec-scale AGN jets as measured by the VLBA at 15 GHz, with their associated gamma-ray properties that are reported in the Fermi LAT 3-month bright source list. We find in our radio selected sample that the gamma-ray photon flux correlates well with the quasi-simultaneously measured compact radio flux density. The LAT-detected jets in our radio-selected complete sample generally have higher compact radio flux densities, and their parsec-scale cores are brighter (i.e., have higher brightness temperature) than the jets in the non-LAT detected objects. This suggests that the jets of bright gamma-ray AGN have preferentially higher Doppler-boosting factors. In addition, jets of the LAT-detected AGN tend to be in a more active radio state, when quasi-simultaneous data are used. This result becomes more pronounced for confirmed gamma-ray flaring sources. We identify the parsec-scale radio core as a likely location for both the gamma-ray and radio flares, which appear within typical timescales of up to a few months of each other.We have compared the radio emission from a sample of parsec-scale active galactic nucleus (AGN) jets as measured by the VLBA at 15 GHz, with their associated γ-ray properties that are reported in the Fermi Large Area Telescope (LAT) three month bright source list. We find in our radio-selected sample that the γ-ray photon flux correlates well with the quasi-simultaneously measured compact radio flux density. The LAT-detected jets in our radio-selected complete sample generally have higher compact radio flux densities, and their parsec-scale cores are brighter (i.e., have higher brightness temperature) than the jets in the LAT nondetected objects. This suggests that the jets of bright γ-ray AGN have preferentially higher Doppler-boosting factors. In addition, AGN jets tend to be found in a more active radio state within several months from LAT-detection of their strong γ-ray emission. This result becomes more pronounced for confirmed γ-ray flaring sources. We identify the parsec-scale radio core as a likely location for both the γ-ray and radio flares, which appear within typical timescales of up to a few months of each other.

Radio Monitoring of OJ 287 and Binary Black Hole Models for Periodic Outbursts

The BL Lac-type active galaxy OJ 287 exhibits a 12 year periodicity with a double-peaked maxima in its optical flux variations. Several models sought to explain this periodicity, the first one firmly established in any active galactic nucleus (AGN), as a result of the orbital motion of a pair of supermassive black holes. In one class of models the orientation of the jets changes in a regular manner, and the optical flaring is due to a consequent increase in the Doppler boosting factor. In another class of models the optical flaring reflects a true increase in luminosity, either due to an enhanced accretion during the pericenter passage or due to a collision between the secondary black hole and the accretion disk of the primary black hole. However, these models have been based solely on the optical data. Here we consider the full radio flux density monitoring data between 8 and 90 GHz from the Michigan, Metsahovi, and Swedish-ESO Submillimeter Telescope AGN monitoring programs. We find that the radio flux density and polarization data, as well as the optical polarization data, indicate that the first of the two optical peaks is a thermal flare occurring in the vicinity of the black hole and the accretion disk, while the second one is a synchrotron flare originating in a shocked region down the jet. None of the proposed binary black hole models for OJ 287 offers satisfactory explanations for these observations. We suggest a new scenario, in which a secondary black hole penetrates the accretion disk of the primary during the pericenter passage, causing a thermal flare visible only in the optical regime. The pericenter passage enhances accretion into the primary black hole, leading to increased jet flow and formation of shocks down the jet. These become visible as standard radio and optical synchrotron flares roughly a year after the pericenter passage and are identified with the second optical peaks. In addition to explaining the radio and the optical data, our model eliminates the need for a strong precession of the binary and for an ultramassive (≥1010 M☉) primary black hole. If our interpretation is correct, the next periodic optical flare, a thermal one, should occur around 2006 September 25. Nonthermal, simultaneous optical and radio flares should follow about a year later.

Multifrequency monitoring of the blazar 0716+714 during the GASP-WEBT-AGILE campaign of 2007

Aims. Since the CGRO operation in 1991–2000, one of the primary unresolved questions about the blazar γ -ray emission has been its possible correlation with the low-energy (in particular optical) emission. To help answer this problem, the Whole Earth Blazar Telescope (WEBT) consortium has organized the GLAST-AGILE Support Program (GASP) to provide the optical-to-radio monitoring data to be compared with the γ -ray detections by the AGILE and GLAST satellites. This new WEBT project started in early September 2007, just before a strong γ -ray detection of 0716+714 by AGILE. Methods. We present the GASP-WEBT optical and radio light curves of this blazar obtained in July–November 2007, about various AGILE pointings at the source. We construct NIR-to-UV spectral energy distributions (SEDs), by assembling GASP-WEBT data together with UV data from the Swift ToO observations of late October. Results. We observe a contemporaneous optical-radio outburst, which is a rare and interesting phenomenon in blazars. The shape of the SEDs during the outburst appears peculiarly wavy because of an optical excess and a UV drop-and-rise. The optical light curve is well sampled during the AGILE pointings, showing prominent and sharp flares. A future cross-correlation analysis of the optical and AGILE data will shed light on the expected relationship between these flares and the γ -ray events.

Results of WEBT, VLBA and RXTE monitoring of 3C 279 during 2006-2007 ⋆

Context. The quasar 3C 279 is among the most extreme blazars in terms of luminosity and variability of flux at all wavebands. Its vari ations in flux and polarization are quite complex and therefore require intensive monitoring observations at multiple wavebands to characterise and interpret the observed changes. Aims. In this paper, we present radio-to-optical data taken by the WEBT, supplemented by our VLBA and RXTE observations, of 3C 279. Our goal is to use this extensive database to draw inferences regarding the physics of the relativistic jet. Methods. We assemble multifrequency light curves with data from 30 ground-based observatories and the space-based instruments SWIFT (UVOT) and RXTE, along with linear polarization vs. time in the optical R band. In addition, we present a sequence of 22 images (with polarization vectors) at 43 GHz at resolution 0.15 milliarcse c, obtained with the VLBA. We analyse the light curves and polarization, as well as the spectral energy distributions at different epochs, corresponding to different brightness states. Results. We find that the IR-optical-UV continuum spectrum of the vari able component corresponds to a power law with a constant slope of−1.6, while in the 2.4‐10 keV X-ray band it varies in slope from−1.1 to−1.6. The steepest X-ray spectrum occurs at a flux minimum. Durin g a decline in flux from maximum in late 2006, the optical and 43 GHz core po larization vectors rotate by∼ 300 ◦ . Conclusions. The continuum spectrum agrees with steady injection of relativistic electrons with a power-law energy distribution of slope−3.2 that is steepened to−4.2 at high energies by radiative losses. The X-ray emission at flux minimum comes most likely from a new component that starts in an upstream section of the jet where inverse Compton scattering of seed photons from outside the jet is important. The rotation of the polarization vector implies that the jet contains a helical magnetic field that extends ∼ 20 pc past the 43 GHz core.

Pearson-readhead survey sources: properties of the centimeter-wavelength flux and polarization of a complete radio sample

We have carried out a monitoring program to observe systematically the integrated total flux density and linear polarization of sources in the Pearson-Readhead sample of extragalactic objects at 4.8, 8.0, and 14.5 GHz. The large data base obtained has been used to investigate the statistical behavior of the flux and polarization in this well-defined sample and to search for dependence of these properties on other radio and optical properties

Pearson-Readhead Survey Sources. II. The Long-Term Centimeter-Band Total Flux and Linear Polarization Properties of a Complete Radio Sample

Using UMRAO centimeter-band total flux density and linear polarization monitoring observations of the complete Pearson-Readhead extragalactic source sample obtained between 1984 August and 2001 March, we identify the range of variability in extragalactic objects as functions of optical and radio morphological classification and relate total flux density variations to structural changes in published coeval VLBI maps in selected objects. As expected, variability is common in flat- or inverted-spectrum ( ≤ 0.5) core-dominated QSOs and BL Lac objects. Unexpectedly, we find flux variations in several steep-spectrum sample members, including the commonly adopted flux standard 3C 147. Such variations are characteristically several-year rises or declines or infrequent outbursts, requiring long-term observations for detection: we attribute them to the brightening of weak core components, a change that is suppressed by contributions from extended structure in all but the strongest events, and identify a wavelength dependence for the amplitude of this variability consistent with the presence of opacity in some portions of the jet flow. One morphological class of steep-spectrum objects, the compact symmetric objects (CSOs), characteristically shows only low-level variability. We examine the statistical relation between fractional polarization and radio class based on the data at 14.5 and 4.8 GHz. The blazars typically exhibit flat-to-inverted polarization spectra, a behavior attributed to opacity effects. Among the steep-spectrum objects, the lobe-dominated FR I galaxies have steep fractional polarization spectra, while the FR II galaxies exhibit fractional polarization spectra ranging from inverted to steep, with no identifiable common property that accounts for the range in behavior. For the CSO/gigahertz-peaked spectrum sources, we verify that the fractional polarizations at 4.8 GHz are only of the order of a few tenths of a percent, but at 14.5 GHz we find significantly higher polarizations, ranging from 1% to 3%; this frequency dependence supports a scenario invoking Faraday depolarization by a circumnuclear torus. We have identified preferred orientations of the electric vector of the polarized emission (EVPA) at 14.5 and 4.8 GHz in roughly half of the objects and compared these with orientations of the flow direction indicated by VLBI morphology. When comparing the distributions of the orientation offsets for the BL Lac objects and the QSOs, we find differences in both range and mean value, in support of intrinsic class differences. In the shock-in-jet scenario, we attribute this to the allowed range of obliquities of shocks developing in the flow relative to the flow direction: in the BL Lac objects the shocks are nearly transverse to the flow direction, while in the QSOs they include a broader range of obliquities and can be at large angles to it. The fact that we find long-term stability in EVPA over many events implies that a dominant magnetic field orientation persists; in the core-dominated objects, with small contribution from the underlying quiescent jet, this plausibly suggests that the magnetic field has a long-term memory, with subsequent shock events exhibiting similar EVPA orientation, or, alternatively, the presence of a standing shock in the core. We have looked for systematic, monotonic changes in EVPA, which might be expected in the emission from a precessing jet, a model currently invoked for some AGNs; none were identified. Further, we carried out a Scargle periodogram analysis of the total flux density observations, but found no strong evidence for periodicity in any of the sample sources. The only well-established case in support of both jet precession and periodic variability remains the non-sample member OJ 287.

The WEBT campaigns on BL Lacertae: Time and cross-correlation analysis of optical and radio light curves 1968-2003

The Whole Earth Blazar Telescope (WEBT) collaboration has collected a large amount of optical and radio data on BL Lacertae in the last years, which, when added to literature data, allow to construct well-sampled light curves of the source from 1968 to the end of 2003. These optical and radio data are here analysed with three statistical methods designed for unevenly-sampled data trains in order to search for possible periodicities. While the main radio outbursts repeat every ∼8 years, with a possible progressive stretching of the period, the evidence of an optical periodicity is much less clear. Radio light curves from 4.8 to 37 GHz are well correlated, with variations at the higher frequencies leading the lower-frequency ones by a few weeks for contiguous bands, up to a few months when considering the largest frequency separations. The radio behaviour reveals the presence of two different components, the softer-spectrum one constituting the bulk of the radio emission. On the other hand, the harder component shows itself as radio events which appear enhanced at the higher frequencies and seem to have optical counterparts. Cross-correlation between the optical light curve and radio hardness ratios indicates a radio time delay of more than 3 months. Thus, our analysis suggests a scenario where flux variations propagate towards less and less opaque regions, giving rise to related optical and hard radio events and, in more extended zones, to soft events apparently uncorrelated with the former ones.