We report on a nine-month monitoring campaign of the blazar PG 1553+113, relying on three observations carried out in 2025 with the Imaging X-ray Polarimetry Explorer (IXPE) and supported by multiwavelength facilities. The source displayed pronounced variability across the electromagnetic spectrum, with X-ray flux changes by up to a factor of similar to 5 and complex evolution of the optical polarization properties, including one of the largest (exceeding 150 degrees) and fastest rotations in the electric vector position angle (EVPA) ever recorded. This swing of the EVPA was also accompanied by a temporary drop of the optical polarization degree to nearly zero. Significant X-ray polarization was observed during the third IXPE pointing, with a polarization degree Pi X = (18.4 +/- 5.8)% and Psi X = 74 degrees +/- 9 degrees in the 2-8 keV band, while only upper limits were obtained in the first two epochs. The optical data show that the second IXPE observation occurred shortly after a dramatic optical polarization event characterized by a rapid EVPA swing and strong depolarization. Two possible scenarios may explain the broadband polarimetric behavior: (i) the superposition of two emitting regions with nearly orthogonal magnetic field configurations and variable relative contributions, and (ii) the interaction of a single emitting region with a shock that temporarily reorders the magnetic field. In both cases, the data support a picture in which the X-ray and optical emissions arise from closely related but not strictly cospatial regions within a dynamically evolving, magnetically structured jet.

Millimeter-band polarimetry offers a powerful probe of active galactic nuclei jets, accessing regions less affected by opacity and Faraday rotation than at longer radio wavelengths. As part of the POLAMI program, we have conducted 14 yr of 1 and 3 mm polarization monitoring of a sample of gamma-ray-bright blazars with the IRAM 30 m telescope, complemented here with long-term optical polarimetric observations from multiple facilities. We aim to test whether current models of parsec-scale jet physics are consistent with observed multiband polarization behavior. Using a Bayesian framework, we derive intrinsic mean flux densities and modulation indices for total flux and fractional polarization, and characterize electric vector position angle (EVPA) variability using circular statistics. We then examine how these quantities reflecting variability properties across millimeter and optical bands relate to synchrotron peak frequency, jet orientation, and radio/gamma-ray luminosities. BL Lacertae (BL Lac) objects exhibit, on average, higher fractional polarization and lower EVPA variability than flat-spectrum radio quasars (FSRQs) at all wavelengths. Fractional polarization increases with frequency, consistent with increasingly ordered magnetic fields at shorter wavelengths. BL Lacs also show more frequent alignment of EVPAs between optical and millimeter bands, whereas FSRQs display weaker coherence. EVPA variability correlates positively with radio and gamma-ray luminosities and negatively with synchrotron peak frequency, most strongly in the optical. We further find a positive correlation between EVPA spread and fractional polarization variability, suggesting a direct link between magnetic-field structure and polarization dynamics.

We analyzed 19 yr of R-band data of the blazar 3C 454.3 from the Whole Earth Blazar Telescope archive, along with new data from its members and from public archives such as those provided by the Small and Moderate Aperture Research Telescope System and the Steward Observatory projects to search for quasiperiodic oscillations (QPOs). We detected a QPO of similar to 433 days using Lomb-Scargle periodogram, which lasted from MJD 54980-58450 as detected by the weighted wavelet Z-transform technique, making it one of the most persistent QPOs ever detected in the optical regime. The phase dispersion minimization technique was also performed to further validate this QPO claim. We detected this signal at a global significance of 2.53 sigma across all methodologies. To explain the observed QPO, we have considered both models focused on the accretion disk around a supermassive black hole (SMBH), and those based purely on jet emissions. Plausible jet-based models involve a shock moving down the jet in a helical magnetic field, whereas the SMBH models could involve Lense-Thirring-effect-induced jet precession or dual jets in a binary SMBH system. We introduce a novel approach to distinguish genuine QPOs from spurious signals arising from annual seasonal gaps, a common limitation of ground-based observations.

Context. The BL Lac object 3C 371 is one of the targets regularly monitored by the Whole-Earth Blazar Telescope (WEBT), a collaboration of observers studying blazar variability on both short and long timescales. Aims. We aim to evaluate the long-term multi-wavelength (MWL) behaviour of 3C 371, comparing it with results derived from its optical emission in our previous study. For this, we make use of the multi-band campaigns organised by the WEBT collaboration in optical and radio between January 2018 and December 2020, and of public data from Swift and Fermi satellites and the MOJAVE Very Large Interferometry programme. Methods. We evaluated the variability shown by the source in each band by quantifying the amplitude variability parameter, and also looked for a possible inter-band correlation using the z-discrete correlation function. We also present a deep analysis of the optical-UV, X-ray, and gamma-ray spectral variability. With the MOJAVE data, we performed a kinematics analysis, looking for components propagating along the jet and calculating its kinematics parameters. We then used this set of parameters to interpret the source MWL behaviour, modelling its broadband spectral energy distribution (SED) with theoretical blazar emission scenarios. Results. The MWL variability of the source in the UV, X-ray, and gamma-ray bands is comparable to that in optical, especially considering the lower coverage of the first two wavebands. On the other hand, the radio bands show variability of much lower magnitude. Moreover, this MWL emission shows a high degree of correlation, which is compatible with zero lag, again with the exception of the radio emission. The radio VLBI images reveal super-luminal motion of one of the identified components, which we used to set constraints on the jet kinematics and parameters, and to estimate a viewing angle of theta = (9.6 +/- 1.6)degrees, a Doppler factor of delta = 6.0 +/- 1.1, and a Lorentz factor of Gamma = 6.0 +/- 1.8. The polarised radio emission was found to be anti-correlated with the total flux, and to follow the same behaviour as the polarised optical radiation. The optical-UV spectral behaviour shows a mild harder-when-brighter trend on long timescales, and other trends such as redder-when-brighter on shorter timescales. We successfully modelled the broadband emission with a leptonic scenario, where we compared the low and high emission states during the period of complete MWL coverage. The difference between these two states can be ascribed mainly to a hardening of the distribution of particles. The derived features of the source confirm that 3C 371 is a BL Lac whose jet is not well aligned with the line of sight.

We report the Imaging X-ray Polarimetry Explorer (IXPE) polarimetric and simultaneous multiwavelength observations of the high-energy-peaked BL Lacertae object (HBL) 1ES 1959+650, performed in 2022 October and 2023 August. In 2022 October, IXPE measured an average polarization degree Π _X = 9.4% ± 1.6% and an electric-vector position angle ψ _X = 53° ± 5°. The polarized X-ray emission can be decomposed into a constant component, plus a rotating component, with the rotation velocity ω _EVPA = (−117 ± 12) deg day ^−1 . In 2023 August, during a period of pronounced activity of the source, IXPE measured an average Π _X = 12.4% ± 0.7% and ψ _X = 20° ± 2°, with evidence (∼0.4% chance probability) for a rapidly rotating component with ω _EVPA = 1864 ± 34 deg day ^−1 . These findings suggest the presence of a helical magnetic field in the jet of 1ES 1959+650 or stochastic processes governing the field in turbulent plasma. Our multiwavelength campaigns from radio to X-ray reveal variability in both polarization and flux from optical to X-rays. We interpret the results in terms of a relatively slowly varying component dominating the radio and optical emission, while rapidly variable polarized components dominate the X-ray and provide minor contribution at optical wavelengths. The radio and optical data indicate that on parsec scales the magnetic field is primarily orthogonal to the jet direction. On the contrary, X-ray measurements show a magnetic field almost aligned with the parsec jet direction. Confronting with other IXPE observations, we guess that the magnetic field of HBLs on subparsec scale should be rather unstable, often changing its direction with respect to the Very Long Baseline Array jet.

OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source, which includes very high energy (VHE, E > 100 GeV) gamma-ray data taken by the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov telescopes) and H.E.S.S. (High Energy Stereoscopic System) imaging Cherenkov telescopes. The discovery of VHE gamma-ray emission happened during a high state of gamma-ray activity in July 2016, observed by many instruments from radio to VHE gamma-rays. We identify four states of activity of the source, one of which includes VHE gamma-ray emission. Variability in the VHE domain is found on daily time-scales. The intrinsic VHE spectrum can be described by a power law with index 3 . 27 +/- 0.44(stat) +/- 0.15(sys) (MAGIC) and 3 . 39 +/- 0.58(stat) +/- 0.64(sys) (H.E.S.S.) in the energy range of 55-300 and 120-500 GeV, respectively. The broadband emission cannot be successfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the data set well and a proton-synchrotron-dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the broad-line region to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be a flat spectrum radio quasar (FSRQ), in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL Lac and FSRQ objects.

The X-ray polarization observations, made possible with the Imaging X-ray Polarimetry Explorer (IXPE), offer new ways of probing high-energy emission processes in astrophysical jets from blazars. Here, we report the first X-ray polarization observation of the blazar S4 0954 +65 in a high optical and X-ray state. During our multi-wavelength (MWL) campaign of the source, we detected an optical flare whose peak coincided with the peak of an X-ray flare. This optical-X-ray flare most likely took place in a feature moving along the parsec-scale jet, imaged at 43 GHz by the Very Long Baseline Array (VLBA). The 43 GHz polarization angle of the moving component underwent a rotation near the time of the flare. In the optical band, prior to the IXPE observation, we measured the polarization angle to be aligned with the jet axis. In contrast, during the optical flare, the optical polarization angle was perpendicular to the jet axis; after the flare, it reverted to being parallel to the jet axis. Due to the smooth behavior of the optical polarization angle during the flare, we favor shocks as the main acceleration mechanism. We also infer that the ambient magnetic field lines in the jet were parallel to the jet position angle. The average degree of optical polarization during the IXPE observation was (14.3 +/- 4.1)%. Despite the flare, we only detected an upper limit of 14% (at 3 sigma level) on the X-ray polarization degree; however, a reasonable assumption on the X-ray polarization angle results in an upper limit of 8.8% (3 sigma). We modeled the spectral energy distribution (SED) and spectral polarization distribution (SPD) of S4 0954 +65 with leptonic (synchrotron self-Compton) and hadronic (proton and pair synchrotron) models. Our combined MWL polarization observations and SED modeling tentatively disfavor the use of hadronic models for the X-ray emission in S4 0954 +65.

The polarimetric properties of blazars enable us to place constraints on the acceleration mechanisms that fuel their powerful jets. By studying the multiwavelength polarimetric behaviour of high-synchrotron peaked (HSP) and low-synchrotron peaked (LSP) blazars, we aim to explore differences in their emission mechanisms and magnetic field structure in the acceleration region. In this study, we take advantage of several X-ray polarisation observations of HSP by the IXPE, including four new observations of Mrk 501, along with optical polarisation observations of LSP from RoboPol and other instruments. We find that the polarisation degree (PD) distribution of HSP in X-rays is systematically higher than in optical and mm-radio wavelengths, as reported in previous IXPE publications. The distribution of the X-ray electric vector position angles (PA) is centred around the jet axis with most of the observations consistent with zero difference within uncertainties. In fact, the distribution of the offset of the PA from the jet axis is consistent between the LSP and HSP populations (with PA measured in optical for the first, X-ray for the latter), suggesting a common magnetic field structure close to the acceleration region. These results offer strong support for the emerging energy stratified scenario of particle acceleration followed by energy loss in blazar jets.

3C 84 is the brightest cluster galaxy in the Perseus Cluster. It is among the closest radio-loud active galaxies and among the very few that can be detected from low-frequency radio up to TeV γ -rays. Here we report on the first X-ray polarization observation of 3C 84 with the Imaging X-ray Polarimetry Explorer, for a total of 2.2 Ms coinciding with a flare in γ -rays. This is the longest observation for a radio-loud active galaxy, which allowed us to reach unprecedented sensitivity, leading to the detection of an X-ray polarization degree of Π _X = 4.2% ± 1.3% (∼3.2 σ confidence) at an X-ray electric vector polarization angle of ψ _X = 163° ± 9°, which is aligned with the radio jet direction on the sky. Optical polarization observations show fast variability about the jet axis as well. Our results strongly favor models in which X-rays are produced by Compton scattering from relativistic electrons—specifically synchrotron self-Compton—that takes place downstream, away from the supermassive black hole.

Polarimetric properties of blazars allow us to put constraints on the acceleration mechanisms that fuel their powerful jets. By studying the multiwavelength polarimetric behaviour of high synchrotron peaked (HSP) and low synchrotron peaked (LSP) blazars, we aim to explore differences in their emission mechanisms and magnetic field structure in the acceleration region. In this study, we take advantage of several X-ray polarisation observations of HSP by the IXPE, including four new observations of Mrk 501, and optical polarisation observations of LSP from RoboPol and many others. We find that the polarisation degree (PD) distribution of HSP in X-rays is systematically higher than in optical and mm-radio wavelengths, as reported in previous IXPE publications. The distribution of the X-ray electric vector position angles (PA) is centered around the jet axis with most of the observations consistent with zero difference within uncertainties. In fact, the distribution of the offset of the PA from the jet axis is consistent between the LSP and HSP populations (with PA measured in optical for the first, X-ray for the latter), suggesting a common magnetic field structure close to the acceleration region, in strong support of the emerging energy stratified picture of particle acceleration followed by energy loss in blazar jets.

The origin of the high-energy emission in astrophysical jets from black holes is a highly debated issue. This is particularly true for jets from supermassive black holes, which are among the most powerful particle accelerators in the Universe. So far, the addition of new observations and new messengers have only managed to create more questions than answers. However, the newly available X-ray polarization observations promise to finally distinguish between emission models. We use extensive multiwavelength and polarization campaigns as well as state-of-the-art polarized spectral energy distribution models to attack this problem by focusing on two X-ray polarization observations of blazar BL Lacertae in flaring and quiescent γ-ray states. We find that, regardless of the jet composition and underlying emission model, inverse-Compton scattering from relativistic electrons dominates at X-ray energies.

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including an ∼90° angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets.

Blazars, supermassive black hole systems with highly relativistic jets aligned with the line of sight, are the most powerful long-lived emitters of electromagnetic emission in the Universe. We report here on a radio-to-gamma-ray multiwavelength campaign on the blazar BL Lacertae with unprecedented polarimetric coverage from radio to X-ray wavelengths. The observations caught an extraordinary event on 2023 November 10–18, when the degree of linear polarization of optical synchrotron radiation reached a record value of 47.5%. In stark contrast, the Imaging X-ray Polarimetry Explorer found that the X-ray (Compton scattering or hadron-induced) emission was polarized at less than 7.4% (3 σ confidence level). We argue here that this observational result rules out a hadronic origin of the high-energy emission and strongly favors a leptonic (Compton scattering) origin, thereby breaking the degeneracy between hadronic and leptonic emission models for BL Lacertae and demonstrating the power of multiwavelength polarimetry to address this question. Furthermore, the multiwavelength flux and polarization variability, featuring an extremely prominent rise and decay of the optical polarization degree, is interpreted for the first time by the relaxation of a magnetic “spring” embedded in the newly injected plasma. This suggests that the plasma jet can maintain a predominant toroidal magnetic field component parsecs away from the central engine.

Blazars are a subclass of active galactic nuclei (AGNs) with a high optical linear polarization that originates in relativistic jets. Polarization parameters such as the degree of polarization (PD) and the electric vector position angle (EVPA) are directly related to the properties of the magnetic field in the jets. A study of the optical polarization of blazars allows conclusions to be drawn about the field geometry, its evolution, and its relation to the emission properties of the blazars. The periods of ordered changes in the electric vector position angle, so-called rotations, are of particular interest. We used a new method to determine EVPA rotations and to estimate their statistical significance with the aim to analyze long-term polarimetric observations of five blazars: OJ 287, S5 0716+71, 3C 454.3, CTA 102, and PG 1553+113. This resultes in the identification of 256 EVPA rotations. We found possible tendencies for the EVPA rotations to occur in a preferred direction in each of these sources: clockwise for OJ 287 and CTA 102, and counterclockwise for the others. The EVPA rotations can be explained by the spiral structure of the magnetic field in the jet. In this case, the observed preferred direction of rotations reflects the global structure of the magnetic field, which can be associated with the direction of rotation of either the black hole ergosphere or the accretion disk.

Observations of quasars show that the polarization position angle of the emission coming from them varies greatly over time, including periods called rotations during which the angle changes in an orderly manner. The study proposes a method for identifying such events and assessing their statistical significance. The operation of the method is demonstrated using the example of long-term polarimetric observations of the blazars CTA 102, 3C 454.3, and OT 081. During the analysis of light curves, 51 rotations of the polarization position angle were found and it was shown that for CTA 102 and 3C 454.3 the rotations are predominantly oriented in one direction.

The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $\nu_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $\nu_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $\nu_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.

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

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

PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.

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.