Search

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 ± 63) × 10−8 ph cm−2 s−1, corresponding to a luminosity of (2.2 ± 0.1) × 1050 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. © 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, The data collected by the WEBT collaboration are stored in the WEBT archive at the Osservatorio Astrofisico di Torino INAF (http://www.oato.inaf.it/blazars/webt/); for questions regarding their availability, contact the WEBT President Massimo Villata(massimo.villata@inaf.it).We acknowledge financial contribution from the agreementASI-INAFn. 2017-14-H.0 and from the contract PRIN-SKA-CTA-INAF 2016. This research was partially supported by the Bulgarian National Science Fund of the Ministry of Education and Science under grants DN 08-1/2016, DN 18-13/2017, and KP-06-H28/3 (2018). The Skinakas Observatory is a collaborative project of the University of Crete, the Foundation for Research and Technology -Hellas, and the Max-Planck-Institut fur Extraterrestrische Physik. The Abastumani team acknowledges financial support by the Shota Rustaveli National Science Foundation under contract FR/217554/16. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. GD and OV gratefully acknowledge the observing grant support from the Institute of Astronomy and Rozhen National Astronomical Observatory, Bulgarian Academy of Sciences, via bilateral joint research project 'Study of ICRF radio-sources and fast variable astronomical objects' (head -G.Damljanovic). This work is a part of the Projects No. 176011 ('Dynamics and Kinematics of Celestial Bodies and Systems'), No. 176004 ('Stellar Physics'), and No. 176021 ('Visible and Invisible Matter in Nearby Galaxies: Theory and Observations') supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia. JE is indebted to DGAPA (Universidad Nacional Autonoma de M ' exico) for financial support, PAPIIT project IN114917. 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. Data from the Steward Observatory blazar monitoring project were used. This program is supported by NASA/Fermi Guest Investigator grants NNX12AO93G and NNX15AU81G. We acknowledge support by 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. This publication makes use of data obtained at Mets ahovi Radio Observatory, operated by Aalto University in Finland. The Astronomical Observatory of the Autonomous Region of theAostaValley (OAVdA) is managed by the Fondazione Clment Fillietroz-ONLUS, which is supported by the Regional Government of the Aosta Valley, the TownMunicipality of Nus and the 'Unit des Communes valdtaines Mont-milius'. The research at the OAVdA was partially funded by two 'Research and Education' grants from Fondazione CRT. RR acknowledges support from CONICYT project Basal AFB-170002. MM and TM acknowledge support through the Russian Government Program of Competitive Growth of Kazan Federal University. The Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat a l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes Spatiales in France. This work performed in part under DOE Contract DE-AC02-76SF00515. 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. We thank the Swift team for making these observations possible, the duty scientists, and science planners. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. FD thanks S. Covino for his help with the REM data reduction. This research has made use of data obtained from the high-energy Astrophysics Science Archive Research Center (HEASARC) provided by NASA's Goddard Space Flight Center.

Blazar emission is dominated by non-thermal radiation from a relativistic jet pointing toward us, therefore undergoing Doppler beaming. This is responsible for flux enhancement and contraction of the variability time scales, so that most blazars appear as luminous sources characterized by noticeable and fast flux changes at all frequencies. The mechanisms producing their unpredictable variability are debated and include injection, acceleration and cooling of particles, with possible intervention of shock waves or turbulence. Changes in the viewing angle of the emitting knots or jet regions have also been suggested to explain flaring events or specific properties such as intraday variability, quasi-periodicities, or the delay of radio flux variations relative to optical changes. However, such a geometric interpretation has not been universally accepted because alternative explanations based on changes of physical conditions can also work in many cases. Here we report the results of optical-to-radio monitoring of the blazar CTA 102 by the Whole Earth Blazar Telescope Collaboration and show that the observed long-term flux and spectral variability is best explained by an inhomogeneous, curved jet that undergoes orientation changes. We propose that magnetohydrodynamic instabilities or rotation of a twisted jet cause different jet regions to change their orientation and hence their relative Doppler factors. In particular, the recent extreme optical outburst (six magnitudes) occurred when the corresponding jet emitting region acquired a minimum viewing angle., Comment: 31 pages, 11 figures. This is the accepted version (before Nature editing) of the paper published by Nature, online version at http://dx.doi.org/10.1038/nature24623, Nature 2017/12/04/online

We have been performing multi-wavelength monitoring of a sample of γ -ray blazars since the launch of the Fermi Gamma-ray Space Telescope in 2008. We present γ -ray and optical light curves for several quasars and BL Lac objects from the sample to illustrate different patterns of variability. We investigate correlations between γ -ray and R-band light curves and, if these are statistically significant, determine delays between variations at the two wavebands. Such time delays can reveal the relative locations of the emitting regions in AGN jets and the origin of the high-energy photons. We present preliminary results of this analysis. Of the 29 blazars with sufficient time coverage, 17 display a significant, singular, correlated time lag when tested over the entire 7-year period. Of these sources, the six that exhibit a consistent time lag across a majority of epochs of high activity have lags of 0 ± 7 days; the 11 without consistency across epochs of high activity generally display longer mean lags, with γ -ray leading optical. Eleven sources display no significant singular correlation over either the entire 7-year period or across shorter intervals. No significant difference is apparent between the BL Lac objects and FSRQs. Even after 7 years of monitoring, our correlation analysis remains plagued with uncertainties due to insufficient data.

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ 287 as revealed by multi-epoch Very Long Baseline Array images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave >core> or, at other epochs, coincident with a bright, quasi-stationary emission feature ∼0.2 mas (∼0.9 pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and γ-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies ≥10 pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies.© 2018. The American Astronomical Society. All rights reserved.., The authors thank Dr. Ioannis Myserlis for a critical review of a draft of this manuscript. The authors acknowledge financial support by MEXT/JSPS Grant-in-Aid for Scientific Research on Innovative Areas (No. 25120007). M.S. was supported during this study by a JSPS Postdoctoral Fellowship for Research Abroad. The research at Boston University was supported in part by National Science Foundation grant AST-1615796 and NASA Fermi Guest Investigator Program grants NNX14AQ58G and 80NSSC17K0649. The St. Petersburg University team acknowledges support from Russian Science Foundation grant 17-12-01029. I.A. acknowledges support by a Ramon y Cajal grant of the Ministerio de Economia y Competitividad (MINECO) of Spain. The research at the IAA-CSIC was supported in part by the MINECO through grants AYA2016-80889-P, AYA2013-40825-P, and AYA2010-14844; and by the regional government of Andalucia through grant P09-FQM-4784. Calar Alto Observatory is jointly operated by the Max-Planck-Institut fur Astronomie and the IAA-CSIC. The VLBA is an instrument of the Long Baseline Observatory (LBO). The LBO is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. Data from the Steward Observatory spectro-polarimetric monitoring project were used. This program was supported by Fermi Guest Investigator grants NNX08AW56G, NNX09AU10G, NNX12AO93G, and NNX15AU81G.

We present $\gamma$-ray, X-ray, ultraviolet, optical, and near-infrared light curves of 33 $\gamma$-ray bright blazars over four years that we have been monitoring since 2008 August with multiple optical, ground-based telescopes and the Swift satellite, and augmented by data from the Fermi Gamma-ray Space Telescope and other publicly available data from Swift. The sample consists of 21 flat-spectrum radio quasars (FSRQs) and 12 BL Lac objects (BL Lacs). We identify quiescent and active states of the sources based on their $\gamma$-ray behavior. We derive $\gamma$-ray, X-ray, and optical spectral indices, $\alpha_\gamma$, $\alpha_X$, and $\alpha_o$, respectively ($F_\nu\propto\nu^\alpha$), and construct spectral energy distributions (SEDs) during quiescent and active states. We analyze the relationships between different spectral indices, blazar classes, and activity states. We find (i) significantly steeper $\gamma$-ray spectra of FSRQs than for BL Lacs during quiescent states, but a flattening of the spectra for FSRQs during active states while the BL Lacs show no significant change; (ii) a small difference of $\alpha_X$ within each class between states, with BL Lac X-ray spectra significantly steeper than in FSRQs; (iii) a highly peaked distribution of X-ray spectral slopes of FSRQs at $\sim-$0.60, but a very broad distribution of $\alpha_X$ of BL Lacs during active states; (iv) flattening of the optical spectra of FSRQs during quiescent states, but no statistically significant change of $\alpha_o$ of BL Lacs between states; and (v) a positive correlation between optical and $\gamma$-ray spectral slopes of BL Lacs, with similar values of the slopes. We discuss the findings with respect to the relative prominence of different components of high-energy and optical emission as the flux state changes., Comment: 40 pages, 28 figures, 10 Tables. Additional figures at http://www.bu.edu/blazars/VLBAproject.html. Accepted by ApJ

A possible source of $\gamma$-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four $\gamma$-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510$-$089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature $\sim1200$ K, plus a much weaker optically thin component at $\sim660$ K. The total luminosity of the thermal dust emission is $7.9\pm0.2 \times 10^{45}$ erg s$^{-1}$. If the dust lies in an equatorial torus, the density of IR photons from the torus is sufficient to explain the $\gamma$-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is $7 \pm 2 \times 10^{45}$ erg s$^{-1}$. However, in CTA102 the far-IR spectra are too noisy to detect the $10 \mu$m silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510-089, and ON231, are, $2.3\times10^{45}$, and $6.6\times10^{43}$ erg s$^{-1}$, respectively. These upper limits do not rule out the possibility of inverse Compton up-scattering of IR photons to $\gamma$-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust., Comment: 23 Pages, 5 Figures, 2 Tables, 1 Machine Readable Table. Accepted to ApJ