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.

We have analysed all the X-ray images centred on Gamma Ray Bursts generated by Swift over the last 15 years using automatic tools that do not require any expertise in X-ray astronomy, producing results in excellent agreement with previous findings. This work, besides presenting the largest medium-deep survey of the X-ray sky and a complete sample of blazars, wishes to be a step in the direction of achieving the ultimate goal of the Open Universe Initiative, that is to enable non expert people to fully benefit of space science data, possibly extending the potential for scientific discovery, currently confined within a small number of highly specialised teams, to a much larger population. We have used the Swift_deepsky Docker container encapsulated pipeline to build the largest existing flux-limited and unbiased sample of serendipitous X-ray sources. Swift_deepsky runs on any laptop or desktop computer with a modern operating system. The tool automatically downloads the data and the calibration files from the archives, runs the official Swift analysis software and produces a number of results including images, the list of detected sources, X-ray fluxes, SED data, and spectral slope estimations. We used our source list to build the LogN-LogS of extra-galactic sources, which perfectly matches that estimated by other satellites. Combining our survey with multi-frequency data we selected a complete radio flux-density limited sample of High Energy Peaked (HBL) blazars.

The origin of the multiwavelength emission from the high-synchrotron-peaked BL Lac 1ES 1218+304 is studied using the data from {\it Swift UVOT/XRT}, {\it NuSTAR} and {\it Fermi-LAT}. A detailed temporal and spectral analysis of the data observed during 2008-2020 in the $\gamma$-ray ($>100$ MeV), X-ray (0.3-70 keV), and optical/UV bands is performed. The $\gamma$-ray spectrum is hard with a photon index of $1.71\pm0.02$ above $100$ MeV. The {\it Swift UVOT/XRT} data show a flux increase in the UV/optical and X-ray bands; the highest $0.3-3$ keV X-ray flux was $(1.13\pm0.02)\times10^{-10}{\rm erg\:cm^{-2}\:s^{-1}}$. In the 0.3-10 keV range the averaged X-ray photon index is $>2.0$ which softens to $2.56 \pm 0.028$ in the 3-50 keV band. However, in some periods, the X-ray photon index became extremely hard ($<1.8$), indicating that the peak of the synchrotron component was above $1$ keV, and so 1ES 1218+304 behaved like an extreme synchrotron BL Lac. The hardest X-ray photon index of 1ES 1218+304 was $1.60 \pm 0.05 $ on MJD 58489. The time-averaged multiwavelength spectral energy distribution is modeled within a one-zone synchrotron self-Compton leptonic model using a broken power-law and power-law with an exponential cutoff electron energy distributions. The data are well explained when the electron energy distribution is $E_{\rm e}^{-2.1}$ extending up to $\gamma_{\rm br/cut}\simeq(1.7-4.3)\times10^{5}$, and the magnetic field is weak ($B\sim1.5\times10^{-2}$ G). By solving the kinetic equation for electron evolution in the emitting region, the obtained electron energy distributions are discussed considering particle injection, cooling, and escape.

Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

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.

Aims: $\gamma$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using new software, we derive instrument response and background models, enabling us to study the diffuse emission in the region. We use pre-existing data of the gas distribution in the GC region to derive the underlying distribution of CRs. Results: We obtain a significant detection for all four model components used to fit our data (Sgr~A*, ``Arc'', G0.9+0.1, and an extended component for the Galactic Ridge). We find that the diffuse component is best described as a power-law with index 2 and an exponential cut-off at around 20~TeV with the significance of the cut-off being only 2~$\sigma$. The derived cosmic-ray profile hints to a peak at the GC position, with a measured profile index of $1.2 \pm 0.3$, supporting the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity.

We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematics of Triangulum II are affected by large uncertainties leading to a bias in the determination of the properties of its dark matter halo. Using a scaling relation between the annihilation J-factor and heliocentric distance of well-known dwarf spheroidal galaxies, we estimate an annihilation J-factor for Triangulum II for WIMP dark matter of $\log[J_{\text{ann}}({0.5^{\circ}})/$ GeV$^{2}$ cm$^{-5}] = 19.35 \pm 0.37$. We also derive a dark matter density profile for the object relying on results from resolved simulations of Milky Way sized dark matter halos. We obtain 95% confidence-level limits on the thermally averaged annihilation cross section for WIMP annihilation into various Standard Model channels. The most stringent limits are obtained in the $\tau^{+}\tau^{-}$ final state, where a cross section for annihilation down to $\langle \sigma_{\text{ann}} v \rangle = 3.05 \times 10^{-24}$ cm$^{3}$ s$^{-1}$ is excluded.

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E >100 GeV) gamma-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ~ 3 times the Crab Nebula flux was measured by MAGIC on the 13th, 14th of June and 1st July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to high values of the Doppler factors (>=30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L_{jet}~10^{46} erg/s and under high values of the magnetic field strength (~100 G) and maximum proton energy (~few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016.

M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3\sigma$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region.

Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording of optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements. Results: We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~56 hr of observations, we detect the gamma-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV gamma-ray signal.

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.

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.

The flat spectrum radio quasar CTA 102 entered an extended period of activity from 2016 to 2017 during which several strong $\gamma$-ray flares were observed. Using Fermi large area telescope data a detailed investigation of \gray spectra of CTA 102 during the flaring period is performed. In several periods the \gray spectrum is not consistent with a simple power-law, having a hard photon index with an index of $\sim(1.8-2.0)$ that shows a spectral cutoff around an observed photon energy of $\sim(9-16)$ GeV. The internal $\gamma$-ray absorption via photon-photon pair production on the broad-line-region-reflected photons cannot account for the observed cut-off/break even if the emitting region is very close to the central source. This cut-off/break is likely due to a similar intrinsic break in the energy distribution of emitting particles. The origin of the spectral break is investigated through the multiwavelength modeling of the spectral energy distribution, considering a different location for the emitting region. The observed X-ray and $\gamma$-ray data is modeled as inverse Compton scattering of synchrotron and/or external photons on the electron population that produce the radio-to-optical emission which allowed to constrain the power-law index and cut-off energy in the electron energy distribution. The obtained results are discussed in the context of a diffusive acceleration of electrons in the CTA 102 jet.

We propose that the "inner engine" of a type I binary-driven hypernova (BdHN) is composed of a Kerr black hole (BH) in a non-stationary state, embedded in a uniform magnetic field $B_0$ aligned with the BH rotation axis, and surrounded by an ionized plasma of extremely low density of $10^{-14}$~g~cm$^{-3}$. Using GRB 130427A as a prototype we show that this "inner engine" acts in a sequence of "elementary impulses". Electrons are accelerated to ultra-relativistic energy near the BH horizon and, propagating along the polar axis, $\theta =0$, they can reach energies of $\sim 10^{18}$ eV, and partially contribute to ultra-high energy cosmic rays (UHECRs). When propagating with ${\theta \neq 0}$ through the magnetic field $B_0$ they give origin by synchrotron emission to GeV and TeV radiation. The mass of BH, $M=2.3 M_\odot$, its spin, $\alpha = 0.47$, and the value of magnetic field $B_0= 3.48 \times 10^{10}$ G, are determined self-consistently in order to fulfill the energetic and the transparency requirement. The repetition time of each elementary impulse of energy ${\cal E} \sim 10^{37}$ erg, is $\sim 10^{-14}$ s at the beginning of the process, then slowly increasing with time evolution. In principle, this "\textit{inner engine}" can operate in a GRB for thousands of years. By scaling the BH mass and the magnetic field the same "inner engine" can describe active galactic nuclei (AGN).

Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multi-wavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within a MWL context, including an analysis of the Fermi-LAT data over ten years of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED.

Extreme high-frequency peaked BL Lac objects (EHBLs) are blazars which exhibit extremely energetic synchrotron emission. They also feature non-thermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1TeV: this is the case of hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, ten targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 262 h of good quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all the sources but two are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma rays, three new sources were detected and a hint of signal was found for another new source. In each case the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone purely leptonic model. The VHE gamma-ray detected sources were also interpreted with a spine-layer and a proton synchrotron models. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents a first mini-catalog of VHE gamma-ray and multi-wavelength observations of EHBLs.

The gamma-ray blazar TXS 0506+056, was discovered in VHE gamma-rays by the MAGIC telescopes in 2017 in a follow-up campaign of a high energy neutrino event IceCube-170922A (IC+Fermi+MAGIC++, Science 361, eaat1378 (2018)). Subsequent multivawelenght (MWL) observations and theoretical modeling in a frame of hadro-leptonic emission confirmed that this source could be a potential cosmic ray and neutrino emitter (MAGIC Collaboration, Ansoldi et al., (2018)). This is, by far, the most significant association between a high-energy neutrino and an astrophysical source emitting gamma rays and X-rays. TXS 0506+056 is a key object to help the astrophysics community to establish connections between high-energy neutrinos and astrophysical sources. Accurate and contemporaneous MWL spectral measurements are essential ingredients to achieve this goal. In the conference, we present the measurements from the MAGIC and MWL monitoring of this source, spanning the time period from November 2017 till February 2019. These include the lowest VHE gamma-ray emission state measured from this source so far as well as a flaring episode in December 2018.

Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-13} cm^{-2} s ^{-1} . This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr^{-1}.

Open Universe for blazars is a set of high-transparency data products for blazar science, and the tools designed to generate them. Blazar astrophysics is becoming increasingly data driven, depending on the integration and combined analysis of large quantities of data from the entire span of observational astrophysics techniques. The project was therefore chosen as one of the pilot activities within the United Nations Open Universe Initiative. In this work we developed a data analysis pipeline called Swift_deepsky, based on the Swift XRTDAS software and the XIMAGE package, encapsulated into a Docker container. Swift_deepsky, downloads and reads low-level data, generates higher-level products, detects X-ray sources and estimates several intensity and spectral parameters for each detection, thus facilitating the generation of complete and up-to-date science-ready catalogues from an entire space-mission dataset. The Docker version of the pipeline and its derived products is publicly available from the Open Universe Website at openuniverse.asi.it. We present the results of a detailed X-ray image analysis based on Swift_deepsky on all Swift XRT observations including a known blazar, carried out during the first 14 years of operations of the Swift Observatory. The resulting database includes over 27,000 images integrated in different X-ray bands, and a catalogue, called 1OUSXB, that provides intensity and spectral information for 33,396 X-ray sources, 8,896 of which are single or multiple detections of 2,308 distinct blazars. All the results can be accessed on-line in a variety of ways: e.g., from the Open Universe portal at openuniverse.asi.it, through Virtual Observatory services, via the VOU-Blazar tool and the SSDC SED builder. One of the most innovative aspects of this work is that the results can be safely reproduced and extended by anyone.

We motivate how the most recent progress in the understanding the nature of the GeV radiation in most energetic gamma-ray bursts (GRBs), the binary-driven hypernovae (BdHNe), has led to the solution of a forty years unsolved problem in relativistic astrophysics: how to extract the rotational energy from a Kerr black hole for powering synchrotron emission and ultra high-energy cosmic rays. The "inner engine" is identified in the proper use of a classical solution introduced by Wald in 1974 duly extended to the most extreme conditions found around the newborn black hole in a BdHN. The energy extraction process occurs in a sequence impulsive processes each accelerating protons to $10^{21}$ eV in a timescale of $10^{-6}$ s and in presence of an external magnetic field of $10^{14}$ G. Specific example is given for a black hole of initial angular momentum $J=0.3\,M^2$ and mass $M\approx 3\,M_\odot$ leading to the GeV radiation of $10^{49}$ erg$\cdot$s$^{-1}$. The process can energetically continue for thousands of years.

Following Fermi and NOT observations, Ruffini et al. (2019b) soon identified GRB 190114C as BdHN I at z=0.424, it has been observed since, with unprecedented accuracy, [...] all the way to the successful optical observation of our predicted supernova (SN). This GRB is a twin of GRB 130427A. Here we take advantage of the GBM data and identify in it three different Episodes. Episode 1 represents the precursor which includes the SN breakout and the creation of the new neutron star ($\nu$NS), the hypercritical accretion of the SN ejecta onto the NS binary companion, exceeding the NS critical mass at $t_{rf}=1.9$s. Episode 2 starting at $t_{rf}=1.9$s includes three major events: the formation of the BH, the onset of the GeV emission and the onset of the ultra-relativistic prompt emission (UPE), which extends all the way up to $t_{rf}=3.99$s. Episode 3 which occurs at times following $t_{rf}=3.99$s reveals the presence of a cavity carved out in the SN ejecta by the BH formation. We perform an in depth time-resolved spectral analysis on the entire UPE with the corresponding determination of the spectra best fit by a cut-off power-law and a black body (CPL+BB) model, and then we repeat the spectral analysis in 5 successive time iterations in increasingly shorter time bins: we find a similarity in the spectra in each stage of the iteration revealing clearly a self-similar structure. We find a power-law dependence of the BB temperature with index $-1.56\pm0.38$, a dependence with index $-1.20\pm0.26$ for the gamma-ray luminosity confirming a similar dependence with index $-1.20\pm0.36$ which we find as well in the GeV luminosity, both expressed in the rest-frame. We thus discover in the realm of relativistic astrophysics the existence of a self-similar physical process and power-law dependencies, extensively described in the micro-physical world by the classical works of Heisenberg-Landau-Wilson.

We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03 to 0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the universe, and the fraction of it which was absorbed by dust in galaxies and re-emitted at longer wavelengths. The EBL can be studied indirectly through its effect on very-high energy photons that are emitted by cosmic sources and absorbed via photon-photon interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1-sigma upper bounds, including systematic uncertainties, are between 13% and 23% above the nominal EBL density in the models. No anomaly in the expected transparency of the universe to gamma rays is observed in any range of optical depth.We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18 - 0.62 $\mu$m range relative to the studied models, yet compatible with them within systematics.

It has been recently proposed that the ejected matter from white dwarf (WD) binary mergers can produce transient, optical and infrared emission similar to the "kilonovae" of neutron star (NS) binary mergers. To confirm this we calculate the electromagnetic emission from WD-WD mergers and compare with kilonova observations. We simulate WD-WD mergers leading to a massive, fast rotating, highly magnetized WD with an adapted version of the smoothed-particle-hydrodynamics (SPH) code Phantom. We thus obtain initial conditions for the ejecta such as escape velocity, mass and initial position and distribution. The subsequent thermal and dynamical evolution of the ejecta is obtained by integrating the energy-conservation equation accounting for expansion cooling and a heating source given by the fallback accretion onto the newly-formed WD and its magneto-dipole radiation. We show that magnetospheric processes in the merger can lead to a prompt, short gamma-ray emission of up to $\approx 10^{46}$ erg in a timescale of $0.1$-$1$ s. The bulk of the ejecta initially expands non-relativistically with velocity $0.01 c$ and then it accelerates to $0.1 c$ due to the injection of fallback accretion energy. The ejecta become transparent at optical wavelengths around $\sim 7$ days post-merger with a luminosity $10^{41}$-$10^{42}$ erg s$^{-1}$. The X-ray emission from the fallback accretion becomes visible around $\sim 150$-$200$ day post-merger with a luminosity of $10^{39}$ erg s$^{-1}$. We also predict the post-merger time at which the central WD should appear as a pulsar depending on the value of the magnetic field and rotation period.

A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $\gamma$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $\gamma$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the globular cluster M15 with the MAGIC telescopes and used 165 hrs in order to search for $\gamma$-ray emission. A strong upper limit on the TeV $\gamma$-ray flux $<3.2\times 10^{-13}\mathrm{cm^{-2}s^{-1}}$ above 300 GeV ($<0.26\%$ of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, $\eta_{\rm e}$, from the MSPs magnetospheres and their surrounding. We conclude that $\eta_{\rm e}$ must be lower than expected from the modelling of high energy processes in MSP inner magnetospheres. For leptons accelerated with the power law spectrum in the MSP wind regions, $\eta_{\rm e}$ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the globular cluster very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars.

The origin of the multiwavelength emission from PKS 0502+049 neighboring the first cosmic neutrino source TXS 0506+056 is studied using the data observed by Fermi-LAT and Swift UVOT/XRT. This source was in a flaring state in the considered bands before and after the neutrino observations in 2014-2015, characterized by hard emission spectra in the X-ray and $\gamma$-ray bands, $1.5-1.8$ and $\leq2.0$, respectively. During the neutrino observations, the $\gamma$-ray spectrum shows a deviation from a simple power-law shape, indicating a spectral cutoff at $E_c =8.50\pm2.06$ GeV. The spectral energy distributions of PKS 0502+049 are modeled within a one-zone leptonic scenario assuming that high energy $\gamma$-ray emission is produced either by IC scattering of synchrotron or dusty torus photons by the electron population that produce the radio-to-optical emission. Alternatively, the observed $\gamma$-rays are modeled considering inelastic interaction of protons, when the jet interacts with a dense gaseous target. During the neutrino observations, the $\gamma$-ray data are best described when the proton energy distribution is $E_p^{-2.61}$ and if the protons are effectively accelerated up to 10 PeV, the expected neutrino rate is $\sim1.1$ events within 110 days. In principle, if the $\gamma$-ray emission with a hard photon index observed during the flaring periods extends up to TeV, the expected rate can be somewhat higher, but such conditions are hardly possible. Within the hadronic interpretation, the $\gamma$-ray data can be reproduced only when the accretion rate of PKS 0502+049 is in the supper-Eddington regime, as opposed to the leptonic scenario. From the point of view of the necessary energetics as well as considering that the required parameters are physically reasonable, when the neutrinos were observed, the broadband emission from PKS 0502+049 is most likely of a leptonic origin.

We describe the afterglows of the long gamma-ray-burst (GRB) 130427A within the context of a binary-driven hypernova (BdHN). The afterglows originate from the interaction between a newly born neutron star ($\nu$NS), created by an Ic supernova (SN), and a mildly relativistic ejecta of a hypernova (HN). Such a HN in turn results from the impact of the GRB on the original SN Ic. The mildly relativistic expansion velocity of the afterglow ($\Gamma \sim 3$) is determined, using our model independent approach, from the thermal emission between $196$~s and $461$~s. The power-law in the optical and X-ray bands of the afterglow is shown to arise from the synchrotron emission of relativistic electrons in the expanding magnetized HN ejecta. Two components contribute to the injected energy: the kinetic energy of the mildly relativistic expanding HN and the rotational energy of the fast rotating highly magnetized $\nu$NS. We reproduce the afterglow in all wavelengths from the optical ($10^{14}$~Hz) to the X-ray band ($10^{19}$~Hz) over times from $604$~s to $5.18\times 10^6$~s relative to the Fermi-GBM trigger. Initially, the emission is dominated by the loss of kinetic energy of the HN component. After $10^5$~s the emission is dominated by the loss of rotational energy of the $\nu$NS, for which we adopt an initial rotation period of $2$~ms and a dipole plus quadrupole magnetic field of $\lesssim \! 7\times 10^{12}$~G or $\sim \! 10^{14}$~G. This scenario with a progenitor composed of a CO$_{\rm core}$ and a NS companion differs from the traditional ultra-relativistic-jetted treatments of the afterglows originating from a single black hole.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

The observation of IceCube-170922A event from the direction of TXS 0506+056 when it was in its enhanced $\gamma$-ray emission state offers a unique opportunity to investigate the lepto-hadronic processes in blazar jets. Here, the observed broadband emission of TXS 0506+056 is explained by boosted synchrotron/synchrotron self Compton emission from the jet whereas the $\gamma$-ray data observed during the neutrino emission- by inelastic interactions of the jet-accelerated protons in a dense gaseous target. The proton energy distribution is $\sim E_{p}^{-2.50}$, calculated straightforwardly from the data obtained by Fermi-LAT and MAGIC and if such distribution continues up to $E_{c,p}=10$ PeV, the expected neutrino rate is as high as $\sim0.46$ events during the long active phase of the source or $\sim0.15$ if the activity lasts 60 days. In this interpretation, the energy content of the protons above $>$ GeV in blazar jets can be estimated as well: the required proton injection luminosity is $\simeq2.0\times10^{48}\:{\rm erg\:s^{-1}}$ exceeding $10^3$ times that of electrons $\simeq10^{45}\:{\rm erg\:s^{-1}}$ which are in equipartition with the magnetic field. As the required parameters are physically realistic, this can be an acceptable model for explanation of the neutrino and $\gamma$-ray emission from TXS 0506+056.

We present the dissection in space, time, and energy of the region around the IceCube-170922A neutrino alert. This study is motivated by: (1) the first association between a neutrino alert and a blazar in a flaring state, TXS 0506+056; (2) the evidence of a neutrino flaring activity during 2014 - 2015 from the same direction; (3) the lack of an accompanying simultaneous $\gamma$-ray enhancement from the same counterpart; (4) the contrasting flaring activity of a neighbouring bright $\gamma$-ray source, the blazar PKS 0502+049, during 2014 - 2015. Our study makes use of multi-wavelength archival data accessed through Open Universe tools and includes a new analysis of Fermi-LAT data. We find that PKS 0502+049 contaminates the $\gamma$-ray emission region at low energies but TXS 0506+056 dominates the sky above a few GeV. TXS 0506+056, which is a very strong (top percent) radio and $\gamma$-ray source, is in a high $\gamma$-ray state during the neutrino alert but in a low though hard $\gamma$-ray state in coincidence with the neutrino flare. Both states can be reconciled with the energy associated with the neutrino emission and, in particular during the low/hard state, there is evidence that TXS 0506+056 has undergone a hadronic flare with very important implications for blazar modelling. All multi-messenger diagnostics reported here support a single coherent picture in which TXS 0506+056, a very high energy $\gamma$-ray blazar, is the only counterpart of all the neutrino emissions in the region and therefore the most plausible first non-stellar neutrino and, hence, cosmic ray source.

We report on broadband observations of CTA 102 ($z=1.037$) during the active states in 2016-2017. In the $\gamma$-ray band, Fermi LAT observed several prominent flares which followed a harder-when-brighter behavior: the hardest photon index $\Gamma=1.61\pm 0.10$ being unusual for FSRQs. The peak $\gamma$-ray flux above 100 MeV $(3.55\pm0.55)\times10^{-5}\:{\rm photon\:cm^{-2}\:s^{-1}}$ observed on MJD 57738.47 within 4.31 minutes, corresponds to an isotropic $\gamma$-ray luminosity of $L_{\gamma}=3.25\times10^{50}\:{\rm erg\:s^{-1}}$, comparable with the highest values observed from blazars so far. The analyses of the Swift UVOT/XRT data show an increase in the UV/optical and X-ray bands which is contemporaneous with the bright $\gamma$-ray periods. The X-ray spectrum observed by Swift XRT and NuSTAR during the $\gamma$-ray flaring period is characterized by a hard photon index of $\sim1.30$. The shortest e-folding time was $4.08\pm1.44$ hours, suggesting a very compact emission region $R\leq\delta\times2.16\times10^{14}$ cm. We modeled the spectral energy distribution of CTA 102 in several periods (having different properties in UV/optical, X-ray and $\gamma$-ray bands) assuming a compact blob inside and outside the BLR. We found that the high-energy data are better described when the infrared thermal radiation of the dusty torus is considered. In the flaring periods when the correlation between the $\gamma$-ray and UV/optical/X-ray bands is lacking, the $\gamma$-ray emission can be produced from the interaction of fresh electrons in a different blob, which does not make a dominant contribution at lower energies.

We report on a detailed investigation of the $\gamma$-ray emission from 26 non-blazar AGNs based on the Fermi LAT data accumulated for 7 years. The photon index of non-blazar AGNs changes in the range of 1.84-2.86 and the flux varies from a few times $10^{-9} photon\: cm^{-2} s^{-1}$ to $10^{-7} photon\: cm^{-2}s^{-1}$. Over long time periods, power-law provides an adequate description of the $\gamma$-ray spectra of almost all sources. Significant curvature is observed in the $\gamma$-ray spectra of NGC 1275, NGC 6251, SBS 0846+513 and PMN J0948+0022 and their spectra are better described by log-parabola or power-law with exponential cut-off models. The $\gamma$-ray spectra of PKS 0625-25 and 3C 380 show a possible deviation from a simple power-law shape, indicating a spectral cutoff around the observed photon energy of $E_{cut}=131.2\pm88.04$ GeV and $E_{cut}=55.57\pm50.74$ GeV, respectively. Our analysis confirms the previous finding of an unusual spectral turnover in the $\gamma$-ray spectrum of Cen A: the photon index changes from $2.75\pm0.02$ to $2.31\pm0.1$ at $2.35\pm0.08$ GeV. In the $\Gamma-L_{\gamma}$ plane, the luminosity of non-blazar AGNs is spread in the range of $10^{41}-10^{47}\: erg\: s^{-1}$, where the lowest luminosity have FRI radio galaxies (but typically appear with a harder photon index) and the highest- SSRQs/NLSY1s (with softer photon indexes). We confirm the previously reported short-timescale flux variability of NGC 1275 and 3C 120. The $\gamma$-ray emission from NLSY1s, 1H 0323+342, SBS 0846+513 and PMN J0948+0022, is variable, showing flares in short scales sometimes accompanied by a moderate hardening of their spectra (e.g., on MJD 56146.8 the $\gamma$-ray photon index of SBS 0846+513 was $1.73\pm0.14$). 3C 111, Cen A core, 3C 207, 3C 275.1, 3C 380, 4C+39.23B, PKS 1502+036 and PKS 2004-447 show a long-timescale flux variability in the $\gamma$-ray band.

Context. Important information on the evolution of the jet can be obtained by comparing the physical state of the plasma at its propagation through the broad-line region (where the jet is most likely formed) into the intergalactic medium, where it starts to significantly decelerate. Aims. We compare the constraints on the physical parameters in the innermost ($\leq$ pc) and outer ($\geq$ kpc) regions of the 3C 120 jet by means of a detailed multiwavelength analysis and theoretical modeling of their broadband spectra. Methods.The data collected by Fermi LAT, Swift and Chandra are analyzed together and the spectral energy distributions are modeled using a leptonic synchrotron and inverse Compton model, taking into account the seed photons originating inside and outside of the jet. The model parameters are estimated using the MCMC method. Results. The $\gamma$-ray flux from the inner jet of 3C 120 was characterized by rapid variation from MJD 56900 to MJD 57300. Two strong flares were observed on April 24, 2015 when, within 19.0 minutes and 3.15 hours the flux was as high as $(7.46\pm1.56)\times10^{-6}photon\:cm^{-2}\:s^{-1}$ and $(4.71\pm0.92)\times10^{-6}photon\:cm^{-2}\:s^{-1}$ respectively. The broadband emission in the quiet and flaring states can be described as SSC emission while IC scattering of dusty torus photons cannot be excluded for the flaring states. The X-ray emission from the knots can be well reproduced by IC scattering of CMB photons only if the jet is highly relativistic (since even when $\delta=10$ still $U_{\rm e}/U_B\geq80$). These extreme requirements can be somewhat softened assuming the X-rays are from the synchrotron emission of a second population of very-high-energy electrons. Conclusions. We found that the jet power estimated at two scales is consistent, suggesting that the jet does not suffer severe dissipation, it simply becomes radiatively inefficient.

We report on a detailed analysis of the $\gamma$-ray light curve of NGC 1275 using the Fermi large area telescope data accumulated in 2008-2017. Major $\gamma$-ray flares were observed in October 2015 and December 2016/January 2017 when the source reached a daily peak flux of $(2.21\pm0.26)\times10^{-6}\:{\rm photon\:cm^{-2}\:s^{-1}}$, achieving a flux of $(3.48\pm0.87)\times10^{-6}\:{\rm photon\:cm^{-2}\:s^{-1}}$ within 3 hours, which corresponds to an apparent isotropic $\gamma$-ray luminosity of $\simeq3.84\times10^{45}\:{\rm erg\:s^{-1}}$. The most rapid flare had e-folding time as short as $1.21\pm0.22$ hours which had never been previously observed for any radio galaxy in $\gamma$-ray band. Also $\gamma$-ray spectral changes were observed during these flares: in the flux versus photon index plane the spectral evolution follows correspondingly a counter clockwise and a clockwise loop inferred from the light curve generated by an adaptive binning method. On December 30, 2016 and January 01, 2017 the X-ray photon index softened ($\Gamma_{\rm X}\simeq 1.75-1.77$) and the flux increased nearly $\sim3$ times as compared with the quiet state. The observed hour-scale variability suggests a very compact emission region ($R_\gamma\leq5.22\times10^{14}\:(\delta/4)\:{\rm cm}$) implying that the observed emission is most likely produced in the subparsec-scale jet if the entire jet width is responsible for the emission. During the active periods the $\gamma$-ray photon index hardened, shifting the peak of the high energy spectral component to $>{\rm GeV}$, making it difficult to explain the observed X-ray and $\gamma$-ray data in the standard one-zone synchrotron self-Compton model.

We present the $\gamma$-ray observations of the flat-spectrum radio quasar PKS 1441+25 (z=0.939), using the {\it Fermi} large Area Telescope data accumulated during January - December 2015. A $\gamma$-ray flare was observed in January 24, when the flux increased up to $(2.22\pm0.38)\times10^{-6}\;{\rm photon\:cm^{-2}\:s^{-1}}$ with the flux-doubling time scale being as short as $\sim1.44$ days. The spectral analysis shows that from April 13 to April 28, 2015 the MeV-to-GeV photon index has hardened and changes in the range of $\Gamma=(1.73-1.79)$ for most of the time. The hardest photon index of $\Gamma=1.54\pm0.16$ has been observed on MJD 57131.46 with $11.8\sigma$ which is not common for flat-spectrum radio quasars. For the same period the \gray spectrum shows a possible deviation from a simple power-law shape, indicating a spectral cutoff at $E_{\rm cut}=17.7\pm8.9$ GeV. The spectral energy distributions during quiescent and flaring states are modeled using one-zone leptonic models that include the synchrotron, synchrotron self Compton and external inverse Compton processes; the model parameters are estimated using the Markov Chain Monte Carlo method. The emission in the flaring states can be modeled assuming that either the bulk Lorentz factor or the magnetic field has increased. The modeling shows that there is a hint of hardening of the low-energy index ($\sim1.98$) of the underlying non-thermal distribution of electrons responsible for the emission in April 2015. Such hardening agrees with the $\gamma$-ray data, which pointed out a significant $\gamma$-ray photon index hardening on April 13 to 28, 2015.

The recent results from ground based $\gamma$-ray detectors (HESS, MAGIC, VERITAS) provide a population of TeV galactic $\gamma$-ray sources which are potential sources of High Energy (HE) neutrinos. Since the $\gamma$-rays and $\nu$ -s are produced from decays of neutral and charged pions, the flux of TeV $\gamma$-rays can be used to estimate the upper limit of $\nu$ flux and vice versa; the detectability of $\nu$ flux implies a minimum flux of the accompanying $\gamma$-rays (assuming the internal and the external absorption of $\gamma$-rays is negligible). Using this minimum flux, it is possible to find the sources which can be detected with cubic-kilometer telescopes. I will discuss the possibility to detect HE neutrinos from powerful galactic accelerators, such as Supernova Remnants (SNRs) and Pulsar Wind Nebulae (PWNe) and show that likely only RX J1713.7-3946 , RX J0852.0-4622 and Vela X can be detected by current generation of instruments (IceCube and Km3Net). It will be shown also, that galactic binary systems could be promising sources of HE $\nu$ -s. In particular, $\nu$-s and $\gamma$-rays from Cygnus X-3 will be discussed during recent gamma-ray activity, showing that in the future such kind of activities could produce detectable flux of HE $\nu$-s

We report the analysis of Fermi Large Area Telescope data from five years of observations of the broad line radio galaxy 3C 120. The accumulation of larger data set results in the detection of high-energy $\gamma$-rays up to 10 GeV, with a detection significance of about $8.7\sigma$. A power-law spectrum with a photon index of $2.72\pm0.1$ and integrated flux of $F_{\gamma}=(2.35\pm0.5)\times10^{-8}\:\mathrm{photon\:cm}^{-2}s^{-1}$ above 100 MeV well describe the data averaged over five year observations. The variability analysis of the light curve with 180-, and 365- day bins reveals flux increase (nearly twice from its average level) during the last year of observation. This variability on month timescales indicates the compactness of the emitting region. The $\gamma$-ray spectrum can be described as synchrotron self-Compton (SSC) emission from the electron population producing the radio-to-X-ray emission in the jet. The required electron energy density exceeds the one of magnetic field only by a factor of 2 meaning no significant deviation from equipartition.

Cygnus X-3 (Cyg X-3) is a remarkable Galactic microquasar (X-ray binary) emitting from radio to $\gamma$-ray energies. In this paper, we consider hadronic model of emission of $\gamma$-rays above 100 MeV and their implications. We focus here on the joint $\gamma$-ray and neutrino production resulting from proton-proton interactions within the binary system. We find that the required proton injection kinetic power, necessary to explain the $\gamma$-ray flux observed by AGILE and Fermi-LAT, is $L_p \sim 10^{38}\:\rm{erg\:s^{-1}}$, a value in agreement with the average bolometric luminosity of the hypersoft state (when Cygnus X-3 was repeatedly observed to produce transient $\gamma$-ray activity). If we assume an increase of the wind density at the superior conjunction, the asymmetric production of $\gamma$-rays along the orbit can reproduce the observed modulation. According to observational constraints and our modelling, a maximal flux of high-energy neutrinos would be produced for an initial proton distribution with a power-law index $\alpha=2.4$. The predicted neutrino flux is almost two orders of magnitude less than the 2-month IceCube sensitivity at $\sim$1 TeV. If the protons are accelerated up to PeV energies, the predicted neutrino flux for a prolonged "soft X-ray state" would be a factor of about 3 lower than the 1-year IceCube sensitivity at $\sim$10 TeV. This study shows that, for a prolonged soft state (as observed in 2006) possibly related with $\gamma$-ray activity and a hard distribution of injected protons, Cygnus X-3 might be close to being detectable by cubic-kilometer neutrino telescopes such as IceCube.

We report on an analysis of \fermi data from four year of observations of the nearby radio galaxy Centaurus A (Cen A). The increased photon statistics results in a detection of high-energy ($>\:100$ MeV) $\gamma$-rays up to 50 GeV from the core of Cen A, with a detection significance of about 44$\sigma$. The average gamma-ray spectrum of the core reveals evidence for a possible deviation from a simple power-law. A likelihood analysis with a broken power-law model shows that the photon index becomes harder above $E_b \simeq 4$ GeV, changing from $\Gamma_1=2.74\pm0.03$ below to $\Gamma_2=2.09\pm0.20$ above. This hardening could be caused by the contribution of an additional high-energy component beyond the common synchrotron-self Compton jet emission. A variability analysis of the light curve with 15-, 30-, and 60-day bins does not provide evidence for variability for any of the components. Indications for a possible variability of the observed flux are found on 45-day time scale, but the statistics do not allow us to make a definite conclusion in this regards. We compare our results with the spectrum reported by H.E.S.S. in the TeV energy range and discuss possible origins for the hardening observed.

The detection of high energy (HE) {\gamma}-ray emission up to about 3 GeV from the giant lobes of the radio galaxy Centaurus A has been recently reported by the Fermi-LAT Collaboration based on ten months of all-sky survey observations. A data set more than three times larger is used here to study the morphology and photon spectrum of the lobes with higher statistics. The larger data set results in the detection of HE {\gamma}-ray emission (up to about 6 GeV) from the lobes with a significance of more than 10 and 20 {\sigma} for the North and the South lobe, respectively. Based on a detailed spatial analysis and comparison with the associated radio lobes, we report evidence for a substantial extension of the HE {\gamma}-ray emission beyond the WMAP radio image in the case of the Northern lobe of Cen A. We reconstruct the spectral energy distribution (SED) of the lobes using radio (WMAP) and Fermi-LAT data from the same integration region. The implications are discussed in the context of hadronic and leptonic scenarios.

With the arrival of km**3 volume scale neutrino detectors the chances to detect the first astronomical sources of TeV neutrinos will be dramatically increased. While the theoretical estimates of the neutrino fluxes contain large uncertainties, we can formulate the conditions for the detectability of certain neutrino sources phenomenologically. In fact, since most galactic neutrino sources are transparent for TeV gamma-rays, their detectability implies a minimum flux of the accompanying gamma-rays. For a typical energy-dependence of detection areas of km**3 volume neutrino detectors, we obtain the quantitative condition I_gamma(20 TeV)>2*10^-15 ph/(cm**2 s), that thanks to the normalization of the gamma-ray spectrum at 20 TeV appears to be quite robust, i.e. almost independent of the shape of energy spectrum of neutrinos. We remark that this condition is satisfied by the young supernova remnants RX J1713.7-3946 and RX J0852.0-4622 (Vela Jr) - two of the strongest galactic gamma-ray sources. The preliminary condition for the detectability of high energy neutrinos is that the bulk of gamma-rays has a hadronic origin: A new way to test this hypothesis for RX J1713.7-3946 is proposed. Finally, we assess the relevance of a neutrino detector located in the Northern Hemisphere for the search for galactic neutrino sources. In particular, we argue that if the TeV neutrino sources correlate with the galactic mass distribution, the probability that some of them will be observed by a detector in the Mediterranean Sea is larger by a factor of 1.4-2.9 compared to the one of IceCube.