I. M. McHardy

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.

Active galactic nuclei as scaled-up Galactic black holes.

A long-standing question is whether active galactic nuclei (AGN) vary like Galactic black hole systems when appropriately scaled up by mass. If so, we can then determine how AGN should behave on cosmological timescales by studying the brighter and much faster varying Galactic systems. As X-ray emission is produced very close to the black holes, it provides one of the best diagnostics of their behaviour. A characteristic timescale—which potentially could tell us about the mass of the black hole—is found in the X-ray variations from both AGN and Galactic black holes, but whether it is physically meaningful to compare the two has been questioned. Here we report that, after correcting for variations in the accretion rate, the timescales can be physically linked, revealing that the accretion process is exactly the same for small and large black holes. Strong support for this linkage comes, perhaps surprisingly, from the permitted optical emission lines in AGN whose widths (in both broad-line AGN and narrow-emission-line Seyfert 1 galaxies) correlate strongly with the characteristic X-ray timescale, exactly as expected from the AGN black hole masses and accretion rates. So AGN really are just scaled-up Galactic black holes.

Non-linear X-ray variability in X-ray binaries and active galaxies

We show that the rms-flux relation recently discovered in the X-ray light curves of Active Galactic Nuclei (AGN) and X-ray binaries (XRBs) implies that the light curves have a formally non-linear, exponential form, provided the rms-flux relation applies to variations on all time-scales (as it appears to). This phenomenological model implies that stationary data will have a log normal flux distribution. We confirm this result using an observation of Cyg X-1, and further demonstrate that our model predicts the existence of the powerful millisecond flares observed in Cyg X-1 in the low/hardstate, and explains the general shape and amplitude of the bicoherence spectrum in that source. Our model predicts that the most variable light curves will show the most extreme non-linearity. This result can naturally explain the apparent non-linear variability observed in some highly variable Narrow Line Seyfert 1 (NLS1) galaxies, as well as the low states observed on long time-scales in the NLS1 NGC 4051, as being nothing more than extreme manifestations of the same variability process that is observed in XRBs and less variable AGN. That variability process must be multiplicative (with variations coupled together on all time-scales) and cannot be additive (such as shot-noise), or related to self-organised criticality, or result from completely independent variations in many separate emitting regions. Successful models for variability must reproduce the observed rms-flux relation and non-linear behaviour, which are more fundamental characteristics of the variability process than the power spectrum or spectral-timing properties. Models where X-ray variability is driven by accretion rate variations produced at different radii remain the most promising.

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.

The flux-dependent amplitude of broadband noise variability in X-ray binaries and active galaxies

Standard shot-noise models, which seek to explain the broadband noise variability that characterizes the X-ray light curves of X-ray binaries and active galaxies, predict that the power spectrum of the X-ray light curve is stationary (i.e. constant amplitude and shape) on short time-scales. We show that the broadband noise power spectra of the black hole candidate Cyg X-1 and the accreting millisecond pulsar SAX J1808.4-3658 are intrinsically non-stationary, in that rms variability scales linearly with flux. Flux-selected power spectra confirm that this effect is due to changes in power-spectral amplitude and not shape. The light curves of three Seyfert galaxies are also consistent with a linear relationship between rms variability and flux, suggesting that it is an intrinsic feature of the broadband noise variability in compact accreting systems over more than six decades of central object mass. The rms variability responds to flux variations on all measured time-scales, raising fundamental difficulties for shot-noise models which seek to explain this result by invoking variations in the shot parameters. We suggest that models should be explored where the longest time-scale variations are fundamental and precede the variations on shorter time-scales. Possible models which can explain the linear rms-flux relation include the fractal break-up of large coronal flares, or the propagation of fluctuations in mass accretion rate through the accretion disc. The linear relationship between rms variability and flux in Cyg X-1 and SAX J1808.4-3658 is offset on the flux axis, suggesting the presence of a second, constant-flux component to the light curve which contributes ~25 per cent of the total flux. The spectrum of this constant component is similar to the total spectrum, suggesting that it may correspond to quiet, non-varying regions in the X-ray emitting corona.

X-RAY AND OPTICAL VARIABILITY IN NGC 4051 AND THE NATURE OF NARROW-LINE SEYFERT 1 GALAXIES

We report on the results of a three-year program of coordinated X-ray and optical monitoring of the narrow-line Seyfert 1 galaxy NGC 4051. The rapid continuum variations observed in the X-ray spectra are not detected in the optical, although the time-averaged X-ray and optical continuum fluxes are well correlated. Variations in the flux of the broad Hβ line are found to lag behind the optical continuum variations by 6 days (with an uncertainty of 2-3 days), and combining this with the line width yields a virial mass estimate of ~1.1 × 106 M☉, at the very low end of the distribution of active galactic nucleus masses measured by line reverberation. Strong variability of He II λ4686 is also detected, and the response time measured is similar to that of Hβ but with a much larger uncertainty. The He II λ4686 line is almost 5 times broader than Hβ, and it is strongly blueward asymmetric, as are the high-ionization UV lines recorded in archival spectra of NGC 4051. The data are consistent with the Balmer lines arising in a low-to-moderate-inclination disklike configuration and the high-ionization lines arising in an outflowing wind, of which we observe preferentially the near side. Previous observations of the narrow-line region morphology of this source suggest that the system is inclined by ~50°, and if this is applicable to the broad Hβ-emitting region, a central mass of ~1.4 × 106 M☉ can be inferred. During the third year of monitoring, both the X-ray continuum and the He II λ4686 line went into extremely low states, although the optical continuum and the Hβ broad line were both still present and variable. We suggest that the inner part of the accretion disk may have gone into an advection-dominated state, yielding little radiation from the hotter inner disk.

The star formation history of the Universe as revealed by deep radio observations

Discerning the exact nature of the sub-mJy radio population has been historically difficult due to the low luminosity of these sources at most wavelengths. Using deep ground based optical follow-up and observations from the Spitzer Space Telescope we are able to disentangle the radio-selected active galactic nuclei (AGN) and star-forming galaxy (SFG) populations for the first time in a deep multifrequency VLA/MERLIN Survey of the 13^H XMM Newton/Chandra Deep Field. The discrimination diagnostics include radio morphology, radio spectral index, radio/near-infrared (near-IR) and mid-IR/radio flux density ratios. We are nowable to calculate the extragalactic Euclidean normalized source counts separately for AGN and SFGs. We find that while SFGs dominate at the faintest flux densities and account for the majority of the upturn in the counts, AGN still make up around one quarter of the counts at ∼50 μJy (1.4 GHz). Using radio luminosity as an unobscured star formation rate (SFR) measure we are then able to examine the comoving SFR density of the Universe up to z = 3 which agrees well with measures at other wavelengths. We find a rough correlation of SFR with stellar mass for both the sample presented here and a sample of local radio-selected SFGs from the 6df-NVSS survey. This work also confirms the existence of, and provides alternative evidence for, the evolution of distribution of star formation by galaxy mass: ‘downsizing’. As both these samples are SFR-selected, this result suggests that there is a maximum SFR for a given galaxy that depends linearly on its stellar mass. The low ‘characteristic times’ (inverse specific SFR) of the SFGs in our sample are similar to those of the 6dF-NVSS sample, implying that most of these sources are in a current phase of enhanced star formation.

Flaring Behavior of the Quasar 3C 454.3 Across the Electromagnetic Spectrum

We analyze the behavior of the parsec-scale jet of the quasar 3C 454.3 during pronounced flaring in 2005-2008. Three major disturbances propagated down the jet along different trajectories with Lorentz factors Γ > 10. The disturbances show a clear connection with millimeter-wave outbursts, in 2005 May/June, 2007 July, and 2007 December. High-amplitude optical events in the R-band light curve precede peaks of the millimeter-wave outbursts by 15-50 days. Each optical outburst is accompanied by an increase in X-ray activity. We associate the optical outbursts with propagation of the superluminal knots and derive the location of sites of energy dissipation in the form of radiation. The most prominent and long lasting of these, in 2005 May, occurred closer to the black hole, while the outbursts with a shorter duration in 2005 autumn and in 2007 might be connected with the passage of a disturbance through the millimeter-wave core of the jet. The optical outbursts, which coincide with the passage of superluminal radio knots through the core, are accompanied by systematic rotation of the position angle of optical linear polarization. Such rotation appears to be a common feature during the early stages of flares in blazars. We find correlations between optical variations and those at X-ray and γ-ray energies. We conclude that the emergence of a superluminal knot from the core yields a series of optical and high-energy outbursts, and that the millimeter-wave core lies at the end of the jets acceleration and collimation zone. We infer that the X-ray emission is produced via inverse Compton scattering by relativistic electrons of photons both from within the jet (synchrotron self-Compton) and external to the jet (external Compton, or EC); which one dominates depends on the physical parameters of the jet. A broken power-law model of the γ-ray spectrum reflects a steepening of the synchrotron emission spectrum from near-IR to soft UV wavelengths. We propose that the γ-ray emission is dominated by the EC mechanism, with the sheath of the jet supplying seed photons for γ-ray events that occur near the millimeter-wave core.

Multiepoch multiwavelength spectra and models for blazar 3C 279

Of the blazars detected by EGRET in GeV γ-rays, 3C 279 is not only the best observed by EGRET but also one of the best monitored at lower frequencies. We have assembled 11 spectra, from GHz radio through GeV γ-rays, from the time intervals of EGRET observations. Although some of the data have appeared in previous publications, most are new, including data taken during the high states in early 1999 and early 2000. All of the spectra show substantial γ-ray contribution to the total luminosity of the object; in a high state, the γ-ray luminosity dominates over that at all other frequencies by a factor of more than 10. There is no clear pattern of time correlation; different bands do not always rise and fall together, even in the optical, X-ray, and γ-ray bands. The spectra are modeled using a leptonic jet, with combined synchrotron self-Compton plus external Compton γ-ray production. Spectral variability of 3C 279 is consistent with variations of the bulk Lorentz factor of the jet, accompanied by changes in the spectral shape of the electron distribution. Our modeling results are consistent with the UV spectrum of 3C 279 being dominated by accretion disk radiation during times of low γ-ray intensity.

Synchrotron Self-Compton Model for Rapid Nonthermal Flares in Blazars with Frequency-dependent Time Lags

We model rapid variability of multifrequency emission from blazars occurring across the electromagnetic spectrum (from radio to gamma rays). Lower energy emission is produced by the synchrotron mechanism, whereas higher energy emission is due to inverse Compton scattering of the synchrotron emission. We take into account energy stratification established by particle acceleration at shock fronts and energy losses due to synchrotron emission. We also consider the effect of light-travel delays for the synchrotron emission that supplies the seed photons for inverse Compton scattering. The production of a flare is caused by the collision between a relativistic shock wave and a stationary feature in the jet (e.g., a Mach disk). The collision leads to the formation of forward and reverse shocks, which confine two contiguous emission regions resulting in complex profiles of simulated flares. Simulations of multifrequency flares indicate that relative delays between the inverse Compton flares and their synchrotron counterparts are dominated by energy stratification and geometry of the emitting regions, resulting in both negative and positive time delays depending on the frequency of observation. Light-travel effects of the seed photons may lead to a noticeable delay of the inverse Compton emission with respect to synchrotron variability if the line of sight is almost perfectly aligned with the jet. We apply the model to a flare in 3C 273 and derive the properties of shocked plasma responsible for the flare. We show that the pronounced negative time delay between the X-ray and IR light curves (X-rays peak after the maximum in the synchrotron emission) can be accounted for if both forward and reverse shocks are considered.