• ### Scale invariant jets: from blazars to microquasars(1711.03979)

Nov. 10, 2017 astro-ph.HE
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.
• ### Full-Stokes polarimetry with circularly polarized feeds - Sources with stable linear and circular polarization in the GHz regime(1706.04200)

June 13, 2017 astro-ph.IM
We present a pipeline that allows recovering reliable information for all four Stokes parameters with high accuracy. Its novelty relies on the treatment of the instrumental effects already prior to the computation of the Stokes parameters contrary to conventional methods, such as the M\"uller matrix one. The instrumental linear polarization is corrected across the whole telescope beam and significant Stokes $Q$ and $U$ can be recovered even when the recorded signals are severely corrupted. The accuracy we reach in terms of polarization degree is of the order of 0.1-0.2 %. The polarization angles are determined with an accuracy of almost 1$^{\circ}$. The presented methodology was applied to recover the linear and circular polarization of around 150 Active Galactic Nuclei. The sources were monitored from July 2010 to April 2016 with the Effelsberg 100-m telescope at 4.85 GHz and 8.35 GHz with a cadence of around 1.2 months. The polarized emission of the Moon was used to calibrate the polarization angle. Our analysis showed a small system-induced rotation of about 1$^{\circ}$ at both observing frequencies. Finally, we identify five sources with significant and stable linear polarization; three sources remain constantly linearly unpolarized over the period we examined; a total of 11 sources have stable circular polarization degree $m_\mathrm{c}$ and four of them with non-zero $m_\mathrm{c}$. We also identify eight sources that maintain a stable polarization angle over the examined period. All this is provided to the community for polarization observations reference. We finally show that our analysis method is conceptually different from the traditionally used ones and performs better than the M\"uller matrix method. Although it was developed for a system equipped with circularly polarized feeds it can easily be modified for systems with linearly polarized feeds as well.
• ### Properties of flat-spectrum radio-loud Narrow-Line Seyfert 1 Galaxies(1409.3716)

Nov. 8, 2019 astro-ph.HE
We have conducted a multiwavelength survey of 42 radio loud narrow-1ine Seyfert 1 galaxies (RLNLS1s), selected by searching among all the known sources of this type and omitting those with steep radio spectra. We analyse data from radio frequencies to X-rays, and supplement these with information available from online catalogs and the literature in order to cover the full electromagnetic spectrum. This is the largest known multiwavelength survey for this type of source. We detected 90% of the sources in X-rays and found 17% at gamma rays. Extreme variability at high energies was also found, down to timescales as short as hours. In some sources, dramatic spectral and flux changes suggest interplay between a relativistic jet and the accretion disk. The estimated masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than those of blazars, while the accretion luminosities span a range from $\sim 0.01$ to $\sim 0.49$ times the Eddington limit, similar to those of quasars. The distribution of the calculated jet power spans a range from $\sim 10^{42.6}$ to $\sim 10^{45.6}$ erg s$^{-1}$, generally lower than quasars and BL Lac objects, but partially overlapping with the latter. Once normalised by the mass of the central black holes, the jet power of the three types of active galactic nuclei are consistent with each other, indicating the scalability of the jet. Despite the observational differences, the central engine of RLNLS1s is apparently quite similar to that of blazars. The historical difficulties in finding radio-loud narrow-line Seyfert 1 galaxies might be due to their low power and to intermittent jet activity.
• ### F-GAMMA: Variability Doppler factors of blazars from multiwavelength monitoring(1701.01452)

Jan. 5, 2017 astro-ph.HE
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.
• ### Exploring the magnetic field configuration in BL Lac using GMVA(1609.04031)

Sept. 13, 2016 astro-ph.HE
The high radio frequency polarization imaging of non-thermal emission from active galactic nuclei (AGN) is a direct way to probe the magnetic field strength and structure in the immediate vicinity of supermassive black holes (SMBHs) and is crucial in testing the jet-launching scenario. To explore the the magnetic field configuration at the base of jets in blazars, we took advantage of the full polarization capabilities of the Global Millimeter VLBI Array (GMVA). With an angular resolution of $\sim$50 micro-arcseconds ($\mu$as) at 86 GHz, one could resolve scales up to $\sim$450 gravitational radii (for a 10$^9$ solar mass black hole at a redshift of 0.1). We present here the preliminary results of our study on the blazar BL~Lac. Our results suggest that on sub-mas scales the core and the central jet of BL Lac are significantly polarized with two distinct regions of polarized intensity. We also noted a great morphological similarity between the 7mm/3mm VLBI images at very similar angular resolution.
• ### Inner jet kinematics and the viewing angle towards the {\gamma}-ray narrow-line Seyfert 1 galaxy 1H 0323+342(1608.03232)

Aug. 10, 2016 astro-ph.HE
Near-Eddington accretion rates onto low-mass black holes are thought to be a prime driver of the multi-wavelength properties of the narrow-line Seyfert 1 (NLS1) population of active galactic nuclei (AGN). Orientation effects have repeatedly been considered as another important factor involved, but detailed studies have been hampered by the lack of measured viewing angles towards this type of AGN. Here we present multi-epoch, 15 GHz VLBA images (MOJAVE program) of the radio-loud and Fermi/LAT-detected NLS1 galaxy 1H 323+342. These are combined with single-dish, multi-frequency radio monitoring of the source's variability, obtained with the Effelsberg 100-m and IRAM 30-m telescopes, in the course of the F-GAMMA program. The VLBA images reveal 6 components with apparent peeds of ~1 to ~7 c, and one quasi-stationary feature. Combining the obtained apparent jet speed ($\beta_{app}$) and variability Doppler factor ($D_{var}$) estimates together with other methods, we constrain the viewing angle towards 1H 0323+342 to $\theta \leq 4 - 13$ deg. Using literature values of $\beta_{app}$ and $D_{var}$, we also deduce a viewing angle of $\leq$ 8-9 deg towards another radio- and {\gamma}-ray-loud NLS1, namely SBS 0846+513.
• ### 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(1608.02580)

Aug. 7, 2016 astro-ph.CO, astro-ph.HE
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.
• ### Location of Gamma-ray emission and magnetic field strengths in OJ 287(1607.00725)

July 4, 2016 astro-ph.HE
The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.
• ### Exploring the nature of broadband variability in the FSRQ 3C 273(1603.06989)

March 22, 2016 astro-ph.HE
Detailed investigation of broadband flux variability in the blazar 3C 273 allows us to probe the location and size of emission regions and their physical conditions. We report the results on correlation studies of the flaring activity observed between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function, and power spectral density analysis (PSD) methods. The PSD analysis showed that the optical/IR light curve slopes are consistent with the slope of white noise processes, while, the PSD slopes at radio, X-ray and gamma-ray energies are consistent with red-noise processes. The flux variations at gamma-ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110\pm27) days between gamma-ray and radio light curves, we constrained the location of the gamma-ray emission region at a de-projected distance of 1.2\pm0.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 \gamma-rays energies. The correlation between X-rays and gamma-rays light curves suggests presence of two components responsible for the X-ray emission. A negative time lag of -(50\pm20) days, where the X-rays are leading the emission, suggests X-rays are emitted closer to the jet apex from a compact region at a distance of ~(0.5\pm0.4) pc from the jet apex. A positive time lag of (110\pm20) days suggests jet-base origin of the other X-ray component at ~(4--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 could be expected in the standard shock-in-jet model.
• ### RoboPol: First season rotations of optical polarization plane in blazars(1505.07467)

March 22, 2016 astro-ph.HE
We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.
• ### What can the 2008/10 broadband flare of PKS 1502+106 tell us? Nuclear opacity, magnetic fields, and the location of gamma rays(1603.04220)

March 14, 2016 astro-ph.HE
Context. The origin of blazar variability, seen from radio up to gamma rays, is still a heavily debated matter and broadband flares offer a unique testbed towards a better understanding of these extreme objects. Such an energetic outburst was detected by Fermi/LAT in 2008 from the blazar PKS 1502+106. The outburst was observed from gamma rays down to radio frequencies. Aims. Through the delay between flare maxima at different radio frequencies, we study the frequency-dependent position of the unit-opacity surface and infer its absolute position with respect to the jet base. This nuclear opacity profile enables the magnetic field tomography of the jet. We also localize the gamma-ray emission region and explore the mechanism producing the flare. Methods. The radio flare of PKS 1502+106 is studied through single-dish flux density measurements at 12 frequencies in the range 2.64 to 226.5 GHz. To quantify it, we employ both a Gaussian process regression and a discrete cross-correlation function analysis. Results. We find that the light curve parameters (flare amplitude and cross-band delays) show a power-law dependence on frequency. Delays decrease with frequency, and the flare amplitudes increase up to about 43 GHz and then decay. This behavior is consistent with a shock propagating downstream the jet. The self-absorbed radio cores are located between about 10 and 4 pc from the jet base and their magnetic field strengths range between 14 and 176 mG, at the frequencies 2.64 to 86.24 GHz. Finally, the gamma-ray active region is located at (1.9 +/- 1.1) pc away from the jet base.

Feb. 26, 2016 astro-ph.GA, astro-ph.HE
• 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.
• ### Multiwavelength survey of a sample of flat-spectrum radio-loud narrow-line Seyfert 1 galaxies(1512.00192)

Dec. 1, 2015 astro-ph.GA, astro-ph.HE
We report on a multiwavelength survey of a sample of 42 flat-spectrum radio-loud narrow-line Seyfert 1 galaxies (RLNLS1s). This is the largest known sample of this type of active galactic nucleus (AGN) to date. We found that 17% of sources were detected at high-energy gamma rays (E>100 MeV), and 90% at X-rays (0.3-10 keV). The masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than the values of blazars. The disk luminosities are about 1-49% of the Eddington value, with one outlier at 0.3%, comparable with the luminosities observed in flat-spectrum radio quasars (FSRQs). The jet powers are $\sim 10^{42-46}$ erg s$^{-1}$, comparable with BL Lac Objects, yet relatively smaller than FSRQs. However, once renormalized by the mass of the central black hole, the jet powers of RLNLS1s, BL Lacs, and FSRQs are consistent each other, indicating the scalability of the jets. We found episodes of extreme variability at high energies on time scales of hours. In some cases, dramatic spectral and flux changes are interpreted as the interplay between the relativistic jet and the accretion disk. We conclude that, despite the distinct observational properties, the central engines of RLNLS1s are similar to those of blazars.
• ### PKS 1502+106: A high-redshift Fermi blazar at extreme angular resolution. Structural dynamics with VLBI imaging up to 86 GHz(1511.01085)

Nov. 3, 2015 astro-ph.CO, astro-ph.HE
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.
• ### Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113(1509.02063)

Oct. 12, 2015 astro-ph.HE
We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.
• ### The Third Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope(1501.06054)

Aug. 26, 2015 astro-ph.HE
The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

July 3, 2015 astro-ph.HE
Aims: Our aim is to explore the broad-band radio continuum spectrum of LSI+61303 during its outbursts by employing the available set of secondary focus receivers of the Effelsberg 100 m telescope. Methods: The clear periodicity of the system LSI+61303 allowed observations to be scheduled covering the large radio outburst in March-April 2012. We observed LSI+61303 on 14 consecutive days at 2.6, 4.85, 8.35, 10.45, 14.3, 23, and 32 GHz with a cadence of about 12 hours followed by two additional observations several days later. Based on these observations we obtained a total of 24 quasi-simultaneous broad-band radio spectra. Results: During onset, the main flare shows an almost flat broad-band spectrum, most prominently seen on March 27, 2012, where - for the first time - a flat spectrum (alpha=0.00+/-0.07, S nu^alpha) is observed up to 32 GHz (9 mm wavelength). The flare decay phase shows superimposed 'sub-flares' with the spectral index oscillating between -0.4 and -0.1 in a quasi-regular fashion. Finally, the spectral index steepens during the decay phase, showing optically thin emission with values alpha $\sim$ -0.5 to -0.7. Conclusions: The radio characteristics of LSI+61303 compare well with those of the microquasars XTE J1752-223 and Cygnus X-3. In these systems the flaring phase is actually also composed of a sequence of outbursts with clearly different spectral characteristics: a first outburst with a flat/inverted spectrum followed by a bursting phase of optically thin emission.
• ### The Gamma-Ray Detected Narrow-Line Seyfert 1 Galaxy 1H 0323+342: Swift Monitoring and Suzaku Spectroscopy(1505.01269)

May 6, 2015 astro-ph.HE
As a radio-loud narrow-line Seyfert 1 galaxy (NLS1) detected by Fermi/LAT in GeV $\gamma$-rays, 1H 0323+342 is a remarkable Active Galactic Nucleus (AGN) showing properties characteristic of both NLS1s and blazars. Here we present results of simultaneous X-ray and UV/optical monitoring observations on 1H 0323+342 taken with the UV/Optical Telescope (UVOT) and X-ray Telescope (XRT) onboard the Swift satellite over six years from 2006. Overall, the object showed statistically correlated variations in both the UV and X-ray bands on timescales of years as well as on timescales of days. A deep Suzaku observation reveals X-ray variability on timescales as short as a few tens of thousand seconds, and an X-ray spectrum typical of Seyfert galaxies. The broad-band spectral energy distribution, for which the data of UV and X-ray observations taken on 2009 July 26-27 were used, can be well modeled with a simple one-zone leptonic jet model plus accretion disk/corona emission. The latter is predominantly responsible for the UV/optical and X-ray (0.3-10 keV) emission and their observed variations. The correlated UV-X-ray variability on the timescale of days is consistent with reprocessing of the X-ray radiation by the accretion disk. The shortest timescale and large normalized excess variance of the X-ray variability detected with Suzaku suggest a relatively small black hole mass of the order of $10^7M_{\odot}$, consistent with the estimation based on the broad H$\beta$ line in the optical.
• ### Intra-day variability observations and the VLBI structure analysis of quasar S4 0917+624(1504.04141)

April 28, 2015 astro-ph.GA, astro-ph.HE
The IDV observations of S4 0917+624 were carried out monthly, from August 2005 to January 2010, with the Urumqi 25m radio telescope at 4.8 GHz. The quasar S4 0917+624 exhibits only very weak or no IDV during our 4.5 year observing interval. Prior to the year 2000, the source S4 0917+624 was one of the most prominent IDV sources. Our new data indicate that the previous strong IDV has ceased. We analyzed the long-term VLBI structural variability using Gaussian model-fitting. From this we obtained the flux densities and the deconvolved sizes of core and inner-jet components of the source. We studied the properties such as core fraction, angular size, spectral index, and brightness temperature of VLBI core for S4 0917+624, as well as the time delay between 5 and 15 GHz variations, and compared them with the IDV properties of S4 0917+624. The source shows ejection of several jet components that are suspected to have partially reduced the IDV amplitude of S4 0917+624. However, during 2005-2006, the VLBI core size was comparable to the size before the year 2000, but no strong IDV was detected in the period, suggesting that the quenching effect due to source size changes may not be responsible for the lack of strong IDV after the year 2000. The refractive scattering properties for the strong IDV phase of S4 0917+624 before the year 2000 are discussed. The disappearance of strong IDV in S4 0917+624 after the year 2000 is a mystery and cannot be explained via the quenching effect by changes in the observable VLBI structure. However, it may be caused by changes in the interstellar medium, i.e. by interstellar weather, which induces changes in the scintillation pattern on timescales of several years. Further coordinated multi-frequency observations will be required to distinguish between the effect of source-intrinsic variability and changing properties of the interstellar medium.
• ### Localizing the $\gamma$ rays from blazar PKS 1502+106(1504.02314)

April 9, 2015 astro-ph.GA, astro-ph.HE
Blazars are among the most variable objects in the universe. They feature energetic jets of plasma that launch from the cores of these active galactic nuclei (AGN), triggering activity from radio up to gamma-ray energies. Spatial localization of the region of their MeV/GeV emission is a key question in understanding the blazar phenomenon. The flat spectrum radio quasar (FSRQ) PKS 1502+106 has exhibited extreme and correlated, radio and high-energy activity that triggered intense monitoring by the Fermi-GST AGN Multi-frequency Monitoring Alliance (F-GAMMA) program and the Global Millimeter VLBI Array (GMVA) down to $\lambda$3 mm (or 86 GHz), enabling the sharpest view to date towards this extreme object. Here, we report on preliminary results of our study of the gamma-ray loud blazar PKS 1502+106, combining VLBI and single dish data. We deduce the critical aspect angle towards the source to be $\theta_{\rm c} = 2.6^{\circ}$, calculate the apparent and intrinsic opening angles and constrain the distance of the 86 GHz core from the base of the conical jet, directly from mm-VLBI but also through a single dish relative timing analysis. Finally, we conclude that gamma rays from PKS 1502+106 originate from a region between ~1-16 pc away from the base of the hypothesized conical jet, well beyond the bulk of broad-line region (BLR) material of the source.
• ### Parsec scale polarization properties of the TeV blazar Markarian 421(1504.02014)

April 8, 2015 astro-ph.HE
In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.
• ### Radio jet emission from GeV-emitting narrow-line Seyfert 1 galaxies(1501.02158)

Jan. 6, 2015 astro-ph.GA, astro-ph.HE
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.
• ### Multi-frequency properties of synthetic blazar radio light curves within the shock-in-jet scenario(1412.7194)

Dec. 22, 2014 astro-ph.GA, astro-ph.HE
Blazars are among the most powerful extragalactic objects, as a sub-class of active galactic nuclei. They launch relativistic jets and their emitted radiation shows strong variability across the entire electro-magnetic spectrum. The mechanisms producing the variability are still controversial and different models have been proposed to explain the observed variations in multi-frequency blazar light curves.We investigate the capabilities of the classical shock-in-jet model to explain and reconstruct the observed evolution of flares in the turnover frequency turnover flux density plane and their frequency-dependent light curve parameters. With a detailed parameter space study we provide the framework for future, detailed comparisons of observed flare signatures with the shock-in-jet scenario. Based on the shock model we compute synthetic single-dish light curves at different radio frequencies (2.6 to 345 GHz) and for different physical conditions in a conical jet (e.g. magnetic field geometry and Doppler factor). From those we extract the slopes of the different energy loss stages within the $\nu_\mathrm{m}$-$S_\mathrm{m}$ plane and deduce the frequency-dependence of different light curve parameters such as flare amplitude, time scale and cross-band delays. The evolution of the Doppler factor along the jet has the largest influence on the evolution of the flare and on the frequency-dependent light curve parameters. The synchrotron stage can be hidden in the Compton or in the adiabatic stage, depending mainly on the evolution of the Doppler factor, which makes it difficult to detect its signature in observations. In addition, we show that the time lags between different frequencies can be used as an efficient tool to better constrain the physical properties of these objects.
• ### Very Long Baseline polarimetry and the gamma-ray connection in Markarian 421 during the broadband campaign in 2011(1410.0884)

Oct. 3, 2014 astro-ph.HE
We analyzed data in polarized intensity obtained with the Very Long Baseline Array (VLBA) at twelve epochs (one observation per month from January to December 2011) at 15, 24, and 43 GHz. For the absolute orientation of the electric vector position angles (EVPA) we used the D-terms method. We also used gamma-ray data from the Fermi Large Area Telescope on weekly time bins throughout 2011. The source shows polarized emission, and its properties vary with time, frequency, and location along the jet. The core mean polarization fraction is generally between 1% and 2%, with a 4% peak at 43 GHz in March; the polarization angle is variable, mainly at 15 GHz, where it changes frequently, and less so at 43 GHz, where it oscillates in the range 114 - 173 deg. The jet polarization properties are more stable, with a fractional polarization of around 16% and a polarization angle nearly perpendicular to the jet axis. The average flux and photon index at gamma-ray energies are (17.7+-0.5)x10^8 ph cm^-2 s^-1 and the photon index is 1.77+-0.02. The gamma-ray light curve shows variability, with a main peak that appears to be associated with the peak in the core polarized emission at 43 GHz, as well as with the total intensity light curve. A discrete correlation function analysis yields a correlation coefficient of 0.54 at zero delay, with a significance level above 99.7%. We accurately determine the polarization properties of Mrk 421, both in the core and in the jet region. The radio and gamma-ray light curves are correlated. The observed EVPA variability at 15 GHz is partly due to opacity and partly to a variable Faraday rotation effect. To explain the residual variability of the intrinsic polarization angle and the low degree of polarization in the core region, we invoke a blend of variable cross-polarized subcomponents with different polarization properties within the beam.