We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F>1TeV < 1.7× 10−12 ph cm−2 s−1) and in the 2-10 keV (X-ray) band (F2 − 10 keV < 3.6 × 10−11 erg cm−2 s−1), during which the Swift-BAT data suggests an additional spectral component beyond the regular synchrotron emission.The highest flux variability occurs in X-rays and very-high-energy (E>0.1 TeV) γ-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HRkeV and HRTeV show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 days centered at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 days centered at a time lag of $43^{+9}_{-6}$ days.This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a LogNormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process.

Markarian 421 and Markarian 501 are Active Galactic Nuclei and two of the most prominent sources of very high energetic gamma rays. In contrast to the predominantly stable Crab Nebula, the two blazars have extremely variable gamma-ray fluxes. Those could be caused by erratic accretion flow instabilities, but also by processes leading to quasi periodic oscillations (QPOs). Recently, hints of such a QPO of Mrk 501 with a period of about 322 days have been observed with the Large Area Telescope on-board of the Fermi satellite at GeV energies featuring global and local significances of the signal above 99%. The First G-APD Cherenkov Telescope (FACT) has monitored a small sample of known TeV gamma-ray sources including Mrk 421 and Mrk 501 over the last seven years. Therefore, it is ideally suited for periodicity studies of Active Galactic Nuclei, as they require unbiased and continuous long-term observations. Based on measurements of the well-known gamma-ray flux originating from the Crab Nebula, corrections to the measured light curves are applied. In total, 2500 hours of data for each of the two blazars have been analysed using the Lomb-Scargle algorithm. The resulting periodogram of Mrk 501 shows a hint of QPO with a similar period as the one observed by Fermi/LAT.

Aims. Radio-to-TeV observations of the bright nearby (z = 0.034) blazar Markarian 501 (Mrk 501), performed from December 2012 to April 2018, are used to study the emission mechanisms in its relativistic jet. Methods. We examined the multi-wavelength variability and the correlations of the light curves obtained by eight different instruments, including the First G-APD Cherenkov Telescope (FACT), observing Mrk 501 in very high-energy (VHE) gamma-rays at TeV energies. We identified individual TeV and X-ray flares and found a sub-day lag between variability in these two bands. Results. Simultaneous TeV and X-ray variations with almost zero lag are consistent with synchrotron self-Compton (SSC) emission, where TeV photons are produced through inverse Compton scattering. The characteristic time interval of 5−25 days between TeV flares is consistent with them being driven by Lense-Thirring precession.

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.

Aims. We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with a special focus on the multi-band flux correlations. Methods. The dataset has been collected through an extensive multi-wavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina, and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. Results. The very-high-energy (VHE; E >  100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance for a few days without a simultaneous increase in the X-ray flux (i.e., orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. The peak frequency of the synchrotron bump varies by two orders of magnitude throughout the campaign. Our multi-band correlation study also hints at an anti-correlation between UV-optical and X-ray at a significance higher than 3σ. A VHE flare observed on MJD 57788 (2017 February 4) shows gamma-ray variability on multi-hour timescales, with a factor ten increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors, from γ′min=2×104 to γ′max=6×105.