Multiwavelength polarimetry offers a powerful tool to probe magnetic field and flow geometries in the relativistic jets of blazars. Sources with synchrotron emission that spans a broad frequency range, from radio to X-rays, such as High Synchrotron Peak (HSP) type BL Lac objects, are particularly interesting. Previous measurements including radio, optical, and X-ray data show a clear trend, with the degree of polarization increasing with frequency. Here we report radio, optical, and X-ray observations (Swift, Nustar, and IXPE) of 1ES 1101-232 and RGB J0710+591, two blazars belonging to the puzzling subclass of extreme BL Lacs (EHBL). For 1ES 1101-232, we find a strong frequency dependence of the degree of polarization, with a ratio ΠX/ΠO ≃ 5.2. For RGB J0710+591, IXPE derives a 1σ upper limit ΠX < 11.6%, comparable to the measured optical degree of polarization (average ΠO ∼ 12%). We discuss the results in the framework of current interpretations and, in particular, we report an improved version of the stratified shock model that reproduces the observed data of both sources.

Active galactic nuclei (AGNs) are powered by supermassive black holes and exhibit different observational types, primarily due to orientation effects according to the unified model of AGN. Changing-look AGNs (CL AGNs), however, challenge this paradigm by transitioning between different types overtime. A hallmark of CL AGNs is the variability of the broad emission line (BEL) flux relative to the continuum. In this study, we investigate the optical spectral evolution of Ton 599 and its correlation with the gamma-ray flux in the GeV band. The optical continuum flux density exhibits strong variability, with a normalized excess variance of F-var= 1 . 089. In contrast, the MgII flux varies with a much smaller amplitude (F-var= 0 . 194), while the equivalent width of Mg II shows intermediate variability (F-var= 0 . 686). A strong correlation is observed between the optical continuum and the gamma-ray flux, suggesting a common origin. In addition, an amplitude-modulated quasi-periodic oscillation (QPO) with a period of similar to 176 d is identified in the gamma-ray light curve during MJD 57000-58400, with modulation on a similar to 400-d time-scale. A similar QPO with a period of similar to 181 d is also detected in the optical g-band light curve during MJD 59500-60100. The presence of these QPOs in both bands likely reflect the influence of geometrical structures in the jet as well as long-term bursts. Our results suggest that jet dynamics may drive CL events in some AGNs, contributing to our understanding of the physical processes governing AGN jets and offering critical refinements to the unified model.

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.

BL Lacertae is not only archetypical of an entire class of jet-dominated active galactic nuclei, blazars, but also one of the most active and rapidly changing objects in this class. In the fall of 2024 (September–November), BL Lacertae underwent another episode of strong optical activity, reaching an R-band magnitude of about 12 and showing extremely rapid and large-amplitude inter- and intra-night flux and polarization variations. During this period, the object was monitored over 40 nights using telescopes with an aperture of up to 2 m at three observatories: Rozhen and Belogradchik in Bulgaria and Skinakas in Greece. The results from this study include some of the most spectacular intra-night variability episodes detected in a blazar. These rapid variations, combined with high photometric accuracy and high time resolution, allowed for confirmation of consistency between different optical bands with zero time delays, down to a minute scale. Unlike previous activity reports, polarization was relatively stable on these short time-scales. Possible connections between polarization, flux, and intra-night variability were explored in order to better model or constrain the physical processes and emission mechanisms in the relativistic 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.

Due to its peculiar and highly variable nature, the blazar 3C 454.3 has been extensively monitored by the WEBT team. Here, we present for the first time these long-term optical flux and color variability results using data acquired in B, V, R, and I bands over a time span of about two decades. We include data from WEBT collaborators and public archives such as SMARTS, Steward Observatory, and Zwicky Transient Facility. The data are binned and segmented to study the source over this long term when more regular sampling was available. During our study, the long-term spectral variability reveals a redder-when-brighter trend, which, however, stabilizes at a particular brightness cutoff of similar to 14.5 mag in the I band, after which it saturates and evolves into a complex state. This trend indicates increasing dominance of jet emission over accretion disk (AD) emission until jet emission completely dominates. Plots of the variation in spectral index (following F nu proportional to nu -alpha ) reveal a bimodal distribution using a one-day binning. These correlate with two extreme phases of 3C 454.3, an outburst or high-flux state and a quiescent or low-flux state, which are respectively jet- and AD-dominated. We have also conducted intraday variability studies of nine light curves and found that six of them are variable. Discrete correlation function analysis between different pairs of optical wave bands peaks at zero lags, indicating cospatial emission in different optical bands.

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.

Context. Blazars are beamed active galactic nuclei (AGNs) known for their strong multi-wavelength variability on timescales ranging from years down to minutes. Many different models have been proposed to explain this variability. Aims. We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. Methods. We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019-2022, together with radio data from the WEBT and other teams, and gamma-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and gamma-ray brightness maxima. Results. The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, and gamma-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. The gamma-ray emitting region is found to be co-spatial with the optical one, and the analysis of the gamma-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. Conclusions. We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability.

Blazars are beamed active galactic nuclei known for their strong multi-wavelength variability on timescales from years down to minutes. We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019-2022, together with radio data from the WEBT and other teams, and gamma-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and gamma-ray brightness maxima. The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, and gamma-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. The gamma-ray emitting region is found to be co-spatial with the optical one, and the analysis of the gamma-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability.

ABSTRACT We present an extensive recent multiband optical photometric observations of the blazar S5 0716+714 carried out over 53 nights with two telescopes in India, two in Bulgaria, one in Serbia, and one in Egypt during 2019 November – 2022 December. We collected 1401, 689, 14726, and 165 photometric image frames in B, V, R, and I bands, respectively. We monitored the blazar quasi-simultaneously during three nights in B, V, R, and I bands; four nights in B, V, and R; two nights in V, R, and I; five nights in B and R; and two nights in V and R bands. We also took 37 nights of data only in R band. Single band data are used to study intraday flux variability and two or more bands quasi-simultaneous observations allow us to search for colour variation in the source. We employ the power-enhanced F-test and the nested ANOVA test to search for genuine flux and colour variations in the light curves of the blazar on intraday time-scales. Out of 12, 11, 53, and 5 nights observations, intraday variations with amplitudes between ∼3 and ∼20 per cent are detected in 9, 8, 31 and 3 nights in B, V, R, and I bands, respectively, corresponding to duty cycles of 75, 73, 58, and 60 per cent. These duty cycles are lower than those typically measured at earlier times. On these time-scales colour variations with both bluer-when-brighter and redder-when-brighter are seen, though nights with no measurable colour variation are also present. We briefly discuss possible explanations for this observed intraday variability.

We report the first extensive optical flux and spectral variability study of the TeV blazar TXS0506 + 056 on intranight to long-term time-scales using BVRI data collected over 220 nights between 2017 January 21 to 2022 April 9 using eight optical ground-based telescopes. In our search for intraday variability (IDV), we have employed two statistical analysis techniques, the nested ANOVA test and the power enhanced F-test. We found the source was variable in 8 nights out of 35 in the R-band and in 2 of 14 in the V-band yielding duty cycles (DC) of 22.8 percent and 14.3 percent, respectively. Clear colour variation in V - R was seen in only 1 out of 14 observing nights, but no IDV was found in the more limited B, I, and B - I data. During our monitoring period the source showed a 1.18 mag variation in the R-band and similar variations are clearly seen at all optical wavelengths. We extracted the optical (BVRI) SEDs of the blazar for 44 nights when observations were carried out in all four of those wavebands. The mean spectral index (alpha) was determined to be 0.897 +/- 0.171.

We study the optical flux and polarization variability of the binary black hole blazar OJ 287 using quasi-simultaneous observations from 2015 to 2023 carried out using telescopes in the USA, Japan, Russia, Crimea, and Bulgaria. This is one of the most extensive quasi-simultaneous optical flux and polarization variability studies of OJ 287. OJ 287 showed large amplitude, ∼3.0 mag flux variability, large changes of ∼37% in degree of polarization, and a large swing of ∼215° in the angle of the electric vector of polarization. During the period of observation, several flares in flux were detected. Those flares are correlated with a rapid increase in the degree of polarization and swings in electric vector of polarization angle. A peculiar behavior of anticorrelation between flux and polarization degree, accompanied by a nearly constant polarization angle, was detected from JD 2,458,156 to JD 2,458,292. We briefly discuss some explanations for the flux and polarization variations observed in OJ 287.

The “blazar sequence” has been proposed for more than 20 yr, yet its nature is still unclear. In this work, for the first time we expand this topic to the TeV band by using a sample of 58 TeV blazars, including 48 blazars in the quiescent state and 21 blazars in the flaring state (11 blazars show in both quiescent and flaring state). We investigate the correlation between the TeV luminosity, which has been compensated for attenuation from extragalactic background light, and the synchrotron peak frequency. We note that there is no correlation between TeV luminosity and peak frequency in the quiescent state and a strong anticorrelation in the flaring state for the observed value. However, there is a strong positive correlation in both the quiescent state and the flaring state for the intrinsic value. This indicates that the blazar sequence is shown in the flaring state rather than in the quiescent state for the observed value and the blazar sequence is not present in both states after removing the beaming effect. In addition, to confirm whether the beaming effect results in the blazar sequence, we compare the Fermi γ -ray luminosity between the quiescent state and the flaring state. We find the Fermi γ -ray luminosity in the flaring state is greater than that in the quiescent state, and the Doppler factor in the flaring state is greater. We suggest the blazar sequence in the flaring state may be due to a stronger beaming effect.

In 2021 BL Lacertae underwent an extraordinary activity phase, which was intensively followed by the Whole Earth Blazar Telescope (WEBT) Collaboration. We present the WEBT optical data in the BVRI bands acquired at 36 observatories around the world. In mid 2021 the source showed its historical maximum, with R = 11.14. The light curves display many episodes of intraday variability, whose amplitude increases with source brightness, in agreement with a geometrical interpretation of the long-term flux behaviour. This is also supported by the long-term spectral variability, with an almost achromatic trend with brightness. In contrast, short-term variations are found to be strongly chromatic and are ascribed to energetic processes in the jet. We also analyse the optical polarimetric behaviour, finding evidence of a strong correlation between the intrinsic fast variations in flux density and those in polarisation degree, with a time delay of about 13 h. This suggests a common physical origin. The overall behaviour of the source can be interpreted as the result of two mechanisms: variability on time scales greater than several days is likely produced by orientation effects, while either shock waves propagating in the jet, or magnetic reconnection, possibly induced by kink instabilities in the jet, can explain variability on shorter time scales. The latter scenario could also account for the appearance of quasi-periodic oscillations, with periods from a few days to a few hours, during outbursts, when the jet is more closely aligned with our line of sight and the time scales are shortened by relativistic effects.

Studies and constraints on the emission region are crucial to the blazar radiation mechanism. Yet previous works have mainly focused on individual sources. In this work, we make use of the largest and the latest spectral energy distribution fitting results in the literature to statistically study the blazar emission region properties in the framework of a one-zone leptonic model. Our results reveal: (1) that flat-spectrum radio quasars (FSRQs) show lower electron energy ( γ _p ≲ 1.6 × 10 ^3 ) than BL Lacertae objects (BL Lacs) and tend to have a stronger magnetic field ( B ) and smaller electron-to-magnetic energy ratio ( U _e / U _B ) than BL Lacs; (2) we find that the electromagnetic equipartition would rather happen in the jets of BL Lacs than happen in the jets of FSRQs; (3) there are 682 blazars with a magnetic field weaker than the critical value for generating the Kelvin–Helmholtz instability, thus one-third of the blazars in our sample are able to produce this instability; and (4) the distance ( d _em ) between the emission region and the central black hole is on the scale of ∼0.1 pc, so the location of the emission region may be evenly distributed inside and outside the broad-line region.

Aims. We aim to investigate the extreme variability properties of the TeV blazar S4 0954+65 using optical photometric and polarisation observations carried out between 2017 and 2023 using three ground-based telescopes. Methods. We examined an extensive dataset comprised of 138 intraday (observing duration shorter than a day) light curves (LCs) of S4 0954+65 for flux, spectral, and polarisation variations on diverse timescales. For the variable LCs, we computed the minimum variability timescales. We investigated flux–flux correlations and colour variations to look for spectral variations on long (several weeks to years) timescales. Additionally, we looked for connections between optical R-band flux and polarisation degree. Results. We found significant variations in 59 out of 138 intraday LCs. We detected a maximum change of 0.58 ± 0.11 in V-band magnitude within ∼2.64 h and a corresponding minimum variability timescale of 18.21 ± 4.87 mins on 2017 March 25. During the course of our observing campaign, the source brightness changed by ∼4 magnitudes in V and R bands; however, we did not find any strong spectral variations. The slope of the relative spectral energy distribution was 1.37 ± 0.04. The degree of polarisation varied from ∼3% to 39% during our monitoring. We observed a change of ∼120 degrees in polarisation angle (PA) within ∼3 h on 2022 April 13. No clear correlation was found between optical flux and the degree of polarisation. Conclusions. The results of our optical flux, colour, and polarisation study provide hints that turbulence in the relativistic jet could be responsible for the intraday optical variations in the blazar S4 0954+65. However, the long-term flux variations may be caused by changes in the Doppler factor.

ABSTRACT In this work, we report the presence of rapid intra-night optical variations in both flux and polarization of the blazar BL Lacertae during its unprecedented 2020–2021 high state of brightness. The object showed significant flux variability and some colour changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both polarization degree and angle (electric vector polarization angle). The object was observed from several observatories throughout the world, covering a total of almost 300 h during 66 nights. Based on our results, we suggest that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight, can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night time-scales.

We report the first extensive optical flux and spectral variability study of the TeV blazar TXS 0506+056 on intra-night to long-term timescales using BVRI data collected over 220 nights between January 21, 2017 to April 9, 2022 using 8 optical ground-based telescopes. In our search for intraday variability (IDV), we have employed two statistical analysis techniques, the nested ANOVA test and the power enhanced F-test. We found the source was variable in 8 nights out of 35 in the R-band and in 2 of 14 in the V-band yielding Duty Cycles (DC) of 22.8 was seen in only 1 out of 14 observing nights, but no IDV was found in the more limited B, I, and B - I data. During our monitoring period the source showed a 1.18 mag variation in the R-band and similar variations are clearly seen at all optical wavelengths. We extracted the optical (BVRI) SEDs of the blazar for 44 nights when observations were carried out in all four of those wavebands. The mean spectral index (α) was determined to be 0.897+-0.171

We study the broadband emission of Mrk 501 using multiwavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi's Large Area Telescope (LAT), NuSTAR, Swift, GASP-WEBT, and the Owens Valley Radio Observatory. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wave bands, 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 the Swift X-Ray Telescope 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 yr 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 toward 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.