Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://jet.sao.ru/hq/marat/publ/articles/aa2567-04.pdf Äàòà èçìåíåíèÿ: Wed Feb 24 17:32:51 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 06:00:17 2012 Êîäèðîâêà: koi8-r Ïîèñêîâûå ñëîâà: ôîìàëüãàóò

A&A 438, 39­53 (2005) DOI: 10.1051/0004-6361:20042567
c ESO 2005

Astronomy & Astrophysics

The WEBT campaign to observe AO 0235+16 in the 2003­2004 observing season
Results from radio-to-optical monitoring and XMM-Newton observations
C. M. Raiteri1 , M. Villata1 ,M. A.Ibrahimov2 , V.M.Larionov3,4 ,M.Kadler5 ,H .D.Aller6 , M. F. Aller6 , Y. Y. Kovalev7,8 , L. Lanteri1 ,K.Nilsson9 , I. E. Papadakis10,11 , T. Pursimo12 ,G.E.Romero13 , H. TerÄsranta14 , M. Tornikoski14 ,A. A.Arkharov15 , D. Barnaby16 , A. Berdyugin9 , M. BÆttcher17 , K. Byckling18 ,M. T.Carini16 , D. Carosati19 , S.A.Cellone20 , S. Ciprini9 ,J. A.Combi13,21 , S. Crapanzano1 , R. Crowe22 ,A.DiPaola23 ,M. Dolci24 , L. Fuhrmann1,5,25 ,M.Gu26 ,V. A.Hagen-Thorn3,4 ,P.Hakala18 , V. Impellizzeri5 ,S.Jorstad27 , J. Kerp5 , G. N. Kimeridze28 , Yu. A. Kovalev8 ,A.Kraus5 ,T.P.Krichbaum5 ,O.M.Kurtanidze28 , A. LÄhteenmÄki14 , E. Lindfors9 , M. G. Mingaliev29 , R. Nesci30 ,M.G.Nikolashvili28 ,J. Ohlert31 ,M.Orio1 ,L.Ostorero32 , M. Pasanen9 , A. Pati33 ,C.Poteet16 , E. Ros5 , J. A. Ros34 ,P.Shastri33 , L.A.Sigua28 , A. SillanpÄÄ9 ,N.Smith35 ,L.O. Takalo9 , G. Tosti25 ,A.Vasileva3 ,S.J. Wagner32 , R. Walters16 ,J. R.Webb36 , W. Wills16 ,A.Witzel5 , and E. Xilouris37
1 INAF - Osservatorio Astronomico di Torino, Italy, e-mail: raiteri@to.astro.it 2 Ulugh Beg Astronomical Inst., Uzbekistan 3 Astronomical Inst., St.-Petersburg State Univ., Russia 4 Isaac Newton Inst. of Chile, St. Petersburg Branch 5 Max-Planck-Institut fÝr Radioastronomie, Germany 6 Dept. of Astronomy, Univ. of Michigan, USA 7 National Radio Astronomy Observatory, Green Bank, USA 8 Astro Space Center of Lebedev Physical Inst., Russia 9 Tuorla Astronomical Observatory, Finland 10 IESL, FORTH, Greece 11 Physics Dept., Univ. of Crete, Greece 12 Nordic Optical Telescope, Santa Cruz de La Palma, Spain 13 Inst. Argentino de RadioastronomÌa, Argentina 14 MetsÄhovi Radio Observatory, Helsinki Univ. of Technology, Finland 15 Pulkovo Observatory, Russia 16 Dept. of Physics & Astronomy, Western Kentucky Univ., USA 17 Dept. of Physics and Astronomy, Ohio Univ., USA 18 Observatory, Univ. of Helsinki, Finland 19 Armenzano Astronomical Observatory, Italy 20 Facultad de Ciencias AstronÑmicas y GeofÌsicas, Argentina 21 Dept. de FÌsica, Escuela PolitÈcnica Superior, Univ. de JaÈn, Spain 22 Dept. of Physics & Astronomy, Univ. of Hawaii, USA 23 INAF - Osservatorio Astronomico di Roma, Italy 24 INAF - Osservatorio Astronomico di Collurania Teramo, Italy 25 Osservatorio Astronomico, Univ. di Perugia, Italy 26 Korea Astronomy Observatory, Korea 27 Inst. for Astrophysical Research, Boston Univ., USA 28 Abastumani Observatory, Georgia 29 Special Astrophysical Observatory, Karachaevo-Cherkessia, Russia 30 Dipartimento di Fisica, Univ. di Roma "La Sapienza", Italy 31 Michael Adrian Observatory, Trebur, Germany 32 Landessternwarte Heidelberg, Germany 33 Indian Inst. of Astrophysics, Bangalore, India 34 AgrupaciÑ AstronÐmica de Sabadell, Spain 35 Dept. of Applied Physics & Instrumentation, Cork Inst. of Technology, Ireland 36 SARA Observatory, Florida International Univ., USA 37 Inst. of Astronomy and Astrophysics, National Observatory of Athens, Greece

Received 17 December 2004 / Accepted 17 February 2005
Abstract. A multiwavelength campaign to observe the BL Lac object AO 0235+16 (z = 0.94) was set up by the Whole Earth Blazar Telescope (WEBT) collaboration during the observing seasons 2003­2004 and 2004­2005, involving radio, near-IR and optical photometric monitoring, VLBA monitoring, optical spectral monitoring, and three pointings by the XMM-Newton satellite. Here we report on the results of the first season, which involved the participation of 24 optical and near-IR telescopes and 4 radio telescopes, as well as the first XMM-Newton pointing, which occurred on January 18­19, 2004. Unpublished data from previous epochs were also collected (from 5 optical-NIR and 3 radio telescopes), in order to fill the gap between the end of the period presented in Raiteri et al. (2001) and the start of the WEBT campaign. The contribution of the southern AGN, 2 arcsec distant from the source, is taken into account. It is found to especially affect the blue part of the optical spectrum when the source is faint. In the optical and near-IR the source has been very active in the last 3 years, although it has been rather faint most of the time, with noticeable variations of more than a magnitude over a few days. In contrast, in the radio bands it appears to have been "quiescent" since early 2000. The major radio (and optical) outburst predicted to peak around February­March 2004 (with a six month uncertainty) has not occurred yet. When comparing our results with the historical light

For questions regarding the availability of the data presented in this paper, please contact the WEBT President Massimo Villata (villata@to.astro.it). Jansky Fellow.

76524002:1636-4000/1501.01/gro.iod.xd//:ptth ro aa/gro.secneicspde.www//:ptth ta elbaliava dna secneicS PDE yb dehsilbup elcitrA


40

C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

curves, two different behaviours seem to characterize the optical outbursts: only the major events present a radio counterpart. The X-ray spectra obtained by the three EPIC detectors are well fitted by a power law with extra-absorption at z = 0.524; the energy index in the 0.2­10 keV range is well constrained: = 0.645 ± 0.028 and the 1 keV flux density is 0.311 ± 0.008 µJy. The analysis of the X-ray light curves reveals that no significant variations occurred during the pointing. In contrast, simultaneous dense radio monitoring with the 100 m telescope at Effelsberg shows a 2­3% flux decrease in 6­7 h, which, if intrinsic, would imply a brightness temperature well above the Compton limit and hence a lower limit to the Doppler factor > 46. We construct the broad-band spectral energy distribution of January 18­19, 2004 with simultaneous radio data from Effelsberg, optical data from the Nordic Optical Telescope (NOT), optical­UV data from the Optical Monitor onboard XMMNewton, and X-ray data by the EPIC instruments. Particular care is taken to correct data for extinction due to both the Milky Way and the z = 0.524 absorber. The resulting SED suggests the existence of a bump in the UV spectral region.
Key words. galaxies: active ­ BL Lacertae objects: general ­ BL Lacertae objects: individual: AO 0235+16

1. Introduction
BL Lac objects belong to the blazar class of active galactic nuclei (AGNs). They show variability at all wavelengths, from the radio to the -ray band, on a variety of time scales. In general, slow oscillations of the flux mean level on time scales of weeks to years have been observed, based on which fast flares lasting as little as a few hours are superposed. The commonly accepted scenario for blazar emission assumes a black hole surrounded by an accretion disc, feeding a plasma jet which is oriented close to the line of sight. This geometrical configuration can explain the apparent superluminal motion of relativistically moving plasma down the jet which is often derived from the analysis of radio maps. Moreover, relativistic Doppler effects can, at least in part, explain the extremely high brightness temperature which otherwise would result in the so-called Compton catastrophe. The best way to study the source flux variability and to understand the corresponding physical processes is the simultaneous observation across the whole accessible electromagnetic spectrum. Indeed, it is commonly accepted that the observed lower-energy flux (from the radio band to the UV­X-ray one) is due to synchrotron emission from the plasma jet, while the higher-energy one (up to GeV­TeV energies) is thought to be produced by inverse-Compton scattering of soft photons off the same relativistic electrons producing the synchrotron radiation. It goes without saying that the high-energy flux should thus be correlated with the low-energy one, but the relative variability amplitude and the amount of time lag expected depend on the particular model adopted. Analysis of the long-term radio and optical light curves of AO 0235+16 led to the conclusion that the main radio outbursts have optical counterparts, and that they repeat every 5.7 ± 0.5 years (Raiteri et al. 2001). This was interpreted in terms of geometrical effects due to the possible presence of a binary black hole system (Romero et al. 2003; Ostorero et al. 2004). The next outburst was predicted to peak around February­March 2004, and a multiwavelength radio-NIR-optical campaign has been organized by the Whole Earth Blazar Telescope (WEBT; http://www.to.astro.it/blazars/webt/; e.g. Villata et al. 2000, 2002, 2004a,b) collaboration to follow the source behaviour. The WEBT campaign began in summer 2003, and was planned to continue until spring 2005. During this longterm campaign three pointings by the XMM-Newton satellite were granted in order to compare the behaviour of the

low-energy emission with the high-energy one. The XMM-Newton observations were scheduled on January 18­19 and August 2, 2004, and January 28, 2005. We know from previous observations by X-ray satellites (Einstein et al. 1990; EXOSAT, Ghosh & Soundararajaperumal 1995; ROSAT, Madejski et al. 1996; Comastri et al. 1997; ASCA, Madejski et al. 1996; Junkkarinen et al. 2004; RXTE, Webb et al. 2000; BeppoSAX, Padovani et al. 2004) that both the X-ray flux and spectral index can be very variable. The variability of AO 0235+16 on short time scales is also of special interest. At radio wavelengths several flatspectrum radio quasars and BL Lacs are known to exhibit rapid variability (Witzel et al. 1986; Heeschen et al. 1987; Kedziora-Chudczer et al. 2001; Lovell et al. 2003; Kraus et al. 2003) on time scales of days to hours (intraday variability, IDV). The cause of the variations seen in these sources is currently controversial, with claims being made for either a source-intrinsic (e.g. shocks in the emitting jet) or extrinsic origin (interstellar scintillation or gravitational microlensing). Interstellar scintillation (ISS) may play an important role in the cm-radio regime (Jauncey & Macquart 2001; Rickett et al. 2001; Dennet-Thorpe & de Bruyn 2002), otherwise intrinsic variations may require extremely high Doppler boosting ( 100) or special source geometries in order to prevent the inverse-Compton catastrophe (see e.g. Wagner & Witzel 1995; Wagner et al. 1996). Since ISS cannot explain correlated IDV over a wide range of the electromagnetic spectrum, the detection of broad-band correlations would directly rule out ISS as sole explanation for radio IDV and would favour a source intrinsic contribution to the IDV pattern. AO 0235+16 has shown IDV over a wide range of the electromagnetic spectrum, from the radio to the optical band (Takalo et al. 1992; Heidt & Wagner 1996; Noble & Miller 1996; Romero et al. 1997; Kraus et al. 1999; Romero et al. 2000; Raiteri et al. 2001). Moreover, hints for the existence of very high Doppler factors come from a variety of other observational studies (Fujisawa et al. 1999; Frey et al. 2000; Jorstad et al. 2001), thus making this blazar an ideal candidate for the detection of correlated broad-band IDV. XMM-Newton has already revealed fast X-ray variability of AO 0235+16 (Kadler et al. 2005). Consequently, one of the aims of the WEBT campaign was to search for radio-NIR-optical-X-ray correlated fast variations during the XMM-Newton pointings, with the support of dense radio monitoring by the 100 m radio telescope at Effelsberg.


C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

41

Side by side with the WEBT campaign, spectroscopic monitoring of the source in the optical band with both the Telescopio Nazionale Galileo (TNG) and the ESO's New Technology Telescope (NTT) has been carried out (Raiteri et al. 2005). 15 VLBA epochs have been awarded to monitor the source structure variability during the WEBT campaign (Wiik et al. 2005). In this paper we present results obtained during the first observing season of the campaign, from July 1, 2003 to April 2004, including the first XMM-Newton pointing of January 18­19, 2004. We also collected data taken before the start of the campaign, back to the last data published by Raiteri et al. (2001), in order to reconstruct the source behaviour in between. Other papers will follow, containing new results and more detailed analyses. The long-term optical, NIR and radio monitoring data obtained by the WEBT during the 2003­2004 observing season, along with the data filling the gap between this WEBT campaign and the data published by Raiteri et al. (2001), are presented in Sect. 2. The multiwavelength observations performed during the XMM-Newton pointing of January 18­19, 2004 are analysed in Sect. 3. In Sect. 4 we discuss the problem of extinction and construct the broad-band spectral energy distribution of AO 0235+16 during the XMM-Newton pointing. Section 5 contains some results of the cross-correlation analysis. Conclusions are drawn in Sect. 6.

The source was calibrated preferentially with respect to Stars 123(UBVRI photometry by Smith et al. 1985). However, since the source was rather faint, in some cases long exposure times produced saturation of Star 2 (and Star 3), so that it was necessary to include other reference stars for the calibration. Hence, we sometimes also used Stars 6 and C1 calibrated in VRI by Fiorucci et al. (1998). We made several checks to test the variation of the source magnitude value when adopting different choices of reference stars for the calibration: in general, we saw that the variation is within a few hundredths of a mag, which is usually inside the errors. In order to construct meaningful light curves, the observations reported in Table 1 were carefully analysed and unreliable data discarded. This was done by zooming into subsequent very short periods of time and inspecting all the light curves at the same time. The same trend in all bands was required within errors; bad points were eliminated. In the case of noisy intranight datasets from the same telescope, some binning over time intervals not exceeding 20 min was performed. In general, data affected by errors larger than 0.2 mag were not included. Moreover, we investigated the possibility that our photometric measurements of AO 0235+16 were affected by the contribution of objects very close to the source. The complex environment of the source was studied by various authors (Yanny et al. 1989; Burbidge et al. 1996; Nilsson et al. 1996); they noticed that there are a few galaxies at redshift of 0.524 in the source surroundings, and two of them are very close to AO 0235+16. One galaxy, probably a normal spiral, is 1.3 arcsec to the east, but it is very faint and thus does not affect the source photometry, even when the source is faint. The object about 2 arcsec to the south is more problematic, since it is brighter than the former one, and it is also known to be an AGN, possibly interacting with the eastern galaxy. The photometric contribution of this southern object (hereafter named ELISA: elusive intervening southern AGN) can affect the source magnitude when it is very faint, especially in the bluer part of the spectrum. This increasing contribution towards the ultraviolet is mainly due to the fact that the AO 0235+16 spectrum is much steeper than that of ELISA, since AO 0235+16 is strongly absorbed in these bands (see Sect. 4). We thus estimated the magnitude of ELISA from the Mt. Maidanak and Hanle frames, where it is often clearly resolved. There is some dispersion in the data, and it is not possible to recognize any clear variability trend in the period 2000­2004. The best estimates are: U = 20.8, B = 21.4, V = 20.9, R = 20.5, I = 19.9. These values are in fair agreement with R = 20.4 ± 0.1, R - I = 0.7 ± 0.2 given by Nilsson et al. (1996), and with m450 = 21.40 ± 0.01 estimated by Burbidge et al. (1996). We then converted the AO 0235+16 magnitudes into fluxes, subtracted the ELISA contribution, and then converted again into magnitudes. In Fig. 1 we present the UBVRI light curves obtained during the WEBT campaign, after applying both the cleaning and ELISA-subtraction procedures described above. These are the light curves that will be used for the subsequent analysis. Data from different telescopes are plotted with different colours and symbols in order to be able to distinguish the single

2. Long-term observations by the WEBT
The WEBT campaign was started on July 1, 2003; the last data presented here were taken in April 2004. However, for the present work we collected data from previous epochs in order to fill the gap between the start of the WEBT campaign and the last data published in Raiteri et al. (2001), and to enrich the historical light curves.

2.1. Optical and near-IR data
Tables 1 and 2 show the lists of the optical and near-IR telescopes participating in the WEBT campaign (and prior data collection) ordered by longitude1. In the tables Col. 1 contains the name of the observatory and its location; Col. 2 gives the telescope diameter in cm; Cols. 3­7 report the number of observations done in the UBVRI and JHK bands; Col. 8 shows the total number of observations performed by each telescope in those bands. Data in brackets indicate the number of points remaining after the "cleaning" procedure described below and which are shown in the light curves. As one can see in the last row, the total number of observations collected is 2509 in the optical and 385 in the near-IR. Data were collected as instrumental magnitudes of the source and comparison stars in the same field.
1 One of the main advantages of WEBT is to have observers spread out at different longitudes. Hence, the observing task can move from east to west allowing us in principle to obtain continuous monitoring over 24 h.


42

C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

Table 1. List of participating optical observatories by longitude and number of observations performed in the various bands. Observatory IAO Hanle, India Mt. Maidanak, Uzbekistan Abastumani, Georgia (FSU) Crimean, Ukraine Skinakas, Crete Tuorla, Finland Vallinfreda, Italy Armenzano, Italy Perugia, Italy Michael Adrian, Germany Greve, Italy Torino, Italy Heidelberg, Germany Calar Alto, Spain Sabadell, Spain Roque (KVA), La Palma Roque (NOT), La Palma Roque (TNG), La Palma CASLEO, Argentina Bell, Kentucky Kitt Peak (SARA), Arizona Mt. Lemmon, Arizona Lowell (Perkins), Arizona Mauna Kea, Hawaii Total d (cm) 200 150 70 70 130 103 50 40 40 120 32 105 70 220 50 35 256 358 215 60 90 100 180 60 N
U

N

B

N

V

N

R

N

I

N

tot

8 (6) 59 (42) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 (7) 0 0 0 0 0 0 0 79 (55)

7 (6) 336 (311) 0 0 24 (17) 0 0 14 (12) 0 0 0 6 (2) 0 0 0 0 14 (13) 0 0 0 0 2 (0) 1 (1) 3 (0) 407 (362)

12 (10) 112 (107) 0 17 (16) 27 (23) 0 0 28 (18) 2(0) 0 0 12 (11) 0 0 0 0 14 (13) 0 55 (55) 0 3 (2) 5 (2) 3 (2) 4 (1) 294 (260)

8 (8) 450 (417) 41 (32) 19 (18) 28 (27) 44 (17) 1 (1) 41 (32) 4(1) 7 (5) 1 (1) 132 (102) 4 (2) 101 (15) 90 (1) 180 (95) 23 (23) 2 (2) 59 (59) 10 (10) 30 (19) 8 (5) 4 (4) 0 (0) 1287 (896)

14 (13) 237 (219) 0 19 (19) 26 (25) 0 0 29 (26) 4(4) 0 0 83 (62) 1 (0) 0 0 0 17 (15) 0 0 0 3 (3) 5 (5) 3 (3) 1 (1) 442 (395)

49 (43) 1194 (1096) 41 (32) 55 (53) 105 (92) 44 (17) 1(1) 112 (88) 10 (5) 7(5) 1(1) 233 (177) 5(2) 101 (15) 90 (1) 180 (95) 80 (71) 2(2) 114 (114) 10 (10) 36 (24) 20 (12) 11 (10) 8(2) 2509 (1968)

Table 2. List of participating near-IR observatories by longitude and number of observations performed in the various bands. Observatory Campo Imperatore, Italy TIRGO, Switzerland Roque (NOT), La Palma Total d (cm) 110 150 256 N
J

N

H

N

K

N

tot

82 (82) 59 (1) 6 (6) 147 (89)

81 (81) 1 (1) 5 (5) 87 (87)

88 (87) 55 (1) 8 (8) 151 (96)

251 (250) 115 (3) 19 (19) 385 (272)

contributions as well as to appreciate how well the datasets agree with one another. The source was in a rather faint state, but was undergoing noticeable variability, with a total variation of about 2 mag in the best sampled R band, and changes of about 1 mag on time scales of about a week. Average colour indices (and standard deviations) obtained by coupling data acquired by the same telescope within 20 min are: B-R = 1.72 (0.12), V -R = 0.76 (0.10), R- I = 0.99 (0.09). These averages are obtained from 223, 151 and 296 indices, respectively. Near-IR observations during the WEBT campaign have been carried out with the 110 cm telescope at Campo Imperatore and the 256 cm Nordic Optical Telescope (NOT). The source magnitude was calibrated with the

GonzÀlez-PÈrez et al. (2001) photometry. The contribution from ELISA in this frequency regime is negligible. Figure 2 shows the JHK light curves compared with the R-band one. The near-IR light curves are generally undersampled to allow a detailed statement about a correspondence with the optical data, but where sufficient sampling has been achieved, they indicate that a close correspondence exists, as expected. Average colour indices (and standard deviations) obtained by coupling data acquired by the same telescope within 20 min are: J - K = 2.12 (0.09), H - K = 1.00 (0.07). A rough estimate of optical-NIR colour indices can be obtained by coupling data within 60 min: R - J = 2.50 (0.23), R - H = 3.60 (0.29), R - K = 4.58 (0.34). JHK data before the campaign were obtained on December 18, 2002 (JD = 2 452 627) with the 150 cm Italian


C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

43

Fig. 1. UBVRI light curves of AO 0235+16 during the observing season 2003­2004; the XMM-Newton pointing of January 18­19 is marked by the vertical (cyan) line at JD = 2 453 023. The ELISA contribution has been subtracted as explained in the text.

national telescope TIRGO, at Gornergrat, Switzerland. These data yield average values of J = 15.28 ± 0.06, H = 14.36 ± 0.10, K = 13.17 ± 0.08.

(Raiteri et al. 2001), nonetheless it reveals that the source can be very active even during faint states.

2.3. Radio data 2.2. Intraday variability
On very short time scales, the source displayed clear variability in various occasions. As an example, well-sampled intranight light curves in V and R bands were obtained from November 10 to 14, 2001, with the 215 cm telescope of the CASLEO, Argentina. The R-band light curve is shown in Fig. 3, with the same calibration and ELISA flux subtraction adopted for the data in Fig. 1. Clear intranight trends are recognizable on November 10, 12, and 13, with smooth variations of about 0.1 mag in a few hours. On internight time scales, variations as fast as 0.5 mag in 1­2 days were registered. Although this behaviour is not as extreme as observed in other occasions, such as the 1.2 mag brightening in one day observed by Romero et al. (2000) in November 1999, or similar fast variations detected in the period 1997­1999 Long-term radio observations during the WEBT campaign were performed with the 14 m antenna of the MetsÄhovi Radio Observatory in Finland, the 26 m paraboloid of the University of Michigan Radio Astronomy Observatory (UMRAO) in the United States, and the 600 m ring transit radio telescope RATAN-600 of the Russian Academy of Sciences in Russia. Details on the observing and data reduction procedures can be found e.g. in TerÄsranta et al. (1998), Aller et al. (1999), and Kovalev et al. (1999a). A comparison among optical magnitudes in the bestsampled R band and radio fluxes at different frequencies is shown in Fig. 4. Data at 37 and 22 GHz are from the MetsÄhovi Radio Observatory; data at 14.5, 8.0, and 4.8 GHz are from UMRAO; those at 21.7, 11.1, 7.7, 3.9, and 2.3 GHz are from RATAN. The two datasets at 37 GHz from MetsÄhovi have been obtained by two different teams following slightly


44

C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

seems to characterize the radio light curves starting from JD 2 452 970, the slope of which decreases with increasing wavelength, from 4.0 â 10-3 at 37 GHz to -0.3 â 10-3 at 5 GHz. A correlation with the optical data, if it exists, is not obvious. The analysis of the long-term radio and optical light curves led Raiteri et al. (2001) to the conclusion that the major radio outbursts repeat quasi-periodically every 5.7 years and, when optical data are available, an optical counterpart is found immediately preceding the radio event. The next major outburst was predicted to peak in February­March 2004, with a possible window of ±6 months around that date. As Fig. 4 shows, no major radio or optical outburst was detected during the WEBT campaign. To make the comparison between the recent and past behaviour of the source possible, Fig. 5 displays the optical and radio light curves over a more extended time period. The period shown overlaps with that presented in Fig. 8 of Raiteri et al. (2001), and includes the last two major radio outbursts which occurred in 1992­1993 and 1998, both of them presenting an optical counterpart. The first (red) vertical line marks the end of the Raiteri et al. (2001) data; the second (green) line indicates the start of the WEBT campaign on July 1, 2003, and the third (cyan) line corresponds to the XMM-Newton pointing of January 18­19, 2004. The R-band data from Raiteri et al. (2001) are here complemented by data from the Hamburg Quasar Monitoring program (partially published in Schramm et al. 1994) and by data from CASLEO (Romero et al. 2002); similarly, past data from RATAN (partially published by Kovalev et al. 1999b; Kiikov et al. 2002) and data by Venturi et al. (2001) are added to the Raiteri et al. (2001) radio light curves. There are only sparse optical data between the last data presented in the Raiteri et al. (2001) paper and the beginning of the WEBT campaign, which show a noticeable mid-term variability, with a maximum variation of 3.3 mag in 260 days. Variations of smaller amplitude, but still noticeable, are found on shorter time scales (see also Sect. 2.2). These sparse optical data do not allow one to rule out completely the possibility that a major outburst took place in the 2000­2003.5 period, while the impressive sampling obtained during the WEBT campaign assures that no major optical outburst occured in the 2003­2004 observing season. In contrast, the continuous radio monitoring clearly indicates that there have been no outbursts since 1998. The new data collected starting from late 2000 show that the radio flux has been low at all frequencies, even if not at the minimum levels registered in 1996­1997.

Fig. 2. JHK light curves of AO 0235+16 during the 2003­2004 observing season compared to the ELISA-subtracted R-band one (top). NIR data from Campo Imperatore are shown by (blue) circles, while (red) triangles indicate NIR data from the NOT.

Fig. 3. R-band light curve of AO 0235+16 during November 10­14, 2001; data acquired with the 215 cm telescope of the CASLEO, Argentina.

different procedures. Further data at 4.85 and 10.45 GHz were taken with the 100 m radio telescope at Effelsberg during the XMM-Newton pointing of January 18­19, 2004 (see Sect. 3.2, for more details). ELISA is not a radio emitting object, so no correction for its contribution is required. The figure clearly shows a decreasing variability amplitude accompanied by an increasing base level when going from the higher to the lower radio frequencies. In general, after an initial small bump and a subsequent lower state, a quasi-linear trend

3. The XMM-Newton pointing of January 18­19, 2004
The X-ray Multi-Mirror Mission (XMM) ­ Newton satellite observed AO 0235+16 during revolution number 0753, from January 18, 2004 at 19:06:29 UT to January 19 at 03:28:24 UT (JD = 2 453 023.29617­2 453 023.64472).


C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

45

Fig. 4. ELISA-subtracted R-band magnitudes of AO 0235+16 (top) and radio fluxes (Jy) at different frequencies during the WEBT campaign. The vertical (cyan) line marks the XMM-Newton pointing of January 18­19, 2004. The dotted horizontal line at F = 1.5 Jy is drawn to guide the eye through the flux variations.

3.1. Results from EPIC
The European Photon Imaging Camera (EPIC) includes two types of CCD cameras: pn and MOS. There are two MOS cameras, MOS1 and MOS2, so that EPIC involves 3 detectors in total. Since an outburst was expected, a medium filter-small window configuration was chosen for the EPIC detectors in order to avoid both possible contamination by lower-energy photons and photon pile-up. Moreover, a small window is also able

to reduce the fraction of out-of-time events recorded by the EPIC-pn detector. Data were reduced with the Science Analysis System (SAS) software version 6.0. When extracting the spectrum of both the source and the background, we included the strings "FLAG==0"and "PATTERN<=4" in the selection expression for all the three EPIC detectors. The first string excludes events next to the edges of the CCDs and next to bad pixels which


46

C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

Fig. 5. ELISA-subtracted R-band magnitudes of AO 0235+16 (top) compared with radio fluxes (Jy) at different frequencies; in each panel the first (red) vertical line marks the last data published by Raiteri et al. (2001); the second (green) line indicates the start of the WEBT campaign, and the third (cyan) one corresponds to the XMM-Newton pointing of January 18­19, 2004.

may have incorrect energies; the second string selects only single and double events, which have the best energy calibration. Source spectra were extracted from a circular region with a30 arcsec radius2 . Notwithstanding the small-window configuration, we were able to extract the background spectrum on
2 We examined an archive image of the AO 0235+16 field taken in August 2000 by Chandra, which has a space resolution of 0.492 arcsec, and verified that no X-ray emission was detected from

the same CCD of the source within circles of 20 arcsec radius in the case of MOS1 and MOS2, and 30 arcsec for pn. These background spectra, when compared with others extracted in different regions, produce the most satisfactory background subtraction.
ELISA. This could be due to intrinsic X-ray faintness, or to Compton absorption, which would favour a type 2 AGN classification of ELISA.


C. M. Raiteri et al.: The WEBT campaign on AO 0235+16 in the 2003­2004 observing season

47

The source spectra were grouped with the task "grppha" of the FTOOL package in order to have a minimum of 25 counts in each bin. The detected net count rates are 0.1876 ± 0.0032 cts/s for MOS1, 0.1907 ± 0.0033 cts/s for MOS2, 0.6261 ± 0.0063 cts/s for pn. In the following we present the results of simultaneous spectral fitting of the data from the three EPIC detectors performed with the Xspec package version 11.3.0. Only data in the channels corresponding to energy greater than 0.2 keV were considered. Once bad channels had been ignored, the channels retained for the subsequent analysis were: 3­148 for MOS1, 6­153 for MOS2, 7­420 for pn, for a total of 708 PHA bins. We report the results in terms of H I column density NH and energy index , which is defined as F E = F0 E - . The first approach (Model 1) was