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

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 study the optical flux and polarization variability of the binary black hole blazar OJ 287 using quasi-simultaneous observations from 2015 to 2023 carried out using telescopes in the USA, Japan, Russia, Crimea, and Bulgaria. This is one of the most extensive quasi-simultaneous optical flux and polarization variability studies of OJ 287. OJ 287 showed large amplitude, ~3.0 mag flux variability, large changes of ~37% in degree of polarization, and a large swing of ~215 degrees in the angle of the electric vector of polarization. During the period of observation, several flares in flux were detected. Those flares are correlated with a rapid increase in the degree of polarization and swings in electric vector of polarization angle. A peculiar behavior of anticorrelation between flux and polarization degree, accompanied by a nearly constant polarization angle, was detected from JD 2,458,156 to JD 2,458,292. We briefly discuss some explanations for the flux and polarization variations observed in OJ 287., Comment: 14 pages, 4 figures, 2 tables, Published in Astrophysical Journal Letters, 2023 November 1

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

Aims. We aim to investigate the extreme variability properties of the TeV blazar S4 0954+65 using optical photometric and polarisation observations carried out between 2017 and 2023 using three ground-based telescopes. Methods. We examined an extensive dataset comprised of 138 intraday (observing duration shorter than a day) light curves (LCs) of S4 0954+65 for flux, spectral, and polarisation variations on diverse timescales. For the variable LCs, we computed the minimum variability timescales. We investigated flux-flux correlations and colour variations to look for spectral variations on long (several weeks to years) timescales. Additionally, we looked for connections between optical R-band flux and polarisation degree. Results. We found significant variations in 59 out of 138 intraday LCs. We detected a maximum change of 0.58$\pm$0.11 in V-band magnitude within $\sim$2.64 h and a corresponding minimum variability timescale of 18.21$\pm$4.87 mins on 2017 March 25. During the course of our observing campaign, the source brightness changed by $\sim$4 magnitudes in V and R bands; however, we did not find any strong spectral variations. The slope of the relative spectral energy distribution was 1.37$\pm$0.04. The degree of polarisation varied from $\sim$ 3% to 39% during our monitoring. We observed a change of $\sim$120 degrees in polarisation angle (PA) within $\sim$3 h on 2022 April 13. No clear correlation was found between optical flux and the degree of polarisation. Conclusions. The results of our optical flux, colour, and polarisation study provide hints that turbulence in the relativistic jet could be responsible for the intraday optical variations in the blazar S4 0954+65. However, the long-term flux variations may be caused by changes in the Doppler factor., Comment: 9 pages, 10 figures, 4 tables, 4 appendix, Astronomy & Astrophysics journal (in press)

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

The studies and constraints on the emission region are crucial to the blazar radiation mechanism. Yet the previous works mainly focus on individual sources. In this work, we make use of the largest and the latest spectral energy distribution (SED) fitting results in the literature to statistically study the blazar emission region property in the framework of leptonic one-zone. Our results reveal (1) FSRQs show lower electron energy ($\gamma_{\rm p} \lesssim 1.6 \times 10^{3}$) than BL Lacs and tend to have a stronger magnetic field ($B$) and smaller electron-to-magnetic energy ratio ($U_{\rm e}/U_{\rm B}$) than BL Lacs; (2) we find the electro-magnetic equipartition would rather happen in the jets of BL Lacs than happen in the jets of FSRQs; (3) there are 682 blazars with a magnetic field weaker critical value of generating the Kelvin-Helmholtz instability, thus one-third of the blazars in our sample are able to produce this instability; (4) the distance ($d_{\rm em}$) between the emission region and the central black hole (BH) is in the scale of $\sim$0.1 pc, the location of the emission region may be evenly distributed inside and outside the broad line region (BLR)., Comment: 13 pages, 5 figures, 2 tables. Accepted to ApJS

In this work, we report the presence of rapid intra-night optical variations in both -- flux and polarization of the blazar BL Lacertae during its unprecedented 2020--2021 high state of brightness. The object showed significant flux variability and some color changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both -- polarization degree and angle (EVPA). The object was observed from several observatories throughout the world, covering in a total of almost 300 hours during 66 nights. Based on our results, we suggest, that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night timescales., Comment: 23 pages,7 figures, 5 Tables (2 as appendix). Accepted for publication in MNRAS

The `blazar sequence' has been proposed for more than 20 years, yet its nature is still unclear. In this work, for the first time, we expand this topic to the TeV band by using a sample of 58 TeV blazars including 48 blazars in the quiescent state and 21 blazars in the flaring state. We investigate the correlation between the TeV luminosity, which has been compensated for attenuation from extragalactic background light, and the synchrotron peak frequency. We note that there is no correlation between TeV luminosity and peak frequency in the quiescent state and a strong anti-correlation in the flaring state for the observed value. However, there is a strong positive correlation in both the quiescent state and the flaring state for the intrinsic value. This indicates that the blazar sequence is shown in the flaring state rather than in the quiescent state for the observed value and the blazar sequence is not present in both two states after removing the beaming effect. In addition, to confirm whether the beaming effect results in the blazar sequence, we compare the \textit{Fermi} $\gamma$-ray luminosity between the quiescent state and the flaring state. We find the \textit{Fermi} $\gamma$-ray luminosity in the flaring state is greater than that in the quiescent state and the Doppler factor in the flaring state is greater. We suggest the blazar sequence in the flaring state may be due to the stronger beaming effect., Comment: 19 pages, 4 figures and 4 tables. Accepted to APJ

In 2021 BL Lacertae underwent an extraordinary activity phase, which was intensively followed by the Whole Earth Blazar Telescope (WEBT) Collaboration. We present the WEBT optical data in the BVRI bands acquired at 36 observatories around the world. In mid 2021 the source showed its historical maximum, with R = 11.14. The light curves display many episodes of intraday variability, whose amplitude increases with source brightness, in agreement with a geometrical interpretation of the long-term flux behaviour. This is also supported by the long-term spectral variability, with an almost achromatic trend with brightness. In contrast, short-term variations are found to be strongly chromatic and are ascribed to energetic processes in the jet. We also analyse the optical polarimetric behaviour, finding evidence of a strong correlation between the intrinsic fast variations in flux density and those in polarisation degree, with a time delay of about 13 h. This suggests a common physical origin. The overall behaviour of the source can be interpreted as the result of two mechanisms: variability on time scales greater than several days is likely produced by orientation effects, while either shock waves propagating in the jet, or magnetic reconnection, possibly induced by kink instabilities in the jet, can explain variability on shorter time scales. The latter scenario could also account for the appearance of quasi-periodic oscillations, with periods from a few days to a few hours, during outbursts, when the jet is more closely aligned with our line of sight and the time scales are shortened by relativistic effects., Comment: 15 pages, 16 figures, submitted to MNRAS

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

BL Lacertae had undergone a series of historical high flux activity over a year, from August 2020 in the optical to VHE $\gamma$-rays. In this paper, we report on optical flux and spectral variability of the first historical maxima outburst event during October-November in g, r and i bands with the 1.26m telescope at Xinglong observatory, China. We detected significant intranight variations with amplitude rising up to $\sim 30$%, when the fastest variability timescale is found to be a few tens of minutes, giving an emitting region size of the order $10^{-3}$ pc, which corresponds to $\sim 100$ Schwarzschild radius of the central black hole, likely coming from some jet mini-structures. Unlike on intranight timescale, a clear frequency dependent pattern along with symmetric timescales ($\sim$ 11d) of flux variation are detected on long timescale. The spectral evolution was predominated by flattening of the spectra with increasing brightness i.e., a bluer-when-brighter trend in 96% of the cases. On the night before the outburst peak, the color indices clustered in two distinct branches in color--magnitude diagram within a period of $\sim$ 6 hours that is connected to a hard-soft-hard spectral evolution trend extracted from time-resolved spectra. Such trend has never seen in BL Lac or any other blazars before to the best of our knowledge. The results obtained in this study can be explained in the context of shock induced particle acceleration or magnetic re-connection in the jet where turbulent processes most likely resulted the asymmetric flux variation on nightly timescale., Comment: 15 pages, 3 tables, 8 figures. Accepted for publication in ApJ

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

The study of short-term variability properties in AGN jets has the potential to shed light on their particle acceleration and emission mechanisms. We report results from a four-day coordinated multi-wavelength campaign on the highly-peaked blazar (HBL) Mkn 421 in 2019 January. We obtained X-ray data from AstroSAT, BVRI photometry with the Whole Earth Blazar Telescope (WEBT), and TeV data from FACT to explore short-term multi-wavelength variability in this HBL. The X-ray continuum is rapidly variable on time-scales of tens of ks. Fractional variability amplitude increases with energy across the synchrotron hump, consistent with previous studies; we interpret this observation in the context of a model with multiple cells whose emission spectra contain cutoffs that follow a power-law distribution. We also performed time-averaged and time-resolved (time-scales of 6 ks) spectral fits; a broken power-law model fits all spectra well; time-resolved spectral fitting reveals the usual hardening when brightening behaviour. Intra-X-ray cross correlations yield evidence for the 0.6-0.8 keV band to likely lead the other bands by an average of 4.6 +- 2.6 ks, but only during the first half of the observation. The source displayed minimal night-to-night variability at all wavebands thus precluding significant interband correlations during our campaign. The broadband SED is modeled well with a standard one-zone leptonic model, yielding jet parameters consistent with those obtained from previous SEDs of this source., Comment: This article has been accepted for publication in The Monthly Notices of the Royal Astronomical Society (2022), and is published in Volume 513, Issue 2, pp.1662-1679. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved. 19 pages; 12 figures

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

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

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

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

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays., Comment: 25 pages, 9 figures, accepted for publication in A&A. Corresponding authors: Gareth Hughes (gareth.hughes@cfa.harvard.edu), David Paneque (dpaneque@mppmu.mpg.de), Amit Shukla (amit.shukla@astro.uni-wuerzburg.de)

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 the AGILE detection and the results of the multifrequency follow-up observations of a bright $\gamma$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $\gamma$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution during the $\gamma$-ray flare. The $\gamma$-ray flaring data, compared with as yet unpublished simultaneous optical data which allow to set constraints on the big blue bump disk luminosity, show very high Compton dominance values of $\sim 100$, with a ratio of $\gamma$-ray to optical emission rising by a factor of three in a few hours. The multi-wavelength behavior of the source during the flare challenges one-zone leptonic theoretical models. The new observations during the June 2015 flare are also compared with already published data and non-simultaneous historical 3C 279 archival data., Comment: 13 pages, 5 figures,5 Tables. Accepted for publication in ApJ

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 examine the 2008-2016 $\gamma$-ray and optical light curves of three bright BL Lac objects, 0716+714, MRK 421, BL Lac, which exhibit large structured variability. We searched for periodicities by using a fully Bayesian approach. For two out of three sources investigated no significant periodic variability was found. In the case of BL Lac we detected a periodicity of ~ 680 days. Although the signal related to this is modest, the coincidence of the periods in both gamma and optical bands is indicative of a physical relevance. Considering previous literature results, possibly related $\gamma$-ray and optical periodicities of about one year time scale are proposed in 4 bright $\gamma$-ray blazars out of the 10 examined in detail. Comparing with results from periodicity search of optical archives of quasars, the presence of quasi-periodicities in blazars might be more frequent by a large factor. This suggests the intriguing possibility that the basic conditions for their observability are related to the relativistic jet in the observer direction, but the overall picture remains uncertain., Comment: 7 pages, 8 figure; minor language corrections included. Accepted for publication in Astronomy and Astrophysics

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

We present results of optical observations of the extremely violently variable blazar S4 0954+65 on intra-night time scales. The object showed flux changes of up to 100% within a few hours. Time delays between optical bands, color changes and "rms-flux" relations are investigated and the results are discussed in terms of existing models of blazar variability., Comment: 5 pages, accepted for publication in the special issue of Galaxies, "Blazars through Sharp Multi-wavelength Eyes", edited by J. L. G\'omez, A. P. Marscher, and S. G. Jorstad

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models., Comment: 16 pages, 17 Figures; ApJ accepted

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 search for low-dimensional chaotic signatures in the optical lightcurve of the $Kepler$ field blazar W2R 1926+42. The frequently used correlation integral method is employed in our analysis. We find no apparent evidence for the presence of low-dimensional chaos in the lightcurve. If further confirmed, these results could be of importance for modeling the blazar emission mechanisms., Comment: 6 pages including 3 figures; to appear in the Proceedings of 14th Marcel Grossman Meeting (MG14), Rome, Italy, 12-18 July 2015; based on the talk given in the meeting in the session GN2

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 study the long term Kepler light curve of the blazar W2R 1926+42 ($\sim$ 1.6 years) which indicates a variety of variability properties during different intervals of observation. The normalized excess variance, $F_{\rm var}$ ranges from 1.8 % in the quiescent phase and 43.3 % in the outburst phase. We find no significant deviation from linearity in the $F_{\rm var}$-flux relation. Time series analysis is conducted using the Fourier power spectrum and the wavelet analysis methods to study the power spectral density (PSD) shape, infer characteristic timescales and statistically significant quasi-periodic oscillations (QPOs). A bending power law with an associated timescale of $T_B = 6.2^{+6.4}_{-3.1}$ hours is inferred in the PSD analysis. We obtain a black hole mass of $M_\bullet = (1.5 - 5.9) \times 10^7 M_\odot$ for the first time using $F_{\rm var}$ and the bend timescale for this source. From a mean outburst lifetime of days, we infer a distance from the jet base $r \leq 1.75$ pc indicating that the outburst originates due to a shock. A possible QPO peaked at 9.1 days and lasting 3.4 cycles is inferred from the wavelet analysis. Assuming that the QPO is a true feature, $r = (152 - 378)~ G M_\bullet/c^2$ and supported by the other timing analysis products such as a weighted mean PSD slope of $-1.5 \pm 0.2$ from the PSD analysis, we argue that the observed variability and the weak and short duration QPO could be due to jet based processes including orbital features in a relativistic helical jet and others such as shocks and turbulence., Comment: 14 pages, 6 figures, 2 tables; Accepted for publication in MNRAS

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