Gaussian processes for blazar variability studies

This article shortly introduces Gaussian processes (GP) as a new approach for modelling time series in the field of blazar physics. In the second part of the paper, recent results from an application of GP modelling to the multi-wavelength light curves of the blazar PKS 1502+106 are discussed.Comment: 7 pages, 2 figures. Accepted for publication in Galaxies. This article belongs to the Special Issue "Blazars through Sharp Multi-wavelength Eyes

Scale invariant jets: from blazars to microquasars

Black holes, anywhere in the stellar-mass to supermassive range, are often associated with relativistic jets. Models suggest that jet production may be a universal process common in all black hole systems regardless of their mass. Although in many cases observations support such hypotheses for microquasars and Seyfert galaxies, little is known on whether boosted blazar jets also comply with such universal scaling laws. We use uniquely rich multiwavelength radio light curves from the F-GAMMA program and the most accurate Doppler factors available to date to probe blazar jets in their emission rest frame with unprecedented accuracy. We identify for the first time a strong correlation between the blazar intrinsic broad-band radio luminosity and black hole mass, which extends over $\sim$ 9 orders of magnitude down to microquasars scales. Our results reveal the presence of a universal scaling law that bridges the observing and emission rest frames in beamed sources and allows us to effectively constrain jet models. They consequently provide an independent method for estimating the Doppler factor, and for predicting expected radio luminosities of boosted jets operating in systems of intermediate or tens-of-solar mass black holes, immediately applicable to cases as those recently observed by LIGO.Comment: 13 pages, 4 figures, accepted for publication in AP

F-GAMMA: Multi-frequency radio monitoring of Fermi blazars. The 2.64 to 43 GHz Effelsberg light curves from 2007-2015

The advent of the Fermi-GST with its unprecedented capability to monitor the entire 4 pi sky within less than 2-3 hours, introduced new standard in time domain gamma-ray astronomy. To explore this new avenue of extragalactic physics the F-GAMMA programme undertook the task of conducting nearly monthly, broadband radio monitoring of selected blazars from January 2007 to January 2015. In this work we release all the light curves at 2.64, 4.85, 8.35, 10.45, 14.6, 23.05, 32, and 43 GHz and present first order derivative data products after all necessary post-measurement corrections and quality checks; that is flux density moments and spectral indices. The release includes 155 sources. The effective cadence after the quality flagging is around one radio SED every 1.3 months. The coherence of each radio SED is around 40 minutes. The released dataset includes more than $4\times10^4$ measurements. The median fractional error at the lowest frequencies (2.64-10.45 GHz) is below 2%. At the highest frequencies (14.6-43 GHz) with limiting factor of the atmospheric conditions, the errors range from 3% to 9%, respectively.Comment: Accepted for publication in Section: Catalogs and data of Astronomy & Astrophysic

The TeV-emitting radio galaxy 3C 264. VLBI kinematics and SED modeling

Context. In March 2018, the detection by VERITAS of very-high-energy emission (VHE > 100 GeV) from 3C 264 was reported. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. Aims: In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this γ-ray emitting jet. Methods: The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and covers a time period of about two years. Through a segmented wavelet decomposition method (WISE code), we estimated the apparent displacement of individual plasma features; we then performed a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the spectral energy distribution (SED). Results: Proper motion is mostly detected along the edges of the flow, which appears strongly limb brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ̃11.5 c. This constrains the jet viewing angle to assume relatively small values (θ ≲ 10°). In the acceleration region, extending up to a de-projected distance of ̃4.8 × 104 Schwarzschild radii (̃11 pc), the jet is collimating (r ∝ z0.40 ± 0.04), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated (UB/Ue > 1), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling

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.Comment: Accepted for publication in 4 - Extragalactic astronomy of Astronomy and Astrophysic

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 $\gamma$-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.Comment: 9 pages, 7 figures, accepted for publication in MNRA

High cadence, linear and circular polarization monitoring of OJ 287 - Helical magnetic field in a bent jet

We present a multi-frequency, dense radio monitoring program of the blazar OJ287 using the 100m Effelsberg radio telescope. We analyze the evolution in total flux density, linear and circular polarization to study the jet structure and its magnetic field geometry. The total flux density is measured at nine bands from 2.64 GHz to 43 GHz, the linear polarization parameters between 2.64 GHz and 10.45 GHz, and the circular polarization at 4.85 GHz and 8.35 GHz. The mean cadence is 10 days. Between MJD 57370 and 57785, OJ287 showed flaring activity and complex linear and circular polarization behavior. The radio EVPA showed a large clockwise (CW) rotation by ~340$^{\circ}$ with a mean rate of -1.04 $^{\circ}$/day. Based on concurrent VLBI data, the rotation seems to originate within the jet core at 43 GHz (projected size $\le$ 0.15 mas or 0.67 pc). Moreover, optical data show a similar monotonic CW EVPA rotation with a rate of about -1.1 $^{\circ}$/day which is superposed with shorter and faster rotations of about 7.8 $^{\circ}$/day. The observed variability is consistent with a polarized emission component propagating on a helical trajectory within a bent jet. We constrained the helix arc length to 0.26 pc and radius to $\le$ 0.04 pc as well as the jet bending arc length projected on the plane of the sky to $\le$ 1.9-7.6 pc. A similar bending is observed in high angular resolution VLBI images at the innermost jet regions. Our results indicate also the presence of a stable polarized emission component with EVPA (-10$^{\circ}$) perpendicular to the large scale jet, suggesting dominance of the poloidal magnetic field component. Finally, the EVPA rotation begins simultaneously with an optical flare and hence the two might be physically connected. That optical flare has been linked to the interaction of a secondary SMBH with the inner accretion disk or originating in the jet of the primary.Comment: 11 pages, 9 figures, 1 table, accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics on August 21, 201

The TeV-emitting radio galaxy 3C 264. VLBI kinematics and SED modeling

Context. In March 2018, the detection by VERITAS of very-high-energy emission (VHE > 100 GeV) from 3C 264 was reported. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. Aims: In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this γ-ray emitting jet. Methods: The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and covers a time period of about two years. Through a segmented wavelet decomposition method (WISE code), we estimated the apparent displacement of individual plasma features; we then performed a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the spectral energy distribution (SED). Results: Proper motion is mostly detected along the edges of the flow, which appears strongly limb brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ̃11.5 c. This constrains the jet viewing angle to assume relatively small values (θ ≲ 10°). In the acceleration region, extending up to a de-projected distance of ̃4.8 × 104 Schwarzschild radii (̃11 pc), the jet is collimating (r ∝ z0.40 ± 0.04), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated (UB/Ue > 1), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling

The F-GAMMA program: Multi-frequency study of Active Galactic Nuclei in the Fermi era. Program description and the first 2.5 years of monitoring

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly 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-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. 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 FSRQs, 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 as well as 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 3sigma, implying that gamma rays are produced very close to the mm-band emission region.Comment: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics (18 pages, 9 figures