Flaring radio lanterns along the ridge line: long-term oscillatory motion in the jet of S5 1803+784

We present a detailed analysis of 30 very long baseline interferometric (VLBI) observations of the BL Lac object S5 1803+784 (z= 0.679), obtained between mean observational time 1994.67 and 2012.91 at observational frequency 15 GHz. The long-term behaviour of the jet ridge line reveals the jet experiences an oscillatory motion superposed on its helical jet kinematics on a time-scale of about 6 yr. The excess variance of the positional variability indicates the jet components being farther from the VLBI core have larger amplitude in their position variations. The fractional variability amplitude shows slight changes in 3 yrbins of the component's position. The temporal variability in the Doppler boosting of the ridge line results in jet regions behaving as flaring `radio lanterns'. We offer a qualitative scenario leading to the oscillation of the jet ridge line that utilizes the orbital motion of the jet emitter black hole due to a binary black hole companion. A correlation analysis implies composite origin of the flux variability of the jet components, emerging due to possibly both the evolving jet structure and its intrinsic variability

The WEBT BL Lacertae Campaign 2001 and its extension : Optical light curves and colour analysis 1994–2002

BL Lacertae has been the target of four observing campaigns by the Whole Earth Blazar Telescope (WEBT) collaboration. In this paper we present UBVRI light curves obtained by theWEBT from 1994 to 2002, including the last, extended BL Lac 2001 campaign. A total of about 7500 optical observations performed by 31 telescopes from Japan to Mexico have been collected, to be added to the ∼15 600 observations of the BL Lac Campaign 2000. All these data allow one to follow the source optical emission behaviour with unprecedented detail. The analysis of the colour indices reveals that the flux variability can be interpreted in terms of two components: longer-term variations occurring on a fewday time scale appear as mildly-chromatic events, while a strong bluer-when-brighter chromatism characterizes very fast (intraday) flares. By decoupling the two components, we quantify the degree of chromatism inferring that longer-term flux changes imply moving along a ∼0.1 bluerwhen- brighter slope in the B − R versus R plane; a steeper slope of ∼0.4 would distinguish the shorter-term variations. This means that, when considering the long-term trend, the B-band flux level is related to the R-band one according to a power law of index ∼1.1. Doppler factor variations on a “convex” spectrum could be the mechanism accounting for both the long-term variations and their slight chromatism.Reig Torres, Pablo, [email protected]

Multi-frequency monitoring of γ-ray loud blazars I. Light curves and spectral energy distributions

Context: Being dominated by non-thermal emission from aligned relativistic jets, blazars allow us to elucidate the physics of extragalactic jets, and, ltimately, how the energy is extracted from the central black hole in radio-loud active galactic nuclei. Aims: Crucial information is provided by broad-band spectral energy distributions (SEDs), their trends with luminosity and correlated multi-frequency variability. With this study we plan to obtain a database of contemporaneous radio-to-optical spectra of a sample of blazars, which are and will be observed by current and future high-energy satellites. Methods: Since December 2004 we are performing a monthly multi-frequency radio monitoring of a sample of 35 blazars at the antennas in Medicina and Noto. Contemporaneous near-IR and optical observations for all our observing epochs are organised. Results: Until June 2006 about 4000 radio measurements and 5500 near-IR and optical measurements were obtained. Most of the sources show significant variability in all observing bands. Here we present the multi-frequency data acquired during the first eighteen months of the project, and construct the SEDs for the best-sampled sources

Multiwavelength behaviour of the blazar 3C 279: Decade-long study from γ -ray to radio

We report the results of decade-long (2008-2018) γ -ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ -ray light curves correlate well, with no delay ≳ 3 h, implying general cospatiality of the emission regions. The γ -ray-optical flux-flux relation changes with activity state, ranging from a linear to amore complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz Very Long Baseline Array 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 γ -ray variability on very short time-scales. The MgII 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 τ= 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 ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. © 2020 The Author(s).We thank the referee for attentive reading and comments that helped to improve presentation of the manuscript. The data collected by the WEBT collaboration are stored in the WEBT archive at the Osservatorio Astrofisico di Torino -INAF (ht tp://www.oato.inaf.it/blazars/webt/); for questions regarding their availability, please contact the WEBT President Massimo Villata([email protected]).TheSt.Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. The research at BU was supported in part by National Science Foundation grant AST-1615796 and NASA Fermi Guest Investigator grants 80NSSC17K0649, 80NSSC19K1504, and 80NSSC19K1505. The PRISM camera at Lowell Observatory was developed by K. Janes et al. at BU and Lowell Observatory, with funding from the NSF, BU, and Lowell Observatory. The emission-line observations made use of the DCT at Lowell Observatory, supported by Discovery Communications, Inc., BU, the University of Maryland, the University of Toledo, and Northern Arizona University. The VLBA is an instrument of the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the US NSF, operated under cooperative agreement by Associated Universities, Inc. This research has used data from the UMRAO which was supported by the University of Michigan; research at this facility was supported by NASA under awards NNX09AU16G, NNX10AP16G, NNX11AO13G, and NNX13AP18G, and by the NSF under award AST-0607523. The Steward Observatory spectropolarimetric monitoring project was supported by NASA Fermi Guest Investigator grants NNX08AW56G, NNX09AU10G, NNX12AO93G, and NNX15AU81G. The Torino group acknowledges financial contribution from agreement ASI-INAF n.2017-14-H.0 and from contract PRIN-SKA-CTA-INAF 2016. I.A. acknowledges support by a Ramon y Cajal grant (RYC-2013-14511) of the 'Ministerio de Ciencia, Innovacion, y Universidades (MICIU)' of Spain and from MCIU through the 'Center of Excellence Severo Ochoa' award for the Instituto de Astrofisica de Andalucia-CSIC (SEV-20170709). Acquisition and reduction of the POLAMI and MAPCAT data were supported by MICIU through grant AYA2016-80889-P. The POLAMI observations were carried out at the IRAM 30-m Telescope, supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). The MAPCAT observations were carried out at theGerman-Spanish Calar Alto Observatory, jointly operated by the Max-Plank-Institut fur Astronomie and the Instituto de Astrofisica de Andalucia-CSIC. The study is based partly on data obtained with the STELLA robotic telescopes in Tenerife, an AIP facility jointly operated by AIP and IAC. The OVRO 40-m monitoring program is supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G, and NSF grants AST-0808050 and AST-1109911. TH was supported by the Academy of Finland projects 317383 and 320085. AZT-24 observations were made within an agreement between Pulkovo, Rome and Teramo observatories. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR/217950/16. r This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 081/2016, DN 18-13/2017, KP-06-H28/3 (2018), and KP-06-PN38/1 (2019), Bulgarian National Science Programme 'Young Scientists and Postdoctoral Students 2019', Bulgarian National Science Fund under grant DN18-10/2017 and National RI Roadmap Projects DO1-157/28.08.2018 and DO1-153/28.08.2018 of the Ministry of Education and Science of the Republic of Bulgaria. GD and OV gratefully acknowledge observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory via bilateral joint research project `Study of ICRF radio-sources and fast variable astronomical objects' (head -G. Damljanovic). This work was partly supported by the National Science Fund of the Ministry of Education and Science of Bulgaria under grant DN 08-20/2016, and by project RD-08-37/2019 of the University of Shumen. This work is a part of projects nos 176011, 176004, and 176021, supported by theMinistry of Education, Science and Technological Development of the Republic of Serbia. MGM acknowledges support through the Russian Government Program of Competitive Growth of Kazan Federal University. The Astronomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA) is managed by the Fondazione Clement Fillietroz-ONLUS, which is supported by the Regional Government of the Aosta Valley, the Town Municipality of Nus and the 'Unite des Communes vald 'otainesMont-Emilius'. The research at the OAVdA was partially funded by several `Research and Education' annual grants from Fondazione CRT. This article is partly based on observations made with the IAC80 and TCS telescopes operated by the Instituto de Astrofisica de Canarias in the Spanish Observatorio del Teide on the island of Tenerife. A part of the observations were carried out using theRATAN-600 scientific equipment (SAO of the Russian Academy of Sciences)

Multiwavelength variability of BL Lacertae measured with high time resolution

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 Transiting Exoplanet Survey Satellite (TESS) at 6000–10000 Å with 2 minute cadence; with the Neil Gehrels Swift satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the Fermi Large Area Telescope at γ-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 ~0.5 hr. The most common timescale is 13 ± 1 hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multiwavelength 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 γ-ray and optical light curves by up to ~0.4 day. 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 ~3 G.Accepted manuscrip

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES 2344+514

Aims. The BL Lac 1ES 2344+514 is known for temporary extreme properties characterised by a shift of the synchrotron spectral energy distribution (SED) peak energy νsynch, p above 1 keV. While those extreme states have only been observed during high flux levels thus far, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to very high energy (VHE) performed so far, focussing on a systematic characterisation of the intermittent extreme states. Methods. We organised a monitoring campaign covering a 3-year period from 2019 to 2021. More than ten instruments participated in the observations in order to cover the emission from radio to VHE. In particular, sensitive X-ray measurements by XMM-Newton, NuSTAR, and AstroSat took place simultaneously with multi-hour MAGIC observations, providing an unprecedented constraint of the two SED components for this blazar. Results. While our results confirm that 1ES 2344+514 typically exhibits νsynch, p > 1 keV 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 in the electron acceleration efficiency without a change in the electron injection luminosity. On the other hand, we also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of ν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. By combining Swift-XRT and Swift-UVOT measurements 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 significantly contribute 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−2 keV band. Using time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.publishedVersio

The variability patterns of the TeV blazar PG 1553 + 113 from a decade of MAGIC and multiband observations

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.publishedVersio

The complex variability of blazars: time-scales and periodicity analysis in S4 0954+65

Among active galactic nuclei, blazars show extreme variability properties. We here investigate the case of the BL Lac object S4 0954+65 with data acquired in 2019–2020 by the Transiting Exoplanet Survey Satellite (TESS) and by the Whole Earth Blazar Telescope (WEBT) Collaboration. The 2-min cadence optical light curves provided by TESS during three observing sectors of nearly 1 month each allow us to study the fast variability in great detail. We identify several characteristic short-term time-scales, ranging from a few hours to a few days. However, these are not persistent, as they differ in the various TESS sectors. The long-term photometric and polarimetric optical and radio monitoring undertaken by the WEBT brings significant additional information, revealing that (i) in the optical, long-term flux changes are almost achromatic, while the short-term ones are strongly chromatic; (ii) the radio flux variations at 37 GHz follow those in the optical with a delay of about 3 weeks; (iii) the range of variation of the polarization degree and angle is much larger in the optical than in the radio band, but the mean polarization angles are similar; (iv) the optical long-term variability is characterized by a quasi-periodicity of about 1 month. We explain the source behaviour in terms of a rotating inhomogeneous helical jet, whose pitch angle can change in time.Accepted manuscrip

The correlated optical and radio variability of BL Lacertae - WEBT data analysis 1994-2005

Since 1997, BL Lacertae has undergone a phase of high optical activity, with the occurrence of several prominent outbursts. Starting from 1999, the Whole Earth Blazar Telescope (WEBT) consortium has organized various multifrequency campaigns on this blazar, collecting tens of thousands of data points. One of the main issues in the study of this huge dataset has been the search for correlations between the optical and radio flux variations, and for possible periodicities in the light curves. The analysis of the data assembled during the first four campaigns (comprising also archival data to cover the period 1968-2003) revealed a fair optical-radio correlation in 1994-2003, with a delay of the hard radio events of ~100 days. Moreover, various statistical methods suggested the existence of a radio periodicity of ~8 years. In 2004 the WEBT started a new campaign to extend the dataset to the most recent observing seasons, in order to possibly confirm and better understand the previous results. In this campaign we have collected and assembled about 11000 new optical observations from twenty telescopes, plus near-IR and radio data at various frequencies. Here, we perform a correlation analysis on the long-term R-band and radio light curves. In general, we confirm the ~100-day delay of the hard radio events with respect to the optical ones, even if longer (~200-300 days) time lags are also found in particular periods. The radio quasi-periodicity is confirmed too, but the "period" seems to progressively lengthen from 7.4 to 9.3 years in the last three cycles. The optical and radio behaviour in the last forty years suggests a scenario where geometric effects play a major role. In particular, the alternation of enhanced and suppressed optical activity (accompanied by hard and soft radio events, respectively) ca

Multiwavelength behaviour of the blazar 3C 279: decade-long study from γ-ray to radio

We report the results of decade-long (2008–2018) γ-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ-ray light curves correlate well, with no delay ≳ 3 h, implying general cospatiality of the emission regions. The γ-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 Very Long Baseline Array 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 γ-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 τ = 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 ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.First author draf