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

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

Multiwavelength variability of BL Lacertae measured with high time resolution

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

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

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

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

We report the results of decade-long (2008–2018) γ-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ-ray light curves correlate well, with no delay ≳ 3 h, implying general cospatiality of the emission regions. The γ-ray–optical flux–flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz Very Long Baseline Array images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain γ-ray variability on very short time-scales. The Mg ii emission line flux in the ‘blue’ and ‘red’ wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands, we find progressive delays of the most prominent light-curve maxima with decreasing frequency, as expected from the frequency dependence of the τ = 1 surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.First author draf

AGILE detection of extreme γ -ray activity from the blazar PKS 1510-089 during March 2009: Multifrequency analysis

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak

The Great Markarian 421 Flare of 2010 February: Multiwavelength Variability and Correlation Studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of 2010 February, when an extraordinary flare reaching a level of ∼27 Crab Units above 1 TeV was measured in very high energy (VHE) γ-rays with the Very Energetic Radiation Imaging Telescope Array System (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, Major Atmospheric Gamma-ray Imaging Cherenkov), high-energy γ-ray (Fermi-LAT), X-ray (Swift, Rossi X-ray Timing Experiment, MAXI), optical (including the GASP-WEBT collaboration and polarization data), and radio (Metsahovi, Owens Valley Radio Observatory, University of Michigan Radio Astronomy Observatory). 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 ∼25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ ⪆ 33) and the size of the emission region (δ-1RB≲ 3.8 × 1013cm) 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 anticorrelation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain with the classic single-zone synchrotron self-Compton model.</p

Insights into the high-energy γ-ray emission of Markarian 501 from extensive multifrequency observations in the Fermi era

We report on the γ-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) γ-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 ± 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size ≲0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (≃1044 erg s-1) constitutes only a small fraction (∼10-3) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude. © 2011. The American Astronomical Society

Optical photometry of the quasar 3C 454.3 during the period 2006–2018 and the long-term periodicity analysis

In this work, we present 8523 pairs of R-band optical photometry observations for the quasar 3C 454.3 made during the period of 2006 October–2018 February on the 70 cm meniscus telescope at Abastumani Observatory, Georgia, to study its intraday variabilities (IDVs) and long-term variations, and we have come to the following results. (1) We detected 10 outbursts, a ΔR = 3.825 mag variation, and some IDVs. The IDV timescales are from 4.1 to 285 minutes, with the corresponding variability amplitude being A = 2.9%–43.67%. The amplitude increases with IDV timescale. (2) The largest variation over a 1 day timescale is ΔR = 1.38 mag. (3) The IDV timescales suggest that the emission sizes are from 8.9 × 1013 cm to 6.20 × 1015 cm, and the magnetic field strengths are B = 0.18–0.79 G. (4) Period analysis results show three possible long-term periods, p = 3.04 ± 0.02 yr, p = 1.66 ± 0.06 yr, and p = 1.20 ± 0.03 yr in the optical light curve. We adopted the accretion disk models and the lighthouse models to period p = 3.04 ± 0.02 yr: in the accretion disk models, the binary black holes have masses M = 1.17 × 109Me; in the lighthouse models, we used two boosted jet flux densities to fit the observational light curve. (5) WWZ analysis gives some short-period (high-frequency) signals associated with strong bursts (JD 2,454,302 and JD 2,454,521) with variable frequencies and lasting for the entire observation time span (11.3 yr)

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

© 2020. The American Astronomical Society. All rights reserved. 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 18 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 ± 0.002 during our monitoring period by fitting a power law (F ν ∝ ν-α) in 23 optical (VRI) spectral energy distributions of 1ES 0806+524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales

Optical and radio variability of BL Lacertae

VK: BIBCODE: 2015A&A...582A.103G; DOI: 10.1051/0004-6361/201526536; eprintid: arXiv:1508.00959Peer reviewe