S5 1803+784 revisited

1. S5 1803+784 revisited R.Nesci1, A. Maselli2, F. Montagni3, S.Sclavi1 1) Universita’ La Sapienza, Roma; 2) INAF IASF-Palermo; 3) Osservatorio Greve in Chianti

2. What is S5 1803+784? S5 1803+784 is a BL Lac object, a special class of Active Galactic Nuclei (AGN). BL Lacs take their names form the prototype source, discovered at Sonnenberg in 1929 and then catalogued as the variable star BL Lacertae. They are charcterized by: • featureless optical continuum • large and fast variability • Appreciable polarization • Flat radio spectrum Most BL Lacs are X-ray sources, and are a substantial part of the extragalactic sources detected in Gamma-rays. S5 1803+784

3. Current unified model The current model to explain the emission of BL Lacs is based on a supermassive black-hole, hosted in the galaxy centre which accrete matter from the surroundings. Part of the matter in the accretion disk is funnelled in narrow relativistic jet, oriented at small angle with the observer line of sight: the existence of such jets has been proved by VLBI radio observations Most of the observed emission is believed to be produced in this jet, by two different processes: • Sinchrotron radiation of relativistic electrons in a strong magnetic field, typically peaking at NIR frequencies; • Inverse Compton radiation from these electrons on the ambient photons, typically peaking in the hard X-rays. S5 1803+784

4. Spectral energy distribution (SED) classification • Depending on the energy of the electrons and intensity of the magnetic field, the peak of the synchrotron emission, in the Log(nu)-Log(nu*Fnu) plane, may range between far IR and X-rays. • Accordingly, BL Lacs are classified as Low-frequency (LBL), intermediate frequency (IBL) and high frequency (HBL) peaked sources. • During a flare, the SED shape is deformed and generally shifted towards higher frequencies. S5 1803+784

5. Unified scheme • A schematic view seen edge-on is in this picture S5 1803+784

6. Time variability • The light curves of several BL Lacs have been studied since their discovery, also using historical plate archives for long term variability studies. • Their behaviour has proven to be generally quite erratic, save the case of OJ 287 (Sillampaa et al 1998) which has a periodicity of about 11 years for the major flares. • In a number of cases time scales of a few years have been found, even though not strict periodicities (e.g. S5 0716+714, Nesci et al. 2005; AO 0235+174, Raiteri et al. 2008; GB6 J1058+5628, Nesci 2010). • Long term trends have also been found in a few sources, superposed to the short term ones, extending for tens of years and possibly being just part of longer recurrency time-scale (e.g. OQ 530, Massaro et al. 2004; ON 231, Massaro et al. 2001; S5 0716+714, Nesci et al. 2005). S5 1803+784

7. Possible models • The long term variability may be due to intrinsic changes in the source or to geometrical effects: in this case the best explanation could be a precession of the relativistic jet, producing a variation of the Doppler boosting and therefore an achromatic variation of the apparent source luminosity. S5 1803+784

8. The source in short • S5 1803+78 is a bright radio source discovered in 1981 (Kuhr et al.) . • Due to its circumpolar position it can be followed for a large part of the year. • It has a large optical polarization and a redshift z=0.680, based on a single emission line assumed to be MgII (2900 A). • It is well monitored in the radio band because it is a source of the ICRS reference frame and is used also as a geodetic reference source. • It was observed by several X-ray satellites, but was not pointed by EGRET. • Its first optical systematic monitoring was made by Nesci et al. (2002) in the years 1996-2001. • The overall SED is of the LBL type S5 1803+784

9. Optical follow-up • Most of the well monitored LBL sources (e.g. BL Lacertae itself), show a correlation between optical luminosity and optical spectral slope, being “bluer when brighter”. • On the contrary, S5 1803+78 showed a very limited variation in the optical spectral slope, despite it varied by more than 2 mag in flux. • We decided therefore to continue its monitoring, to check if: • A) the color variation remained small in time; • B) the source showed any recurrent time scale in flux variations; • C) the source showed monotonic long term trend. The monitoring was performed with the telescopes of Asiago (183 cm), Loiano (152cm), Vallinfreda (50cm) and Greve in Chianti (30cm), mainly in the R band, but with several observations also in the B, V, and I bands. S5 1803+784

10. The light curve from 1996 to 2010 S5 1803+784

11. expanded S5 1803+784

12. Main characteristics • The source underwent a number of flares: • We derived by eye estimates of the starting and ending point for the raising and falling branches of each major well sampled flare; • The luminosity variation rates range between 0.03 and 0.07 mag/day. • There is no strong evidence that raising and falling rates are statistically different. • Inspection of the light curve (by eye and by FFT analysis) shows no strict periodicity for the flares. • A possible time scale of 1300 days between major flares may be present. S5 1803+784

13. Table of major flares parameters S5 1803+784

14. Optical color index • The (V-I) color index showed small variations around the average value of <V-I>=1.1 during all our monitoring, a typical value for LBL sources, corresponding to a spectral index -1.6. • The most accurate measures, made with the larger Asiago and Loiano telescopes during the strong flare at JD 3550, showed a bluer color during the raising phase and a redder one in the falling part. A similar behaviour was detected in a few other BL Lacs (S5 0716+714, 3C 66A, Ciprini et al. 2004). • The UV spectral slope observed by Swift-UVOT was consistent with the optical one from our BVRI data, indicating a common origin (tail of the Synchrotron emission) for the optical and UV light. S5 1803+784

15. Color variations during major flares: case 1 S5 1803+784

16. Case 2 S5 1803+784

17. Case 3 S5 1803+784

18. The X-ray spectrum • The X-ray spectrum was observed by SWIFT in four occasions between 2007 and 2009, • It always showed a flux level comparable to that observed by BeppoSAX in 1998, with small variations. • Also the X-ray spectral slope (photon index) was stable, around +1.6, suggesting an Inverse Compton origin. • During these X-ray observations the R magnitude was not very different, ranging between 15.4 to 16.3: no clear correlation is evident between the X-ray and optical bands. S5 1803+784

19. X-ray spectral slope S5 1803+784

20. FermiGST observations • S5 1803+78 was not pointed by EGRET. • It was detected by Fermi-GST, which performs an all-sky monitoring, and is present in the 11-month catalogue (1FGL, Abdo et al. 2010a). • The Gamma-rays light curve, sampled weekly, shows an average level around 2x10-8 phot/cm2/s with oscillations within a factor 2. • Our optical light curve has only a few points around the end of this 11 months time interval, and shows a short increase of 0.3 mag two weeks before a short Gamma-ray bump. S5 1803+784

21. The Gamma ray flare • A strong flare, at 9x10-7 phot/cm2/s (E>100 MeV) was reported by Donato et al. (2010), 40 times brighter than the average value. • We could observe the source in the optical 4 days after the Gamma burst at R=14.9 with color index (V-I)=1.2. • It was already in the decreasing phase at a rate of 0.04 mag/day, a typical value for this source. • Extrapolating backwards this trend to the epoch of the Gamma-ray flare gives R=14.7 for the peak value, rather lower than the maximum recorded flares of this source (R=14.0-14.2). S5 1803+784

22. Spectral Energy Distribution (SED) We used literature data, as well as our own observations and data analysis, to build the SED of the source. The hard X-ray flux in the 15-30 keV band has been derived from the 54-months BAT all-sky map with the BATIMAGER software (Segreto et al. 2010). Fermi-GST data have been retrieved from the ASI Science Data Center (ASDC) website. For clarity, not all available data are reported in the figure. S5 1803+784

23. The SED of S5 1803+784 S5 1803+784

24. SED description • The SED is typical of an LBL source, with two large bell-shaped parts: • the synchrotron component peaking between 1013 and 1014 Hz, is well fitted by a log-parabola (Maselli et al. 2008); • the Inverse Compton branch is well defined by the Swift-XRT, BAT and Fermi-GST instruments: a log-parabolic fit gives a peak at 3.8x1021 Hz, similar to that of other LBL objects observed by Fermi-GST. • The Gamma-ray photon spectral index is 2.33, within the range of values for LBL objects (2.21+-0.16, see Abdo et al. 2010b). S5 1803+784

25. Conclusions • The optical light curve shows several large (2 mag) flares with a possible time scale of about 1300 days. • Flux variations are slow, with typical rate of 0.04 mag/day, without marked difference between raising and falling branches. • No monotonic long term trend is apparent on a 14-years time span. • The optical spectral slope shows little (10%) variations even in large flares. Accuracy of 0.01 mag are necessary to study these variations. • In one strong flare we had enough accuracy to detect a “bluer when rising” behaviour. • The overall SED is well described by two log-parabola, with peak positions and slopes typical of LBL sources. • Up to now, correlation between Gamma-ray and optical flares is based on one case only: if true we statistically expect another Gamma-ray flare within the lifetime of FermiGST. S5 1803+784

26. References • Abdo et al. 2010a, arXiv:1002.0150v2; ApJ in press. • Abdo et al. 2010b, ApJ 710, 1271. • Donato D. et al. 2010, Atel #2386. • Fiorucci et al. 2004, A&A 419, 25. • Maselli A et al. 2008, Blazar Variability Across EM Spectrum, at http://pos.sissa.it, p.77 . • Massaro et al. 2001, A&A 374, 435 • Massaro et al. 2004, A&A 423, 935 • Nesci R., et al. 2002, AJ 124, 53. • Nesci R. et al. 2005, AJ 130, 1466 • Nesci R, 2010, AJ 139, 2425 • Raiteri et al. 2008, A&A 485, L17 • Segreto et al. 2010, A&A 510, 47 • Sillampaa A et al. 1988 ApJ 325, 628 S5 1803+784