Hugh D. Aller

The inner jet of an active galactic nucleus as revealed by a radio-to-γ-ray outburst

Blazars are the most extreme active galactic nuclei. They possess oppositely directed plasma jets emanating at near light speeds from accreting supermassive black holes. According to theoretical models, such jets are propelled by magnetic fields twisted by differential rotation of the black hole’s accretion disk or inertial-frame-dragging ergosphere. The flow velocity increases outward along the jet in an acceleration and collimation zone containing a coiled magnetic field. Detailed observations of outbursts of electromagnetic radiation, for which blazars are famous, can potentially probe the zone. It has hitherto not been possible to either specify the location of the outbursts or verify the general picture of jet formation. Here we report sequences of high-resolution radio images and optical polarization measurements of the blazar BL Lacertae. The data reveal a bright feature in the jet that causes a double flare of radiation from optical frequencies to TeV γ-ray energies, as well as a delayed outburst at radio wavelengths. We conclude that the event starts in a region with a helical magnetic field that we identify with the acceleration and collimation zone predicted by the theories. The feature brightens again when it crosses a standing shock wave corresponding to the bright ‘core’ seen on the images.

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.

Multiwavelength Observations of a Dramatic High-Energy Flare in the Blazar 3C 279

The blazar 3C 279, one of the brightest identified extragalactic objects in the γ-ray sky, underwent a large (factor of ~10 in amplitude) flare in γ-rays toward the end of a 3 week pointing by Compton Gamma Ray Observatory (CGRO), in 1996 January-February. The flare peak represents the highest γ-ray intensity ever recorded for this object. During the high state, extremely rapid γ-ray variability was seen, including an increase of a factor of 2.6 in ~8 hr, which strengthens the case for relativistic beaming. Coordinated multifrequency observations were carried out with Rossi X-Ray Timing Explorer (RXTE), Advanced Satellite for Cosmology and Astrophysics (ASCA; or, Astro-D), Roentgen Satellite (ROSAT), and International Ultraviolet Explorer (IUE) and from many ground-based observatories, covering most accessible wavelengths. The well-sampled, simultaneous RXTE light curve shows an outburst of lower amplitude (factor of 3) well correlated with the γ-ray flare without any lag larger than the temporal resolution of ~1 day. The optical-UV light curves, which are not well sampled during the high-energy flare, exhibit more modest variations (factor of ~2) and a lower degree of correlation. The flux at millimetric wavelengths was near a historical maximum during the γ-ray flare peak, and there is a suggestion of a correlated decay. We present simultaneous spectral energy distributions of 3C 279 prior to and near to the flare peak. The γ-rays vary by more than the square of the observed IR-optical flux change, which poses some problems for specific blazar emission models. The synchrotron self-Compton (SSC) model would require that the largest synchrotron variability occurred in the mostly unobserved submillimeter/far-infrared region. Alternatively, a large variation in the external photon field could occur over a timescale of a few days. This occurs naturally in the mirror model, wherein the flaring region in the jet photoionizes nearby broad emission line clouds, which, in turn, provide soft external photons that are Comptonized to γ-ray energies.

Sub-Milliarcsecond Imaging of Quasars and Active Galactic Nuclei. IV. Fine-Scale Structure

We have examined the compact structure in 250 flat-spectrum extragalactic radio sources using interferometric fringe visibilities obtained with the Very Long Baseline Array (VLBA) at 15 GHz. With projected baselines out to 440 Mλ, we are able to investigate source structure on typical angular scales as small as 0.05 mas. This scale is similar to the resolution of the VLBI Space Observatory Programme data obtained on longer baselines at a lower frequency and with somewhat poorer accuracy. For 171 sources in our sample, more than half of the total flux density seen by the VLBA remains unresolved on the longest baselines. There are 163 sources in our list with a median correlated flux density at 15 GHz in excess of 0.5 Jy on the longest baselines; these will be useful as fringe finders for short-wavelength VLBA observations. The total flux densities recovered in the VLBA images at 15 GHz are generally close to the values measured around the same epoch at the same frequency with the RATAN-600 and University of Michigan Radio Astronomy Observatory telescopes. We have modeled the core of each source with an elliptical Gaussian component. For about 60% of the sources, we have at least one observation in which the core component appears unresolved (generally smaller than 0.05 mas) in one direction, usually transverse to the direction into which the jet extends. BL Lac objects are on average more compact than quasars, while active galaxies are on average less compact. Also, in an active galaxy the sub-milliarcsecond core component tends to be less dominant. Intraday variability (IDV) sources typically have a more compact, more core-dominated structure on sub-milliarcsecond scales than non-IDV sources, and sources with a greater amplitude of intraday variations tend to have a greater unresolved VLBA flux density. The objects known to be GeV gamma-ray-loud appear to have a more compact VLBA structure than the other sources in our sample. This suggests that the mechanisms for the production of gamma-ray emission and for the generation of compact radio synchrotron–emitting features are related. The brightness temperature estimates and lower limits for the cores in our sample typically range between 1011 and 1013 K, but they extend up to 5 × 1013 K, apparently in excess of the equipartition brightness temperature or the inverse Compton limit for stationary synchrotron sources. The largest component speeds are observed in radio sources with high observed brightness temperatures, as would be expected from relativistic beaming. Longer baselines, which can be obtained by space VLBI observations, will be needed to resolve the most compact high brightness temperature regions in these sources.

A MULTIWAVELENGTH VIEW OF THE TeV BLAZAR MARKARIAN 421: CORRELATED VARIABILITY, FLARING, AND SPECTRAL EVOLUTION

We report results from an intensive multiwavelength monitoring campaign on the TeV blazar Mrk 421 over the period of 2003-2004. The source was observed simultaneously at TeV energies with the Whipple 10 m telescope and at X-ray energies with the Rossi X-Ray Timing Explorer (RXTE) during each clear night within the Whipple observing windows. Supporting observations were also frequently carried out at optical and radio wavelengths to provide simultaneous or contemporaneous coverages. The large amount of simultaneous data has allowed us to examine the variability of Mrk 421 in detail, including cross-band correlation and broadband spectral variability, over a wide range of flux. The variabilities are generally correlated between the X-ray and gamma-ray bands, although the correlation appears to be fairly loose. The light curves show the presence of flares with varying amplitudes on a wide range of timescales at both X-ray and TeV energies. Of particular interest is the presence of TeV flares that have no coincident counterparts at longer wavelengths, because the phenomenon seems difficult to understand in the context of the proposed emission models for TeV blazars. We have also found that the TeV flux reached its peak days before the X-ray flux did during a giant flare (or outburst) in 2004 (with the peak flux reaching ~135 mcrab in X-rays, as seen by the RXTE ASM, and ~3 crab in gamma rays). Such a difference in the development of the flare presents a further challenge to both the leptonic and hadronic emission models. Mrk 421 varied much less at optical and radio wavelengths. Surprisingly, the normalized variability amplitude in the optical seems to be comparable to that in the radio, perhaps suggesting the presence of different populations of emitting electrons in the jet. The spectral energy distribution of Mrk 421 is seen to vary with flux, with the two characteristic peaks moving toward higher energies at higher fluxes. We have failed to fit the measured spectral energy distributions (SEDs) with a one-zone synchrotron self-Compton model; introducing additional zones greatly improves the fits. We have derived constraints on the physical properties of the X-ray/gamma-ray flaring regions from the observed variability (and SED) of the source. The implications of the results are discussed.

Optical and radio behaviour of the BL Lacertae object 0716+714

Eight optical and four radio observatories have been intensively monitoring the BL Lac object 0716+714 in the last years: 4854 data points have been collected in the UBVRI bands since 1994, while radio light curves extend back to 1978. Many of these data, which all together constitute the widest optical and radio database available on this object, are presented here for the first time. Four major optical outbursts were observed at the beginning of 1995, in late 1997, at the end of 2000, and in fall 2001. In particular, an exceptional brightening of 2.3 mag in 9 days was detected in the R band just before the BeppoSAX pointing of October 30, 2000. A big radio outburst lasted from early 1998 to the end of 1999. The long-term trend shown by the optical light curves seems to vary with a characteristic time scale of about 3.3 years, while a longer period of 5.5–6 years seems to characterize the radio long-term variations. In general, optical colour indices are only weakly correlated with brightness; a clear spectral steepening trend was observed during at least one long-lasting dimming phase. Moreover, the optical spectrum became steeper after

Multiepoch Very Long Baseline Array Observations of EGRET-detected Quasars and BL Lacertae Objects: Connection between Superluminal Ejections and Gamma-Ray Flares in Blazars

\rm JD \sim 2\,451\,000

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

, the change occurring in the decaying phase of the late-1997 outburst. The radio flux behaviour at different frequencies is similar, but the flux variation amplitude decreases with increasing wavelength. The radio spectral index varies with brightness (harder when brighter), but the radio fluxes seem to be the sum of two different-spectrum contributions: a steady base level and a harder-spectrum variable component. Once the base level is removed, the radio variations appear as essentially achromatic, similarly to the optical behaviour. Flux variations at the higher radio frequencies lead the lower-frequency ones with week–month time scales. The behaviour of the optical and radio light curves is quite different, the broad radio outbursts not corresponding in time to the faster optical ones and the cross-correlation analysis indicating only weak correlation with long time lags. However, minor radio flux enhancements simultaneous with the major optical flares can be recognized, which may imply that the mechanism producing the strong flux increases in the optical band also marginally affects the radio one. On the contrary, the process responsible for the big radio outbursts does not seem to affect the optical emission.

Multiepoch multiwavelength spectra and models for blazar 3C 279

We examine the coincidence of times of high γ-ray flux and ejections of superluminal components from the core in EGRET blazars based on a Very Long Baseline Array (VLBA) monitoring program at 22 and 43 GHz from 1993 November to 1997 July. In 23 cases of γ-ray flares for which sufficient VLBA data exist, 10 of the flares (in eight objects) fall within 1 σ uncertainties of the extrapolated epoch of zero separation from the core of a superluminal radio component. In each of two sources (0528+134 and 1730-130), two successive γ-ray flares were followed by the appearance of new superluminal components. We carried out statistical simulations that show that if the number of coincidences is ≥10, the radio and γ-ray events are associated with each other at greater than 99.999% confidence. Our analysis of the observed behavior, including variability of the polarized radio flux, of the sources before, during, and after the γ-ray flares suggests that the γ-ray events occur in the superluminal radio knots. This implies that the γ-ray flares are caused by inverse Compton scattering by relativistic electrons in the parsec-scale regions of the jet rather than closer to the central engine.

Optical and radio variability of the BL Lacertae object AO 0235+16: A possible 5-6 year periodicity

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.