Search

Blazars are beamed active galactic nuclei known for their strong multi-wavelength variability on timescales from years down to minutes. We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019-2022, together with radio data from the WEBT and other teams, and gamma-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and gamma-ray brightness maxima. The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, and gamma-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. The gamma-ray emitting region is found to be co-spatial with the optical one, and the analysis of the gamma-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability., Comment: In press for A&A

The BL Lac object 3C 371 has been observed by the Transiting Exoplanet Survey Satellite (\textit{TESS}) for approximately a year, between July 2019 and July 2020, with an unmatched 2-minute observing cadence. In parallel, the Whole Earth Blazar Telescope (WEBT) Collaboration organized an extensive observing campaign, providing three years of continuous optical monitoring between 2018 and 2020. These datasets allow for a thorough investigation of the variability of the source. The goal of this study is evaluating the optical variability of 3C 371. Taking advantage of the remarkable cadence of \textit{TESS} data, we aim to characterize the intra-day variability (IDV) displayed by the source and identify its shortest variability timescale. With this estimate, constraints on the size of the emitting region and black hole mass can be calculated. Moreover, WEBT data is used to investigate long-term variability (LTV), including understanding spectral behaviour of the source and the polarization variability. Based on the derived characteristics, information on the origin of the variability on different timescales is extracted. We evaluated the variability applying the variability amplitude tool that quantifies how variable the emission is. Moreover, we employed common tools like ANOVA (ANalysis Of VAariance) tests, wavelet and power spectral density (PSD) analyses to characterize the shortest variability timescales present in the emission and the underlying noise affecting the data. Short- and long-term colour behaviours have been evaluated to understand the spectral behaviour of the source. The polarized emission was analyzed, studying its variability and possible rotation patterns of the electric vector position angle (EVPA). Flux distributions of IDV and LTV were also studied with the aim of linking the flux variations to turbulent and/or accretion disc related processes., Comment: Accepted in A&A, 20 pages, 19 figures

PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Context: T CrB (NOVA CrB 1946) is a famous recurrent nova with a recurrence timescale of 80 years. Aims: We aim to estimate the colours, luminosity, and mass-accretion rate for T CrB (NOVA CrB 1946) during and after the superactive state. Methods and Results: We performed and analysed $UBV$ photometry of the recurrent nova T~CrB. For the hot component of T~CrB, we find average dereddened colours of $(U-B)_0 = -0.70 \pm 0.08$ and $(B-V)_0 = 0.23 \pm 0.06$, which correspond to an effective temperature of $9400 \pm 500$~K and an optical luminosity of $40-110~L_\odot$ during the superactive state (2016-2022). After the end of the superactive state, the hot component became significantly redder, $(U-B)_0 \approx -0.3$ and $(B-V)_0 \approx 0.6$ in August 2023, and its luminosity decreased markedly to $20-25$~$L_\odot$ in April-May 2023, and to $8-9~L_\odot$ in August 2023. The total mass accreted during the superactive state from 2014 to 2023 is $\sim 2 \times 10^{-7}$~M$_\odot$. This is a significant fraction of the mass required to cause a thermonuclear runaway (TNR). Overall our results support a model in which a large accretion disc acts as a reservoir with increased accretion rate onto the central white dwarf during disc high states, ultimately leading to a TNR explosion, which now seems to be imminent., Comment: Astronomy & Astrophysics (accepted)

We present an extensive recent multi-band optical photometric observations of the blazar S5 0716+714 carried out over 53 nights with two telescopes in India, two in Bulgaria, one in Serbia, and one in Egypt during 2019 November -- 2022 December. We collected 1401, 689, 14726, and 165 photometric image frames in B, V, R, and I bands, respectively. We montiored the blazar quasi-simultaneously during 3 nights in B, V, R, and I bands; 4 nights in B, V, and R; 2 nights in V, R, and I; 5 nights in B and R; and 2 nights in V and R bands. We also took 37 nights of data only in R band. Single band data are used to study intraday flux variability and two or more bands quasi-simultaneous observations allow us to search for colour variation in the source. We employ the power-enhanced F-test and the nested ANOVA test to search for genuine flux and color variations in the light curves of the blazar on intraday timescales. Out of 12, 11, 53, and 5 nights observations, intraday variations with amplitudes between ~3% and ~20% are detected in 9, 8, 31 and 3 nights in B, V, R, and I bands, respectively, corresponding to duty cycles of 75%, 73%, 58% and 60%. These duty cycles are lower than those typically measured at earlier times. On these timescales color variations with both bluer-when-brighter and redder-when-brighter are seen, though nights with no measurable colour variation are also present. We briefly discuss possible explanations for this observed intraday variability., Comment: 19 pages, 5 figures, 4 tables, Accepted for Publication in MNRAS

We report the first extensive optical flux and spectral variability study of the TeV blazar TXS 0506+056 on intra-night to long-term timescales using BVRI data collected over 220 nights between January 21, 2017 to April 9, 2022 using 8 optical ground-based telescopes. In our search for intraday variability (IDV), we have employed two statistical analysis techniques, the nested ANOVA test and the power enhanced F-test. We found the source was variable in 8 nights out of 35 in the R-band and in 2 of 14 in the V-band yielding Duty Cycles (DC) of 22.8% and 14.3%, respectively. Clear colour variation in V - R was seen in only 1 out of 14 observing nights, but no IDV was found in the more limited B, I, and B - I data. During our monitoring period the source showed a 1.18 mag variation in the R-band and similar variations are clearly seen at all optical wavelengths. We extracted the optical (BVRI) SEDs of the blazar for 44 nights when observations were carried out in all four of those wavebands. The mean spectral index (\alpha) was determined to be 0.897+-0.171, Comment: 13 pages, 8 figures, 8 tables, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2212.09016

In 2022 the BL Lac object S4 0954+65 underwent a major variability phase, reaching its historical maximum brightness in the optical and $\gamma$-ray bands. We present optical photometric and polarimetric data acquired by the Whole Earth Blazar Telescope (WEBT) Collaboration from 2022 April 6 to July 6. Many episodes of unprecedented fast variability were detected, implying an upper limit to the size of the emitting region as low as $10^{-4}$ parsec. The WEBT data show rapid variability in both the degree and angle of polarization. We analyse different models to explain the polarization behaviour in the framework of a twisting jet model, which assumes that the long-term trend of the flux is produced by variations in the emitting region viewing angle. All the models can reproduce the average trend of the polarization degree, and can account for its general anticorrelation with the flux, but the dispersion of the data requires the presence of intrinsic mechanisms, such as turbulence, shocks, or magnetic reconnection. The WEBT optical data are compared to $\gamma$-ray data from the Fermi satellite. These are analysed with both fixed and adaptive binning procedures. We show that the strong correlation between optical and $\gamma$-ray data without measurable delay assumes different slopes in faint and high brightness states, and this is compatible with a scenario where in faint states we mainly see the imprint of the geometrical effects, while in bright states the synchrotron self-Compton process dominates., Comment: 12 pages, 13 figures, in press for MNRAS

The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $\nu_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $\nu_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $\nu_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model., Comment: Accepted for publication in Astronomy & Astrophysics

We report the flux and spectral variability of PG 1553+113 on intra-night (IDV) to short-term timescales using BVRI data collected over 91 nights from 28 February to 8 November 2019 employing ten optical telescopes: three in Bulgaria, two each in India and Serbia, and one each in Greece, Georgia, and Latvia. We monitored the blazar quasi-simultaneously for 16 nights in the V and R bands and 8 nights in the V, R, I bands and examined the light curves (LCs) for intra-day flux and colour variations using two powerful tests: the power-enhanced F-test and the nested ANOVA test. The source was found to be significantly (> 99%) variable in 4 nights out of 27 in R-band, 1 out of 16 in V-band, and 1 out of 6 nights in I-band. No temporal variations in the colours were observed on IDV timescale. During the course of these observations the total variation in R-band was 0.89 mag observed. We also investigated the spectral energy distribution (SED) using B, V, R, and I band data. We found optical spectral indices in the range of 0.878+-0.029 to 1.106+-0.065 by fitting a power law to these SEDs of PG 1553+113. We found that the source follows a bluer-when-brighter trend on IDV timescales. We discuss possible physical causes of the observed spectral variability., Comment: 13 pages,8 figures, 7 tables, Accepted for publication in MNRAS

We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $\gamma$-rays. A significant correlation ($>$3$\sigma$) between X-rays and VHE $\gamma$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $\gamma$-rays and radio, with the radio lagging by more than 100 days, placing the $\gamma$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $\gamma$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock., Comment: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina and Mets\"ahovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. The very-high-energy (VHE; E > 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance during a few days without a simultaneous increase in the X-ray flux (i.e. orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. Our multi-band correlation study also hints at an anti-correlation between UV/optical and X-ray at a significance higher than 3 sigmas. A VHE flare observed on 2017 February 4 shows gamma-ray variability on multi-hour timescales, with a factor 10 increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors., Comment: 38 pages, 23 figures, accepted for publication in Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, David Paneque

Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to investigate the spectral variability during the TESS pointing. The spectral analysis is further extended in frequency to the UV and X-ray bands with data from the Neil Gehrels Swift Observatory. We develop a new method to unveil the shortest optical variability time-scales. This is based on progressive de-trending of the TESS light curve by means of cubic spline interpolations through the binned fluxes, with decreasing time bins. The de-trended light curves are then analysed with classical tools for time-series analysis (periodogram, auto-correlation and structure functions). The results show that below 3 d there are significant characteristic variability time-scales of about 1.7, 0.5, and 0.2 d. Variability on time-scales < 0.2 d is strongly chromatic and must be ascribed to intrinsic energetic processes involving emitting regions, likely jet sub-structures, with dimension less than about milliparsec. In contrast, flux changes on time-scales > 0.5 d are quasi-achromatic and are probably due to Doppler factor changes of geometric origin., Comment: 17 pages, 18 figures, accepted by MNRAS

In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the \emph{Transiting Exoplanet Survey Satellite} (TESS) at 6,000-10,000 \AA\ with 2-minute cadence; with the Neil Gehrels \emph{Swift} satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the \emph{Fermi} Large Area Telescope at $\gamma$-ray energies; and with the Whole Earth Blazar Telescope for measurement of the optical flux density and polarization. All light curves are correlated, with similar structure on timescales from hours to days. The shortest timescale of variability at optical frequencies observed with TESS is $\sim 0.5$ hr. The most common timescale is $13\pm1$~hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multi-wavelength variability properties cannot be explained by a change solely in the Doppler factor of the emitting plasma. The polarization behavior implies that there are both ordered and turbulent components to the magnetic field in the jet. Correlation analysis indicates that the X-ray variations lag behind the $\gamma$-ray and optical light curves by up to $\sim 0.4$ days. The timescales of variability, cross-frequency lags, and polarization properties can be explained by turbulent plasma that is energized by a shock in the jet and subsequently loses energy to synchrotron and inverse Compton radiation in a magnetic field of strength $\sim3$ G, Comment: 33 pages, 25 figures, 14 tables. Accepted to ApJ

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer., Comment: 46 pages, 17 figures, 9 tables, accepted for publication in ApJ Supplements

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet., Comment: 21 pages, 18 figures, accepted for publication in MNRAS

We report the results of our optical (VRI) photometric observations of the TeV blazar 1ES 0806$+$524 on 153 nights during 2011-2019 using seven optical telescopes in Europe and Asia. We investigated the variability of the blazar on intraday as well as on long-term timescales. We examined eighteen intraday light curves for flux and color variations using the most reliable power-enhanced F-test and the nested ANOVA test. Only on one night was a small, but significant, variation found, in both $V$ band and $R$ band light curves. The $V-R$ color index was constant on every one of those nights. Flux density changes of around 80 % were seen over the course of these eight years in multiple bands. We found a weighted mean optical spectral index of 0.639$\pm$0.002 during our monitoring period by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales., Comment: 14 pages, 6 figures, 7 tables, accepted for publication in ApJ

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission., Comment: 24 pages, 18 figures, and 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We present the results of photometric observations of three TeV blazars, 3C 66A, S5 0954+658 and BL Lacertae, during the period 2013--2017. Our extensive observations were performed in a total of 360 nights which produced $\sim$6820 image frames in BVRI bands. We study flux and spectral variability of these blazars on these lengthy timescales. We also examine the optical Spectral Energy Distributions of these blazars, which are crucial in understanding the emission mechanism of long-term variability in blazars. All three TeV blazars exhibited strong flux variability during our observations. The colour variations are mildly chromatic on long timescales for two of them. The nature of the long-term variability of 3C 66A and S5 0954+658 is consistent with a model of a non-thermal variable component that has a continuous injection of relativistic electrons with power law distributions around 4.3 and 4.6, respectively. However, the long-term flux and colour variability of BL Lac suggests that these can arise from modest changes in velocities or viewing angle toward the emission region, leading to variations in the Doppler boosting of the radiation by a factor ~1.2 over the period of these observations., Comment: 13 pages, 15 figures, Accepted for publication in MNRAS

The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $\gamma$-loud blazar whose optical emission is dominated by the thermal radiation from accretion disc. 4C $+$71.07 has been detected in outburst twice by the AGILE $\gamma$-ray satellite during the period end of October - mid November 2015, when it reached a $\gamma$-ray flux of the order of $F_{\rm E>100\,MeV} = (1.2 \pm 0.3)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ and $F_{\rm E>100\,MeV} = (3.1 \pm 0.6)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, respectively, allowing us to investigate the properties of the jet and of the emission region. We investigated its spectral energy distribution by means of almost simultaneous observations covering the cm, mm, near-infrared, optical, ultra-violet, X-ray and $\gamma$-ray energy bands obtained by the GASP-WEBT Consortium, the Swift and the AGILE and Fermi satellites. The spectral energy distribution of the second $\gamma$-ray flare (the one whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about $4\times10^{47}$ erg s$^{-1}$. During the most prominent $\gamma$-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, large-mass black-hole flat-spectrum radio quasars might be good targets for future $\gamma$-ray satellites such as e-ASTROGAM., Comment: Accepted for publication in A&A. 9 pages, 4 Figures, 3 Tables

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence., Comment: 13 pages, 4 figures, 2 tables, MNRAS accepted version

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability., Comment: 16 pages, 8 figures, 7 tables; Accepted for Publication to AJ

We monitored BL Lacertae frequently during 2014 - 2016 when it was generally in a high state. We searched for intra-day variability for 43 nights using quasi-simultaneous measurements in the B, V, R, and I bands (totaling 143 light curves); the typical sampling interval was about eight minutes. On hour-like timescales, BL Lac exhibited significant variations during 13 nights in various optical bands. Significant spectral variations are seen during most of these nights such that the optical spectrum becomes bluer when brighter. The amplitude of variability is usually greater for longer observations but is lower when BL Lac is brighter. No evidence for periodicities or characteristic variability time-scales in the light curves was found. The color variations are mildly chromatic on long timescales., Comment: 15 pages, 5 Figures, 3 Tables; Accepted for publication in Galaxies; a special issue on Microvariability of Blazars

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump., Comment: 24 pages, 6 figures, 2 tables, MNRAS in press

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous., Comment: 11 pages, 4 figures, 2 tables; MNRAS in press

We obtained and analyzed more than 100 hours of multicolour optical time series of the blazar CTA 102 during its 2012 and 2016 outbursts. The object reached almost 11-th mag at the end of 2016, which is perhaps the brightest blazar state ever observed! During both outbursts, CTA 102 showed significant and rapid variability on intra-night time scales, reaching up to 0.2 mag for 30 min on some occasions. The "rms-flux" relation, built for all datasets, shows a large scatter and no apparent saturation on the magnitude scale. The ensemble structure function of the light curves can be fitted well with a straight line of a slope of $\sim$0.4. The time lags between the different optical bands appear to be consistent with zero, taking into account our time resolution. We discuss different variability scenarios and favor the changing Doppler factor of the emitting blobs as the most plausible one to account for the observed intra-night variability., Comment: MNRAS, accepted

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability., Comment: 11 pages, 6 figures, 4 tables; Accepted for publication in MNRAS

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet., Comment: 11 pages, 14 figures, accepted in MNRAS

We present our optical photometric observations of three TeV blazars, PKS 1510-089, PG 1553+113 and Mrk 501 taken using two telescopes in India, one in Bulgaria, one in Greece and one in Serbia during 2012 - 2014. These observations covered a total of 95 nights with a total of 202 B filter frames, 247 images in V band, 817 in R band while 229 images were taken in the I filter. This work is focused on multi-band flux and colour variability studies of these blazars on diverse timescales which are useful in understanding the emission mechanisms. We studied the variability characteristics of above three blazars and found all to be active over our entire observational campaigns. We also searched for any correlation between the brightness of the sources and their colour indices. During the times of variability, no significant evidence for the sources to display spectral changes correlated with magnitude was found on timescales of a few months. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux variability., Comment: 11 pages, 7 figures, 9 tables. Accepted for publication in MNRAS

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states., Comment: 32 pages, 14 figures; accepted for publication in ApJ

We analyzed the multi-band optical behaviour of the BL Lacertae object, S5 0716+714, during its outburst state from 2014 November - 2015 March. We took data on 23 nights at three observatories, one in India and two in Bulgaria, making quasi-simultaneous observations in B, V, R, and I bands. We measured multi-band optical fluxes, colour and spectral variations for this blazar on intraday and short timescales. The source was in a flaring state during the period analyzed and displayed intense variability in all wavelengths. R band magnitude of 11.6 was attained by the target on 18 Jan 2015, which is the brightest value ever recorded for S5 0716+714. The discrete correlation function method yielded good correlation between the bands with no measurable time lags, implying that radiation in these bands originate from the same region and by the same mechanism. We also used the structure function technique to look for characteristic timescales in the light curves. During the times of rapid variability, no evidence for the source to display spectral changes with magnitude was found on either of the timescales. The amplitude of variations tends to increase with increasing frequency with a maximum of $\sim$ 22% seen during flaring states in B band. A mild trend of larger variability amplitude as the source brightens was also found. We found the duty cycle of our source during the analyzed period to be $\sim$ 90%. We also investigated the optical spectral energy distribution of S5 0716+714 using B, V, R, and I data points for 21 nights. We briefly discuss physical mechanisms most likely responsible for its flux and spectral variations., Comment: 11 pages, 8 figures, 5 tables. Accepted for publication in MNRAS

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission., Comment: Accepted for publication in ApJ

We observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source., Comment: 10 pages, 9 figures, Accepted for publication in A&A

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times., Comment: 12 Pages, 5 figures, 3 Tables, Accepted for Publication in MNRAS

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability., Comment: Accepted for publication in MNRAS, 16 pages, 11 figures, 5 tables, plus supplementary material containing additional figures and tables (please contact authors for it)

We report the results of optical monitoring for a sample of 11 blazars including 10 BL Lacs and 1 Flat Spectrum Radio Quasar (FSRQ). We have measured the multiband optical flux and colour variations in these blazars on intra-day and short-term timescales of months and have limited data for 2 more blazars. These photometric observations were made during 2009 to 2011, using six optical telescopes, four in Bulgaria, one in Greece and one in India. On short-term timescales we found significant flux variations in 9 of the sources and colour variations in 3 of them. Intra-day variability was detected on 6 nights for 2 sources out of the 18 nights and 4 sources for which we collected such data. These new optical observations of these blazars plus data from our previous published papers (for 3 more blazars) were used to analyze their spectral flux distributions in the optical frequency range. Our full sample for this purpose includes 6 high-synchrotron-frequency-peaked BL Lacs (HSPs), 3 intermediate-synchrotron-frequency-peaked BL Lacs (ISPs) and 6 low-synchrotron-frequency-peaked BL Lacs (LSPs; including both BL Lacs and FSRQs). We also investigated the spectral slope variability and found that the average spectral slopes of LSPs show a good accordance with the Synchrotron Self-Compton (SSC) loss dominated model. Our analysis supports previous studies that found that the spectra of the HSPs and FSRQs have significant additional emission components. The spectra of all these HSPs and LSPs get flatter when they become brighter, while for FSRQs the opposite appears to hold. This supports the hypothesis that there is a significant thermal contribution to the optical spectrum for FSRQs., Comment: 26 pages, 7 figures, 7 tables, Accepted for Publication in MNRAS

In this paper we present results of a short-term optical monitoring of 13 blazars. The objects were monitored mostly in the R-band for a total of ~ 160 hours between 2006 and 2011. We study the nature of the short-term variations and show that most of them could be described as slow, smooth, and (almost) linear changes of up to ~ 0.1 mag/hour, but many objects show no short-term variations at all. In fact, we found only ~ 2 per cent chance to observe variability of more than 0.1 mag/hour for the sample we observed. Hints for quasi-periodic oscillations at very low amplitude levels are also found for some objects. We briefly discuss some of the possible mechanisms to generate the intra-night variability and the quasi-periodic oscillations., Comment: 10 pages, 13 figures, 1 table, Accepted for Publication in MNRAS

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high., Comment: 16 pages, 15 figures, table 2; Accepted for Publication in MNRAS

Context. 3C 454.3 is a very active flat spectrum radio quasar (blazar) that has undergone a recent outburst in all observed bands, including the optical. Aims. In this work we explore the short-term optical variability of 3C 454.3 during its outburst by searching for time delays between different optical bands. Finding one would be important for understanding the evolution of the spectrum of the relativistic electrons, which generate the synchrotron jet emission. Methods. We performed photometric monitoring of the object by repeating exposures in different optical bands (BVRI). Occasionally, different telescopes were used to monitor the object in the same band to verify the reliability of the smallest variations we observed. Results. Except on one occasion, where we found indications of a lag of the blue wavelengths behind the red ones, the results are inconclusive for most of the other cases. There were either no structures in the light curves to be able to search for patterns, or else different approaches led to different conclusions., Comment: 4 pages, 4 figures. Accepted for publication in A&A

We present the results from a monitoring campaign of the Narrow-Line Seyfert~1 galaxy PG 1211+143. The object was monitored with ground-based facilities (UBVRI photometry; from February to July, 2007) and with Swift (X-ray photometry/spectroscopy and UV/Optical photometry; between March and May, 2007). We found PG 1211+143 in a historical low X-ray flux state at the beginning of the Swift monitoring campaign in March 2007. It is seen from the light curves that while violently variable in X-rays, the quasar shows little variations in optical/UV bands. The X-ray spectrum in the low state is similar to other Narrow-Line Seyfert 1 galaxies during their low-states and can be explained by a strong partial covering absorber or by X-ray reflection onto the disk. With the current data set, however, it is not possible to distinguish between both scenarios. The interband cross-correlation functions indicate a possible reprocessing of the X-rays into the longer wavelengths, consistent with the idea of a thin accretion disk, powering the quasar. The time lags between the X-ray and the optical/UV light curves, ranging from ~2 to ~18 days for the different wavebands, scale approximately as ~lambda^(4/3), but appear to be somewhat larger than expected for this object, taking into account its accretion disk parameters. Possible implications for the location of the X-ray irradiating source are discussed., Comment: 14 pages, 11 figures, to appear in MNRAS