• ### 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.
• ### Reconciling inverse-Compton Doppler factors with variability Doppler factors in blazar jets(1503.04780)

June 25, 2017 astro-ph.HE
Blazar population models have shown that the inverse-Compton and variability Doppler factor estimates yield consistent results at the population level for flat spectrum radio quasars (FSRQs). The two methods, however, are inconsistent when compared on a source-by-source basis. In this work, we attempt to understand the source of the discrepancy by tracing the potential sources of systematic and statistical error for the inverse-Compton Doppler factors. By eliminating these sources of error, we provide stronger constrains on the value of the Doppler factor in blazar jets. We re-estimate the inverse-Compton Doppler factor for 11 sources that meet certain criteria for their synchrotron peak frequency and the availability of Doppler factor estimates in the literature. We compare these estimates with the average of two different estimates of the variability Doppler factor obtained using various datasets and methodologies to identify any discrepancies and, in each case, trace their sources in the methodology or assumptions adopted. We identify three significant sources of error for the inverse-Compton Doppler factors: a) contamination of the X-ray flux by non-synchrotron self-Compton emission; b) radio observations at frequencies other than the synchrotron turnover frequency; c) non-simultaneity between radio and X-ray observations. We discuss key aspects in the correct application of the inverse-Compton method in light of these potential errors. We are able to constrain the Doppler factor of 3C273, 3C345, 3C454.3, PKS1510-089, and PKS1633+382 effectively, since all available estimates from both methods converge to the same values for these five sources.
• ### 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.
• ### Radio and gamma-ray loud narrow-line Seyfert 1 galaxies in the spotlight(1701.09111)

Jan. 31, 2017 astro-ph.HE
Narrow-line Seyfert 1 (NLS1) galaxies provide us with unique insights into the drivers of AGN activity under extreme conditions. Given their low black hole (BH) masses and near-Eddington accretion rates, they represent a class of galaxies with rapidly growing supermassive BHs in the local universe. Here, we present the results from our multi-frequency radio monitoring of a sample of {\gamma}-ray loud NLS1 galaxies ({\gamma}NLS1s), including systems discovered only recently, and featuring both the nearest and the most distant {\gamma}NLS1s known to date. We also present high-resolution radio imaging of 1H 0323+342, which is remarkable for its spiral or ring-like host. Finally, we present new radio data of the candidate {\gamma}-emitting NLS1 galaxy RX J2314.9+2243, characterized by a very steep radio spectrum, unlike other {\gamma}NLS1s.
• ### 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.
• ### Detecting the Elusive Blazar Counter-Jets(1610.06561)

Oct. 20, 2016 astro-ph.HE
Detection of blazar pc scale counter-jets is difficult, but it can provide invaluable insight into the relativistic effects, radiative processes and the complex mechanisms of jet production, collimation and accelation in blazars. We build on recent populations models (optimized using the MOJAVE apparent velocity and redshift distributions) in order to derive the distribution of jet-to-counter-jet ratios and the flux densities of the counter-jet at different frequencies, in an effort to set minimum sensitivity limits required for existing and future telescope arrays in order to detect these elusive counter-jets. We find that: for the BL Lacs $5\%$ of their counter-jets have a flux-density higher than 100mJy, $15\%$ are higher than 10 mJy, and $32\%$ have higher flux-density than 1 mJy, whereas for the FSRQs $8\%$ have a flux-density higher than 10mJy, $17\%$ are higher than 1 mJy, and $32\%$ are higher than 0.1 mJy (at 15 GHz). Future telescopes like the SKA and newly operating like e-MERLIN and JVLA may detect up to $99\%$ of the BL Lac and $77\%$ of the FSRQ counter-jets. Sources with both low apparent velocity and a low Doppler factor make prime candidates for counter-jet detection. Combining our findings with literature values we have identified five such counter-jet detection candidates. Finally, we discuss possible effects beyond relativistic deboosting that may complicate the detection of counter-jets and that need to be accounted for in the interpretation of detections.
• ### RoboPol: The optical polarization of gamma-ray--loud and gamma-ray--quiet blazars(1609.00640)

Sept. 1, 2016 astro-ph.CO, astro-ph.HE
We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.
• ### Optical polarization of high-energy BL Lac objects(1608.08440)

Aug. 30, 2016 astro-ph.CO, astro-ph.HE
We investigate the optical polarization properties of high-energy BL Lac objects using data from the RoboPol blazar monitoring program and the Nordic Optical Telescope. We wish to understand if there are differences in the BL Lac objects that are detected with the current-generation TeV instruments compared to those that have not yet been detected. The mean polarization fraction of the TeV-detected BL Lacs is 5% while the non-TeV sources show a higher mean polarization fraction of 7%. This difference in polarization fraction disappears when the dilution by the unpolarized light of the host galaxy is accounted for. The TeV sources show somewhat lower fractional polarization variability amplitudes than the non-TeV sources. Also the fraction of sources with a smaller spread in the Q/I - U/I -plane and a clumped distribution of points away from the origin, possibly indicating a preferred polarization angle, is larger in the TeV than in the non-TeV sources. These differences between TeV and non-TeV samples seems to arise from differences between intermediate and high spectral peaking sources instead of the TeV detection. When the EVPA variations are studied, the rate of EVPA change is similar in both samples. We detect significant EVPA rotations in both TeV and non-TeV sources, showing that rotations can occur in high spectral peaking BL Lac objects when the monitoring cadence is dense enough. Our simulations show that we cannot exclude a random walk origin for these rotations. These results indicate that there are no intrinsic differences in the polarization properties of the TeV-detected and non-TeV-detected high-energy BL Lac objects. This suggests that the polarization properties are not directly related to the TeV-detection, but instead the TeV loudness is connected to the general flaring activity, redshift, and the synchrotron peak location. (Abridged)
• ### 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.
• ### RoboPol: Do optical polarization rotations occur in all blazars?(1607.04292)

July 14, 2016 astro-ph.HE
We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.
• ### 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.
• ### Optical polarization map of the Polaris Flare with RoboPol(1503.03054)

June 20, 2016 astro-ph.SR
The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $\mu m$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.
• ### 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
• ### RoboPol: optical polarization-plane rotations and flaring activity in blazars(1601.03392)

Jan. 15, 2016 astro-ph.HE
We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.
• 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.
• We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.
• ### 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.
• We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.