Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use the longest and densest ongoing unbiased observing campaign obtained at TeV and GeV energies during 5.5 years with the FACT telescope and the Fermi-LAT detector. The contemporaneous multi-wavelength observations were used to characterize the variability of the source and to constrain the underlying physical mechanisms. We study and correlate light curves obtained by nine different instruments from radio to gamma rays and found two significant results. The TeV and X-ray light curves are very well correlated with lag, if any, shorter than a day. The GeV light curve varies independently and accurately leads the variations observed at long wavelengths, in particular in the radio band. We find that the observations match the predictions of leptonic models and suggest that the physical conditions vary along the jet, when the emitting region moves outwards.

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use multi-wavelength light curves of Mrk 421 spanning 5.5 years with FACT (TeV) and Fermi LAT (GeV) in the gamma rays, Swift BAT, Swift XRT and MAXI in the X-rays, together with optical and radio data and investigate the physical processes driving the emission and variability. Observations by FACT are continuous and not triggered, so the source was found in a wide range of flux states and more than 30 flares were identified from X-rays to TeV. The light curves in TeV and X-rays feature very similar flares with rise and decay times of a few days and zero lag, characteristic for electron processes. At least two parameters per flare, the amplitude and the cut-off energy, are required to explain the observed variability. In addition, the GeV light curve leads and is strongly correlated with the optical and radio light curves as expected from SSC emitting shock propagating in a conical jet.