Search

CTA 102 is a $\gamma$-ray bright blazar that exhibited multiple flares in observations by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope during the period of 2016-2018. We present results from the analysis of multi-wavelength light curves aiming at revealing the nature of $\gamma$-ray flares from the relativistic jet in the blazar. We analyse radio, optical, X-ray, and $\gamma$-ray data obtained in a period from 2012 September 29 to 2018 October 8. We identify six flares in the $\gamma$-ray light curve, showing a harder-when-brighter-trend in the $\gamma$-ray spectra. We perform a cross-correlation analysis of the multi-wavelength light curves. We find nearly zero time lags between the $\gamma$-ray and optical and X-ray light curves, implying a common spatial origin for the emission in these bands. We find significant correlations between the $\gamma$-ray and radio light curves as well as negative/positive time lags with the $\gamma$-ray emission lagging/leading the radio during different flaring periods. The time lags between $\gamma$-ray and radio emission propose the presence of multiple $\gamma$-ray emission sites in the source. As seen in 43 GHz images from the Very Long Baseline Array, two moving disturbances (or shocks) were newly ejected from the radio core. The $\gamma$-ray flares from 2016 to 2017 are temporally coincident with the interaction between a traveling shock and a quasi-stationary one at $\sim$0.1 mas from the core. The other shock is found to emerge from the core nearly simultaneous with the $\gamma$-ray flare in 2018. Our results suggest that the $\gamma$-ray flares originated from shock-shock interactions., Comment: 14 pages, 15 figures, 4 tables, accepted for publication in A&A

We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3\sigma$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions., Comment: Accepted for publication in ApJ

Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

Context. The blazar AO 0235+164, located at redshift $z=0.94$, has displayed interesting and repeating flaring activity in the past, the latest episodes occurring in 2008 and 2015. In 2020, the source brightened again, starting a new flaring episode that peaked in 2021. Aims. We study the origin and properties of the 2021 flare in relation to previous studies and the historical behavior of the source, in particular to the 2008 and 2015 flaring episodes. Methods. We analyze the multi-wavelength photo-polarimetric evolution of the source. From Very Long Baseline Array images, we derive the kinematic parameters of new components associated with the 2021 flare. We use this information to constrain a model for the spectral energy distribution of the emission during the flaring period. We propose an analytical geometric model to test whether the observed wobbling of the jet is consistent with precession. Results. We report the appearance of two new components that are ejected in a different direction than previously, confirming the wobbling of the jet. We find that the direction of ejection is consistent with that of a precessing jet.The derived period independently agrees with the values commonly found in the literature. Modeling of the spectral energy distribution further confirm that the differences between flares can be attributed to geometrical effects., Comment: Accepted 15 May 2024

Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $\Pi_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $\psi_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $\Pi_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $\Pi_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day., Comment: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

One of the most well-known extragalactic sources in the sky, quasar 3C 454.3, shows a curved parsec-scale jet that has been exhaustively monitored with very-long-baseline interferometry (VLBI) over the recent years. In this work, we present a comprehensive analysis of four years of high-frequency VLBI observations at 43 GHz and 86 GHz, between 2013-2017, in total intensity and linear polarization. The images obtained from these observations enabled us to study the jet structure and the magnetic field topology of the source on spatial scales down to 4.6 parsec in projected distance. The kinematic analysis reveals the abrupt vanishing of at least four new superluminal jet features in a characteristic jet region (i.e., region C), which is located at an approximate distance of 0.6 milliarcseconds from the VLBI core. Our results support a model in which the jet bends, directing the relativistic plasma flow almost perfectly toward our line of sight, co-spatially with the region where components appear to stop., Comment: 15 pages, 7 figures

Context. Blazar AO 0235+164, located at redshift z = 0.94, has undergone several sharp multi-spectral-range flaring episodes during the last decades. In particular, the episodes peaking in 2008 and 2015, that received extensive multi-wavelength coverage, exhibited interesting behavior. Aims. We study the actual origin of these two observed flares by constraining the properties of the observed photo-polarimetric variability, those of the broad-band spectral energy-distribution and the observed time-evolution behavior of the source as seen by ultra-high resolution total-flux and polarimetric Very-long-baseline interferometry (VLBI) imaging. Methods. The analysis of VLBI images allows us to constrain kinematic and geometrical parameters of the 7 mm jet. We use the Discrete Correlation Function to compute the statistical correlation and the delays between emission at different spectral ranges. Multi-epoch modeling of the spectral energy distributions allows us to propose specific models of emission; in particular for the unusual spectral features observed in this source in the X-ray region of the spectrum during strong multi spectral-range flares. Results. We find that these X-ray spectral features can be explained by an emission component originating in a separate particle distribution than the one responsible for the two standard blazar bumps. This is in agreement with the results of our correlation analysis that do not find a strong correlation between the X-rays and the remaining spectral ranges. We find that both external Compton dominated and synchrotron self-Compton dominated models can explain the observed spectral energy distributions. However, synchrotron self-Compton models are strongly favored by the delays and geometrical parameters inferred from the observations.

We present the multiwavelength flaring activity of the blazar AO 0235+164 during its recent active period from 2013 to 2019. From a discrete correlation function (DCF) analysis, we find a significant (>95%) correlation between radio and $\gamma$-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz VLBA data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of $\delta_{\rm var}=28.5\pm8.4$, a bulk Lorentz factor of $\Gamma=16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta_{\rm v}=1.42^{+1.07}_{-0.52}\textrm{ degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of $10-70$ GHz and the peak flux density varying in the range of $0.7-4$ Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}$ mG at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}$ mG calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235+164 suggests that it did not significantly deviate from equipartition during its recent flaring activity., Comment: 13 pages, 11 figures, 4 tables; Accepted for publication in MNRAS

We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$\sigma$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polarizations favor models where the X-ray band is dominated by unpolarized photons upscattered by relativistic electrons in the jets of blazars, although hadronic models are not completely eliminated. We discuss the X-ray polarization upper limits in the context of our contemporaneous multiwavelength polarization campaigns., Comment: 19 pages, 10 figures, 6 tables; submitted to the Astrophysical Journal

The lower energy peak of the spectral energy distribution of blazars has commonly been ascribed to synchrotron radiation from relativistic particles in the jets. Despite the consensus regarding jet emission processes, the particle acceleration mechanism is still debated. Here, we present the first X-ray polarization observations of PG 1553+113, a high-synchrotron-peak blazar observed by the Imaging X-ray Polarimetry Explorer (IXPE). We detect an X-ray polarization degree of $(10\pm2)\%$ along an electric-vector position angle of $\psi_X=86^{\circ}\pm8^{\circ}$. At the same time, the radio and optical polarization degrees are lower by a factor of $\sim$3. During our IXPE pointing, we observed the first orphan optical polarization swing of the IXPE era, as the optical angle of PG 1553+113 underwent a smooth monotonic rotation by about 125$^\circ$, with a rate of $\sim$17 degrees per day. We do not find evidence of a similar rotation in either radio or X-rays, which suggests that the X-ray and optically emitting regions are separate or, at most, partially co-spatial. Our spectro-polarimetric results provide further evidence that the steady-state X-ray emission in blazars originates in a shock-accelerated and energy-stratified electron population., Comment: Accepted for pubblication in The Astrophysical Journal Letters. 3 Figures, 13 pages

We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$\gamma$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\it NuSTAR} in hard X-rays (3--70~keV), and optical polarization from the Calar Alto, and Perkins Telescope observatories. Our contemporaneous X-ray data suggest that the {\it IXPE} energy band is at the crossover between the low- and high-frequency blazar emission humps. The source displays significant variability during the observation, and we measure polarization in three separate time bins. Contemporaneous X-ray spectra allow us to determine the relative contribution from each emission hump. We find $>99\%$ confidence X-ray polarization $\Pi_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $\psi_{2-4{\rm keV}} = -28.7 \pm 8.7^{\circ}$ in the time bin with highest estimated synchrotron flux contribution. We discuss possible implications of our observations, including previous {\it IXPE} BL Lacertae pointings, tentatively concluding that synchrotron self-Compton emission dominates over hadronic emission processes during the observed epochs., Comment: Accepted to ApJL. 16 pages, 8 figures

The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $\Psi$. Here we report the discovery of a $\Psi_{\mathrm x}$ rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, $\Psi_{\mathrm x}$ rotated in total by $\geq360^\circ$. Over the two respective date ranges, we find constant, within uncertainties, rotation rates ($80 \pm 9$ and $91 \pm 8 ^\circ/\rm day$) and polarization degrees ($\Pi_{\mathrm x}=10\%\pm1\%$). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of $\Psi$ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released.

The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow., Comment: 50 pages, 18 figures, 3 tables, author's version of the paper published in Nature

Blazars are a class of jet-dominated active galactic nuclei with a typical double-humped spectral energy distribution. It is of common consensus the Synchrotron emission to be responsible for the low frequency peak, while the origin of the high frequency hump is still debated. The analysis of X-rays and their polarization can provide a valuable tool to understand the physical mechanisms responsible for the origin of high-energy emission of blazars. We report the first observations of BL Lacertae performed with the Imaging X-ray Polarimetry Explorer ({IXPE}), from which an upper limit to the polarization degree $\Pi_X<$12.6\% was found in the 2-8 keV band. We contemporaneously measured the polarization in radio, infrared, and optical wavelengths. Our multiwavelength polarization analysis disfavors a significant contribution of proton synchrotron radiation to the X-ray emission at these epochs. Instead, it supports a leptonic origin for the X-ray emission in BL Lac., Comment: 17 pages, 5 figures, accepted for publication in ApJL

The nature of dark matter (DM) is still under intense debate. Sub-galactic scales are particularly critical, as different, currently viable DM models make diverse predictions on the expected abundance and density profile of DM haloes on these scales. We investigate the ability of sub-galactic DM haloes to act as strong lenses on background compact sources, producing gravitational lensing events on milli-arcsecond scales (milli-lenses), for different DM models. For each DM scenario, we explore whether a sample of $\sim$ 5000 distant sources is sufficient to detect at least one milli-lens. We develop a semi-analytical model to estimate the milli-lensing optical depth as a function of the source's redshift for various DM models. We employ the Press-Schechter formalism, as well as results from recent N-body simulations to compute the halo mass function, taking into account the appropriate spherically averaged density profile of haloes for each DM model. We treat the lensing system as a point-mass lens and invoke the effective surface mass density threshold to calculate the fraction of a halo that acts as a gravitational lens. We study three classes of dark matter models: cold DM, warm DM, and self-interacting DM. We find that haloes consisting of warm DM turn out to be optically thin for strong gravitational milli-lensing (zero expected lensing events). CDM haloes may produce lensing events depending on the steepness of the concentration-mass relation. Self-interacting DM haloes can efficiently act as gravitational milli-lenses only if haloes experience gravothermal collapse, resulting in highly dense central cores., Comment: 11 pages, 6 figures. Submitted for publication in A&A

Most of the light from blazars, active galactic nuclei with jets of magnetized plasma that point nearly along the line of sight, is produced by high-energy particles, up to $\sim 1$ TeV. Although the jets are known to be ultimately powered by a supermassive black hole, how the particles are accelerated to such high energies has been an unanswered question. The process must be related to the magnetic field, which can be probed by observations of the polarization of light from the jets. Measurements of the radio to optical polarization - the only range available until now - probe extended regions of the jet containing particles that left the acceleration site days to years earlier (Jorstad et al., 2005; Marin et al., 2018; Blinov et al., 2021), and hence do not directly explore the acceleration mechanism, as could X-ray measurements. Here we report the detection of X-ray polarization from the blazar Markarian~501 (Mrk~501). We measure an X-ray linear polarization degree $\Pi_X \sim10\%$, a factor of $\sim2$ higher than the value at optical wavelengths, with a polarization angle parallel to the radio jet. This points to a shock front as the source of particle acceleration, and also implies that the plasma becomes increasingly turbulent with distance from the shock., Comment: 34 pages, 4 figures, 4 tables

The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of $35\times 35$ arcminutes and will measure the Stokes parameters $I$, $q$, and $u$ in a single exposure in the SDSS-r broadband filter. The design goal for the instrument is to achieve an overall polarimetric measurement accuracy of 0.1 % over the entire field of view. We present here the complete polarimetric modeling of the instrument, characterizing the amount and sources of instrumental polarization. To accurately retrieve the real Stokes parameters of a source from the measured values, we have developed a calibration method for the instrument. Using this calibration method and simulated data, we demonstrate how to correct instrumental polarization and obtain 0.1 % accuracy in the degree of polarization, $p$. Additionally, we tested and validated the calibration method by implementing it on a table-top WALOP-like test-bed polarimeter in the laboratory., Comment: 34 pages, 24 figures. Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems

We present multi-epoch optical spectra of the $\gamma$-ray bright blazar 1156+295 (4C +29.45, Ton 599) obtained with the 4.3~m Lowell Discovery Telescope. During a multi-wavelength outburst in late 2017, when the $\gamma$-ray flux increased to $2.5\times 10^{-6} \; \rm phot\; cm^{-2}\; s^{-1}$ and the quasar was first detected at energies $\geq100$ GeV, the flux of the Mg II $\lambda 2798$ emission line changed, as did that of the Fe emission complex at shorter wavelengths. These emission line fluxes increased along with the highly polarized optical continuum flux, which is presumably synchrotron radiation from the relativistic jet, with a relative time delay of $\lesssim2$ weeks. This implies that the line-emitting clouds lie near the jet, which points almost directly toward the line of sight. The emission-line radiation from such clouds, which are located outside the canonical accretion-disk related broad-line region, may be a primary source of seed photons that are up-scattered to $\gamma$-ray energies by relativistic electrons in the jet., Comment: 25 pages, 12 figures

We present the first polarimetric space VLBI observations of OJ 287, observed with RadioAstron at 22 GHz during a perigee session on 2014 April 4 and five near-in-time snapshots, together with contemporaneous ground VLBI observations at 15, 43, and 86 GHz. Ground-space fringes were obtained up to a projected baseline of 3.9 Earth diameters during the perigee session, and at a record 15.1 Earth diameters during the snapshot sessions, allowing us to image the innermost jet at an angular resolution of $\sim50\mu$as, the highest ever achieved at 22 GHz for OJ 287. Comparison with ground-based VLBI observations reveals a progressive jet bending with increasing angular resolution that agrees with predictions from a supermassive binary black hole model, although other models cannot be ruled out. Spectral analyses suggest that the VLBI core is dominated by the internal energy of the emitting particles during the onset of a multi-wavelength flare, while the parsec-scale jet is consistent with being in equipartition between the particles and magnetic field. Estimated minimum brightness temperatures from the visibility amplitudes show a continued rising trend with projected baseline length up to $10^{13}$ K, reconciled with the inverse Compton limit through Doppler boosting for a jet closely oriented to the line of sight. The observed electric vector position angle suggests that the innermost jet has a predominantly toroidal magnetic field, which together with marginal evidence of a gradient in rotation measure across the jet width indicate that the VLBI core is threaded by a helical magnetic field, in agreement with jet formation models.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017)., Comment: Published in Astronomy & Astrophysics

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ~287 as revealed by multi-epoch Very Long Baseline Array (VLBA) images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle (EVPA) of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave "core" or, at other epochs, coincident with a bright, quasi-stationary emission feature $\sim0.2$~milliarcsec ($\sim$0.9~pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and $\gamma$-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies $\geq10$~pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We present total and linearly polarized 3 mm Global mm-VLBI Array images of a sample of blazars and radio galaxies from the VLBA-BU-BLAZAR 7 mm monitoring program designed to probe the innermost regions of active galactic nuclei (AGN) jets and locate the sites of gamma-ray emission observed by the Fermi-LAT. The lower opacity at 3 mm and improved angular resolution, on the order of 50 microarcseconds, allow us to distinguish features in the jet not visible in the 7 mm VLBA data. We also compare two different methods used for the calibration of instrumental polarisation and we analyze the resulting images for some of the sources in the sample., Comment: Polarised Emission from Astrophysical Jets, June 12-16, 2017, Ierapetra, Greece

We analyze the parsec-scale jet kinematics from 2007 June to 2013 January of a sample of $\gamma$-ray bright blazars monitored roughly monthly with the Very Long Baseline Array at 43~GHz. In a total of 1929 images, we measure apparent speeds of 252 emission knots in 21 quasars, 12 BL~Lacertae objects (BLLacs), and 3 radio galaxies, ranging from 0.02$c$ to 78$c$; 21\% of the knots are quasi-stationary. Approximately 1/3 of the moving knots execute non-ballistic motions, with the quasars exhibiting acceleration along the jet within 5~pc (projected) of the core, and knots in the BLLacs tending to decelerate near the core. Using apparent speeds of components and timescales of variability from their light curves, we derive physical parameters of 120 superluminal knots, including variability Doppler factors, Lorentz factors, and viewing angles. We estimate the half-opening angle of each jet based on the projected opening angle and scatter of intrinsic viewing angles of knots. We determine characteristic values of physical parameters for each jet and AGN class based on the range of values obtained for individual features. We calculate intrinsic brightness temperatures of the cores, $T_{\rm b,int}^{\rm core}$, at all epochs, finding that the radio galaxies usually maintain equipartition conditions in the cores, while $\sim$30\% of $T_{\rm b,int}^{\rm core}$ measurements in the quasars and BLLacs deviate from equipartition values by a factor $>$10. This probably occurs during transient events connected with active states. In the Appendix we briefly describe the behavior of each blazar during the period analyzed., Comment: 48 pages, 18 figures, Tables 2,3,5, & 7 in a full length and figures 2 and 4 for each source can be found in the journal version of the paper

We analyse the circular polarisation data accumulated in the first 7 years of the POLAMI project introduced in an accompanying paper (Agudo et al.). In the 3mm wavelength band, we acquired more than 2600 observations, and all but one of our 37 sample sources were detected, most of them several times. For most sources, the observed distribution of the degree of circular polarisation is broader than that of unpolarised calibrators, indicating that weak (<0.5%) circular polarisation is present most of the time. Our detection rate and the maximum degree of polarisation found, 2.0%, are comparable to previous surveys, all made at much longer wavelengths. We argue that the process generating circular polarisation must not be strongly wavelength dependent, and we propose that the widespread presence of circular polarisation in our short wavelength sample dominated by blazars is mostly due to Faraday conversion of the linearly polarised synchrotron radiation in the helical magnetic field of the jet. Circular polarisation is variable, most notably on time scales comparable to or shorter than our median sampling interval <1 month. Longer time scales of about one year are occasionally detected, but severely limited by the weakness of the signal. At variance with some longer wavelength investigations we find that the sign of circular polarisation changes in most sources, while only 7 sources, including 3 already known, have a strong preference for one sign. The degrees of circular and linear polarisation do not show any systematic correlation. We do find however one particular event where the two polarisation degrees vary in synchronism during a time span of 0.9 years. The paper also describes a novel method for calibrating the sign of circular polarisation observations., Comment: Accepted for publication in MNRAS. 16 pages, 8 figures, and 4 tables

We describe the POLAMI program for the monitoring of all four Stokes parameters of a sample of bright radio-loud active galactic nuclei with the IRAM 30m telescope at 3.5 and 1.3mm. The program started in October 2006 and accumulated, until August 2014, 2300 observations at 3.5mm, achieving a median time sampling interval of 22 days for the sample of 37 sources. This first paper explains the source selection, mostly blazars, the observing strategy and data calibration, and gives the details of the instrumental polarisation corrections. The sensitivity (1sigma) reached at 3.5mm is 0.5% (linear polarisation degree), 4.7 deg. (polarisation angle), and 0.23% (circular polarisation), while the corresponding values at 1.3mm are 1.7%, 9.9 deg., and 0.72%, respectively. The data quality is demonstrated by the time sequences of our calibrators Mars and Uranus. For the quasar 3C286, widely used as a linear polarisation calibrator, we give improved estimates of its linear polarisation, and show for the first time occasional detections of its weak circular polarisation, which suggests a small level of variability of the source at millimeter wavelengths., Comment: Accepted for publication in MNRAS. 10 pages, 5 figures, and 5 tables

We report on the first results of the POLAMI program, a simultaneous 3.5 and 1.3mm full-Stokes-polarisation monitoring of a sample of 36 of the brightest active galactic nuclei in the northern sky with the IRAM 30m Telescope. Through a systematic statistical study of data taken from October 2006 (from December 2009 for the case of the 1.3mm observations) to August 2014, we characterise the variability of the total flux density and linear polarisation. We find that all sources in the sample are highly variable in total flux density at both 3.5 and 1.3mm, as well as in spectral index, that is found to be optically thin in general. The total flux-density variability at 1.3mm is found, in general, to be faster, and to have larger amplitude and flatter PSD slopes than 3.5mm. The polarisation degree is on average larger at 1.3mm than at 3.5mm, by a factor of 2.6. The variability of linear polarisation degree is faster and has higher fractional amplitude than for total flux density, with the typical time scales during prominent polarisation peaks being significantly faster at 1.3mm than at 3.5mm. The polarisation angle at both 3.5 and 1.3mm is highly variable. Most of the sources show one or two excursions of >180 deg. on time scales from a few weeks to about a year during the course of our observations. The 3.5 and 1.3mm polarisation angle evolution follow rather well each other, although the 1.3mm data show a clear preference to more prominent variability on the short time scales, i.e. weeks. The data are compatible with multi-zone models of conical jets involving smaller emission regions for the shortest-wavelength emitting sites. Such smaller emitting regions should also be more efficient in energising particle populations. The data also favours the integrated emission at 1.3mm to have better ordered magnetic fields than the one at 3.5mm., Comment: Accepted for publication in MNRAS. 36 pages, 11 figures, and 6 tables

We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $\gamma$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 2014. We find that the $\gamma$-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful $\gamma$-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight ($\theta\sim$1.2$^{\circ}$) during the ejection of the knot and the $\gamma$-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the $\gamma$-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of near-infrared to ultraviolet photons is the probable mechanism for the $\gamma$-ray production., Comment: Accepted for publication in The Astrophysical Journal

We present the analysis of the radio jet evolution of the radio galaxy 3C 120 during a period of prolonged gamma-ray activity detected by the Fermi satellite between December 2012 and October 2014. We find a clear connection between the gamma-ray and radio emission, such that every period of gamma-ray activity is accompanied by the flaring of the mm-VLBI core and subsequent ejection of a new superluminal component. However, not all ejections of components are associated with gamma-ray events detectable by Fermi. Clear gamma-ray detections are obtained only when components are moving in a direction closer to our line of sight.This suggests that the observed gamma-ray emission depends not only on the interaction of moving components with the mm-VLBI core, but also on their orientation with respect to the observer. Timing of the gamma-ray detections and ejection of superluminal components locate the gamma-ray production to within almost 0.13 pc from the mm-VLBI core, which was previously estimated to lie about 0.24 pc from the central black hole. This corresponds to about twice the estimated extension of the broad line region, limiting the external photon field and therefore suggesting synchrotron self Compton as the most probable mechanism for the production of the gamma-ray emission. Alternatively, the interaction of components with the jet sheath can provide the necessary photon field to produced the observed gamma-rays by Compton scattering., Comment: Already accepted for publication in The Astrophysical Journal

The relativistic jet in M87 offers a unique opportunity for understanding the detailed jet structure and emission processes due to its proximity. In particular, the peculiar jet region HST-1 at ~1 arcsecond (or 80 pc, projected) from the nucleus has attracted a great deal of interest in the last decade because of its superluminal motion and broadband radio-to-X-ray outbursts, which may be further connected to the gamma-ray productions up to TeV energies. Over the last five years, we have been doing an intensive monitoring of HST-1 with EVN at 5GHz in order to examine the detailed structural evolution and its possible connection to high-energy activities. While this program already yielded interesting results in terms of the detailed mas-scale structure, proper motion measurements and structural variations, the recent HST-1 brightness is continuously decreasing at this frequency. To counter this, we have shifted our monitoring frequency to 1.7GHz from October 2013. This strategy successfully recovered the fainter emission that was missed in the last 5GHz session. Moreover, we again discovered the sudden emergence of a new component at the upstream edge of HST-1, demonstrating that the use of EVN 1.7GHz is indeed powerful to probe the current weak nature of HST-1. Here we report early results from the 1.7GHz monitoring as well as further progress on the long-term kinematic study., Comment: Proceedings of the 12th European VLBI Network Symposium and Users Meeting, 7-10 October 2014, Cagliari, Italy

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet., Comment: 11 pages, 7 figures and 3 tables. Accepted for publication in MNRAS. Online-only Table 1 is available as ancillary file with this submission

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

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells., Comment: 80 pages, 10 Tables, 31 figures, accepted for publication in ApJ

We present Very Long Baseline Array observations of the radio galaxy 3C120 at 5, 8, 12, and 15 GHz designed to study a peculiar stationary jet feature (hereafter C80) located ~80 mas from the core, which was previously shown to display a brightness temperature ~600 times lager than expected at such distances. The high sensitivity of the images -- obtained between December 2009 and June 2010 -- has revealed that C80 corresponds to the eastern flux density peak of an arc of emission (hereafter A80), downstream of which extends a large (~20 mas in size) bubble-like structure that resembles an inverted bow shock. The linearly polarized emission closely follows that of the total intensity in A80, with the electric vector position angle distributed nearly perpendicular to the arc-shaped structure. Despite the stationary nature of C80/A80, superluminal components with speeds up to ~3 c have been detected downstream from its position, resembling the behavior observed in the HST-1 emission complex in M87. The total and polarized emission of the C80/A80 structure, its lack of motion, and brightness temperature excess are best reproduced by a model based on synchrotron emission from a conical shock with cone opening angle \eta=10 degrees, jet viewing angle \theta=16 degrees, a completely tangled upstream magnetic field, and upstream Lorentz factor \gamma=8.4. The good agreement between our observations and numerical modeling leads us to conclude that the peculiar feature associated with C80/A80 corresponds to a conical recollimation shock in the jet of 3C120 located at a de-projected distance of ~190 pc downstream from the nucleus., Comment: Accepted for publication in ApJ

(Context.) A number of millimeter and submillimeter facilities with linear polarization observing capabilities have started operating during last years. These facilities, as well as other previous millimeter telescopes and interferometers, require bright and stable linear polarization calibrators to calibrate new instruments and to monitor their instrumental polarization. The current limited number of adequate calibrators implies difficulties in the acquisition of these calibration observations. (Aims.) Looking for additional linear polarization calibrators in the millimeter spectral range, in mid-2006 we started monitoring 3C 286, a standard and highly stable polarization calibrator for radio observations. (Methods.) Here we present the 3 and 1 mm monitoring observations obtained between September 2006 and January 2012 with the XPOL polarimeter on the IRAM 30 m Millimeter Telescope. (Results.) Our observations show that 3C 286 is a bright source of constant total flux with 3 mm flux density S_3mm = (0.91 \pm 0.02) Jy. The 3mm linear polarization degree (p_3mm =[13.5\pm0.3]%) and polarization angle (chi_3mm =[37.3\pm0.8]deg.,expressed in the equatorial coordinate system) are also constant during the time span of our observations. Although with poorer time sampling and signal-to-noise ratio, our 1 mm observations of 3C 286 are also reproduced by a constant source of 1 mm flux density (S_1mm = [0.30 \pm 0.03] Jy), polarization fraction (p_1mm = [14.4 \pm 1.8] %), and polarization angle (chi_1mm = [33.1 \pm 5.7]deg.). (Conclusions.) This, together with the previously known compact structure of 3C 286 -extended by ~3.5" in the sky- allow us to propose 3C 286 as a new calibrator for both single dish and interferometric polarization observations at 3 mm, and possibly at shorter wavelengths., Comment: Accepted for publication in A&A. 7 pages, 4 figures, 8 tables. Updated data sets with regard to previous version. New discussion about multi frequency properties of the source. Section 3.3, Figures 3 and 4, and Tables 7 and 8 are new