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.

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 γ-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and γ-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 γ-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 γ-ray emitting region is found to be co-spatial with the optical one, and the analysis of the γ-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.

We present the results of our study of cross-correlations between long-term multi-band observations of the radio variability of the blazar 3C 279. More than a decade (2008−2022) of radio data were collected at seven different frequencies ranging from 2 GHz to 230 GHz. The multi-band radio light curves show variations in flux, with the prominent flare features appearing first at higher-frequency and later in lower-frequency bands. This behavior is quantified by cross-correlation analysis, which finds that the emission at lower-frequency bands lags that at higher-frequency bands. Lag versus frequency plots are well fit by straight lines with negative slope, typically ∼−30 day GHz⁻¹. We discuss these flux variations in conjunction with the evolution of bright moving knots seen in multi-epoch VLBA maps to suggest possible physical changes in the jet that can explain the observational results. Some of the variations are consistent with the predictions of shock models, while others are better explained by a changing Doppler beaming factor as the knot trajectory bends slightly, given a small viewing angle to the jet.

We report the first detection at 4.7 and 8.2 GHz with the RATAN-600 radio telescope of the source PSO J047.4478+27.2992, which has been announced as a blazar at z = 6.1. The average flux densities are 12 ± 3 and 8 ± 4 mJy at 4.7 and 8.2 GHz respectively, and an upper limit is estimated as 3 mJy at 11.2 GHz. The light curve at 4.7 GHz, obtained on a time-scale of four months, exhibits moderate variability of the source (fractional variability Fvar = 0.28 ± 0.02). The new RATAN-600 measurements extend previous literature data to higher frequencies, indicating a single power-law radio spectrum with α0.147–8.2 = −0.51 ± 0.1. The monochromatic radio luminosity at 4.7 GHz Pν ∼ 2 × 1027 W Hz−1 is consistent with values for high-redshift quasars at z ≥ 3.

We present a study of the radio variability of bright, $S_{1.4}\ge 100$ mJy, high-redshift quasars at $z\ge 3$ on timescales of up to 30–40 yrs. The study involved simultaneous RATAN-600 measurements at the frequencies of 2.3, 4.7, 8.2, 11.2, and 22.3 GHz in 2017–2020. In addition, data from the literature were used. We have found that the variability index, $V_S$, which quantifies the normalized difference between the maximum and minimum flux density while accounting for measurement uncertainties, ranges from 0.02 to 0.96 for the quasars. Approximately half of the objects in the sample exhibit a variability index within the range from 0.25 to 0.50, which is comparable to that observed in blazars at lower redshifts. The distribution of $V_S$ at 22.3 GHz is significantly different from that at 2.3–11.2 GHz, which may be attributed to the fact that a compact AGN core dominates at the source’s rest frame frequencies greater than 45 GHz, leading to higher variability indices obtained at 22.3 GHz (the $V_S$ distribution peaks around 0.4) compared to the lower frequencies (the $V_S$ distribution at 2.3 and 4.7 GHz peaks around 0.1–0.2). Several source groups with distinctive variability characteristics were found using the cluster analysis of quasars. We propose seven new candidates for gigahertz-peaked spectrum (GPS) sources and five new megahertz-peaked spectrum (MPS) sources based on their spectrum shape and variability features. Only 6 out of the 23 sources previously reported as GPS demonstrate a low variability level typical of classical GPS sources ($V_S<0.25$) at 4.7–22.3 GHz. When excluding the highly variable peaked-spectrum blazars, we expect no more than 20% of the sources in the sample to be GPS candidates and no more than 10% to be MPS candidates.

We present a study of the multiwavelength (MW) variability of the blazar AO 0235$+$164 based on the radio-to-$\gamma$-ray data covering a long time period from 1997 to 2023. The radio data are represented by the 1–22 GHz measurements from the SAO RAS RATAN-600 radio telescope, the 5 and 8 GHz data from the IAA RAS RT-32 telescopes, and the 37 GHz data from the RT-22 telescope of CrAO RAS. The optical measurements in the R-band were collected with the SAO RAS 1-m Zeiss-1000 and 0.5-m AS-500/2 telescopes. Additionally, we used the archive data at 230 GHz from the Submillimetre Array and the $\gamma$-ray data in the 0.1–100 GeV band from the Fermi-LAT point source 4FGL-DR2 catalogue. The variability properties during four epochs containing major flares and one epoch of relatively low activity were analysed using the fractional variability indices, discrete correlation functions, Lomb–Scargle periodograms, and structure functions. A significant correlation ($\ge \!2\sigma$) between the radio, optical, and $\gamma$-ray bands is found for all these periods with time delays from 0 to 1.7 yr. The relation between time delay and frequency is described by a linear law with a negative slope of $-10$ d GHz$^{-1}$. The discovered properties of MW variability for the low-activity period and for flaring states suggest that the mechanisms dominating the radio–$\gamma$-ray variations are not substantially different. The detected quasi-periodic oscillations of about 6 and 2 yr are tentative, as the time span of the observations includes fewer than four full cycles for the radio and optical data and only about three cycles for the Fermi-LAT data. These results should be interpreted with caution, given the limited number of observed cycles and the influence of red noise. We used cluster analysis to reliably separate the high and low-activity states and determined statistical differences in the main properties of AO 0235$+$164 non-thermal emission. The physical parameters of the radio jet were obtained using the Hedgehog model applied to the average radio spectrum of AO 0235$+$164 in the range 0.1–300 GHz. The effectiveness of replacing electrons with protons in the synchrotron radio emission of relativistic jets is shown for describing the nature of blazars and the generation of high-energy neutrinos.