A. Sillanpää

OJ 287 - binary pair of supermassive black holes

A historical light curve of the BL Lacertae object OJ 287 is constructed in the optical V band using observations between 1890 and the present. The curve exhibits periodic outbursts at intervals of 11.6 yr or 9 yr. The individual outbursts show an initial maximum with decaying submaxima at intervals of about 1 month. It is proposed that these characteristics are caused by a binary in the nucleus of OJ 287, with the binary members having an accretion disk around them. The outbursts at 9-yr intervals would be caused by the tidal action of the companion on the disk of the larger black hole. A computer simulation is used to demonstrate that the inflow into the center of such a disk during the repeated periastron passages of a companion will produce an outburst similar to the ones observed. Based on indirect evidence, probable masses of 5 billion solar and 20 million solar are determined for the two objects. 38 references.

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.

A massive binary black-hole system in OJ 287 and a test of general relativity

Tests of Einstein’s general theory of relativity have mostly been carried out in weak gravitational fields where the space-time curvature effects are first-order deviations from Newton’s theory. Binary pulsars provide a means of probing the strong gravitational field around a neutron star, but strong-field effects may be best tested in systems containing black holes. Here we report such a test in a close binary system of two candidate black holes in the quasar OJ 287. This quasar shows quasi-periodic optical outbursts at 12-year intervals, with two outburst peaks per interval. The latest outburst occurred in September 2007, within a day of the time predicted by the binary black-hole model and general relativity. The observations confirm the binary nature of the system and also provide evidence for the loss of orbital energy in agreement (within 10 per cent) with the emission of gravitational waves from the system. In the absence of gravitational wave emission the outburst would have happened 20 days later.

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.

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.

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

The BL Lacertae object AO 0235+16 is well known for its extreme optical and radio variability. New optical and radio data have been collected in the last four years by a wide international collaboration, which conrm the intense activity of this source: on the long term, overall variations of 5 mag in the R band and up to a factor 18 in the radio fluxes were detected, while short-term variability up to 0:5 mag in a few hours and 1: 3m ag in one day was observed in the optical band. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) rst-light campaign organized in November 1997, involving a dozen optical observatories. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection and Discrete Correlation Function (DCF) analysis. On the long term, the autocorrelation function of the optical data shows a double-peaked maximum at 4100{4200 days (11:2{11:5 years), while a double-peaked maximum at 3900{4200 days (10:7{11:5 years) is visible in the radio autocorrelation functions. The existence of this similar characteristic time scale of variability in the two bands is by itself an indication of optical-radio correlation. A further analysis by means of Discrete Fourier Transform (DFT) technique and folded light curves reveals that the major radio outbursts repeat quasi-regularly with a periodicity of5:7 years, i.e. half the above time scale. This period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them. Visual inspection and DCF analysis of the optical and radio light curves then reveal that in some cases optical outbursts seem to be simultaneous with radio ones, but in other cases they lead the radio events. Moreover, a deep inspection of the radio light curves suggests that in at least two occasions (the 1992{1993 and 1998 outbursts) flux variations at the higher frequencies may have led those at the lower ones.

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.

Detection of the host galaxy of S5 0716+714

We have acquired a deep i-band image of the BL Lacertae object S5 0716+714 while the target was in an low optical state. Due to the faintness of the nucleus, we were able to detect the underlying host galaxy. The host galaxy is measured to have an I-band magnitude of 17.5 ± 0. 5a nd an effective radius of (2.7 ± 0.8) arcsec. Using the host galaxy as a “standard candle”, we derive z = 0.31 ± 0.08 (1σ error) for the host galaxy of S5 0716+714. This redshift is consistent with the redshift z = 0.26 determined by spectroscopy for 3 galaxies close to S5 0716+714. The effective radius at z = 0.31 would be 12 ± 4 kpc, which is consistent with values obtained for BL Lac host galaxies. An optical spectrum acquired during the same epoch shows no identifiable spectral lines.

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.

R-band imaging of the host galaxies of RGB BL Lacertae objects

We present optical R-band images of 100 BL Lacertae objects in the ROSAT-Green Bank (RGB) sample obtained at the Nordic Optical Telescope (NOT). For 62 objects we could not find any previously published images in the literature, of these 35 are new BL Lacs discovered during the RGB BL Lac survey. We were able to resolve the host galaxy in 62% of the objects. We analyze the host galaxies using a two-dimensional fitting procedure and tabulate the best-fit core magnitude and host galaxy parameters (magnitude, effective radius, ellipticity, position angle and shape parameter β) for each object. With two exceptions, all objects are better fit by a model representing an elliptical galaxy (β = 0.25) than by a disk galaxy model (β = 1.0). In the two exceptions the disk host fit is only marginally better than the elliptical host fit, however, and we do not find a single clear example of a disk-dominated host galaxy. The host galaxies have an average brightness MR = −23.9 ± 0.8 and average effective radius reff = 13.2 ± 0.8 kpc. Their bulk properties are indistinguishable from normal elliptical and radio galaxies. The median β is 0.18, well below the value 0.25 that is normally used to describe ellipticals. However, the difference may be caused by as election effect due to the dependence of β on host galaxy luminosity.