Showing 1–50 of 83 results for author: Kurtanidze, O

Optical variability of the blazar 3C 371: from minute to year timescales

Authors: J. Otero-Santos, C. M. Raiteri, J. A. Acosta-Pulido, M. I. Carnerero, M. Villata, S. S. Savchenko, D. Carosati, W. P. Chen, S. O. Kurtanidze, M. D. Joner, E. Semkov, T. Pursimo, E. Benítez, G. Damljanovic, G. Apolonio, G. A. Borman, V. Bozhilov, F. J. Galindo-Guil, T. S. Grishina, V. A. Hagen-Thorn, D. Hiriart, H. Y. Hsiao, S. Ibryamov, R. Z. Ivanidze, G. N. Kimeridze , et al. (19 additional authors not shown)

Abstract: 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. The… ▽ More 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. △ Less

Submitted 4 April, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: Accepted in A&A, 20 pages, 19 figures

Journal ref: A&A 686, A228 (2024)

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

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: 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 a… ▽ More 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. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Extreme photometric and polarimetric variability of blazar S4 0954+65 at its maximum optical and $γ$-ray brightness levels

Authors: C. M. Raiteri, M. Villata, M. I. Carnerero, S. S. Savchenko, S. O. Kurtanidze, V. V. Vlasyuk, A. Marchini, K. Matsumoto, C. Lorey, M. D. Joner, K. Gazeas, D. Carosati, D. O. Mirzaqulov, J. A. Acosta Pulido, I. Agudo, R. Bachev, E. Benítez, G. A. Borman, P. Calcidese, W. P. Chen, G. Damljanovic, S. A. Ehgamberdiev, D. Elsässer, M. Feige, A. Frasca , et al. (42 additional authors not shown)

Abstract: In 2022 the BL Lac object S4 0954+65 underwent a major variability phase, reaching its historical maximum brightness in the optical and $γ$-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… ▽ More In 2022 the BL Lac object S4 0954+65 underwent a major variability phase, reaching its historical maximum brightness in the optical and $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 12 pages, 13 figures, in press for MNRAS

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

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch , et al. (210 additional authors not shown)

Abstract: The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{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, f… ▽ More The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{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 $ν_{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 $ν_{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. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

The optical behaviour of BL Lacertae at its maximum brightness levels: a blend of geometry and energetics

Authors: C. M. Raiteri, M. Villata, S. G. Jorstad, A. P. Marscher, J. A. Acosta Pulido, D. Carosati, W. P. Chen, M. D. Joner, S. O. Kurtanidze, C. Lorey, A. Marchini, K. Matsumoto, D. O. Mirzaqulov, S. S. Savchenko, A. Strigachev, O. Vince

Abstract: 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… ▽ More 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. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 15 pages, 16 figures, submitted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 1, June 2023, Pages 102-116

Multi-band Optical Variability of the TeV Blazar PG 1553+113 in 2019

Authors: Vinit Dhiman, Alok C. Gupta, Sofia O. Kurtanidze, I. Eglitis, A. Strigachev, G. Damljanovic, Paul J. Wiita, Minfeng Gu, Haritma Gaur, Oliver Vince, R. Bachev, D. P. Bisen, S. Ibryamov, R. Z. Ivanidze, Miljana D. Jovanovic, Omar M. Kurtanidze, M. G. Nikolashvili, E. Semkov, B. Spassov, M. Stojanovic, Beatriz Villarroel, Haiguang Xu, Zhongli Zhang

Abstract: 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 nigh… ▽ More 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. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 13 pages,8 figures, 7 tables, Accepted for publication in MNRAS

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (300 additional authors not shown)

Abstract: 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… ▽ More 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) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-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 $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-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 $γ$-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. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, Ž. Bošnjak, G. Busetto, R. Carosi , et al. (263 additional authors not shown)

Abstract: 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, Medicin… ▽ More 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ähovi. 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. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 38 pages, 23 figures, accepted for publication in Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, David Paneque

Journal ref: A&A 655, A89 (2021)

The dual nature of blazar fast variability. Space and ground observations of S5 0716+714

Authors: C. M. Raiteri, M. Villata, D. Carosati, E. Benítez, S. O. Kurtanidze, A. C. Gupta, D. O. Mirzaqulov, F. D'Ammando, V. M. Larionov, T. Pursimo, J. A. Acosta-Pulido, G. V. Baida, B. Balmaverde, G. Bonnoli, G. A. Borman, M. I. Carnerero, W. -P. Chen, V. Dhiman, A. Di Maggio, S. A. Ehgamberdiev, D. Hiriart, G. N. Kimeridze, O. M. Kurtanidze, C. S. Lin, J. M. Lopez , et al. (17 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 17 pages, 18 figures, accepted by MNRAS

Multi-Wavelength Variability of BL Lacertae Measured with High Time Resolution

Authors: Zachary R. Weaver, K. E. Williamson, S. G. Jorstad, A. P. Marscher, V. M. Larionov, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, R. Bachev, G. V. Baida, T. J. Balonek, E. Benitez, G. A. Borman, V. Bozhilov, M. I. Carnerero, D. Carosati, W. P. Chen, G. Damljanovic, V. Dhiman, D. J. Dougherty, S. A. Ehgamberdiev, T. S. Grishina, A. C. Gupta, M. Hart, D. Hiriart , et al. (32 additional authors not shown)

Abstract: 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 Å with 2-minute cadence; with the Neil Gehrels \emph{Swift} satellite at optical, UV, and X-ray ba… ▽ More 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 Å 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 $γ$-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 $γ$-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 △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 33 pages, 25 figures, 14 tables. Accepted to ApJ

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

Authors: A. U. Abeysekara, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, J. Dumm, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, N. Galante, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (234 additional authors not shown)

Abstract: We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a co… ▽ More We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: 41 pages including 3 appendices, 13 figures; version accepted to Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 890, Number 2 (2020)

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

Authors: MAGIC collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babic, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (215 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 46 pages, 17 figures, 9 tables, accepted for publication in ApJ Supplements

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

Authors: V. M. Larionov, S. G. Jorstad, A. P. Marscher, M. Villata, C. M. Raiteri, P. S. Smith, I. Agudo, S. S. Savchenko, D. A. Morozova, J. A. Acosta-Pulido, M. F. Aller, H. D. Aller, T. S. Andreeva, A. A. Arkharov, R. Bachev, G. Bonnoli, G. A. Borman, V. Bozhilov, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, W. -P. Chen, G. Damljanovic, A. V. Dementyev , et al. (62 additional authors not shown)

Abstract: We report the results of decade-long (2008-2018) $γ$-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 $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activi… ▽ More We report the results of decade-long (2008-2018) $γ$-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 $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality 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 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 $γ$-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 $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Journal ref: MNRAS 492 (2020) 3829

Optical Variability of the TeV Blazar 1ES 0806+524 on Diverse Timescales

Authors: Ashwani Pandey, Alok C. Gupta, Sofia O. Kurtanidze, Paul J. Wiita, G. Damljanovic, R. Bachev, Jin Zhang, O. M. Kurtanidze, A. Darriba, R. A. Chigladze, G. Latev, M. G. Nikolashvili, S. Peneva, E. Semkov, A. Strigachev, S. N. Tiwari, O. Vince

Abstract: 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… ▽ More 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_ν \propto ν^{-α}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 14 pages, 6 figures, 7 tables, accepted for publication in ApJ

Investigating the multiwavelength behaviour of the flat spectrum radio quasar CTA 102 during 2013-2017

Authors: F. D'Ammando, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, I. Agudo, A. A. Arkharov, R. Bachev, G. V. Baida, E. Benitez, G. A. Borman, W. Boschin, V. Bozhilov, M. S. Butuzova, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, N. Castro-Segura, W. -P. Chen, G. Damljanovic, A. Di Paola, J. Echevarria, N. V. Efimova, Sh. A. Ehgamberdiev, C. Espinosa , et al. (72 additional authors not shown)

Abstract: 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… ▽ More 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 ($γ$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $γ$-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 $γ$-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 $γ$-ray activity is found. The $γ$-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 $γ$-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. △ Less

Submitted 11 November, 2019; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 24 pages, 18 figures, and 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Optical Variability of TeV Blazars on long time-scales

Authors: Haritma Gaur, Alok C. Gupta, R. Bachev, A. Strigachev, E. Semkov, P. J. Wiita, O. M. Kurtanidze, A. Darriba, G. Damljanovic, R. G. Chanishvili, S. Ibryamov, S. O. Kurtanidze, M. G. Nikolashvili, L. A. Sigua, O. Vince

Abstract: 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 Distributi… ▽ More 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. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 13 pages, 15 figures, Accepted for publication in MNRAS

AGILE, Fermi, Swift, and GASP-WEBT multi-wavelength observations of the high-redshift blazar 4C $+$71.07 in outburst

Authors: S. Vercellone, P. Romano, G. Piano, V. Vittorini, I. Donnarumma, P. Munar-Adrover, C. M. Raiteri, M. Villata, F. Verrecchia, F. Lucarelli, C. Pittori, A. Bulgarelli, V. Fioretti, M. Tavani, J. A. Acosta-Pulido, I. Agudo, A. A. Arkharov, U. Bach, R. Bachev, G. A. Borman, M. S. Butuzova, M. I. Carnerero, C. Casadio, G. Damljanovic, F. D'Ammando , et al. (34 additional authors not shown)

Abstract: The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $γ$-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 $γ$-ray satellite during the period end of October - mid November 2015, when it reached a $γ$-ray flux of the order of… ▽ More The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $γ$-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 $γ$-ray satellite during the period end of October - mid November 2015, when it reached a $γ$-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 $γ$-ray energy bands obtained by the GASP-WEBT Consortium, the Swift and the AGILE and Fermi satellites. The spectral energy distribution of the second $γ$-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 $γ$-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 $γ$-ray satellites such as e-ASTROGAM. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: Accepted for publication in A&A. 9 pages, 4 Figures, 3 Tables

Journal ref: A&A 621, A82 (2019)

The extreme HBL behaviour of Markarian 501 during 2012

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino , et al. (254 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 14 August, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 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)

Journal ref: A&A 620, A181 (2018)

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

Authors: E. Kun, M. Karouzos, K. É. Gabányi, S. Britzen, O. M. Kurtanidze, L. Á. Gergely

Abstract: We present a detailed analysis of 30 very long baseline interferometric 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 years. The… ▽ More We present a detailed analysis of 30 very long baseline interferometric 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 years. 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-year bins 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. △ Less

Submitted 25 April, 2018; originally announced April 2018.

Comments: 12 pages, 9 figures, 5 tables, accepted to publish in MNRAS

Journal ref: MNRAS, 478, 359 (2018)

The ever-surprising blazar OJ 287: multi-wavelength study and appearance of a new component in X-rays

Authors: Pankaj Kushwaha, Alok C. Gupta, Paul J. Wiita, Main Pal, Haritma Gaur, E. M. de Gouveia Dal Pino, O. M. Kurtanidze, E. Semkov, G. Damljanovic, S. M. Hu, M. Uemura, O. Vince, A. Darriba, M. F. Gu, R. Bachev, Xu Chen, R. Itoh, M. Kawabata, S. O. Kurtanidze, T. Nakaoka, M. G. Nikolashvili, L. A. Sigua, A. Strigachev, Z. Zhang

Abstract: 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 $γ$-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 ofte… ▽ More 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 $γ$-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. △ Less

Submitted 14 June, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: 13 pages, 4 figures, 2 tables, MNRAS accepted version

Journal ref: MNRAS, 479, 1672 - 1684 (2018)

The Bright $γ$-ray Flare of 3C 279 in June 2015: AGILE Detection and Multifrequency Follow-up Observations

Authors: C. Pittori, F. Lucarelli, F. Verrecchia, C. M. Raiteri, M. Villata, V. Vittorini, M. Tavani, S. Puccetti, M. Perri, I. Donnarumma, S. Vercellone, J. A. Acosta-Pulido, R. Bachev, E. Benitez, G. A. Borman, M. I. Carnerero, D. Carosati, W. P. Chen, Sh. A. Ehgamberdiev, A. Goded, T. S. Grishina, D. Hiriart, H. Y. Hsiao, S. G. Jorstad, G. N. Kimeridze , et al. (22 additional authors not shown)

Abstract: We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $γ$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $γ$-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 dur… ▽ More We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $γ$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $γ$-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 $γ$-ray flare. The $γ$-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 $γ$-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. △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: 13 pages, 5 figures,5 Tables. Accepted for publication in ApJ

Characterizing optical variability of OJ 287 in 2016 - 2017

Authors: Alok C. Gupta, Haritma Gaur, Paul J. Wiita, A. Pandey, P. Kushwaha, S. M. Hu, O. M. Kurtanidze, E. Semkov, G. Damljanovic, A. Goyal, M. Uemura, A. Darriba, Xu Chen, O. Vince, M. F. Gu, Z. Zhang, R. Bachev, R. Chanishvili, R. Itoh, M. Kawabata, S. O. Kurtanidze, T. Nakaoka, M. G. Nikolashvili, L. Stawarz, A. Strigachev

Abstract: 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 1… ▽ More 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. △ Less

Submitted 13 January, 2019; v1 submitted 11 March, 2018; originally announced March 2018.

Comments: 16 pages, 8 figures, 7 tables; Accepted for Publication to AJ

Journal ref: 2019, AJ, 157, 95

Multi-wavelength temporal and spectral variability of the blazar OJ 287 during and after the December 2015 flare: a major accretion disc contribution

Authors: P. Kushwaha, A. C. Gupta, P. J. Wiita, H. Gaur, E. M. de Gouveia Dal Pino, J. Bhagwan, O. M. Kurtanidze, V. M. Larionov, G. Damljanovic, M. Uemura, E. Semkov, A. Strigachev, R. Bachev, O. Vince, Minfeng Gu, Z. Zhang, T. Abe, A. Agarwal, G. A. Borman, J. H. Fan, T. S. Grishina, J. Hirochi, R. Itoh, M. Kawabata, E. N. Kopatskaya , et al. (11 additional authors not shown)

Abstract: 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 $γ$-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 wit… ▽ More 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 $γ$-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 $γ$-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 $γ$-ray emission that shows a shift in the SED peak is consistent with $γ$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump. △ Less

Submitted 14 September, 2017; originally announced September 2017.

Comments: 24 pages, 6 figures, 2 tables, MNRAS in press

Journal ref: 2018, MNRAS, 473, 1145

Radio and optical intra-day variability observations of five blazars

Authors: X. Liu, P. P. Yang, J. Liu, B. R. Liu, S. M. Hu, O. M. Kurtanidze, S. Zola, A. Kraus, T. P. Krichbaum, R. Z. Su, K. Gazeas, K. Sadakane, K. Nilson, D. E. Reichart, M. Kidger, K. Matsumoto, S. Okano, M. Siwak, J. R. Webb, T. Pursimo, F. Garcia, R. Naves Nogues, A. Erdem, F. Alicavus, T. Balonek , et al. (1 additional authors not shown)

Abstract: We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were… ▽ More We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations. △ Less

Submitted 28 April, 2017; originally announced May 2017.

Comments: 6 pages, 6 figures, accepted by MNRAS

The Variability and Period Analysis for the BL Lac AO 0235+164

Authors: J. H. Fan, O. Kurtanidze, Y. Liu, X. Liu, J. H. Yang, G. M. Richter, M. G. Nikolashvili, S. O. Kurtanidze, H. T. Wang, M. Sasada, A. Y. Zhou, C. Lin, Y. H. Yuan, Y. T. Zhang, D. Constantin

Abstract: Variability is one of the extreme observational properties of BL Lacertae objects. AO 0235+164 is a well studied BL Lac through the whole electro-magnetic wavebands. In the present work, we show its optical R band photometric observations carried out during the period of Nov, 2006 to Dec. 2012 using the Ap6E CCD camera attached to the primary focus of the $\rm 70-cm$ meniscus telescope at Abastuma… ▽ More Variability is one of the extreme observational properties of BL Lacertae objects. AO 0235+164 is a well studied BL Lac through the whole electro-magnetic wavebands. In the present work, we show its optical R band photometric observations carried out during the period of Nov, 2006 to Dec. 2012 using the Ap6E CCD camera attached to the primary focus of the $\rm 70-cm$ meniscus telescope at Abastumani Observatory, Georgia. It shows a large variation of $ΔR$ = 4.88 mag (14.19 - 19.07 mag) and a short time scale of $ΔT_v$ = 73.5 min during our monitoring period. During the period of Dec. 2006 to Nov. 2009, we made radio observations of the source using the 25-m radio telescope at Xinjiang Astronomical Observatory. When a discrete correlation function (DCF) is adopted to the optical and radio observations, we found that the optical variation leads the radio variation by 23.2$\pm$12.9 days. △ Less

Submitted 18 February, 2017; originally announced February 2017.

Comments: Accepted by ApJ

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (268 additional authors not shown)

Abstract: We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux… ▽ More We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties. △ Less

Submitted 30 December, 2016; originally announced December 2016.

Comments: 33 pages, 15 figures, accepted for publication in A&A. Corresponding authors: Marlene Doert (marlene.doert@tu-dortmund.de) and David Paneque (dpaneque@mppmu.mpg.de)

Journal ref: A&A 603, A31 (2017)

Multiband optical variability of the blazar OJ 287 during its outbursts in 2015 -- 2016

Authors: Alok C. Gupta, A. Agarwal, A. Mishra, H. Gaur, P. J. Wiita, M. F. Gu, O. M. Kurtanidze, G. Damljanovic, M. Uemura, E. Semkov, A. Strigachev, R. Bachev, O. Vince, Z. Zhang, B. Villarroel, P. Kushwaha, A. Pandey, T. Abe, R. Chanishvili, R. A. Chigladze, J. H. Fan, J. Hirochi, R. Itoh, Y. Kanda, M. Kawabata , et al. (9 additional authors not shown)

Abstract: 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 bas… ▽ More 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. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: 11 pages, 6 figures, 4 tables; Accepted for publication in MNRAS

Journal ref: MNRAS 465, 4423 (2017)

The Exceptional Flare of Mrk 501 in 2014: Combined Observations with H.E.S.S. and FACT

Authors: Gabriele Cologna, Nachiketa Chakraborty, Agnieszka Jacholkowska, Matthias Lorentz, Mahmoud Mohamed, Cedric Perennes, Carlo Romoli, Stefan J. Wagner, Alicja Wierzcholska, Daniela Dorner, Omar Kurtanidze

Abstract: The BL Lac type object Mrk 501 was observed at very high energies (E > 100 GeV) in 2014 with the upgraded H.E.S.S. (High Energy Stereoscopic System) phase 2 array. The data collected with the central 28m telescope allow for a broader energy range extending to lower energies when compared to the one obtained with the four small telescopes alone. A strong flaring event with a flux level comparable t… ▽ More The BL Lac type object Mrk 501 was observed at very high energies (E > 100 GeV) in 2014 with the upgraded H.E.S.S. (High Energy Stereoscopic System) phase 2 array. The data collected with the central 28m telescope allow for a broader energy range extending to lower energies when compared to the one obtained with the four small telescopes alone. A strong flaring event with a flux level comparable to the 1997 historical maximum has been detected as a consequence of target of opportunity observations triggered by alerts from the FACT collaboration. The First G-APD Cherenkov Telescope (FACT) is continuously monitoring bright blazars at TeV energies providing important pre- and post-flare information. For the first time, the data and lightcurves from H.E.S.S. and FACT are compared. These contemporaneous observations allow for a better characterization of the source emission. In a multiwavelength context, more precise correlation studies between VHE and lower energies are possible thanks to the dense sampling of the FACT observations. The hard intrinsic spectrum detected by H.E.S.S. during the flare allows the derivation of strong constraints on the scale of Lorentz invariance violation via the non-detection of EBL opacity modifications and from time-of-flight studies. △ Less

Submitted 12 November, 2016; originally announced November 2016.

Comments: 6 pages, 5 figures, proceeding Gamma 2016

Multifrequency Photo-polarimetric WEBT Observation Campaign on the Blazar S5 0716+714: Source Microvariability and Search for Characteristic Timescales

Authors: G. Bhatta, Ł. Stawarz, M. Ostrowski, A. Markowitz, H. Akitaya, A. A. Arkharov, R. Bachev, E. Benítez, G. A. Borman, D. Carosati, A. D. Cason, R. Chanishvili, G. Damljanovic, S. Dhalla, A. Frasca, D. Hiriart, S-M. Hu, R. Itoh, D. Jableka, S. Jorstad, M. D. Jovanovic, K. S. Kawabata, S. A. Klimanov, O. Kurtanidze, V. M. Larionov , et al. (21 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2016; originally announced August 2016.

Comments: 16 pages, 17 Figures; ApJ accepted

Long-term multi-wavelength variability and correlation study of Markarian 421 from 2007 to 2009

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (154 additional authors not shown)

Abstract: We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux v… ▽ More We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Metsähovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Metsähovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades. △ Less

Submitted 29 May, 2016; originally announced May 2016.

Comments: Corresponding authors: Ann-Kristin Overkemping (ann-kristin.overkemping@tu-dortmund.de), Marina Manganaro (manganaro@iac.es), Diego Tescaro (diego.tescaro@gmail.com), To be published in Astronomy&Astrophysics (A&A), 12 pages, 9 figures

Journal ref: A&A 593, A91 (2016)

Polarization angle swings in blazars: The case of 3C 279

Authors: S. Kiehlmann, T. Savolainen, S. G. Jorstad, K. V. Sokolovsky, F. K. Schinzel, A. P. Marscher, V. M. Larionov, I. Agudo, H. Akitaya, E. Benítez, A. Berdyugin, D. A. Blinov, N. G. Bochkarev, G. A. Borman, A. N. Burenkov, C. Casadio, V. T. Doroshenko, N. V. Efimova, Y. Fukazawa, J. L. Gómez, T. S. Grishina, V. A. Hagen-Thorn, J. Heidt, D. Hiriart, R. Itoh , et al. (39 additional authors not shown)

Abstract: Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting regio… ▽ More Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity. △ Less

Submitted 21 June, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 19 pages, 15 figures, 6 tables

Journal ref: A&A 592, C1 (2016)

Multiwavelength Study of Quiescent States of Mrk 421 with Unprecedented Hard X-Ray Coverage Provided by NuSTAR in 2013

Authors: M. Baloković, D. Paneque, G. Madejski, A. Furniss, J. Chiang, the NuSTAR team, :, M. Ajello, D. M. Alexander, D. Barret, R. Blandford, S. E. Boggs, F. E. Christensen, W. W. Craig, K. Forster, P. Giommi, B. W. Grefenstette, C. J. Hailey, F. A. Harrison, A. Hornstrup, T. Kitaguchi, J. E. Koglin, K. K. Madsen, P. H. Mao, H. Miyasaka , et al. (286 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Comments: 32 pages, 14 figures; accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 819, 156 (2016)

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

Authors: A. Furniss, K. Noda, S. Boggs, J. Chiang, F. Christensen, W. Craig, P . Giommi, C. Hailey, F. Harisson, G. Madejski, K. Nalewajko, M. Perri, D. Stern, M. Urry, F. Verrecchia, W. Zhang, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio , et al. (251 additional authors not shown)

Abstract: 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 F… ▽ More 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ähovi 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. △ Less

Submitted 24 September, 2015; v1 submitted 16 September, 2015; originally announced September 2015.

Comments: Accepted for publication in ApJ

Long term lightcurve of the BL Lac object 1ES 0229+200 at TeV energies

Authors: Gabriele Cologna, Mahmoud Mohamed, Stefan J. Wagner, Alicja Wierzcholska, Carlo Romoli, Omar Kurtanidze

Abstract: The high-frequency peaked BL Lac object 1ES 0229+200 (z = 0.14) was first detected in very high energy (VHE, E > 100 GeV) $γ$-rays by the H.E.S.S. (High Energy Stereoscopic System) collaboration in 2006. No flux variability was reported in the initial study and its spectral characteristics have been used to derive constraints on the extragalactic background light (EBL) and on the intergalactic mag… ▽ More The high-frequency peaked BL Lac object 1ES 0229+200 (z = 0.14) was first detected in very high energy (VHE, E > 100 GeV) $γ$-rays by the H.E.S.S. (High Energy Stereoscopic System) collaboration in 2006. No flux variability was reported in the initial study and its spectral characteristics have been used to derive constraints on the extragalactic background light (EBL) and on the intergalactic magnetic field (IGMF). 1ES 0229+200 has been observed with H.E.S.S. for $\sim$130 hours from 2004 to 2013: the full dataset analysed with a more sensitive method will be presented here. The results indicate that the source is not constant and displays flux variability on yearly and monthly timescales. The existence of flux variability affects the derivation of the constraints on the IGMF. The H.E.S.S. observations cover several simultaneous multi-frequency campaigns and the VHE variations are compared with those reported in different bands. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: ICRC 2015 proceeding

Optical and Radio Variability of BL Lacertae

Authors: Haritma Gaur, Alok C. Gupta, R. Bachev, A. Strigachev, E. Semkov, Paul J. Wiita, A. E. Volvach, Minfeng Gu, A. Agarwal, I. Agudo, M. F. Aller, H. D. Aller, O. M. Kurtanidze, S. O. Kurtanidze, A. Lahteenmaki, S. Peneva, M. G. Nikolashvili, L. A. Sigua, M. Tornikoski, L. N. Volvach

Abstract: 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 Lac… ▽ More 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. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 10 pages, 9 figures, Accepted for publication in A&A

Discovery of a Highly Polarized Optical Microflare in the Blazar S5 0716+714 During 2014 WEBT Campaign

Authors: Gopal Bhatta, Arti Goyal, Michal Ostrowski, Łukasz Stawarz, A. A. Arkharov, Hiroshi Akitaya, Rumen Bachev, Erika Benítez, G. A. Borman, Daniele Carosati, Andy Cason, Goran Damljanovic, Sarah Dhalla, Antonio Frasca, David Hiriart, Shao Ming Hu, Ryosuke Itoh, Damian Jableka, Svetlana Jorstad, Koji Kawabata, Sergey Klimanov, Omar Kurtanidze, Valeri Larionov, Douglas Laurence, Giuseppe Leto , et al. (20 additional authors not shown)

Abstract: The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we pr… ▽ More The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings. △ Less

Submitted 31 July, 2015; v1 submitted 30 July, 2015; originally announced July 2015.

Comments: 9 pages, 4 Figures, ApJ Letter accepted

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

Authors: MAGIC Collaboration, J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, D. Carreto Fidalgo, P. Colin, E. Colombo , et al. (249 additional authors not shown)

Abstract: We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show… ▽ More We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies. △ Less

Submitted 10 February, 2015; v1 submitted 9 February, 2015; originally announced February 2015.

Comments: Accepted for publication in A&A, 18 pages, 14 figures (v2 has a small modification in the acknowledgments, and also corrects a typo in the field "author" in the metadata)

Report number: MPP-2015-17

Journal ref: A&A 576, A126 (2015)

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Authors: The MAGIC Collaboration, J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, A. Boller, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, P. Colin, E. Colombo, J. L. Contreras , et al. (230 additional authors not shown)

Abstract: A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameter… ▽ More A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($δ$, $B$ and $R$), in addition to the parameters related to the particle population. △ Less

Submitted 14 June, 2015; v1 submitted 11 December, 2014; originally announced December 2014.

Comments: Published in A&A, 30 pages, 15 figures, 6 tables. Online-data: multi-wavelength light curves (data in Fig. 1) and broadband spectral energy distributions (the data in Figs. 6, and B1-B4) are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/578/A22. Corresponding authors: David Paneque (dpaneque@mpp.mpg.de), Shangyu Sun (sysun@mpp.mpg.de), Hajime Takami (takami@post.kek.jp)

Journal ref: A&A 578, A22 (2015)

Multiwavelength observations of Mrk 501 in 2008

Authors: MAGIC Collaboration, J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, D. Carreto Fidalgo, P. Colin, E. Colombo , et al. (237 additional authors not shown)

Abstract: Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on… ▽ More Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity. △ Less

Submitted 23 October, 2014; originally announced October 2014.

Comments: Accepted for publication in A&A

Report number: MPP-2014-381

Future mmVLBI Research with ALMA: A European vision

Authors: R. P. J. Tilanus, T. P. Krichbaum, J. A. Zensus, A. Baudry, M. Bremer, H. Falcke, G. Giovannini, R. Laing, H. J. van Langevelde, W. Vlemmings, Z. Abraham, J. Afonso, I. Agudo, A. Alberdi, J. Alcolea, D. Altamirano, S. Asadi, K. Assaf, P. Augusto, A-K. Baczko, M. Boeck, T. Boller, M. Bondi, F. Boone, G. Bourda , et al. (143 additional authors not shown)

Abstract: Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the… ▽ More Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants △ Less

Submitted 1 July, 2014; v1 submitted 18 June, 2014; originally announced June 2014.

Comments: Replaced figures 2 and 3: corrected position SRT. Corrected minor typo in 5.1

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

Authors: M. Ackermann, M. Ajello, A. Allafort, E. Antolini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. Bonamente, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, C. Cecchi, R. C. G. Chaves, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprini, R. Claus , et al. (266 additional authors not shown)

Abstract: The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifr… ▽ More The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole. △ Less

Submitted 27 June, 2014; v1 submitted 28 March, 2014; originally announced March 2014.

Comments: 46 pages, 10 figures, 4 tables. The Astrophysical Journal, Volume 786, Issue 2, article id. 157, 17 pp. Contact Authors: Filippo D'Ammando, Justin Finke, Davide Donato, Josefa Becerra Gonzalez, Tomislav Terzic. Figure 1 updated

MAGIC gamma-ray and multifrequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

Authors: MAGIC Collaboration, J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, D. Carreto-Fidalgo, P. Colin, E. Colombo, J. L. Contreras, J. Cortina , et al. (196 additional authors not shown)

Abstract: Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of… ▽ More Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock. △ Less

Submitted 3 July, 2014; v1 submitted 22 January, 2014; originally announced January 2014.

Comments: 21 pages, 15 figures, A&A in press

Analyzing polarization swings in 3C 279

Authors: S. Kiehlmann, T. Savolainen, S. G. Jorstad, K. V. Sokolovsky, F. K. Schinzel, I. Agudo, A. A. Arkharov, E. Benitez, A. Berdyugin, D. A. Blinov, N. G. Bochkarev, G. A. Borman, A. N. Burenkov, C. Casadio, V. T. Doroshenko, N. V. Efimova, Y. Fukazawa, J. L. Gomez, V. A. Hagen-Thorn, J. Heidt, D. Hiriart, R. Itoh, M. Joshi, G. N. Kimeridze, T. S. Konstantinova , et al. (31 additional authors not shown)

Abstract: Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of b… ▽ More Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state. △ Less

Submitted 15 November, 2013; v1 submitted 13 November, 2013; originally announced November 2013.

Comments: 6 pages, 2 figures, 4 tables, proceeding to appear on EPJ Web of Conferences "The Innermost Regions of Relativistic Jets and their Magnetic Fields (Granada, Spain,10-14 June 2013)"

The 72-Hour WEBT Microvariability Observation of Blazar S5 0716+714 in 2009

Authors: G. Bhatta, J. R. Webb, H. Hollingsworth, S. Dhalla, A. Khanuja, R. Bachev, D. A. Blinov, M. Bottcher, O. J. A. Bravo Calle, P. Calcidese, D. Capezzali, D. Carosati, R. Chigladze, A. Collins, J. M. Coloma, Y. Efimov, A. C. Gupta, S-M. Hu, O. Kurtanidze, A. Lamerato, V. M. Larionov, C. -U. Lee, E. Lindfors, B. Murphy, K. Nilsson , et al. (20 additional authors not shown)

Abstract: Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous op… ▽ More Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model. △ Less

Submitted 17 October, 2013; originally announced October 2013.

Journal ref: 2013A&A...558A..92B

A Tight Connection between Gamma-Ray Outbursts and Parsec-Scale Jet Activity in the Quasar 3C 454.3

Authors: Svetlana G. Jorstad, Alan P. Marscher, Paul S. Smith, Valeri M. Larionov, Iván Agudo, Mark Gurwell, Ann E. Wehrle, Anne Lähteenmäki, Maria G. Nikolashvili, Gary D. Schmidt, Arkady A. Arkharov, Dmitry A. Blinov, Kelly Blumenthal, Carolina Casadio, Revaz A. Chigladze, Natatia V. Efimova, Joseph R. Eggen, José L. Gómez, Dirk Grupe, Vladimir A. Hagen-Thorn, Manasvita Joshi, Givi N. Kimeridze, Tatiana S. Konstantinova, Evgenia N. Kopatskaya, Omar M. Kurtanidze , et al. (20 additional authors not shown)

Abstract: We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent γ-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each γ-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events… ▽ More We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent γ-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each γ-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the γ-ray outbursts, although the γ-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the γ-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the γ-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that γ-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells. △ Less

Submitted 9 July, 2013; originally announced July 2013.

Comments: 80 pages, 10 Tables, 31 figures, accepted for publication in ApJ

Long-Term Multiwavelength Studies of High-Redshift Blazar 0836+710

Authors: A. Akyuz, D. J. Thompson, D. Donato, J. S. Perkins, L. Fuhrmann, E. Angelakis, J. A. Zensus, S. Larsson, K. Sokolovsky, O. Kurtanidze

Abstract: Aims. The observation of gamma -ray flares from blazar 0836+710 in 2011, following a period of quiescence, offered an opportunity to study correlated activity at different wavelengths for a high-redshift (z=2.218) active galactic nucleus. Methods. Optical and radio monitoring, plus Fermi-LAT gamma-ray monitoring provided 2008-2012 coverage, while Swift offered auxiliary optical, ultraviolet, and X… ▽ More Aims. The observation of gamma -ray flares from blazar 0836+710 in 2011, following a period of quiescence, offered an opportunity to study correlated activity at different wavelengths for a high-redshift (z=2.218) active galactic nucleus. Methods. Optical and radio monitoring, plus Fermi-LAT gamma-ray monitoring provided 2008-2012 coverage, while Swift offered auxiliary optical, ultraviolet, and X-ray information. Other contemporaneous observations were used to construct a broad-band spectral energy distribution. Results. There is evidence of correlation but not a measurable lag between the optical and gamma-ray flaring emission. On the contrary, there is no clear correlation between radio and gamma-ray activity, indicating radio emission regions that are unrelated to the parts of the jet that produce the gamma-rays. The gamma-ray energy spectrum is unusual in showing a change of shape from a power law to a curved spectrum when going from the quiescent state to the active state. △ Less

Submitted 1 July, 2013; originally announced July 2013.

Comments: 11 pages, 10 figures, Accepted for publication in A&A

Locating the Gamma-ray Flaring Emission of Blazar AO 0235+164 in the Jet at Parsec Scales Through Multi Spectral Range Monitoring

Authors: Ivan Agudo, Alan P. Marscher, Svetlana G. Jorstad, Valeri M. Larionov, Jose L. Gomez, Anne Lahteenmaki, Paul S. Smith, Kari Nilsson, Anthony C. S. Readhead, Margo F. Aller, Jochen Heidt, Mark Gurwell, Clemens Thum, Ann E. Wehrle, Omar M. Kurtanidze

Abstract: We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wav… ▽ More We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary core and in the superluminal knot, both at >12 parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes. △ Less

Submitted 8 March, 2013; originally announced March 2013.

Comments: 2012 Fermi Symposium proceedings - eConf C121028

Multiwavelength Intraday Variability of the BL Lac S5 0716+714

Authors: Alok C. Gupta, T. P. Krichbaum, P. J. Wiita, B. Rani, K. V. Sokolovsky, P. Mohan, A. Mangalam, N. Marchili, L. Fuhrmann, I. Agudo, U. Bach, R. Bachev, M. Böttcher, K. E. Gabanyi, H. Gaur, K. Hawkins, G. N. Kimeridze, O. M. Kurtanidze, S. O. Kurtanidze, C. -U. Lee, X. Liu, B. McBreen, R. Nesci, G. Nestoras, M. G. Nikolashvili , et al. (9 additional authors not shown)

Abstract: 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 contributio… ▽ More 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. △ Less

Submitted 25 July, 2012; originally announced July 2012.

Comments: 16 pages, 15 figures, table 2; Accepted for Publication in MNRAS

Journal ref: MNRAS, 425, 1357 (2012)

Variability of the blazar 4C 38.41 (B3 1633+382) from GHz frequencies to GeV energies

Authors: C. M. Raiteri, M. Villata, P. S. Smith, V. M. Larionov, J. A. Acosta-Pulido, M. F. Aller, F. D'Ammando, M. A. Gurwell, S. G. Jorstad, M. Joshi, O. M. Kurtanidze, A. Lähteenmäki, D. O. Mirzaqulov, I. Agudo, H. D. Aller, M. J. Arévalo, A. A. Arkharov, U. Bach, E. Benítez, A. Berdyugin, D. A. Blinov, K. Blumenthal, C. S. Buemi, A. Bueno, T. M. Carleton , et al. (52 additional authors not shown)

Abstract: The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energ… ▽ More The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and γ-ray flux apparently fades in time, likely because of an increasing optical to γ-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor δ. Under the hypothesis that in the period 2008-2011 all the γ-ray and optical variability on a one-week timescale were due to changes in δ, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle θ only, then θ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model. △ Less

Submitted 17 July, 2012; originally announced July 2012.

Comments: 19 pages, 23 figures, in press for Astronomy and Astrophysics

Journal ref: A&A 545, A48 (2012)

Multi-wavelength observations of blazar AO 0235+164 in the 2008-2009 flaring state

Authors: M. Ackermann, M. Ajello, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, R. D. Blandford, E. D. Bloom, E. Bonamente, A. W. Borgland, E. Bottacini, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang , et al. (186 additional authors not shown)

Abstract: The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the… ▽ More The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the γ -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus. △ Less

Submitted 12 July, 2012; originally announced July 2012.

Comments: 48 pages, 8 figures

Journal ref: ApJ 751, 159 (2012)