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.

The detection of mJy/sub-mJy point sources is a significant challenge for single-dish radio telescopes. Detection or upper limits on the faint afterglow from GRBs or other sources at cosmological distances are important means of constraining the source modeling. Using the Sardinia Radio Telescope (SRT), we compare the sensitivity and robustness of three methods applied to the detection of faint radio sources from raster maps around a known source position: the smart quick-look method, the source extraction method (typical of high-energy astronomy), and the fit with a 2-D Gaussian. We developed a Python code specific for the analysis of point-like radio sources applied to the SRT C-band (6.9 GHz) observations of both undetected sources (GRB afterglows of 181201A and 190114C) and the detected Galactic X-ray binary GRS 1915+105. Our comparative analysis of the different detection methods made extensive use of simulations as a useful complement to actual radio observations. The best method for the SRT data analysis is the fit with a 2-D Gaussian, as it pushes down the sensitivity limits of single-dish observations -- with respect to more traditional techniques -- to ~ 1.8 mJy, improving by ~ 40 % compared with the initial value. This analysis shows that -- especially for faint sources -- good maps of the scanned region pre- or post-outburst are essential.

We report multi-color optical imaging and polarimetry observations of the afterglow of the first TeV- detected gamma-ray burst, GRB 190114C, using RINGO3 and MASTER II polarimeters. Observations begin 31 s after the onset of the GRB and continue until $\sim 7000\,$s post-burst. The light curves reveal a chromatic break at $\sim 400- 500\,$s, with initial temporal decay $\alpha = 1.669 \pm 0.013$ flattening to $\alpha \sim 1$ post-break, which we model as a combination of reverse and forward-shock components, with magnetization parameter $R_{\rm B} \sim 70$. The observed polarization degree decreases from $7.7 \pm 1.1\%$ to $2-4\%$ during $52-109\,$s post-burst and remains steady at this level for the subsequent $\sim 2000$-s, at constant position angle. Broadband spectral energy distribution modeling of the afterglow confirms GRB 190114C is highly obscured (A$_{\rm v, HG} = 1.49 \pm 0.12 \,$mag; N$_{\rm H, HG}= (9.0 \pm 0.3) \times 10^{22}\,$cm$^{-2}$). We interpret the measured afterglow polarization as intrinsically low and dominated by dust, in contrast to ${\rm P} >10\%$ measured previously for other GRB reverse shocks, with a small contribution from polarized prompt photons in the first minute. We test whether 1st and higher-order inverse Compton scattering in a magnetized reverse shock can explain the low optical polarization and the sub-TeV emission but conclude neither is explained in the reverse shock Inverse Compton model. Instead, the unexpectedly low intrinsic polarization degree in GRB 190114C can be explained if large-scale jet magnetic fields are distorted on timescales prior to reverse shock emission.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We present radio and optical afterglow observations of the TeV-bright long Gamma Ray Burst (GRB) 190114C, which was detected by the MAGIC telescope at a redshift of z=0.425. Our observations with ALMA, ATCA, and uGMRT were obtained by our low frequency observing campaign and range from ~1 to ~140 days after the burst and the optical observations were done with the 0.7-m GROWTH-India telescope up to ~25 days after the burst. Long term radio/mm observations reveal the complex nature of the afterglow, which does not conform to the predictions of the standard afterglow model. We find that the microphysical parameters of the external forward shock, representing the share of shock-created energy in the non-thermal electron population and magnetic field, are evolving with time. The kinetic energy in the blast-wave is almost an order of magnitude higher than that measured in the prompt emission.

Context. Fast Radio Bursts (FRBs) are millisecond-long bursts uniquely detected at radio frequencies. FRB 131104 is the only case for which a $\gamma$-ray transient positionally and temporally consistent was claimed. This high-energy transient had a duration of $\sim400$~s and a 15-150~keV fluence $S_{\gamma}\sim4\times10^{-6}$ erg $\mathrm{cm}^{-2}$. However, the association with the FRB is still debated. Aims. We aim at testing the systematic presence of an associated transient high-energy counterpart throughout a sample of the FRB population. Methods. We used an approach like that used in machine learning methodologies to accurately model the highly-variable Fermi/GBM instrumental background on a time interval comparable to the duration of the proposed $\gamma$-ray counterpart of FRB 131104. A possible $\gamma$-ray signal is then constrained considering sample average lightcurves. Results. We constrain the fluence of the possible $\gamma$-ray signal in the 8-1000 keV band down to $6.4 \times 10^{-7}$ ($7.1 \times 10^{-8}$) erg cm$^{-2}$ for a 200-s (1-s) integration time. Furthermore, we found the radio-to-gamma fluence ratio to be $\eta>10^{8}$ Jy ms erg$^{-1}$ cm$^2$. Conclusions. Our fluence limits exclude $\sim 94\%$ of Fermi/GBM detected long gamma-ray bursts and $\sim 96\%$ of Fermi/GBM detected short gamma-ray bursts. In addition, our limits on the radio-to-gamma fluence ratio point to a different emission mechanism from that of magnetar giant flares. Finally, we exclude a $\gamma$-ray counterpart as fluent as the one possibly associated with FRB 131104 to be a common feature of FRBs.

No transient electromagnetic emission has yet been found in association to fast radio bursts (FRBs), the only possible exception (3sigma confidence) being the putative gamma-ray signal detected in Swift/BAT data in the energy band 15-150 keV at the time and position of FRB131104. Systematic searches for hard X/gamma-ray counterparts to other FRBs ended up with just lower limits on the radio/gamma-ray fluence ratios. In 2001, at the time of the earliest discovered FRBs, the BeppoSAX Gamma-Ray Burst Monitor (GRBM) was one of the most sensitive open sky gamma-ray monitors in the 40-700~keV energy band. During its lifetime, one of the FRBs with the highest radio fluence ever recorded, FRB010724 (800 +- 400 Jy ms), also known as the Lorimer burst, was promptly visible to the GRBM. Upon an accurate modeling of the GRBM background, eased by its equatorial orbit, we searched for a possible gamma-ray signal in the first 400 s following the FRB, similar to that claimed for FRB131104 and found no significant emission down to a 5-sigma limit in the range (0.24-4.7)x10^-6 erg cm^-2 (corresponding to 1 and 400 s integration time, respectively), in the energy band 40-700 keV. This corresponds to eta = F_radio/F_gamma>10^{8-9} Jy ms erg^-1 cm^2, i.e. the deepest limit on the ratio between radio and gamma-ray fluence, which rules out a gamma-ray counterpart similar to that of FRB131104. We discuss the implications on the possible mechanisms and progenitors that have been proposed in the literature, also taking into account its relatively low dispersion measure (375 +- 3 pc cm^-3) and an inferred redshift limit of z<0.4.

In the era of multi-messenger astronomy the exploration of the early emission from transients is key for understanding the encoded physics. At the same time, current generation networks of fully-robotic telescopes provide new opportunities in terms of fast followup and sky coverage. This work describes our pipeline designed for robotic optical followup of gamma-ray bursts with the Las Cumbres Observatory network. We designed a Python code to promptly submit observation requests to the Las Cumbres Observatory network within 3 minutes of the receipt of the socket notice. Via Telegram the pipeline keeps the users informed, allowing them to take control upon request. Our group was able to track the early phases of the evolution of the optical output from gamma-ray bursts with a fully-robotic procedure and here we report the case of GRB180720B as an example. The developed pipeline represent a key ingredient for any reliable and rapid (minutes timescale) robotic telescope system. While successfully utilized and adapted for LCO, it can also be adapted to any other robotic facilities.

Broad-line type Ic Supernovae (BL-Ic SNe) are characterized by high ejecta velocity ($\gtrsim 10^4$ km s$^{-1}$) and are sometimes associated with the relativistic jets typical of long duration ($\gtrsim 2$ s) Gamma-Ray Bursts (L-GRBs). The reason why a small fraction of BL-Ic SNe harbor relativistic jets is not known. Here we present deep X-ray and radio observations of the BL-Ic SN2014ad extending from $13$ to $930$ days post explosion. SN2014ad was not detected at either frequency and has no observational evidence of a GRB counterpart. The proximity of SN2014ad ($d\sim 26$ Mpc) enables very deep constraints on the progenitor mass-loss rate $\dot{M}$ and on the total energy of the fast ejecta $E$. We consider two synchrotron emission scenarios for a wind-like circumstellar medium (CSM): (i) uncollimated non-relativistic ejecta, and (ii) off-axis relativistic jet. Within the first scenario our observations are consistent with GRB-less BL-Ic SNe characterized by a modest energy budget of their fast ejecta ($E \lesssim 10^{45}$ erg), like SNe 2002ap and 2010ay. For jetted explosions, we cannot rule out a GRB with $E \lesssim 10^{51}$ erg (beam-corrected) with a narrow opening angle ($\theta_j \sim 5^{\circ}$) observed moderately off-axis ($\theta_{\rm obs} \gtrsim 30^{\circ}$) and expanding in a very low CSM density ($\dot{M}$ $\lesssim 10^{-6}$ M$_{\odot}$ yr$^{-1}$). Our study shows that off-axis low-energy jets expanding in a low-density medium cannot be ruled out even in the most nearby BL-Ic SNe with extensive deep observations, and might be a common feature of BL-Ic SNe.

The main characteristics in the radio continuum spectra of Supernova Remnants (SNRs) result from simple synchrotron emission. In addition, electron acceleration mechanisms can shape the spectra in specific ways, especially at high radio frequencies. These features are connected to the age and the peculiar conditions of the local interstellar medium interacting with the SNR. Whereas the bulk radio emission is expected at up to $20-50$ GHz, sensitive high-resolution images of SNRs above 10 GHz are lacking and are not easily achievable, especially in the confused regions of the Galactic Plane. In the framework of the early science observations with the Sardinia Radio Telescope in February-March 2016, we obtained high-resolution images of SNRs Tycho, W44 and IC443 that provided accurate integrated flux density measurements at 21.4 GHz: 8.8 $\pm$ 0.9 Jy for Tycho, 25 $\pm$ 3 Jy for W44 and 66 $\pm$ 7 Jy for IC443. We coupled the SRT measurements with radio data available in the literature in order to characterise the integrated and spatially-resolved spectra of these SNRs, and to find significant frequency- and region-dependent spectral slope variations. For the first time, we provide direct evidence of a spectral break in the radio spectral energy distribution of W44 at an exponential cutoff frequency of 15 $\pm$ 2 GHz. This result constrains the maximum energy of the accelerated electrons in the range $6-13$ GeV, in agreement with predictions indirectly derived from AGILE and \textit{Fermi}-LAT gamma-ray observations. With regard to IC443, our results confirm the noticeable presence of a bump in the integrated spectrum around $20-70$ GHz that could result from a spinning dust emission mechanism.

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.

Observations of supernova remnants (SNRs) are a powerful tool for investigating the later stages of stellar evolution, the properties of the ambient interstellar medium, and the physics of particle acceleration and shocks. For a fraction of SNRs, multi-wavelength coverage from radio to ultra high-energies has been provided, constraining their contributions to the production of Galactic cosmic rays. Although radio emission is the most common identifier of SNRs and a prime probe for refining models, high-resolution images at frequencies above 5 GHz are surprisingly lacking, even for bright and well-known SNRs such as IC443 and W44. In the frameworks of the Astronomical Validation and Early Science Program with the 64-m single-dish Sardinia Radio Telescope, we provided, for the first time, single-dish deep imaging at 7 GHz of the IC443 and W44 complexes coupled with spatially-resolved spectra in the 1.5-7 GHz frequency range. Our images were obtained through on-the-fly mapping techniques, providing antenna beam oversampling and resulting in accurate continuum flux density measurements. The integrated flux densities associated with IC443 are S_1.5GHz = 134 +/- 4 Jy and S_7GHz = 67 +/- 3 Jy. For W44, we measured total flux densities of S_1.5GHz = 214 +/- 6 Jy and S_7GHz = 94 +/- 4 Jy. Spectral index maps provide evidence of a wide physical parameter scatter among different SNR regions: a flat spectrum is observed from the brightest SNR regions at the shock, while steeper spectral indices (up to 0.7) are observed in fainter cooling regions, disentangling in this way different populations and spectra of radio/gamma-ray-emitting electrons in these SNRs.

We present a long-term spectral monitoring of the unique Double Pulsar binary PSR J0737-3039 corresponding to two "Large Programs" performed by XMM-Newton in 2006 and 2011. Spectral variability of pulsar emission in soft X-rays is not evident over 5 years, despite the significant relativistic spin precession in the considered time span (about 25deg). We provide, for the first time, evidence of hard X-ray emission from the system in the 5-8 keV energy band. The standard spectral analysis was coupled to the energy dependent spatial analysis to confirm this excess, most likely ascribed to iron line emission. The Fe K-alpha emission line at 6.4-6.97 keV was previously unheard-of in non-accreting binary systems and could testify to the presence of a relic disk that survived the supernova explosions that terminated the lives of the Double Pulsar's stellar progenitors. The existence of a relic disk in this system reinforces speculation about the presence of similar structures around other peculiar classes of isolated neutron stars.

In the frame of the Astronomical Validation activities for the 64m Sardinia Radio Telescope, we performed 5-22 GHz imaging observations of the complex-morphology supernova remnants (SNRs) W44 and IC443. We adopted innovative observing and mapping techniques providing unprecedented accuracy for single-dish imaging of SNRs at these frequencies, revealing morphological details typically available only at lower frequencies through interferometry observations. High-frequency studies of SNRs in the radio range are useful to better characterize the spatially-resolved spectra and the physical parameters of different regions of the SNRs interacting with the ISM. Furthermore, synchrotron-emitting electrons in the high-frequency radio band are also responsible for the observed high-energy phenomenology as -e.g.- Inverse Compton and bremsstrahlung emission components observed in gamma-rays, to be disentangled from hadron emission contribution (providing constraints on the origin of cosmic rays).

Supernova Remnants (SNRs) exhibit spectra featured by synchrotron radio emission arising from the relativistic electrons, and high-energy emission from both leptonic (Bremsstrahlung and Inverse Compton) and hadronic processes (${\pi}^0$ mesons decay) which are a direct signature of cosmic rays acceleration. Thanks to radio single-dish imaging observations obtained in three frequency bands (1.6, 7, 22 GHz) with the Sardinia Radio Telescope (www.srt.inaf.it), we can model different SNR regions separately. Indeed, in order to disentangle interesting and peculiar hadron contributions in the high-energy spectra (gamma-ray band) and better constrain SNRs as cosmic rays emitters, it is crucial to fully constrain lepton contributions first through radio-observed parameters. In particular, the Bremsstrahlung and Inverse Compton bumps observed in gamma-rays are bounded to synchrotron spectral slope and cut-off in the radio domain. Since these parameters vary for different SNR regions and electron populations, spatially-resolved radio spectra are then required for accurate multi-wavelength modelling.