Showing 1–50 of 76 results for author: Myserlis, I

X-ray and multiwavelength polarization of Mrk 501 from 2022 to 2023

Authors: Chien-Ting J. Chen, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Laura Di Gesu, Alessandro Di Marco, Steven R. Ehlert, Manel Errando, Michela Negro, Svetlana G. Jorstad, Alan P. Marscher, Kinwah Wu, Iván Agudo, Juri Poutanen, Tsunefumi Mizuno, Pouya M. Kouch, Elina Lindfors, George A. Borman, Tatiana S. Grishina, Evgenia N. Kopatskaya, Elena G. Larionova, Daria A. Morozova, Sergey S. Savchenko, Ivan S. Troitsky, Yulia V. Troitskaya , et al. (121 additional authors not shown)

Abstract: 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 optic… ▽ More 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σ$ 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. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: Accepted for publication in ApJ

IXPE observation of PKS 2155-304 reveals the most highly polarized blazar

Authors: Pouya M. Kouch, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Fabrizio Tavecchio, Alan P. Marscher, Herman L. Marshall, Steven R. Ehlert, Laura Di Gesu, Svetlana G. Jorstad, Iván Agudo, Grzegorz M. Madejski, Roger W. Romani, Manel Errando, Elina Lindfors, Kari Nilsson, Ella Toppari, Stephen B. Potter, Ryo Imazawa, Mahito Sasada, Yasushi Fukazawa, Koji S. Kawabata, Makoto Uemura, Tsunefumi Mizuno, Tatsuya Nakaoka , et al. (111 additional authors not shown)

Abstract: We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the… ▽ More We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

The flaring activity of blazar AO 0235+164 during year 2021

Authors: Juan Escudero Pedrosa, Iván Agudo, Till Moritz, Alan P. Marscher, Svetlana Jorstad, Andrea Tramacere, Carolina Casadio, Clemens Thum, Ioannis Myserlis, Albrecht Sievers, Jorge Otero-Santos, Daniel Morcuende, Rubén López-Coto, Filippo D'Ammando, Giacomo Bonnoli, Mark Gurwell, José Luis Gómez, Ramprasad Rao, Garrett Keating

Abstract: 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 sou… ▽ More 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. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted 15 May 2024

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, J. Bernete, A. Berti, J. Besenrieder , et al. (239 additional authors not shown)

Abstract: We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and… ▽ More We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

Detection of X-ray Polarization from the Blazar 1ES 1959+650 with the Imaging X-ray Polarimetry Explorer

Authors: Manel Errando, Ioannis Liodakis, Alan P. Marscher, Herman L. Marshall, Riccardo Middei, Michela Negro, Abel Lawrence Peirson, Matteo Perri, Simonetta Puccetti, Pazit L. Rabinowitz, Iván Agudo, Svetlana G. Jorstad, Sergey S. Savchenko, Dmitry Blinov, Ioakeim G. Bourbah, Sebastian Kiehlmann, Evangelos Kontopodis, Nikos Mandarakas, Stylianos Romanopoulos, Raphael Skalidis, Anna Vervelaki, Francisco José Aceituno, Maria I. Bernardos, Giacomo Bonnoli, Víctor Casanova , et al. (121 additional authors not shown)

Abstract: 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 part… ▽ More 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 $Π_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $ψ_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $Π_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $Π_\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. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 963 (2024) 5

Lost in the curve: Investigating the disappearing knots in the blazar 3C 454.3

Authors: Efthalia Traianou, Thomas P. Krichbaum, José L. Gómez, Rocco Lico, Georgios Filippos Paraschos, Ilje Cho, Eduardo Ros, Guang-Yao Zhao, Ioannis Liodakis, Rohan Dahale, Teresa Toscano, Antonio Fuentes, Marianna Foschi, Carolina Casadio, Nicholas MacDonald, Jae-Young Kim, Olivier Hervet, Svetlana Jorstad, Andrei P. Lobanov, Jeffrey Hodgson, Ioannis Myserlis, Ivan Agudo, Anton J. Zensus, Alan P. Marscher

Abstract: 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.… ▽ More 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. △ Less

Submitted 24 December, 2023; originally announced December 2023.

Comments: 15 pages, 7 figures

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland , et al. (229 additional authors not shown)

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

Unveiling the Bent Jet Structure and Polarization of OJ 287 at 1.7 GHz with Space VLBI

Authors: Ilje Cho, José L. Gómez, Rocco Lico, Guang-Yao Zhao, Efthalia Traianou, Rohan Dahale, Antonio Fuentes, Teresa Toscano, Marianna Foschi, Yuri Y. Kovalev, Andrei Lobanov, Alexander B. Pushkarev, Leonid I. Gurvits, Jae-Young Kim, Mikhail Lisakov, Petr Voitsik, Ioannis Myserlis, Felix Pötzl, Eduardo Ros

Abstract: We present total intensity and linear polarization images of OJ287 at 1.68GHz, obtained through space-based VLBI observations with RadioAstron on April 16, 2016. The observations were conducted using a ground array consisting of the VLBA and the EVN. Ground-space fringes were detected with a maximum projected baseline length of 5.6 Earth's diameter, resulting in an angular resolution of 530 uas. W… ▽ More We present total intensity and linear polarization images of OJ287 at 1.68GHz, obtained through space-based VLBI observations with RadioAstron on April 16, 2016. The observations were conducted using a ground array consisting of the VLBA and the EVN. Ground-space fringes were detected with a maximum projected baseline length of 5.6 Earth's diameter, resulting in an angular resolution of 530 uas. With this unprecedented resolution at such a low frequency, the progressively bending jet structure of OJ287 has been resolved up to 10 pc of the projected distance from the radio core. In comparison with close-in-time VLBI observations at 15, 43, 86 GHz from MOJAVE and VLBA-BU-BLAZAR monitoring projects, we obtain the spectral index map showing the opaque core and optically thin jet components. The optically thick core has a brightness temperature of 10$^{13}$ K, and is further resolved into two sub-components at higher frequencies labeled C1 and C2. These sub-components exhibit a transition from optically thick to thin, with a SSA turnover frequency estimated to be 33 and 11.5 GHz, and a turnover flux density 4 and 0.7 Jy, respectively. Assuming a Doppler boosting factor of 10, the SSA values provide the estimate of the magnetic field strengths from SSA of 3.4 G for C1 and 1.0 G for C2. The magnetic field strengths assuming equipartition arguments are also estimated as 2.6 G and 1.6 G, respectively. The integrated degree of linear polarization is found to be approximately 2.5 %, with the electric vector position angle being well aligned with the local jet direction at the core region. This alignment suggests a predominant toroidal magnetic field, which is in agreement with the jet formation model that requires a helical magnetic field anchored to either the black hole ergosphere or the accretion disk. Further downstream, the jet seems to be predominantly threaded by a poloidal magnetic field. △ Less

Submitted 25 March, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 15 pages, 11 figures, 3 tables

Journal ref: A&A 2024, 683, A248

The Repeating Flaring Activity of Blazar AO 0235+164

Authors: Juan Escudero Pedrosa, Iván Agudo, Andrea Tramacere, Alan P. Marscher, Svetlana Jorstad, Z. R. Weaver, Carolina Casadio, Clemens Thum, Ioannis Myserlis, Antonio Fuentes, Efthalia Traianou, Jae-Young Kim, Joana Kramer, Rubén López-Coto, Filippo D'Ammando, M. Bernardos, Giacomo Bonnoli, Dmitriy A. Blinov, G. A. Borman, T. S. Grishina, V. A. Hagen-Thorn, E. N. Kopatskaya, E. G. Larionova, V. M. Larionov, L. V. Larionova , et al. (5 additional authors not shown)

Abstract: 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 ph… ▽ More 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. △ Less

Submitted 2 November, 2023; originally announced November 2023.

A Near Magnetic-to-kinetic Energy Equipartition Flare from the Relativistic Jet in AO 0235+164 during 2013-2019

Authors: Whee Yeon Cheong, Sang-Sung Lee, Sang-Hyun Kim, Sincheol Kang, Jae Young Kim, Bindu Rani, Anthony C. S. Readhead, Sebastian Kiehlmann, Anne Lähteenmäki, Merja Tornikoski, Joni Tammi, Venkatessh Ramakrishnan, Iván Agudo, Antonio Fuentes, Efthalia Traianou, Juan Escudero, Clemens Thum, Ioannis Myserlis, Carolina Casadio, Mark Gurwell

Abstract: 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 $γ$-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 ejecte… ▽ More 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 $γ$-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 $δ_{\rm var}=28.5\pm8.4$, a bulk Lorentz factor of $Γ=16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $θ_{\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. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 13 pages, 11 figures, 4 tables; Accepted for publication in MNRAS

Observations of Low and Intermediate Spectral Peak Blazars with the Imaging X-ray Polarimetry Explorer

Authors: Herman L. Marshall, Ioannis Liodakis, Alan P. Marscher, Niccolo Di Lalla, Svetlana G. Jorstad, Dawoon E. Kim, Riccardo Middei, Michela Negro, Nicola Omodei, Abel L. Peirson, Matteo Perri, Simonetta Puccetti, Ivan Agudo, Giacomo Bonnoli, Andrei V. Berdyugin, Elisabetta Cavazzuti, Nicole Rodriguez Cavero, Immacolata Donnarumma, Laura Di Gesu, Jenni Jormanainen, Henric Krawczynski, Elina Lindfors, Frederic Marin, Francesco Massaro, Luigi Pacciani , et al. (133 additional authors not shown)

Abstract: 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$σ$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polari… ▽ More 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$σ$ 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. △ Less

Submitted 19 October, 2023; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures, 6 tables; submitted to the Astrophysical Journal

Magnetic Field Properties inside the Jet of Mrk 421: Multiwavelength Polarimetry Including the Imaging X-ray Polarimetry Explorer

Authors: Dawoon E. Kim, Laura Di Gesu, Ioannis Liodakis, Alan P. Marscher, Svetlana G. Jorstad, Riccardo Midde, Herman L. Marshall, Luigi Pacciani, Iván Agudo, Fabrizio Tavecchio, Nicolò Cibrario, Stefano Tugliani, Raffaella Bonino, Michela Negro, Simonetta Puccetti, Francesco Tombesi, Enrico Costa, Immacolata Donnarumma, Paolo Soffitta, Tsunefumi Mizuno, Yasushi Fukazawa, Koji S. Kawabata, Tatsuya Nakaoka, Makoto Uemura, Ryo Imazawa , et al. (111 additional authors not shown)

Abstract: We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $Π_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $ψ_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time varia… ▽ More We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $Π_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $ψ_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time variability analysis, we find a significant episodic variation in $ψ_{\rm X}$. During 7 months from the first IXPE pointing of Mrk 421 in 2022 May, $ψ_{\rm X}$ varied across the range of 0$^{\circ}$ to 180$^{\circ}$, while $Π_{\rm X}$ maintained similar values within $\sim$10-15$\%$. Furthermore, a swing in $ψ_{\rm X}$ in 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that the X-ray polarization degree was generally $\sim$2-3 times greater than that at longer wavelengths, while the polarization angle fluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that rotation of $ψ$ occurred in the opposite direction with respect to the rotation of $ψ_{\rm X}$ over longer timescales at similar epochs. The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. The accompanying spectral variation during the $ψ_{\rm X}$ rotation can be explained by a fluctuation in the physical conditions, e.g., in the energy distribution of relativistic electrons. The opposite rotation direction of $ψ$ between the X-ray and longer-wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, 4 tables; Accepted for publication in A&A

Journal ref: A&A 681, A12 (2024)

X-ray Polarization of the BL Lac Type Blazar 1ES 0229+200

Authors: Steven R. Ehlert, Ioannis Liodakis, Riccardo Middei, Alan P. Marscher, Fabrizio Tavecchio, Iván Agudo, Pouya M. Kouch, Elina Lindfors, Kari Nilsson, Ioannis Myserlis, Mark Gurwell, Ramprasad Rao, Francisco Jose Aceituno, Giacomo Bonnoli, Victor Casanova, Beatriz Agiz-Gonzalez, Juan Escudero, Jorge Otero Santos, Alfredo Sota, Emmanouil Angelakis, Alexander Kraus, Garrett K. Keating, Lucio A. Antonelli, Matteo Bachetti, Luca Baldini , et al. (88 additional authors not shown)

Abstract: We present polarization measurements in the $2-8 \thinspace \mathrm{keV}$ band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be $Π_{X} = 17.9 \pm 2.8 \%$ at an electric-vector position angle… ▽ More We present polarization measurements in the $2-8 \thinspace \mathrm{keV}$ band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be $Π_{X} = 17.9 \pm 2.8 \%$ at an electric-vector position angle $ψ_X = 25.0 \pm 4.6^{\circ}$ using a spectropolarimetric fit from joint IXPE and XMM-Newton observations. There is no evidence for the polarization degree or angle varying significantly with energy or time on both short time scales (hours) or longer time scales (days). The contemporaneous polarization degree at optical wavelengths was $>$7$\times$ lower, making 1ES 0229+200 the most strongly chromatic blazar yet observed. This high X-ray polarization compared to the optical provides further support that X-ray emission in high-peaked blazars originates in shock-accelerated, energy-stratified electron populations, but is in tension with many recent modeling efforts attempting to reproduce the spectral energy distribution of 1ES 0229+200 which attribute the extremely high energy synchrotron and Compton peaks to Fermi acceleration in the vicinity of strongly turbulent magnetic fields. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 17 Pages, 6 Figures, Resubmitted to ApJ after addressing referee comments

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

IXPE and multi-wavelength observations of blazar PG 1553+113 reveal an orphan optical polarization swing

Authors: Riccardo Middei, Matteo Perri, Simonetta Puccetti, Ioannis Liodakis, Laura Di Gesu, Alan P. Marscher, Nicole Rodriguez Cavero, Fabrizio Tavecchio, Immacolata Donnarumma, Marco Laurenti, Svetlana G. Jorstad, Iván Agudo, Herman L. Marshall, Luigi Pacciani, Dawoon E. Kim, Francisco José Aceituno, Giacomo Bonnoli, Víctor Casanova, Beatriz Agís-González, Alfredo Sota, Carolina Casadio, Juan Escudero, Ioannis Myserlis, Albrecht Sievers, Pouya M. Kouch , et al. (97 additional authors not shown)

Abstract: 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 Imagin… ▽ More 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 $ψ_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. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for pubblication in The Astrophysical Journal Letters. 3 Figures, 13 pages

The RoboPol sample of optical polarimetric standards

Authors: D. Blinov, S. Maharana, F. Bouzelou, C. Casadio, E. Gjerløw, J. Jormanainen, S. Kiehlmann, J. A. Kypriotakis, I. Liodakis, N. Mandarakas, L. Markopoulioti, G. V. Panopoulou, V. Pelgrims, A. Pouliasi, S. Romanopoulos, R. Skalidis, R. M. Anche, E. Angelakis, J. Antoniadis, B. J. Medhi, T. Hovatta, A. Kus, N. Kylafis, A. Mahabal, I. Myserlis , et al. (12 additional authors not shown)

Abstract: Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have po… ▽ More Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have polarization parameters that differ from their cataloged values. Our goal is to establish a sample of stable standards suitable for calibrating linear optical polarimeters with an accuracy down to $10^{-3}$ in fractional polarization. For five years, we have been running a monitoring campaign of a sample of standard candidates comprised of 107 stars distributed across the northern sky. We analyzed the variability of the linear polarization of these stars, taking into account the non-Gaussian nature of fractional polarization measurements. For a subsample of nine stars, we also performed multiband polarization measurements. We created a new catalog of 65 stars (see Table 2) that are stable, have small uncertainties of measured polarimetric parameters, and can be used as calibrators of polarimeters at medium- and large-size telescopes. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Journal ref: A&A 677, A144 (2023)

X-ray Polarization of BL Lacertae in Outburst

Authors: Abel L. Peirson, Michela Negro, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Alan P. Marscher, Herman L. Marshall, Luigi Pacciani, Roger W. Romani, Kinwah Wu, Alessandro Di Marco, Niccolo Di Lalla, Nicola Omodei, Svetlana G. Jorstad, Ivan Agudo, Pouya M. Kouch, Elina Lindfors, Francisco Jose Aceituno, Maria I. Bernardos, Giacomo Bonnoli, Victor Casanova, Maya Garcia-Comas, Beatriz Agis-Gonzalez, Cesar Husillos, Alessandro Marchini , et al. (106 additional authors not shown)

Abstract: We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$γ$-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), {\… ▽ More We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$γ$-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 $Π_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $ψ_{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. △ Less

Submitted 3 July, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: Accepted to ApJL. 16 pages, 8 figures

Discovery of X-ray polarization angle rotation in active galaxy Mrk 421

Authors: Laura Di Gesu, Herman L. Marshall, Steven R. Ehlert, Dawoon E. Kim, Immacolata Donnarumma, Fabrizio Tavecchio, Ioannis Liodakis, Sebastian Kiehlmann, Iván Agudo, Svetlana G. Jorstad, Fabio Muleri, Alan P. Marscher, Simonetta Puccetti, Riccardo Middei, Matteo Perri, Luigi Pacciani, Michela Negro, Roger W. Romani, Alessandro Di Marco, Dmitry Blinov, Ioakeim G. Bourbah, Evangelos Kontopodis, Nikos Mandarakas, Stylianos Romanopoulos, Raphael Skalidis , et al. (118 additional authors not shown)

Abstract: 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 $Ψ$. Here we report the discovery of a $Ψ_{\mathrm x}$ rotation in the X-ray band in th… ▽ More 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 $Ψ$. Here we report the discovery of a $Ψ_{\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, $Ψ_{\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 ($Π_{\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 $Ψ$ 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. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

MOMO VI: Multifrequency radio variability of the blazar OJ 287 from 2015-2022, absence of predicted 2021 precursor-flare activity, and a new binary interpretation of the 2016/2017 outburst

Authors: S. Komossa, A. Kraus, D. Grupe, A. G. Gonzalez, M. A. Gurwell, L. C. Gallo, F. K. Liu, I. Myserlis, T. P. Krichbaum, S. Laine, U. Bach, J. L. Gomez, M. L. Parker, S. Yao, M. Berton

Abstract: Based on our dedicated Swift monitoring program, MOMO, OJ 287 is one of the best-monitored blazars in the X-ray--UV--optical regime. Here, we report results from our accompanying, dense, multi-frequency (1.4--44 GHz) radio monitoring of OJ 287 between 2015 and 2022 covering a broad range of activity states. Fermi gamma-ray observations are added. We characterize the radio flux and spectral variabi… ▽ More Based on our dedicated Swift monitoring program, MOMO, OJ 287 is one of the best-monitored blazars in the X-ray--UV--optical regime. Here, we report results from our accompanying, dense, multi-frequency (1.4--44 GHz) radio monitoring of OJ 287 between 2015 and 2022 covering a broad range of activity states. Fermi gamma-ray observations are added. We characterize the radio flux and spectral variability in detail, including DCF and other variability analyses, and discuss its connection with the multiwavelength emission. Deep fades of radio and optical--UV fluxes are found to occur every 1--2 years. Further, it is shown that a precursor flare of thermal bremsstrahlung predicted by one of the binary supermassive black hole (SMBH) models of OJ 287 was absent. We then focus on the nature of the extraordinary, nonthermal 2016/2017 outburst that we initially discovered with Swift. We interpret it as the latest of the famous optical double-peaked outbursts of OJ 287, favoring binary scenarios that do not require a highly precessing secondary SMBH. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: ApJ, in press. The full radio-data tables will be available in the online material of the journal. Note that in a follow-up publication (MNRAS Letters, in press), the black hole mass of OJ 287 will be revised

X-ray Polarization Observations of BL Lacertae

Authors: Riccardo Middei, Ioannis Liodakis, Matteo Perri, Simonetta Puccetti, Elisabetta Cavazzuti, Laura Di Gesu, Steven R. Ehlert, Grzegorz Madejski, Alan P. Marscher, Herman L. Marshall, Fabio Muleri, Michela Negro, Svetlana G. Jorstad, Beatriz Agís-González, Iván Agudo, Giacomo Bonnoli, Maria I. Bernardos, Víctor Casanova, Maya García-Comas, César Husillos, Alessandro Marchini, Alfredo Sota, Pouya M. Kouch, George A. Borman, Evgenia N. Kopatskaya , et al. (121 additional authors not shown)

Abstract: 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 responsibl… ▽ More 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 $Π_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. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 17 pages, 5 figures, accepted for publication in ApJL

Journal ref: 2023, ApJ Letters, 942, L10

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (300 additional authors not shown)

Abstract: We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the… ▽ More We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

The X-ray Polarization View of Mrk~421 in an Average Flux State as Observed by the Imaging X-ray Polarimetry Explorer

Authors: Laura Di Gesu, Immacolata Donnarumma, Fabrizio Tavecchio, Ivan Agudo, Thibault Barnounin, Nicolò Cibrario, Niccolò Di Lalla, Alessandro Di Marco, Juan Escudero, Manel Errando, Svetlana G. Jorstad, Dawoon Kim, Pouya M. Kouch, Elina Lindfors, Ioannis Liodakis, Grzegorz Madejski, Herman L. Marshall, Alan P. Marscher, Riccardo Middei, Fabio Muleri, Ioannis Myserlis, Michela Negro, Nicola Omodei, Luigi Pacciani, Alessandro Paggi , et al. (78 additional authors not shown)

Abstract: Particle acceleration mechanisms in supermassive black hole jets, such as shock acceleration, magnetic reconnection, and turbulence, are expected to have observable signatures in the multi-wavelength polarization properties of blazars. The recent launch of the Imaging X-ray Polarimetry Explorer (IXPE) enables us, for the first time, to use polarization in the X-ray band (2-8 keV) to probe the prop… ▽ More Particle acceleration mechanisms in supermassive black hole jets, such as shock acceleration, magnetic reconnection, and turbulence, are expected to have observable signatures in the multi-wavelength polarization properties of blazars. The recent launch of the Imaging X-ray Polarimetry Explorer (IXPE) enables us, for the first time, to use polarization in the X-ray band (2-8 keV) to probe the properties of the jet synchrotron emission in high-frequency-peaked BL Lac objects (HSPs). We report the discovery of X-ray linear polarization (degree $Π_{\rm x}=15\pm$2\% and electric-vector position angle $Ψ_{\rm x}=35^\circ\pm4^\circ$) from the jet of the HSP Mrk~421 in an average X-ray flux state. At the same time, the degree of polarization at optical, infrared, and millimeter wavelengths was found to be lower by at least a factor of 3. During the IXPE pointing, the X-ray flux of the source increased by a factor of 2.2, while the polarization behavior was consistent with no variability. The higher level of $Π_{\rm x}$ compared to longer wavelengths, and the absence of significant polarization variability, suggest a shock as the most likely X-ray emission site in the jet of Mrk 421 during the observation. The multiwavelength polarization properties are consistent with an energy-stratified electron population, where the particles emitting at longer wavelengths are located farther from the acceleration site, where they experience a more disordered magnetic field. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: accepted for ApJ Letters

Polarized Blazar X-rays imply particle acceleration in shocks

Authors: Ioannis Liodakis, Alan P. Marscher, Iván Agudo, Andrei V. Berdyugin, Maria I. Bernardos, Giacomo Bonnoli, George A. Borman, Carolina Casadio, Víctor Casanova, Elisabetta Cavazzuti, Nicole Rodriguez Cavero, Laura Di Gesu, Niccoló Di Lalla, Immacolata Donnarumma, Steven R. Ehlert, Manel Errando, Juan Escudero, Maya García-Comas, Beatriz Agís-González, César Husillos, Jenni Jormanainen, Svetlana G. Jorstad, Masato Kagitani, Evgenia N. Kopatskaya, Vadim Kravtsov , et al. (103 additional authors not shown)

Abstract: 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 magnet… ▽ More 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 $Π_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. △ Less

Submitted 3 September, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 34 pages, 4 figures, 4 tables

Journal ref: 2022, Nature, 611, 677-681

MOMO V. Effelsberg, Swift and Fermi study of the blazar and supermassive binary black hole candidate OJ 287 in a period of high activity

Authors: S. Komossa, D. Grupe, A. Kraus, A. Gonzalez, L. C. Gallo, M. J. Valtonen, S. Laine, T. P. Krichbaum, M. A. Gurwell, J. L. Gomez, S. Ciprini, I. Myserlis, U. Bach

Abstract: We report results from our ongoing project MOMO (Multiwavelength Observations and Modelling of OJ 287). In this latest publication of a sequence, we combine our Swift UVOT--XRT and Effelsberg radio data (2.6-44 GHz) between 2019 and 2022.04 with public SMA data and gamma-ray data from the Fermi satellite. The observational epoch covers OJ 287 in a high state of activity from radio to X-rays. The e… ▽ More We report results from our ongoing project MOMO (Multiwavelength Observations and Modelling of OJ 287). In this latest publication of a sequence, we combine our Swift UVOT--XRT and Effelsberg radio data (2.6-44 GHz) between 2019 and 2022.04 with public SMA data and gamma-ray data from the Fermi satellite. The observational epoch covers OJ 287 in a high state of activity from radio to X-rays. The epoch also covers two major events predicted by the binary supermassive black hole (SMBH) model of OJ 287. Spectral and timing analyses clearly establish: a new UV-optical minimum state in 2021 December at an epoch where the secondary SMBH is predicted to cross the disk surrounding the primary SMBH; an overall low level of gamma-ray activity in comparison to pre-2017 epochs; the presence of a remarkable, long-lasting UV--optical flare event of intermediate amplitude in 2020--2021; a high level of activity in the radio band with multiple flares; and particularly a bright, ongoing radio flare peaking in 2021 November that may be associated with a gamma-ray flare, the strongest in 6 years. Several explanations for the UV--optical minimum state are explored, including the possibility that a secondary SMBH launches a temporary jet, but the observations are best explained by variability associated with the main jet. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 16 pages, 11 figures. Published in MNRAS (submitted in January, accepted in March)

New jet feature in the parsec-scale jet of the blazar OJ287 connected to the 2017 teraelectronvolt flaring activity

Authors: R. Lico, C. Casadio, S. G. Jorstad, J. L. Gomez, A. P. Marscher, E. Traianou, J. Y. Kim, G. Y. Zhao, A. Fuentes, I. Cho, T. P. Krichbaum, O. Hervet, S. O'Brien, B. Boccardi, I. Myserlis, I. Agudo, A. Alberdi, Z. R. Weaver, J. A. Zensus

Abstract: In February 2017 the blazar OJ287, one of the best super-massive binary-black-hole-system candidates, was detected for the first time at very high energies (VHEs; E>100GeV) with the ground-based gamma-ray observatory VERITAS. Very high energy gamma rays are thought to be produced in the near vicinity of the central engine in active galactic nuclei. For this reason, and with the main goal of provid… ▽ More In February 2017 the blazar OJ287, one of the best super-massive binary-black-hole-system candidates, was detected for the first time at very high energies (VHEs; E>100GeV) with the ground-based gamma-ray observatory VERITAS. Very high energy gamma rays are thought to be produced in the near vicinity of the central engine in active galactic nuclei. For this reason, and with the main goal of providing useful information for the characterization of the physical mechanisms connected with the observed teraelectronvolt flaring event, we investigate the parsec-scale source properties by means of high-resolution very long baseline interferometry observations. We use 86 GHz Global Millimeter-VLBI Array (GMVA) observations from 2015 to 2017 and combine them with additional multiwavelength radio observations at different frequencies from other monitoring programs. We investigate the source structure by modeling the brightness distribution with two-dimensional Gaussian components in the visibility plane. In the GMVA epoch following the source VHE activity, we find a new jet feature (labeled K) at about 0.2 mas from the core region and located in between two quasi-stationary components (labeled S1 and S2). Multiple periods of enhanced activity are detected at different radio frequencies before and during the VHE flaring state. Based on the findings of this work, we identify as a possible trigger for the VHE flaring emission during the early months of 2017 the passage of a new jet feature through a recollimation shock (represented by the model-fit component S1) in a region of the jet located at a de-projected distance of about 10 pc from the radio core. △ Less

Submitted 11 February, 2022; v1 submitted 5 February, 2022; originally announced February 2022.

Comments: 5 pages, 3 figure, 1 table

Journal ref: A&A 658, L10 (2022)

Magnetic Field Strengths of the Synchrotron Self-Absorption Region in the Jet of CTA 102 During Radio Flares

Authors: Sang-Hyun Kim, Sang-Sung Lee, Jee Won Lee, Jeffrey A. Hodgson, Sincheol Kang, Juan-Carlos Algaba, Jae-Young Kim, Mark Hodges, Ivan Agudo, Antonio Fuentes, Juan Escudero, Ioannis Myserlis, Efthalia Traianou, Anne Lähteenmäki, Merja Tornikoski, Joni Tammi, Venkatessh Ramakrishnan, Emilia Järvelä

Abstract: CTA 102 is a blazar implying that its relativistic jet points towards Earth and emits synchrotron radiation produced by energetic particles gyrating in the magnetic field. This study aims to figure out the physical origins of radio flares in the jet, including the connection between the magnetic field and the radio flares. The dataset in the range 2.6-343.5 GHz was collected over a period of 5.5 y… ▽ More CTA 102 is a blazar implying that its relativistic jet points towards Earth and emits synchrotron radiation produced by energetic particles gyrating in the magnetic field. This study aims to figure out the physical origins of radio flares in the jet, including the connection between the magnetic field and the radio flares. The dataset in the range 2.6-343.5 GHz was collected over a period of 5.5 years (2012 November 20-2018 September 23). During the data collection period, seven flares at 15 GHz with a range of the variability time-scale of roughly 26-171 days were detected. The quasi-simultaneous radio data were used to investigate the synchrotron spectrum of the source. We found that the synchrotron radiation is self-absorbed. The turnover frequency and the peak flux density of the synchrotron self-absorption (SSA) spectra are in the ranges of 42-167 GHz and 0.9-10.2 Jy, respectively. From the SSA spectra, we derived the SSA magnetic field strengths to be 9.20 mG, 12.28 mG, and 50.97 mG on 2013 December 24, 2014 February 28, and 2018 January 13, respectively. We also derived the equipartition magnetic field strengths to be in the range 24-109 mG. The equipartition magnetic field strengths are larger than the SSA magnetic field strengths in most cases, which indicates that particle energy mainly dominates in the jet. Our results suggest that the flares in the jet of CTA 102 originated due to particle acceleration. We propose the possible mechanisms of particle acceleration. △ Less

Submitted 27 November, 2021; originally announced November 2021.

Comments: 19 pages, 8 figures, 6 tables, accepted for publication in MNRAS

Probing the innermost regions of AGN jets and their magnetic fields with RadioAstron. V. Space and ground millimeter-VLBI imaging of OJ 287

Authors: Jose L. Gómez, Efthalia Traianou, Thomas P. Krichbaum, Andrei Lobanov, Antonio Fuentes, Rocco Lico, Guang-Yao Zhao, Gabriele Bruni, Yuri Y. Kovalev, Anne Lahteenmaki, Petr A. Voitsik, Mikhail M. Lisakov, Emmanouil Angelakis, Uwe Bach, Carolina Casadio, Ilje Cho, Lankeswar Dey, Achamveedu Gopakumar, Leonid Gurvits, Svetlana G. Jorstad, Yuri A. Kovalev, Matthew L. Lister, Alan P. Marscher, Ioannis Myserlis, Alexander Pushkarev , et al. (5 additional authors not shown)

Abstract: 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 dia… ▽ More 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μ$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. △ Less

Submitted 28 November, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

Journal ref: ApJ 924 (2022) 122

Project MOMO: Multiwavelength Observations and Modelling of OJ 287

Authors: S. Komossa, D. Grupe, A. Kraus, L. C. Gallo, A. Gonzalez, M. L. Parker, M. J. Valtonen, A. R. Hollett, U. Bach, J. L. Gómez, I. Myserlis, S. Ciprini

Abstract: Our project MOMO (Multiwavelength observations and modelling of OJ 287) consists of dedicated, dense, long-term flux and spectroscopic monitoring and deep follow-up observations of the blazar OJ 287 at >13 frequencies from the radio to the X-ray band since late 2015. In particular, we are using Swift to obtain optical-UV-X-ray spectral energy distributions (SEDs) and the Effelsberg telescope to ob… ▽ More Our project MOMO (Multiwavelength observations and modelling of OJ 287) consists of dedicated, dense, long-term flux and spectroscopic monitoring and deep follow-up observations of the blazar OJ 287 at >13 frequencies from the radio to the X-ray band since late 2015. In particular, we are using Swift to obtain optical-UV-X-ray spectral energy distributions (SEDs) and the Effelsberg telescope to obtain radio measurements between 2 and 40 GHz. MOMO is the densest long-term monitoring of OJ 287 involving X-rays and broad-band SEDs. The theoretical part of the project aims at understanding jet and accretion physics of the blazar central engine in general and the supermassive binary black hole scenario in particular. Results are presented in a sequence of publications and so far included: detection and detailed analysis of the bright 2016/17 and 2020 outbursts and the long-term light curve; Swift, XMM and NuSTAR spectroscopy of the 2020 outburst around maximum; and interpretation of selected events in the context of the binary black hole scenario of OJ 287 (papers I-IV). Here, we provide a description of the project MOMO, a summary of previous results, the latest results, and we discuss future prospects. △ Less

Submitted 30 June, 2021; originally announced July 2021.

Comments: Review. 18 pages, 11 figures. Submitted to Universe Special Issue "Panchromatic View of the Life-Cycle of AGN"

The Time-Dependent Distribution of Optical Polarization Angle Changes in Blazars

Authors: S. Kiehlmann, D. Blinov, I. Liodakis, V. Pavlidou, A. C. S. Readhead, E. Angelakis, C. Casadio, T. Hovatta, N. Kylafis, A. Mahabal, N. Mandarakas, I. Myserlis, G. V. Panopoulou, T. J. Pearson, A. Ramaprakash, P. Reig, R. Skalidis, A. Slowikowska, K. Tassis, J. A. Zensus

Abstract: At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-inte… ▽ More At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-interval between successive observations -- i.e. the cadence -- is short enough to ensure that the correct modulo 180° value is selected. This optimal cadence depends on the maximum intrinsic EVPA rotation speed in blazars, which is currently not known. In this paper we address the following questions for the RoboPol sample: What range of rotation speeds for rotations greater than 90° can we expect? What observation cadence is required to detect such rotations? Have rapid rotations been missed in EVPA rotation studies thus far? What fraction of data is affected by the ambiguity? And how likely are detected rotations affected by the ambiguity? We answer these questions with three seasons of optical polarimetric observations of a statistical sample of blazars sampled weekly with the RoboPol instrument and an additional season with daily observations. We model the distribution of EVPA changes on time scales from 1-30 days and estimate the fraction of changes exceeding 90°. We show that at least daily observations are necessary to measure >96% of optical EVPA variability in the RoboPol sample of blazars correctly and that intra-day observations are needed to measure the fastest rotations that have been seen thus far. △ Less

Submitted 14 July, 2021; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: accepted to MNRAS, 20 pages, 15 figures

An underlying universal pattern in galaxy halo magnetic fields

Authors: Ioannis Myserlis, Ioannis Contopoulos

Abstract: Magnetic fields in galaxy halos are in general very difficult to observe. Most recently, the CHANG-ES collaboration (Continuum HAlos in Nearby Galaxies - an EVLA Survey) investigated in detail the radio halos of 35 nearby edge-on spiral galaxies and detected large scale magnetic fields in 16 of them. We used the CHANG-ES radio polarization data to create Rotation Measure (RM) maps for all galaxies… ▽ More Magnetic fields in galaxy halos are in general very difficult to observe. Most recently, the CHANG-ES collaboration (Continuum HAlos in Nearby Galaxies - an EVLA Survey) investigated in detail the radio halos of 35 nearby edge-on spiral galaxies and detected large scale magnetic fields in 16 of them. We used the CHANG-ES radio polarization data to create Rotation Measure (RM) maps for all galaxies in the sample and stack them with the aim to amplify any underlying universal toroidal magnetic field pattern in the halo above and below the disk of the galaxy. We discovered a large-scale magnetic field in the central region of the stacked galaxy profile, attributable to an axial electric current that universally outflows from the center both above and below the plane of the disk. A similar symmetry-breaking has also been observed in astrophysical jets but never before in galaxy halos. This is an indication that galaxy halo magnetic fields are probably not generated by pure ideal magnetohydrodynamic (MHD) processes in the central regions of galaxies. One such promising physical mechanism is the Cosmic Battery operating in the innermost accretion disk around the central supermassive black hole. We anticipate that our discovery will stimulate a more general discussion on the origin of astrophysical magnetic fields. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 10 pages, 6 figures, accepted for publication in the 4. Extragalactic astronomy section of Astronomy & Astrophysics

Journal ref: A&A 649, A94 (2021)

RoboPol: AGN polarimetric monitoring data

Authors: D. Blinov, S. Kiehlmann, V. Pavlidou, G. V. Panopoulou, R. Skalidis, E. Angelakis, C. Casadio, E. N. Einoder, T. Hovatta, K. Kokolakis, A. Kougentakis, A. Kus, N. Kylafis, E. Kyritsis, A. Lalakos, I. Liodakis, S. Maharana, E. Makrydopoulou, N. Mandarakas, G. M. Maragkakis, I. Myserlis, I. Papadakis, G. Paterakis, T. J. Pearson, A. N. Ramaprakash , et al. (7 additional authors not shown)

Abstract: We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was… ▽ More We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was used to process and calibrate the entire dataset, and we make the data publicly available for use by the astronomical community. Average quantities summarising optopolarimetric behaviour (average degree of polarization, polarization variability index) are also provided for each source we have observed and for the time interval we have followed it. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: Accepted to MNRAS

Multiwavelength behaviour of the blazar 3C279: decade-long study from $γ$-ray to radio

Authors: V. M. Larionov, S. G. Jorstad, A. P. Marscher, M. Villata, C. M. Raiteri, P. S. Smith, I. Agudo, S. S. Savchenko, D. A. Morozova, J. A. Acosta-Pulido, M. F. Aller, H. D. Aller, T. S. Andreeva, A. A. Arkharov, R. Bachev, G. Bonnoli, G. A. Borman, V. Bozhilov, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, W. -P. Chen, G. Damljanovic, A. V. Dementyev , et al. (62 additional authors not shown)

Abstract: We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activi… ▽ More We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $γ$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $τ=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Journal ref: MNRAS 492 (2020) 3829

The relativistic jet of the $γ$-ray emitting narrow-line Seyfert 1 galaxy PKS J1222$+$0413

Authors: Daniel Kynoch, Hermine Landt, Martin J. Ward, Chris Done, Catherine Boisson, Mislav Baloković, Emmanouil Angelakis, Ioannis Myserlis

Abstract: We present a multi-frequency study of PKS J1222$+$0413 (4C$+$04.42), currently the highest redshift $γ$-ray emitting narrow-line Seyfert 1 ($γ$-NLS1). We assemble a broad spectral energy distribution (SED) including previously unpublished datasets: X-ray data obtained with the NuSTAR and Neil Gehrels Swift observatories; near-infrared, optical and UV spectroscopy obtained with VLT X-shooter; and m… ▽ More We present a multi-frequency study of PKS J1222$+$0413 (4C$+$04.42), currently the highest redshift $γ$-ray emitting narrow-line Seyfert 1 ($γ$-NLS1). We assemble a broad spectral energy distribution (SED) including previously unpublished datasets: X-ray data obtained with the NuSTAR and Neil Gehrels Swift observatories; near-infrared, optical and UV spectroscopy obtained with VLT X-shooter; and multiband radio data from the Effelsberg telescope. These new observations are supplemented by archival data from the literature. We apply physical models to the broadband SED, parameterising the accretion flow and jet emission to investigate the disc-jet connection. PKS J1222$+$0413 has a much greater black hole mass than most other NLS1s, $M_\mathrm{BH}\approx2\times10^{8}$ M$_\odot$, similar to those found in flat spectrum radio quasars (FSRQs). Therefore this source provides insight into how the jets of $γ$-NLS1s relate to those of FSRQs. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: 18 pages, 10 figures and 6 tables. Accepted for publication in MNRAS

High-Energy Polarimetry - a new window to probe extreme physics in AGN jets

Authors: B. Rani, H. Zhang, S. D. Hunter, F. Kislat, M. Böttcher, J. E. McEnery, D. J. Thompson, D. Giannios, F. Guo, H. Li, M. Baring, I. Agudo, S. Buson, M. Petropoulou, V. Pavlidou, E. Angelakis, I. Myserlis, Z. Wadiasingh, R. M. Curado da Silva, P. Kilian, S. Guiriec, V. V. Bozhilov, J. Hodgson, S. Antón, D. Kazanas , et al. (9 additional authors not shown)

Abstract: The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is m… ▽ More The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020 (Astronomy and Astrophysics Decadal Survey)

RoboPol: A four-channel optical imaging polarimeter

Authors: A. N. Ramaprakash, C. V. Rajarshi, H. K. Das, P. Khodade, D. Modi, G. Panopoulou, S. Maharana, D. Blinov, E. Angelakis, C. Casadio, L. Fuhrmann, T. Hovatta, S. Kiehlmann, O. G. King, N. Kylafis, A. Kougentakis, A. Kus, A. Mahabal, A. Marecki, I. Myserlis, G. Paterakis, E. Paleologou, I. Liodakis, I. Papadakis, I. Papamastorakis , et al. (8 additional authors not shown)

Abstract: We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters… ▽ More We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters (B, V and I), enabling multi-wavelength imaging polarimetry over a large field of view (13.6' $\times$ 13.6'). We demonstrate the accuracy and stability of the instrument throughout its five years of operation. Best performance is achieved within the central region of the field of view and in the R band. For such measurements the systematic uncertainty is below 0.1% in fractional linear polarization, $p$ (0.05% maximum likelihood). Throughout all observing seasons the instrumental polarization varies within 0.1% in $p$ and within 1$^\circ$ in polarization angle. △ Less

Submitted 22 February, 2019; originally announced February 2019.

Comments: 13 pages, 15 figures, accepted for publication in MNRAS

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

Authors: E. Angelakis, L. Fuhrmann, I. Myserlis, J. A. Zensus, I. Nestoras, V. Karamanavis, N. Marchili, T. P. Krichbaum, A. Kraus, J. P. Rachen

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

Submitted 7 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Section: Catalogs and data of Astronomy & Astrophysics

Search for AGN counterparts of unidentified Fermi-LAT sources with optical polarimetry: Demonstration of the technique

Authors: N. Mandarakas, D. Blinov, I. Liodakis, K. Kouroumpatzakis, A. Zezas, G. V. Panopoulou, I. Myserlis, E. Angelakis, T. Hovatta, S. Kiehlmann, K. Kokolakis, E. Paleologou, A. Pouliasi, R. Skalidis, V. Pavlidou

Abstract: The third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and pra… ▽ More The third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and practical tool in the hunt for counterparts of the unidentified gamma-ray sources (UGSs). Using data from the RoboPol project, we validated that a significant fraction of active galactic nuclei (AGN) associated with 3FGL sources can be identified due to their high optical polarization exceeding that of the field stars. We performed an optical polarimetric survey within $3σ$ uncertainties of four unidentified 3FGL sources. We discovered a previously unknown extragalactic object within the positional uncertainty of 3FGL J0221.2+2518. We obtained its spectrum and measured a redshift of $z=0.0609\pm0.0004$. Using these measurements and archival data we demonstrate that this source is a candidate counterpart for 3FGL J0221.2+2518 and most probably is a composite object: a star-forming galaxy accompanied by AGN. We conclude that polarimetry can be a powerful asset in the search for AGN candidate counterparts for unidentified Fermi sources. Future extensive polarimetric surveys at high galactic latitudes (e.g., PASIPHAE) will allow the association of a significant fraction of currently unidentified gamma-ray sources. △ Less

Submitted 12 February, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: accepted to A&A

Journal ref: A&A 623, A61 (2019)

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

Authors: I. Myserlis, S. Komossa, E. Angelakis, J. L. Gómez, V. Karamanavis, T. P. Krichbaum, U. Bach, D. Grupe

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

Submitted 5 September, 2018; originally announced September 2018.

Comments: 11 pages, 9 figures, 1 table, accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics on August 21, 2018

Journal ref: A&A 619, A88 (2018)

Optical polarisation variability of narrow line Seyfert 1 galaxies

Authors: E. Angelakis, S. Kiehlmann, I. Myserlis, D. Blinov, J. Eggen, R. Itoh, S. Komossa, N. Marchilli, J. A. Zensus

Abstract: We have monitored the R-band optical linear polarisation of ten jetted NLSy1 galaxies with the aim to quantify their variability and search for candidate long rotation of the polarisation plane. In all cases for which adequate datasets are available we observe significant variability of both the polarisation fraction and angle. In the best-sampled cases we identify candidate long rotations of the… ▽ More We have monitored the R-band optical linear polarisation of ten jetted NLSy1 galaxies with the aim to quantify their variability and search for candidate long rotation of the polarisation plane. In all cases for which adequate datasets are available we observe significant variability of both the polarisation fraction and angle. In the best-sampled cases we identify candidate long rotations of the polarisation plane. We present an approach that assesses the probability that the observed phenomenology is the result of pure noise. We conclude that although this possibility cannot be excluded it is much more likely that the EVPA undergoes an intrinsic evolution. We compute the most probable parameters of the intrinsic event which forecasts events consistent with the observations. In one case we find that the EVPA shows a preferred direction which, however, does not imply any dominance of a toroidal or poloidal component of the magnetic field at those scales. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: Presented at the conference "Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe" (9-13 April 2018, Padova Botanical Garden, Italy) and is included in Volume 328 of the PoS as article 019

Optical polarisation variability of radio loud narrow line Seyfert 1 galaxies. Search for long rotations of the polarisation plane

Authors: E. Angelakis, S. Kiehlmann, I. Myserlis, D. Blinov, J. Eggen, R. Itoh, N. Marchili, J. A. Zensus

Abstract: Narrow line Seyfert 1 galaxies (NLSy1s) constitute the AGN subclass associated with systematically smaller black hole masses. A few radio loud ones have been detected in MeV -- GeV energy bands by Fermi and evidence for the presence of blazar-like jets has been accumulated. In this study we wish to quantify the temporal behaviour of the optical polarisation, fraction and angle, for a selected samp… ▽ More Narrow line Seyfert 1 galaxies (NLSy1s) constitute the AGN subclass associated with systematically smaller black hole masses. A few radio loud ones have been detected in MeV -- GeV energy bands by Fermi and evidence for the presence of blazar-like jets has been accumulated. In this study we wish to quantify the temporal behaviour of the optical polarisation, fraction and angle, for a selected sample of radio loud NLSy1s. We also search for rotations of the polarisation plane similar to those commonly observed in blazars. We have conducted R-band optical polarisation monitoring of a sample of 10 RL NLSy1s 5 of which have been previously detected by Fermi. The dataset includes observations with the RoboPol, KANATA, Perkins and Steward polarimeters. In the cases where evidences for long rotations of the polarisation plane are found, we carry out numerical simulations to assess the probability that they are caused by intrinsically evolving EVPAs instead of observational noise. Even our moderately sampled sources show indications of variability, both in polarisation fraction and angle. For the four best sampled objects in our sample we find multiple periods of significant polarisation angle variability. In the two best sampled cases, namely J1505+0326 and J0324+3410, we find indications for three long rotations. We show that although noise can induce the observed behaviour, it is much more likely that the apparent rotation is caused by intrinsic evolution of the EVPA. To our knowledge this is the very first detection of such events in this class of sources. In the case of the largest dataset (J0324+3410) we find that the EVPA concentrates around a direction which is at 49.3\degr to the 15-GHz radio jet implying a projected magnetic field at an angle of 40.7\degr to that axis. △ Less

Submitted 6 July, 2018; originally announced July 2018.

Comments: Accepted for publication in section 2. Astrophysical processes of Astronomy and Astrophysics

Journal ref: A&A 618, A92 (2018)

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (165 additional authors not shown)

Abstract: The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light… ▽ More The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling. △ Less

Submitted 1 July, 2018; originally announced July 2018.

Journal ref: Astronomy & Astrophysics 619, A45 (2018)

Detection of persistent VHE gamma-ray emission from PKS 1510-089 by the MAGIC telescopes during low states between 2012 and 2017

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, G. Ceribella, A. Chatterjee , et al. (143 additional authors not shown)

Abstract: PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize… ▽ More PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5σ) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR. △ Less

Submitted 3 September, 2018; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: 12 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 619, A159 (2018)

The dependence of optical polarisation of blazars on the synchrotron component peak frequency

Authors: Emmanouil Angelakis, Dmitry Blinov, Markus Böttcher, Talvikki Hovatta, Sebastian Kiehlmann, Ioannis Myserlis, Vassiliki Pavlidou, J. Anton Zensus

Abstract: The RoboPol instrument and the relevant program was developed in order to conduct a systematic study of the optical polarisation variability of blazars. Driven by the discovery that long smooth rotations of the optical polarisation plane can be associated with the activity in other bands and especially in gamma rays, the program was meant to investigate the physical mechanisms causing them and qua… ▽ More The RoboPol instrument and the relevant program was developed in order to conduct a systematic study of the optical polarisation variability of blazars. Driven by the discovery that long smooth rotations of the optical polarisation plane can be associated with the activity in other bands and especially in gamma rays, the program was meant to investigate the physical mechanisms causing them and quantify the optical polarisation behaviour in blazars. Over the first three nominal observing seasons (2013, 2014 and 2015) RoboPol detected 40 rotations in 24 blazars by observing a gamma-ray-loud and gamma-ray-quite unbiassed sample of blazars, providing a reliable set of events for exploring the phenomenon. The obtain datasets provided the ground for a systematic quantification of the variability of the optical polarisation in such systems. In the following after a brief review of the discoveries that relate to the gamma-ray loudness of the sources we move on to discuss a simple jet model that explains the observed dichotomy in terms of polarisation between gamma-ray-loud and quite sources and the dependence of polarisation and the stability of the polarisation angle on the synchrotron peak frequency. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: Presented at the 7th International Fermi Symposium, 15-20 October 2017, Garmisch-Partenkirchen

Scale invariant jets: from blazars to microquasars

Authors: I. Liodakis, V. Pavlidou, I. Papadakis, E. Angelakis, N. Marchili, J. A. Zensus, L. Fuhrmann, V. Karamanavis, I. Myserlis, I. Nestoras, E. Palaiologou, A. C. S. Readhead

Abstract: 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 u… ▽ More 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. △ Less

Submitted 10 November, 2017; originally announced November 2017.

Comments: 13 pages, 4 figures, accepted for publication in APJ

RoboPol: Connection between optical polarization plane rotations and gamma-ray flares in blazars

Authors: D. Blinov, V. Pavlidou, I. Papadakis, S. Kiehlmann, I. Liodakis, G. V. Panopoulou, E. Angelakis, M. Baloković, T. Hovatta, O. G. King, A. Kus, N. Kylafis, A. Mahabal, S. Maharana, I. Myserlis, E. Paleologou, I. Papamastorakis, E. Pazderski, T. J. Pearson, A. Ramaprakash, A. C. S. Readhead, P. Reig, K. Tassis, J. A. Zensus

Abstract: We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confir… ▽ More We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time scales of rotations and flares are marginally correlated. △ Less

Submitted 24 October, 2017; originally announced October 2017.

Comments: 12 pages, 16 figures, accepted to MNRAS

Radio QPO in the $γ$-ray-loud X-ray binary LS I +61${^\circ}$303

Authors: Frederic Jaron, Richa Sharma, Maria Massi, Lars Fuhrmann, Emmanouil Angelakis, Ioannis Myserlis, Guang-Xing Li, Xun Shi

Abstract: LS I +61${^\circ}$303 is a $γ$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate thes… ▽ More LS I +61${^\circ}$303 is a $γ$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate these microflaring activity in LS I +61${^\circ}$303 we observed the source with the 100-m Effelsberg radio telescope at 4.85, 8.35, and 10.45 GHz and performed timing analysis on the obtained data. Radio oscillations of 15 hours time scales are detected at all three frequencies. We also compare the spectral index evolution of radio data to that of the photon index of GeV data observed by Fermi-LAT. We conclude that the observed QPO could result from multiple shocks in a jet. △ Less

Submitted 26 July, 2017; originally announced July 2017.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS

Full-Stokes polarimetry with circularly polarized feeds - Sources with stable linear and circular polarization in the GHz regime

Authors: I. Myserlis, E. Angelakis, A. Kraus, C. A. Liontas, N. Marchili, M. F. Aller, H. D. Aller, V. Karamanavis, L. Fuhrmann, T. P. Krichbaum, J. A. Zensus

Abstract: 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üller matrix one. The instrumental linear polarization is corrected across the whole telescope beam and significan… ▽ More 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üller 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üller 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. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 19 pages, 17 figures, accepted for publication in Astronomy & Astrophysics on May 30, 2017

Journal ref: A&A 609, A68 (2018)

Multiwavelength picture of the blazar S5 0716+714 during its brightest outburst

Authors: Marina Manganaro, Giovanna Pedaletti, Marlene Doert, Denis Bastieri, Vandad Fallah Ramazani, Dario Gasparrini, Elina Lindfors, Benoit Lott, Mireia Nievas, Bindu Rani, David J. Thompson, Emmanouil Angelakis, George Borman, Mark Gurwell, Talvikki Hovatta, Ryosuke Itoh, Svetlana Jorstad, Alex Kraus, Thomas P. Krichbaum, Paul Kuin, Anne Lähteenmäki, Valeri Larionov, Amy Yarleen Lien, Ioannis Myserlis, Merja Tornikoski , et al. (2 additional authors not shown)

Abstract: S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13σ significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to… ▽ More S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13σ significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Metsähovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst. △ Less

Submitted 11 April, 2017; originally announced April 2017.

Comments: 7 pages, 3 figures, proceeding of the JETS 2016 conference

Journal ref: Galaxies, vol. 4, issue 4, p. 69 (2016)