On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

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

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

Our team at the INAF-Astronomical Observatory of Padova and the University of Padova is engaged in the design, construction and operations of instruments with very high time accuracy in the optical band for applications to High Time Resolution Astrophysics and Quantum Astronomy. Two instruments were built to perform photon counting with sub-nanosecond temporal accuracy, Aqueye+ and Iqueye. Aqueye+ is regularly mounted at the 1.8m Copernicus telescope in Asiago, while Iqueye was mounted at several 4m class telescopes around the world and is now attached through the Iqueye Fiber Interface at the 1.2m Galileo telescope in Asiago. They are used to perform coordinated high time resolution optical observations and, for the first time ever, experiments of optical intensity interferometry on a baseline of a few kilometers. I will report on recent technological developments and scientific results obtained within the framework of this project.

We report the first detection of an optical millisecond pulsar with the fast photon counter Aqueye+ in Asiago. This is an independent confirmation of the detection of millisecond pulsations from PSR J1023+0038 obtained with SiFAP at the Telescopio Nazionale Galileo. We observed the transitional millisecond pulsar PSR J1023+0038 with Aqueye+ mounted at the Copernicus telescope in January 2018. Highly significant pulsations were detected. The rotational period is in agreement with the value extrapolated from the X-ray ephemeris, while the time of passage at the ascending node is shifted by $11.55 \pm 0.08$ s from the value predicted using the orbital period from the X-rays. An independent optical timing solution is derived over a baseline of a few days, that has an accuracy of $\sim 0.007$ in pulse phase ($\sim 12$ $\mu$s in time). This level of precision is needed to derive an accurate coherent timing solution for the pulsar and to search for possible phase shifts between the optical and X-ray pulses using future simultaneous X-ray and optical observations.

The Vela pulsar is among a number of pulsars which show detectable optical pulsations. We performed optical observations of this pulsar in January and December 2009 with the Iqueye instrument mounted at the ESO 3.5 m New Technology Telescope. Our aim was to perform phase fitting of the Iqueye data, and to measure the optical pulse profile of the Vela pulsar at high time resolution, its absolute phase and rotational period. We calculated for the first time an independent optical timing solution and obtained the most detailed optical pulse profile available to date. Iqueye detected a distinct narrow component on the top of one of the two main optical peaks, which was not resolved in previous observations, and a third statistically significant optical peak not aligned with the radio one. The quality of the Iqueye data allowed us to determine the relative time of arrival of the radio-optical-gamma-ray peaks with an accuracy of a fraction of a millisecond. We compare the shape of the Iqueye pulse profile with that observed in other energy bands and discuss its complex multi-wavelength structure.