A. Sievers

Simultaneous Planck, Swift, and Fermi observations of X-ray and γ-ray selected blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and -ray bands, and we compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi Large Area Telescope (LAT), whereas 30 to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the -ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, withh i 0 up to about 70 GHz, above which it steepens toh i 0:65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency ( S ) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples withh S i = 10 13:1 0:1 Hz, while the mean inverse-Compton peak frequency,h IC i, ranges from 10 21 to 10 22 Hz. The distributions of S and of IC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with -ray selected blazars peaking at 7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection e ect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the -ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and S predicted by the blazar sequence.

Detection of significant cm to sub-mm band radio and γ-ray correlated variability in Fermi bright blazars

The exact location of the γ-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA programme) and γ-ray (0.1–300 GeV) ∼3.5 yr light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multiband radio and γ-ray correlations (radio lagging γ rays) when averaging over the whole sample, (ii) average time delays (source frame: 76 ± 23 to 7 ± 9 d), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence τ_(r, γ)(ν) ∝ ν^(−1), in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk γ-ray production region typically located within/upstream of the 3 mm core region (τ_(3mm, γ) = 12 ± 8 d), (iv) mean distances between the region of γ-ray peak emission and the radio ‘τ = 1 photosphere’ decreasing from 9.8 ± 3.0 pc (11 cm) to 0.9 ± 1.1 pc (2 mm) and 1.4 ± 0.8 pc (0.8 mm), (v) 3 mm/γ-ray correlations in nine individual sources at a significance level where one is expected by chance (probability: 4 × 10^(−6)), (vi) opacity and ‘time lag core shift’ estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk γ-ray production region from the supermassive black hole (SMBH) of ∼0.8–1.6 pc, i.e. at the outer edge of the broad-line region (BLR) or beyond. A 3 mm τ = 1 surface at ∼2–3 pc from the jet base (i.e. well outside the ‘canonical BLR’) finally suggests that BLR material extends to several parsec distances from the SMBH.

Radio and γ-ray follow-up of the exceptionally high-activity state of PKS 1510−089 in 2011

We investigate the radio and γ-ray variability of the flat spectrum radio quasar PKS 1510−089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major γ-ray flares detected in 2011 July and October. During the latter episode both the γ-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380° of the optical polarization angle close in time with the rapid and strong γ-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands are observed about three weeks later, close in time with another γ-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimetre regime followed with some time delay at centimetre down to decimetre wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the γ-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the γ-ray variability is not within the broad line, but a few parsecs downstream along the jet.

Radio jet emission from GeV-emitting narrow-line Seyfert 1 galaxies

We studied the radio emission from four radio-loud and gamma-ray-loud narrow-line Seyfert 1 galaxies. The goal was to investigate whether a relativistic jet is operating at the source, and quantify its characteristics. We relied on the most systematic monitoring of such system in the cm and mm radio bands which is conducted with the Effelsberg 100 m and IRAM 30 m telescopes and covers the longest time-baselines and the most radio frequencies to date. We extract variability parameters and compute variability brightness temperatures and Doppler factors. The jet powers were computed from the light curves to estimate the energy output. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. All the sources display intensive variability that occurs at a pace faster than what is commonly seen in blazars. The flaring events show intensive spectral evolution indicative of shock evolution. The brightness temperatures and Doppler factors are moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case clearly implies a quiescent jet underlying the recursive flaring activity. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets.

The F-GAMMA programme: multi-frequency study of active galactic nuclei in the Fermi era - Programme description and the first 2.5 years of monitoring

Context. To fully exploit the scientific potential of the Fermi mission for the physics of active galactic nuclei (AGN), we initiated the F-GAMMA programme. Between 2007 and 2015 the F-GAMMA was the prime provider of complementary multi-frequency monitoring in the radio regime. Aims. We quantify the radio variability of γ-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the γ-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. Methods. The F-GAMMA performed monthly monitoring of a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-yr data set to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and γ-ray fluxes. Results. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100 m and IRAM 30 m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60–80 GHz as expected by shock-in-jet models. Compared to flat-spectrum radio quasars (FSRQs), BL Lacertae objects (BL Lacs) show systematically lower variability amplitudes, brightness temperatures, and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes, brightness temperatures, and Doppler factors than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3σ, implying that γ rays are produced very close to the mm-band emission region.

PKS 1502+106: A high-redshift Fermi blazar at extreme angular resolution - Structural dynamics with VLBI imaging up to 86 GHz

Context. Blazars are among the most energetic objects in the Universe. In 2008 August, Fermi/LAT detected the blazar PKS 1502+106 showing a rapid and strong gamma-ray outburst followed by high and variable flux over the next months. This activity at high energies triggered an intensive multi-wavelength campaign covering also the radio, optical, UV, and X-ray bands indicating that the flare was accompanied by a simultaneous outburst at optical/UV/X-rays and a delayed outburst at radio bands. Aims: In the current work we explore the phenomenology and physical conditions within the ultra-relativistic jet of the gamma-ray blazar PKS 1502+106. Additionally, we address the question of the spatial localization of the MeV/GeV-emitting region of the source. Methods: We utilize ultra-high angular resolution mm-VLBI observations at 43 and 86 GHz complemented by VLBI observations at 15 GHz. We also employ single-dish radio data from the F-GAMMA program at frequencies matching the VLBI monitoring. Results: PKS 1502+106 shows a compact core-jet morphology and fast superluminal motion with apparent speeds in the range 5--22 c. Estimation of Doppler factors along the jet yield values between ~7 up to ~50. This Doppler factor gradient implies an accelerating jet. The viewing angle towards the source differs between the inner and outer jet, with the former at ~3 degrees and the latter at ~1 degree, after the jet bends towards the observer beyond 1 mas. The de-projected opening angle of the ultra-fast, magnetically-dominated jet is found to be (3.8 +/- 0.5) degrees. A single jet component can be associated with the pronounced flare both at high-energies and in radio bands. Finally, the gamma-ray emission region is localized at less than 5.9 pc away from the jet base.

F-GAMMA: variability Doppler factors of blazars from multiwavelength monitoring

Recent population studies have shown that the variability Doppler factors can adequately describe blazars as a population. We use the flux density variations found within the extensive radio multi-wavelength datasets of the F-GAMMA program, a total of 10 frequencies from 2.64 up to 142.33 GHz, in order to estimate the variability Doppler factors for 58

Exploring the nature of the broadband variability in the flat spectrum radio quasar 3C 273

\gamma

What can the 2008/10 broadband flare of PKS 1502+106 tell us? - Nuclear opacity, magnetic fields, and the location of γ rays

-ray bright sources, for 20 of which no variability Doppler factor has been estimated before. We employ specifically designed algorithms in order to obtain a model for each flare at each frequency. We then identify each event and track its evolution through all the available frequencies for each source. This approach allows us to distinguish significant events producing flares from stochastic variability in blazar jets. It also allows us to effectively constrain the variability brightness temperature and hence the variability Doppler factor as well as provide error estimates. Our method can produce the most accurate (16\% error on average) estimates in the literature to date.

POLAMI: Polarimetric Monitoring of Active Galactic Nuclei at Millimetre Wavelengths – II. Widespread circular polarization

The detailed investigation of the broadband flux variability in the blazar 3C 273 allowed us to probe the location and size of emission regions and their physical conditions. We conducted correlation studies of the flaring activity in 3C 273, which was observed for the period between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function methods. Starting from the commonly used power spectral density (PSD) analysis at X-ray frequencies, we extended our investigation to characterise the nature of variability at radio, optical, and γ -ray frequencies. The PSD analysis showed that the optical and infrared light-curve slopes are consistent with the slope of white-noise processes, while the PSD slopes at radio, X-ray, and γ -ray energies are consistent with red-noise processes. We found that the estimated fractional variability amplitudes strongly depend on the observed frequency. The flux variations at γ -ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110 ± 27) days between γ -ray and radio light-curves, where γ -ray variations lead the radio bands, we constrained the location of the γ -ray emission region at a de-projected distance of 1.2 ± 0.9 pc from the jet apex. Flux variations at X-ray bands were found to have a significant correlation with variations at both radio and γ -ray energies. The correlation between X-ray and γ -ray light curves indicates two possible time lags, which suggests that two components are responsible for the X-ray emission. A negative time lag of −(50 ± 20) days, where the X-rays are leading the emission, suggests that X-rays are emitted closer to the jet apex from a compact region (0.02–0.05 pc in size), most likely from the corona at a distance of (0.5 ± 0.4) pc from the jet apex. A positive time lag of (110 ± 20) days ( γ -rays are leading the emission) suggests a jet-base origin of the other X-ray component at ~4 to 5 pc from the jet apex. The flux variations at radio frequencies were found to be well correlated with each other such that the variations at higher frequencies are leading the lower frequencies, which is expected from the standard shock-in-jet model.