We report on the lowest-frequency detection to date of three bursts from the fast radio burst FRB 180916, observed at 328 MHz with the Sardinia Radio Telescope (SRT). The SRT observed the periodic repeater FRB 180916 for five days from 2020 February 20 to 24 during a time interval of active radio bursting, and detected the three bursts during the first hour of observations; no more bursts were detected during the remaining ~ 30 hours. Simultaneous SRT observations at 1548 MHz did not detect any bursts. Burst fluences are in the range 37 to 13 Jy ms. No relevant scattering is observed for these bursts. We also present the results of the multi-wavelength campaign we performed on FRB 180916, during the five days of the active window. Simultaneously with the SRT observations, others with different time spans were performed with the Northern Cross at 408 MHz, with XMM-Newton, NICER, INTEGRAL, AGILE, and with the TNG and two optical telescopes in Asiago, which are equipped with fast photometers. XMM-Newton obtained data simultaneously with the three bursts detected by the SRT, and determined a luminosity upper limit in the 0.3-10 keV energy range of ~$10^{45}$ erg/s for the burst emission. AGILE obtained data simultaneously with the first burst and determined a fluence upper limit in the MeV range for millisecond timescales of $ 10^{-8}$ erg cm$^{-2}$.Our results show that absorption from the circumburst medium does not significantly affect the emission from FRB 180916, thus limiting the possible presence of a superluminous supernova around the source, and indicate that a cutoff for the bursting mechanism, if present, must be at lower frequencies. Our multi-wavelength campaign sensitively constrains the broadband emission from FRB 180916, and provides the best limits so far for the electromagnetic response to the radio bursting of this remarkable source of fast radio bursts.

We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time.

Fast radio bursts (FRBs) are short (millisecond) radio pulses originating from enigmatic sources at extragalactic distances so far lacking a detection in other energy bands. Magnetized neutron stars (magnetars) have been considered as the sources powering the FRBs, but the connection is controversial because of differing energetics and the lack of radio and X-ray detections with similar characteristics in the two classes. We report here the detection by the AGILE satellite on April 28, 2020 of an X-ray burst in coincidence with the very bright radio burst from the Galactic magnetar SGR 1935+2154. The burst detected by AGILE in the hard X-ray band (18-60 keV) lasts about 0.5 seconds, it is spectrally cutoff above 80 keV, and implies an isotropically emitted energy ~ $10^{40}$ erg. This event is remarkable in many ways: it shows for the first time that a magnetar can produce X-ray bursts in coincidence with FRB-like radio bursts; it also suggests that FRBs associated with magnetars may emit X-ray bursts of both magnetospheric and radio-pulse types that may be discovered in nearby sources. Guided by this detection, we discuss SGR 1935+2154 in the context of FRBs, and especially focus on the class of repeating-FRBs. Based on energetics, magnetars with fields B ~ $10^{15}$ G may power the majority of repeating-FRBs. Nearby repeating-FRBs offer a unique occasion to consolidate the FRB-magnetar connection, and we present new data on the X-ray monitoring of nearby FRBs. Our detection enlightens and constrains the physical process leading to FRBs: contrary to previous expectations, high-brightness temperature radio emission coexists with spectrally-cutoff X-ray radiation.

FRB 180916 is a most intriguing source at 150 Mpc distance capable of producing repeating fast radio bursts with a periodic 16.35 day temporal pattern. We report on the X-ray and $\gamma$-ray observations of FRB 180916 obtained by AGILE and Swift. We focused on the recurrent 5-day time intervals of active radio bursting and present results obtained on Feb. 3 - 8; Feb. 25; Mar. 5 - 10; Mar. 22 - 28, 2020 during a multiwavelength campaign involving high-energy and radio observations. We also searched for temporal coincidences at millisecond timescales between all known radio bursts of FRB 180916 and X-ray and MeV events detectable by AGILE. We do not detect any simultaneous event or any extended X-ray and $\gamma$-ray emission on timescales of hours/days/weeks. Our cumulative X-ray (0.3-10 keV) flux upper limit of $5 \times\,10^{-14} \rm \, erg \, cm^{-2} s^{-1}$ (obtained during 5-day active intervals) translates into an isotropic luminosity upper limit of $L_{X,UL} \sim 1.5 \times\, 10^{41} \rm erg \, s^{-1}$. Observations above 100 MeV over a many-year timescale provide an average luminosity upper limit one order of magnitude larger. These results provide the so-far most stringent limits on high-energy emission from FRB 180916 and constrain the dissipation of magnetic energy from a magnetar-like source of radius $R_m$, internal magnetic field $B_m$ and dissipation timescale $\tau_d$ to satisfy the relation $R_{m,6}^3 B_{m,16}^2 \tau_{d,8}^{-1} \lesssim 1$, where $R_{m,6}$ is $R_m$ in units of $10^6$ cm, $B_{m,16}$ is $B_m$ in units of $10^{16}$ G, and $\tau_{d,8}$ in units of $10^8$ s.

Fast radio bursts remain one of the most enigmatic astrophysical sources. Observations have significantly progressed over the last few years, thanks to the capabilities of new radio telescopes and the refurbishment of existing ones. Here we describe the upgrade of the Northern Cross radio telescope, operating in the 400-416~MHz frequency band, with the ultimate goal of turning the array into a dedicated instrument to survey the sky for fast radio bursts. We present test observations of the pulsar B0329+54 to characterize the system performance and forecast detectability. Observations with the system currently in place are still limited by modest sky coverage ($\sim 9.4$~deg$^2$) and biased by smearing of high dispersion measure events within each frequency channels. In its final, upgraded configuration, however, the telescope will be able to carry out unbiased fast radio burst surveys over a $\sim 350$~deg$^2$ instantaneous field of view up to $z \sim 5$, with a (nearly constant) $\sim 760 \, (\tau/{\rm ms})^{-0.5}$~mJy rms sensitivity.

In this work, we propose a new method of computing the Karhunen-Lo\`{e}ve Transform (KLT) applied to complex voltage data for the detection and noise level reduction in astronomical signals. We compared this method with the standard KLT techniques based on the Toeplitz correlation matrix and we conducted a performance analysis for the detection and extraction of astrophysical and artificial signals via Monte Carlo simulations. We applied our novel method to a real data study-case: the Voyager 1 telemetry signal. We evaluated the KLT performance in an astrophysical context: our technique provides a remarkable improvement in computation time and Monte-Carlo simulations show significant reconstruction results for signal-to-noise ratio (SNR) down to -10 dB and comparable results with standard signal detection techniques. The application to artificial signals, such as the Voyager 1 data, shows a notable gain in SNR after the KLT.

We focus on two repeating fast radio bursts (FRBs) recently detected by the CHIME/FRB experiment in 2018--2019 (Source 1: 180916.J0158+65, and Source 2: 181030.J1054+73). These sources have low excess dispersion measures (DMs) ($ < 100 \rm \, pc \, cm^{-3}$ and $ < 20 \rm \, pc \, cm^{-3}$, respectively), implying relatively small maximal distances. They were repeatedly observed by AGILE in the MeV--GeV energy range. We do not detect prompt emission simultaneously with these repeating events. This search is particularly significant for the submillisecond and millisecond integrations obtainable by AGILE. The sources are constrained to emit a MeV-fluence in the millisecond range below $F'_{MeV} = 10^{-8} \, \rm erg \, cm^{-2}$ corresponding to an isotropic energy near $E_{MeV,UL} \simeq 2 \times 10^{46}\,$erg for a distance of 150 Mpc (applicable to Source 1). We also searched for $\gamma$-ray emission for time intervals up to 100 days, obtaining 3$\,\sigma$ upper limits (ULs) for the average isotropic luminosity above 50 MeV, $L_{\gamma,UL} \simeq \,$(5-10)$\,\times 10^{43} \rm \, erg \, s^{-1}$. For a source distance near 100 kpc (possibly applicable to Source 2), our ULs imply $E_{MeV,UL}\simeq10^{40} \rm erg$, and $L_{\gamma,UL} \simeq \,$2$\,\times 10^{37} \rm \, erg \, s^{-1}$. Our results are significant in constraining the high-energy emission of underlying sources such as magnetars, or other phenomena related to extragalactic compact objects, and show the prompt emission to be lower than the peak of the 2004 magnetar outburst of SGR 1806-20 for source distances less than about 100 Mpc.

The discovery of the ubiquity of habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, have ushered in a renaissance in the millenia-old quest to answer our most profound question about the Universe and our place within it - Are we alone? The Breakthrough Listen Initiative, announced in July 2015 as a 10-year 100M USD program, is the most comprehensive effort in history to quantify the distribution of advanced, technologically capable life in the universe. In this white paper, we outline the status of the on-going observing campaign with our primary observing facilities, as well as planned activities with these instruments over the next few years. We also list collaborative facilities which will conduct searches for technosignatures in either primary observing mode, or commensally. We highlight some of the novel analysis techniques we are bringing to bear on multi-petabyte data sets, including machine learning tools we are deploying to search for a broader range of technosignatures than was previously possible.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We present a study on the gamma-ray emission detected by the AGILE-GRID from the region of the SNR G78.2+2.1 (Gamma Cygni). In order to investigate the possible presence of gamma rays associated with the SNR below 1 GeV, it is necessary to analyze the gamma-ray radiation underlying the strong emission from the pulsar PSR J2021+4026, which totally dominates the field. An "off-pulse" analysis has been carried out, by considering only the emission related to the pulsar off-pulse phase of the AGILE-GRID light curve. We found that the resulting off-pulsed emission in the region of the SNR - detected by the AGILE-GRID above 400 MeV - partially overlaps the radio shell boundary. By analyzing the averaged emission on the whole angular extent of the SNR, we found that a lepton-dominated double population scenario can account for the radio and gamma-ray emission from the source. In particular, the MeV-GeV averaged emission can be explained mostly by Bremsstrahlung processes in a high density medium, whereas the GeV-TeV radiation by both Bremsstrahlung (E < 250 GeV) and inverse Compton processes (E > 250 GeV) in a lower density medium.

Pulsar studies in the recent years have shown, more than others, to have benefited from a multi-wavelength approach. The INAF - Astronomical Observatory in Cagliari (INAF-OAC) is a growing facility with a young group devoted to pulsar and fast transients studies across the electromagnetic spectrum. Taking advantage of this expertise we have worked to provide a suite of multi-wavelength software and databases for the observations of pulsars and compact Galactic objects at the Sardinia Radio Telescope (SRT). In turn, radio pulsar observations at SRT will be made available, in a processed format, to gamma-ray searches using AGILE and Fermi gamma-ray satellite and, in a near future, they will be complementary to polarimetric X-ray observations with IXPE.

Aims. We present the second AGILE-GRID Catalog (2AGL) of {\gamma}-ray sources in the 100 MeV-10 GeV energy range. Methods. With respect to previous AGILE-GRID catalogs, the current 2AGL Catalog is based on the first 2.3 years of science data from the AGILE mission (the so called 'pointing mode') and incorporates more data and several analysis improvements, including better calibrations at the event reconstruction level, an updated model for the Galactic diffuse gamma-ray emission, a refined procedure for point-like source detection, and the inclusion of a search for extended {\gamma}-ray sources. Results. The 2AGL Catalog includes 175 high-confidence sources (above 4{\sigma} significance) with their location regions and spectral properties, and a variability analysis with 4-day light curves for the most significant ones. Relying on the error region of each source position, including systematic uncertainties, 121 sources are considered as positionally associated with known couterparts at different wavelengths or detected by other {\gamma}-ray instruments. Among the identified or associated sources, 62 are Active Galactic Nuclei (AGNs) of the blazar class. Pulsars represent the largest Galactic source class, with 40 associated pulsars, 7 of them with detected pulsation; 8 Supernova Remnants and 4 high-mass X-ray binaries have also been identified. A substantial number of 2AGL sources are unidentified: for 54 sources no known counterpart is found at different wavelengths. Among these sources, we discuss a sub-class of 29 AGILE-GRID-only {\gamma}-ray sources that are not present in 1FGL, 2FGL or 3FGL catalogs; the remaining sources are unidentified in both 2AGL and 3FGL Catalogs. We also present an extension of the analysis of 2AGL sources detected in the 50-100 MeV energy range.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

We determined Faraday rotation measures (RMs) towards 137 pulsars in the northern sky, using Low-Frequency Array (LOFAR) observations at 110-190 MHz. This low-frequency RM catalogue, the largest to date, improves the precision of existing RM measurements on average by a factor of 20 - due to the low frequency and wide bandwidth of the data, aided by the RM synthesis method. We report RMs towards 25 pulsars for the first time. The RMs were corrected for ionospheric Faraday rotation to increase the accuracy of our catalogue to approximately 0.1 rad m$^{\rm -2}$. The ionospheric RM correction is currently the largest contributor to the measurement uncertainty. In addition, we find that the Faraday dispersion functions towards pulsars are extremely Faraday thin - mostly less than 0.001 rad m$^{\rm -2}$. We use these new precise RM measurements (in combination with existing RMs, dispersion measures, and distance estimates) to estimate the scale height of the Galactic halo magnetic field: 2.0$\pm$0.3 kpc for Galactic quadrants I and II above and below the Galactic plane (we also evaluate the scale height for these regions individually). Overall, our initial low-frequency catalogue provides valuable information about the 3-D structure of the Galactic magnetic field.

The origin of cosmic neutrinos is still largely unknown. Using data obtained by the gamma-ray imager on board of the AGILE satellite, we systematically searched for transient gamma-ray sources above 100 MeV that are temporally and spatially coincident with ten recent high-energy neutrino IceCube events. We find three AGILE candidate sources that can be considered possible counterparts to neutrino events. Detecting 3 gamma-ray/neutrino associations out of 10 IceCube events is shown to be unlikely due to a chance coincidence. One of the sources is related to the BL Lac source TXS 0506+056. For the other two AGILE gamma-ray sources there are no obvious known counterparts, and both Galactic and extragalactic origin should be considered.

AGILE is a mission of the Italian Space Agency (ASI) Scientific Program dedicated to gamma-ray astrophysics, operating in a low Earth orbit since April 23, 2007. It is designed to be a very light and compact instrument, capable of simultaneously detecting and imaging photons in the 18 keV to 60 keV X-ray energy band and in the 30 MeV{50 GeV gamma-ray energy with a good angular resolution (< 1 deg at 1 GeV). The core of the instrument is the Silicon Tracker complemented with a CsI calorimeter and a AntiCoincidence system forming the Gamma Ray Imaging Detector (GRID). Before launch, the GRID needed on-ground calibration with a tagged gamma-ray beam to estimate its performance and validate the Monte Carlo simulation. The GRID was calibrated using a tagged gamma-ray beam with energy up to 500 MeV at the Beam Test Facilities at the INFN Laboratori Nazionali di Frascati. These data are used to validate a GEANT3 based simulation by comparing the data and the Monte Carlo simulation by measuring the angular and energy resolutions. The GRID angular and energy resolutions obtained using the beam agree well with the Monte Carlo simulation. Therefore the simulation can be used to simulate the same performance on-light with high reliability.

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.

We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.

The LIGO-Virgo Collaboration (LVC) detected, on 2017 August 17, an exceptional gravitational-wave (GW) event temporally consistent within $\sim\,1.7 \, \rm s$ with the GRB 1708117A observed by Fermi-GBM and INTEGRAL. The event turns out to be compatible with a neutron star-neutron star (NS-NS) coalescence that subsequently produced a radio/optical/X-ray transient detected at later times. We report the main results of the observations by the AGILE satellite of the GW170817 localization region (LR) and its electromagnetic (e.m.) counterpart. At the LVC detection time $T_0$, the GW170817 LR was occulted by the Earth. The AGILE instrument collected useful data before and after the GW-GRB event because in its spinning observation mode it can scan a given source many times per hour. The earliest exposure of the GW170817 LR by the gamma-ray imaging detector (GRID) started about 935 s after $T_0$. No significant X-ray or gamma-ray emission was detected from the LR that was repeatedly exposed over timescales of minutes, hours, and days before and after GW170817, also considering Mini-calorimeter and Super-AGILE data. Our measurements are among the earliest ones obtained by space satellites on GW170817 and provide useful constraints on the precursor and delayed emission properties of the NS-NS coalescence event. We can exclude with high confidence the existence of an X-ray/gamma-ray emitting magnetar-like object with a large magnetic field of $10^{15} \, \rm G$. Our data are particularly significant during the early stage of evolution of the e.m. remnant.

The LIGO/Virgo Collaboration (LVC) detected on 2017 January 4, a significant gravitational-wave (GW) event (now named GW170104). We report in this Letter the main results obtained from the analysis of hard X-ray and gamma-ray data of the AGILE mission that repeatedly observed the GW170104 localization region (LR). At the LVC detection time $T_0$ AGILE observed about 36% of the LR. The gamma-ray imaging detector did not reveal any significant emission in the energy range 50 MeV--30 GeV. Furthermore, no significant gamma-ray transients were detected in the LR that was repeatedly exposed over timescales of minutes, hours, and days. We also searched for transient emission using data near $T_0$ of the omnidirectional detector MCAL operating in the energy band 0.4--100 MeV. A refined analysis of MCAL data shows the existence of a weak event (that we call "E2") with a signal-to-noise ratio of $4.4\,\sigma$ lasting about 32 ms and occurring $0.46\,\pm\,0.05 \,\rm s$ before $T_0$. A study of the MCAL background and of the false-alarm rate of E2 leads to the determinination of a post-trial significance between $2.4\,\sigma$ and $2.7\,\sigma$ for a temporal coincidence with GW170104. We note that E2 has characteristics similar to those detected from the weak precursor of the short GRB 090510. The candidate event E2 is worth consideration for simultaneous detection by other satellites. If associated with GW170104, it shows emission in the MeV band of a short burst preceding the final coalescence by 0.46 sec and involving $\sim 10^{-7}$ of the total rest mass energy of the system.

On July 31st, 2016, the ICECUBE collaboration reported the detection of a high-energy starting event induced by an astrophysical neutrino. We report here about the search for a gamma-ray counterpart of the ICECUBE-160731 event made with the AGILE satellite. No detection was found spanning the time interval of +/- 1 ks around the neutrino event time T0 using the AGILE "burst search" system. Looking for a possible gamma-ray precursor in the results of the AGILE-GRID automatic Quick Look procedure over predefined 48-hours time-bins, we found an excess above 100 MeV between one and two days before T0, positionally consistent with the ICECUBE error circle, having a post-trial significance of about 4 sigma. A refined data analysis of this excess confirms a-posteriori the automatic detection. The new AGILE transient source, named AGL J1418+0008, thus stands as possible ICECUBE-160731 gamma-ray precursor. No other space missions nor ground observatories have reported any detection of transient emission consistent with the ICECUBE event. We show that Fermi-LAT had a low exposure of the ICECUBE region during the AGILE gamma-ray transient. Based on an extensive search for cataloged sources within the error regions of ICECUBE-160731 and AGL J1418+0008, we find a possible common counterpart showing some of the key features associated to the high-energy peaked BL Lac (HBL) class of blazars. Further investigations on the nature of this source using dedicated SWIFT ToO data are presented.

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.

We measure the effects of interstellar scattering on average pulse profiles from 13 radio pulsars with simple pulse shapes. We use data from the LOFAR High Band Antennas, at frequencies between 110 and 190~MHz. We apply a forward fitting technique, and simultaneously determine the intrinsic pulse shape, assuming single Gaussian component profiles. We find that the constant $\tau$, associated with scattering by a single thin screen, has a power-law dependence on frequency $\tau \propto \nu^{-\alpha}$, with indices ranging from $\alpha = 1.50$ to $4.0$, despite simplest theoretical models predicting $\alpha = 4.0$ or $4.4$. Modelling the screen as an isotropic or extremely anisotropic scatterer, we find anisotropic scattering fits lead to larger power-law indices, often in better agreement with theoretically expected values. We compare the scattering models based on the inferred, frequency dependent parameters of the intrinsic pulse, and the resulting correction to the dispersion measure (DM). We highlight the cases in which fits of extreme anisotropic scattering are appealing, while stressing that the data do not strictly favour either model for any of the 13 pulsars. The pulsars show anomalous scattering properties that are consistent with finite scattering screens and/or anisotropy, but these data alone do not provide the means for an unambiguous characterization of the screens. We revisit the empirical $\tau$ versus DM relation and consider how our results support a frequency dependence of $\alpha$. Very long baseline interferometry, and observations of the scattering and scintillation properties of these sources at higher frequencies, will provide further evidence.

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.

[Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015.

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.

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).

A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.

We present new observations of the galaxy cluster 3C 129 obtained with the Sardinia Radio Telescope in the frequency range 6000-7200 MHz, with the aim to image the large-angular-scale emission at high-frequency of the radio sources located in this cluster of galaxies. The data were acquired using the recently-commissioned ROACH2-based backend to produce full-Stokes image cubes of an area of 1 deg x 1 deg centered on the radio source 3C 129. We modeled and deconvolved the telescope beam pattern from the data. We also measured the instrumental polarization beam patterns to correct the polarization images for off-axis instrumental polarization. Total intensity images at an angular resolution of 2.9 arcmin were obtained for the tailed radio galaxy 3C 129 and for 13 more sources in the field, including 3C 129.1 at the galaxy cluster center. These data were used, in combination with literature data at lower frequencies, to derive the variation of the synchrotron spectrum of 3C 129 along the tail of the radio source. If the magnetic field is at the equipartition value, we showed that the lifetimes of radiating electrons result in a radiative age for 3C 129 of t_syn = 267 +/- 26 Myrs. Assuming a linear projected length of 488 kpc for the tail, we deduced that 3C 129 is moving supersonically with a Mach number of M=v_gal/c_s=1.47. Linearly polarized emission was clearly detected for both 3C 129 and 3C 129.1. The linear polarization measured for 3C 129 reaches levels as high as 70% in the faintest region of the source where the magnetic field is aligned with the direction of the tail.

We report the results of an extensive search in the AGILE data for a gamma-ray counterpart of the LIGO gravitational wave event GW150914. Currently in spinning mode, AGILE has the potential of covering with its gamma-ray instrument 80 % of the sky more than 100 times a day. It turns out that AGILE came within a minute from the event time of observing the accessible GW150914 localization region. Interestingly, the gamma-ray detector exposed about 65 % of this region during the 100 s time intervals centered at -100 s and +300 s from the event time. We determine a 2-sigma flux upper limit in the band 50 MeV - 10 GeV, $UL = 1.9 \times 10^{-8} \rm \, erg \, cm^{-2} \, s^{-1}$ obtained about 300 s after the event. The timing of this measurement is the fastest ever obtained for GW150914, and significantly constrains the electromagnetic emission of a possible high-energy counterpart. We also carried out a search for a gamma-ray precursor and delayed emission over timescales ranging from minutes to days: in particular, we obtained an optimal exposure during the interval -150 / -30 s. In all these observations, we do not detect a significant signal associated with GW150914. We do not reveal the weak transient source reported by Fermi-GBM 0.4 s after the event time. However, even though a gamma-ray counterpart of the GW150914 event was not detected, the prospects for future AGILE observations of gravitational wave sources are decidedly promising.

Multi-messenger astronomy is becoming the key to understanding the Universe from a comprehensive perspective. In most cases, the data and the technology are already in place, therefore it is important to provide an easily-accessible package that combines datasets from multiple telescopes at different wavelengths. In order to achieve this, we are working to produce a data analysis pipeline that allows the data reduction from different instruments without needing detailed knowledge of each observation. Ideally, the specifics of each observation are automatically dealt with, while the necessary information on how to handle the data in each case is provided by a tutorial that is included in the program. We first focus our project on the study of pulsars and their wind nebulae (PWNe) at radio and gamma-ray frequencies. In this way, we aim to combine time-domain and imaging datasets at two extremes of the electromagnetic spectrum. In addition, the emission has the same non-thermal origin in pulsars at radio and gamma-ray frequencies, and the population of electrons is believed to be the same at these energies in PWNe. The final goal of the project will be to unveil the properties of these objects by tracking their behaviour using all of the available multi-wavelength data.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing.

We present broadband, low-frequency (25-80 MHz and 110-190 MHz) LOFAR observations of PSR B0943+10, with the goal of better illuminating the nature of its enigmatic mode-switching behaviour. This pulsar shows two relatively stable states: a Bright (B) and Quiet (Q) mode, each with different characteristic brightness, profile morphology, and single-pulse properties. We model the average profile evolution both in frequency and time from the onset of each mode, and highlight the differences between the two modes. In both modes, the profile evolution can be well explained by radius-to-frequency mapping at altitudes within a few hundred kilometres of the stellar surface. If both B and Q-mode emission originate at the same magnetic latitude, then we find that the change of emission height between the modes is less than 6%. We also find that, during B-mode, the average profile is gradually shifting towards later spin phase and then resets its position at the next Q-to-B transition. The observed B-mode profile delay is frequency-independent (at least from 25-80 MHz) and asymptotically changes towards a stable value of about 0.004 in spin phase by the end of mode instance, much too large to be due to changing spin-down rate. Such a delay can be interpreted as a gradual movement of the emission cone against the pulsar's direction of rotation, with different field lines being illuminated over time. Another interesting explanation is a possible variation of accelerating potential inside the polar gap. This explanation connects the observed profile delay to the gradually evolving subpulse drift rate, which depends on the gradient of the potential across the field lines.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

High-energy phenomena in the cosmos, and in particular processes leading to the emission of gamma- rays in the energy range 10 MeV - 100 GeV, play a very special role in the understanding of our Universe. This energy range is indeed associated with non-thermal phenomena and challenging particle acceleration processes. The technology involved in detecting gamma-rays is challenging and drives our ability to develop improved instruments for a large variety of applications. GAMMA-LIGHT is a Small Mission which aims at an unprecedented advance of our knowledge in many sectors of astrophysical and Earth studies research. The Mission will open a new observational window in the low-energy gamma-ray range 10-50 MeV, and is configured to make substantial advances compared with the previous and current gamma-ray experiments (AGILE and Fermi). The improvement is based on an exquisite angular resolution achieved by GAMMA-LIGHT using state-of-the-art Silicon technology with innovative data acquisition. GAMMA-LIGHT will address all astrophysics issues left open by the current generation of instruments. In particular, the breakthrough angular resolution in the energy range 100 MeV - 1 GeV is crucial to resolve patchy and complex features of diffuse sources in the Galaxy as well as increasing the point source sensitivity. This proposal addresses scientific topics of great interest to the community, with particular emphasis on multifrequency correlation studies involving radio, optical, IR, X-ray, soft gamma-ray and TeV emission. At the end of this decade several new observatories will be operational including LOFAR, SKA, ALMA, HAWK, CTA. GAMMA-LIGHT will "fill the vacuum" in the 10 MeV-10 GeV band, and will provide invaluable data for the understanding of cosmic and terrestrial high-energy sources.

We present a variability study of a sample of bright gamma-ray (30 MeV -- 50 GeV) sources. This sample is an extension of the first AGILE catalogue of gamma-ray sources (1AGL), obtained using the complete set of AGILE observations in pointing mode performed during a 2.3 year period from July 9, 2007 until October 30, 2009. The dataset of AGILE pointed observations covers a long time interval and its gamma-ray data archive is useful for monitoring studies of medium-to-high brightness gamma-ray sources. In the analysis reported here, we used data obtained with an improved event filter that covers a wider field of view, on a much larger (about 27.5 months) dataset, integrating data on observation block time scales, which mostly range from a few days to thirty days. The data processing resulted in a better characterized source list than 1AGL was, and includes 54 sources, 7 of which are new high galactic latitude (|BII| >= 5) sources, 8 are new sources on the galactic plane, and 20 sources from the previous catalogue with revised positions. Eight 1AGL sources (2 high-latitude and 6 on the galactic plane) were not detected in the final processing either because of low OB exposure and/or due to their position in complex galactic regions. We report the results in a catalogue of all the detections obtained in each single OB, including the variability results for each of these sources. In particular, we found that 12 sources out of 42 or 11 out of 53 are variable, depending on the variability index used, where 42 and 53 are the number of sources for which these indices could be calculated. Seven of the 11 variable sources are blazars, the others are Crab pulsar+nebula, LS I +61{\deg}303, Cyg X-3, and 1AGLR J2021+4030.

Context: AGILE is a gamma-ray astrophysics mission which has been in orbit since 23 April 2007 and continues to operate reliably. The gamma-ray detector, AGILE-GRID, has observed Galactic and extragalactic sources, many of which were collected in the first AGILE Catalog. Aims: We present the calibration of the AGILE-GRID using in-flight data and Monte Carlo simulations, producing Instrument Response Functions (IRFs) for the effective area A_eff), Energy Dispersion Probability (EDP), and Point Spread Function (PSF), each as a function of incident direction in instrument coordinates and energy. Methods: We performed Monte Carlo simulations at different gamma-ray energies and incident angles, including background rejection filters and Kalman filter-based gamma-ray reconstruction. Long integrations of in-flight observations of the Vela, Crab and Geminga sources in broad and narrow energy bands were used to validate and improve the accuracy of the instrument response functions. Results: The weighted average PSFs as a function of spectra correspond well to the data for all sources and energy bands. Conclusions: Changes in the interpolation of the PSF from Monte Carlo data and in the procedure for construction of the energy-weighted effective areas have improved the correspondence between predicted and observed fluxes and spectra of celestial calibration sources, reducing false positives and obviating the need for post-hoc energy-dependent scaling factors. The new IRFs have been publicly available from the Agile Science Data Centre since November 25, 2011, while the changes in the analysis software will be distributed in an upcoming release.

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

We present the results of multi-year gamma-ray observations by the AGILE satellite of the black hole binary system Cygnus X-1. In a previous investigation we focused on gamma-ray observations of Cygnus X-1 in the hard state during the period mid-2007/2009. Here we present the results of the gamma-ray monitoring of Cygnus X-1 during the period 2010/mid-2012 carried out for which includes a remarkably prolonged `soft state' phase (June 2010 -- May 2011). Previous 1--10 MeV observations of Cyg X-1 in this state hinted at a possible existence of a non-thermal particle component with substantial modifications of the Comptonized emission from the inner accretion disk. Our AGILE data, averaged over the mid-2010/mid-2011 soft state of Cygnus X-1, provide a significant upper limit for gamma-ray emission above 100 MeV of F_soft < 20 x 10^{-8} ph/cm^2/s, excluding the existence of prominent non-thermal emission above 100 MeV during the soft state of Cygnus X-1. We discuss theoretical implications of our findings in the context of high-energy emission models of black hole accretion. We also discuss possible gamma-ray flares detected by AGILE. In addition to a previously reported episode observed by AGILE in October 2009 during the hard state, we report a weak but important candidate for enhanced emission which occurred at the end of June 2010 (2010-06-30 10:00 - 2010-07-02 10:00 UT) exactly in coincidence with a hard-to-soft state transition and before an anomalous radio flare. An appendix summarizes all previous high-energy observations and possible detections of Cygnus X-1 above 1 MeV.

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

Magnetars are an extreme, highly magnetized class of isolated neutron stars whose large X-ray luminosity is believed to be driven by their high magnetic field. In this work we study for the first time the possible very high energy gamma-ray emission above 100 GeV from magnetars, observing the sources 4U 0142+61 and 1E 2259+586. We observed the two sources with atmospheric Cherenkov telescopes in the very high energy range (E > 100 GeV). 4U 0142+61 was observed with the MAGIC I telescope in 2008 for ~25 h and 1E 2259+586 was observed with the MAGIC stereoscopic system in 2010 for ~14 h. The data were analyzed with the standard MAGIC analysis software. Neither magnetar was detected. Upper limits to the differential and integral flux above 200 GeV were computed using the Rolke algorithm. We obtain integral upper limits to the flux of 1.52*10^-12cm^-2 s^-1 and 2.7*10^-12cm^-2 s^-1 with a confidence level of 95% for 4U 0142+61 and 1E 2259+586, respectively. The resulting differential upper limits are presented together with X-ray data and upper limits in the GeV energy range.

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

Context. We present the discovery of very high energy (VHE, E > 100GeV) gamma-ray emission from the BL Lac object 1ES 1215+303 by the MAGIC telescopes and simultaneous multi-wavelength data in a broad energy range from radio to gamma-rays. Aims. We study the VHE gamma-ray emission from 1ES 1215+303 and its relation to the emissions in other wavelengths. Methods. Triggered by an optical outburst, MAGIC observed the source in January-February 2011 for 20.3 hrs. The target was monitored in the optical R-band by the KVA telescope that also performed optical polarization measurements. We triggered target of opportunity observations with the Swift satellite and obtained simultaneous and quasi-simultaneous data from the Fermi Large Area Telescope and from the Mets\"ahovi radio telescope. We also present the analysis of older MAGIC data taken in 2010. Results. The MAGIC observations of 1ES 1215+303 carried out in January-February 2011 resulted in the first detection of the source at VHE with a statistical significance of 9.4 sigma. Simultaneously, the source was observed in a high optical and X-ray state. In 2010 the source was observed in a lower state in optical, X-ray, and VHE, while the GeV gamma-ray flux and the radio flux were comparable in 2010 and 2011. The spectral energy distribution obtained with the 2011 data can be modeled with a simple one zone SSC model, but it requires extreme values for the Doppler factor or the electron energy distribution.

We present the results of a long M87 monitoring campaign in very high energy $\gamma$-rays with the MAGIC-I Cherenkov telescope. We aim to model the persistent non-thermal jet emission by monitoring and characterizing the very high energy $\gamma$-ray emission of M87 during a low state. A total of 150\,h of data were taken between 2005 and 2007 with the single MAGIC-I telescope, out of which 128.6\,h survived the data quality selection. We also collected data in the X-ray and \textit{Fermi}--LAT bands from the literature (partially contemporaneous). No flaring activity was found during the campaign. The source was found to be in a persistent low-emission state, which was at a confidence level of $7\sigma$. We present the spectrum between 100\,GeV and 2\,TeV, which is consistent with a simple power law with a photon index $\Gamma=2.21\pm0.21$ and a flux normalization at 300\,GeV of $(7.7\pm1.3) \times 10^{-8}$ TeV$^{-1}$ s$^{-1}$ m$^{-2}$. The extrapolation of the MAGIC spectrum into the GeV energy range matches the previously published \textit{Fermi}--LAT spectrum well, covering a combined energy range of four orders of magnitude with the same spectral index. We model the broad band energy spectrum with a spine layer model, which can satisfactorily describe our data.

The high frequency peaked BL Lac PKS 2155-304 with a redshift of z=0.116 was discovered in 1997 in the very high energy (VHE, E >100GeV) gamma-ray range by the University of Durham Mark VI gamma-ray Cherenkov telescope in Australia with a flux corresponding to 20% of the Crab Nebula flux. It was later observed and detected with high significance by the Southern Cherenkov observatory H.E.S.S. Detection from the Northern hemisphere is difficult due to challenging observation conditions under large zenith angles. In July 2006, the H.E.S.S. collaboration reported an extraordinary outburst of VHE gamma-emission. During the outburst, the VHE gamma-ray emission was found to be variable on the time scales of minutes and with a mean flux of ~7 times the flux observed from the Crab Nebula. Follow-up observations with the MAGIC-I standalone Cherenkov telescope were triggered by this extraordinary outburst and PKS 2155-304 was observed between 28 July to 2 August 2006 for 15 hours at large zenith angles. Here we present our studies on the behavior of the source after its extraordinary flare and an enhanced analysis method for data taken at high zenith angles. We developed improved methods for event selection that led to a better background suppression. The averaged energy spectrum we derived has a spectral index of -3.5 +/- 0.2 above 400GeV, which is in good agreement with the spectral shape measured by H.E.S.S. during the major flare on MJD 53944. Furthermore, we present the spectral energy distribution modeling of PKS 2155-304. With our observations we increased the duty cycle of the source extending the light curve derived by H.E.S.S. after the outburst. Finally, we find night-by-night variability with a maximal amplitude of a factor three to four and an intranight variability in one of the nights (MJD 53945) with a similar amplitude.

In 2008 AGILE and Fermi detected gamma-ray flaring activity from the unidentified EGRET source 3EG J1236+0457, recently associated with a flat spectrum radio quasar GB6 J1239+0443 at z=1.762. The optical counterpart of the gamma-ray source underwent a flux enhancement of a factor 15-30 in 6 years, and of ~10 in six months. We interpret this flare-up in terms of a transition from an accretion-disk dominated emission to a synchrotron-jet dominated one. We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum taken during a period of low radio and optical activity of the source. We estimated the mass of the central black hole using the width of the CIV emission line. In our work, we have also investigated SDSS archival optical photometric data and UV GALEX observations to estimate the thermal-disk emission contribution of GB6 J1239+0443. Our analysis of the gamma-ray data taken during the flaring episodes indicates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no statistically-relevant sign of absorption from the broad line region, suggesting that the blazar-zone is located beyond the broad line region. This result is confirmed by the modeling of the broad-band spectral energy distribution (well constrained by the available multiwavelength data) of the flaring activity periods and by the accretion disk luminosity and black hole mass estimated by us using archival data.

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) -- which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) -- represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306. Its evolutionary state is investigated by means of deep dedicated observations obtained with XMM-Newton, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi and AGILE missions. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula, has a soft spectrum with absorption features and no magnetospheric emission. The optical counterpart is fainter than V=27 and no gamma-ray emission is significantly detected. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P=18.6 ms. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.

Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii.

The very high energy (VHE) \gamma-ray source HESS J0632+057 has recently been confirmed to be a \gamma-ray binary. The optical counterpart is the Be star MWC 148, and a compact object of unknown nature orbits it every ~321 d with a high eccentricity of ~0.8. We monitored HESS J0632+057 with the stereoscopic MAGIC telescopes from 2010 October to 2011 March and detected significant VHE \gamma-ray emission during 2011 February, when the system exhibited an X-ray outburst. We find no \gamma-ray signal in the other observation periods when the system did not show increased X-ray flux. Thus HESS J0632+057 exhibits \gamma-ray variability on timescales of the order of one to two months possibly linked to the X-ray outburst that takes place about 100 days after the periastron passage. Furthermore our measurements provide for the first time the \gamma-ray spectrum down to about 140 GeV and indicate no turnover of the spectrum at low energies. We compare the properties of HESS J0632+057 with the similar \gamma-ray binary LS I +61 303, and discuss on the possible origin of the multi-wavelength emission of the source

We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 500 GeV can be described by a power law with a steep spectral index of \Gamma=-4.1+/-0.7stat+/-0.3syst, and the average flux above 100 GeV is F_{gamma}=(1.3+/-0.2stat+/-0.3syst) x 10^-11 cm^-2 s^-1. These results, combined with the power-law spectrum measured in the first two years of observations by the Fermi-LAT above 100 MeV, with a spectral index of Gamma ~ -2.1, strongly suggest the presence of a break or cut-off around tens of GeV in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Fermi-LAT and MAGIC. The derived flux upper limit F^{U.L.}_{gamma} (>300 GeV)=1.2 x 10^-12 cm^-2 s^-1 is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and February 2010, thus confirming the year time-scale variability of the source at VHE.

We use 73 h of stereoscopic data taken with the MAGIC telescopes to investigate the very high-energy (VHE) gamma-ray emission of the Crab pulsar. Our data show a highly significant pulsed signal in the energy range from 50 to 400 GeV in both the main pulse (P1) and the interpulse (P2) phase regions. We provide the widest spectra to date of the VHE components of both peaks, and these spectra extend to the energy range of satellite-borne observatories. The good resolution and background rejection of the stereoscopic MAGIC system allows us to cross-check the correctness of each spectral point of the pulsar by comparison with the corresponding (strong and well-known) Crab nebula flux. The spectra of both P1 and P2 are compatible with power laws with photon indices of 4.0 \pm 0.8 (P1) and 3.42 \pm 0.26 (P2), respectively, and the ratio P1/P2 between the photon counts of the two pulses is 0.54 \pm 0.12. The VHE emission can be understood as an additional component produced by the inverse Compton scattering of secondary and tertiary e\pm pairs on IR-UV photons.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

The W51 complex hosts the supernova remnant W51C which is known to interact with the molecular clouds in the star forming region W51B. In addition, a possible pulsar wind nebula CXO J192318.5+140305 was found likely associated with the supernova remnant. Gamma-ray emission from this region was discovered by Fermi/LAT (between 0.2 and 50 GeV) and H.E.S.S. (>1 TeV). The spatial distribution of the events could not be used to pinpoint the location of the emission among the pulsar wind nebula, the supernova remnant shell and/or the molecular cloud. However, the modeling of the spectral energy distribution presented by the Fermi/LAT collaboration suggests a hadronic emission mechanism. We performed observations of the W51 complex with the MAGIC telescopes for more than 50 hours. The good angular resolution in the medium (few hundred GeV) to high (above 1 TeV) energies allow us to perform morphological studies. We detect an extended emission of very-high-energy gamma rays, with a significance of 11 standard deviations. We extend the spectrum from the highest Fermi/LAT energies to \sim 5 TeV and find that it follows a single power law with an index of 2.58 \pm 0.07stat \pm 0.22syst . The main part of the emission coincides with the shocked cloud region, while we find a feature extending towards the pulsar wind nebula. The possible contribution of the pulsar wind nebula, assuming a point-like source, shows no dependence on energy and it is about 20% of the overall emission. The broad band spectral energy distribution can be explained with a hadronic model that implies proton acceleration above 100 TeV. This result, together with the morphology of the source, tentatively suggests that we observe ongoing acceleration of ions in the interaction zone between supernova remnant and cloud. These results shed light on the long-standing problem of the origin of galactic cosmic rays.

During its first 2 years of operation, the gamma-ray AGILE satellite accumulated an extensive dataset for the Galactic plane. The data have been monitored for transient sources and several gamma-ray sources were detected. Their variability and possible association were studied. In this talk we will focus on the results of extensive observations of the Carina Region during the time period 2007 July - 2009 January, for a total livetime of ~130 days. The region is extremely complex, hosting massive star formation, with the remarkable colliding wind binary Eta Carinae, massive star clusters and HII regions (e.g. NGC 3324, RCW49, Westerlund II) and a giant molecular cloud extending over 150 pc (between l=284.7 and l=289). The Carina Nebula itself is the largest and IR highest surface brightness nebula of the Southern emisphere. We monitored several gamma ray sources in the Carina Region. In particular we detect a gamma ray source (1AGL J1043-5931) consistent with the position of Eta Carinae and report a remarkable 2-days gamma-ray flaring episode from this source on 2008 Oct 11-13. If 1AGL J1043-5931 is associated with the Eta Car system, our data provides the long sought first detection above 100 MeV of a colliding wind binary.

At the core of the AGILE scientific instrument, designed to operate on a satellite, there is the Gamma Ray Imaging Detector (GRID) consisting of a Silicon Tracker (ST), a Cesium Iodide Mini-Calorimeter and an Anti-Coincidence system of plastic scintillator bars. The ST needs an on-ground calibration with a $\gamma$-ray beam to validate the simulation used to calculate the energy response function and the effective area versus the energy and the direction of the $\gamma$ rays. A tagged $\gamma$-ray beam line was designed at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali of Frascati (LNF), based on an electron beam generating $\gamma$ rays through bremsstrahlung in a position-sensitive target. The $\gamma$-ray energy is deduced by difference with the post-bremsstrahlung electron energy \cite{prest}-\cite{hasan}. The electron energy is measured by a spectrometer consisting of a dipole magnet and an array of position sensitive silicon strip detectors, the Photon Tagging System (PTS). The use of the combined BTF-PTS system as tagged photon beam requires understanding the efficiency of $\gamma$-ray tagging, the probability of fake tagging, the energy resolution and the relation of the PTS hit position versus the $\gamma$-ray energy. This paper describes this study comparing data taken during the AGILE calibration occurred in 2005 with simulation.

We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 14.2 hours of good quality very high energy (VHE) $\gamma$-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy $\gamma$-ray (HE, 0.1 GeV-100 GeV) data from the Fermi satellite. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE $\gamma$-rays at a statistical significance of 5.6 $\sigma$. A soft VHE spectrum with a photon index of -3.2 $\pm$ 0.6 was determined. No significant short term flux variations were found. We model the spectral energy distribution using a one-zone SSC-model, which can successfully describe our data.

We present the results of five years (2005-2009) of MAGIC observations of the BL Lac object PG 1553+113 at very high energies (VHEs, E > 100 GeV). Power law fits of the individual years are compatible with a steady mean photon index \Gamma = 4.27 $\pm$ 0.14. In the last three years of data, the flux level above 150 GeV shows a clear variability (probability of constant flux < 0.001%). The flux variations are modest, lying in the range from 4% to 11% of the Crab Nebula flux. Simultaneous optical data also show only modest variability that seems to be correlated with VHE gamma ray variability. We also performed a temporal analysis of (all available) simultaneous Fermi/LAT data of PG 1553+113 above 1 GeV, which reveals hints of variability in the 2008-2009 sample. Finally, we present a combination of the mean spectrum measured at very high energies with archival data available for other wavelengths. The mean spectral energy distribution can be modeled with a one-zone Synchrotron Self Compton (SSC) model, which gives the main physical parameters governing the VHE emission in the blazar jet.

We report the AGILE gamma-ray observations and the results of the multiwavelength campaigns on seven flaring blazars detected by the mission: During two multiwavelength campaigns, we observed gamma-ray activity from two Flat Spectrum Radio Quasars of the Virgo region, e.g. 3C 279 and 3C 273 (the latter being the first extragalactic source simultaneously observed with the gamma-ray telescope and the hard X ray imager of the mission). Due to the large FOV of the AGILE/GRID instrument, we achieved an almost continuous coverage of the FSRQ 3C 454.3. The source showed flux above 10E-6 photons/cm2/s (E > 100 MeV) and showed day by day variability during all the AGILE observing periods. In the EGRET era, the source was found in high gamma-ray activity only once. An other blazar, PKS 1510-089 was frequently found in high gamma-ray activity. S5 0716+71, an intermediate BL Lac object, exhibited a very high gamma-ray activity and fast gamma-ray variability during a period of intense optical activity. We observed high gamma-ray activity from W Comae, a BL Lac object, and Mrk 421, an high energy peaked BL Lac object. For this source, a multiwavelength campaign from optical to TeV has been performed.

Since its early phases of operation, the AGILE mission is successfully observing Gamma Ray Bursts (GRBs) in the hard X-ray band with the SuperAGILE imager and in the MeV range with the Mini-Calorimeter. Up to now, three firm GRB detections were obtained above 25 MeV and some bursts were detected with lower statistical confidence in the same energy band. When a GRB is localized, either by SuperAGILE or Swift/BAT or INTEGRAL/IBIS or Fermi/GBM or IPN, inside the field of view of the Gamma Ray Imager of AGILE, a detection is searched for in the gamma ray band or an upper limit is provided. A promising result of AGILE is the detection of very short gamma ray transients, a few ms in duration and possibly identified with Terrestrial Gamma-ray Flashes. In this paper we show the current status of the observation of Gamma Ray Bursts and Terrestrial Gamma-ray Flashes with AGILE.

The AGILE scientific instrument has been calibrated with a tagged $\gamma$-ray beam at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali di Frascati (LNF). The goal of the calibration was the measure of the Point Spread Function (PSF) as a function of the photon energy and incident angle and the validation of the Monte Carlo (MC) simulation of the silicon tracker operation. The calibration setup is described and some preliminary results are presented.

We present very high energy (VHE, E > 100 GeV) {\gamma}-ray observations of the {\gamma}-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3{\sigma} {\gamma}-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 GeV is F(E>300 GeV)=(1.4 +- 0.3stat +- 0.4syst) * 10^{-12} cm^{-2} s^{-1}, which corresponds to about 1.3% of the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61 303 in the orbital phase interval 0.6--0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude compared to our previous measurements between 2005 September and 2008 January. This provides evidence for a new low emission state in LS I +61 303. We find that the change to the low flux state cannot be solely explained by an increase of photon-photon absorption around the compact star.

Cygnus X-3 (Cyg X-3) is a well-known microquasar producing variable emission at all wavelengths. Cyg X-3 is a prominent X-ray binary producing relativistic jets, and studying its high energy emission is crucial for the understanding of the fundamental acceleration processes in accreting compact objects. Aims. Our goal is to study extreme particle acceleration and {\gamma}-ray production above 100 MeV during special spectral states of Cyg X- 3 usually characterized by a low hard X-ray flux and enhanced soft X-ray states. We observed Cyg X-3 with the AGILE satellite in extended time intervals from 2009 Jun.-Jul., and 2009 Nov.-2010 Jul. We report here the results of the AGILE {\gamma}-ray monitoring of Cyg X-3 as well as the results from extensive multiwavelength campaigns involving radio (RATAN-600, AMI-LA and Mets\"{a}hovi Radio Observatories) and X-ray monitoring data (XTE and Swift). We detect a series of repeated {\gamma}-ray flaring activity from Cyg X-3 that correlate with the soft X-ray states and episodes of decreasing or non-detectable hard X-ray emission. Furthermore, we detect {\gamma}-ray enhanced emission that tends to be associated with radio flares greater than 1 Jy at 15 GHz, confirming a trend already detected in previous observations. The source remained active above 100 MeV for an extended period of time (almost 1.5 months in 2009 Jun.-Jul. and 1 month in 2010 May). We study in detail the short timescale {\gamma}-ray flares that occurred before or near the radio peaks. Our results confirm the transient nature of the extreme particle acceleration from the microquasar Cyg X-3. A series of repeated {\gamma}-ray flares shows correlations with radio and X-ray emission confirming a well established trend of emission. We compare our results with Fermi-LAT and MAGIC TeV observations of Cyg X-3.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50-60 GeV. Advanced stereo analysis techniques allow MAGIC to achieve a sensitivity as good as (0.76 +/- 0.03)% of the Crab Nebula flux in 50 h of observations above 290 GeV. The angular resolution at those energies is better than ~0.07 degree. We also perform a detailed study of possible systematic effects which may influence the analysis of the data taken with the MAGIC telescopes.

Compilation of the papers contributed by the MAGIC collaboration to the 32nd International Cosmic Ray Conference (ICRC 2011), which took place between August 11 and 18, 2011 in Beijing, China. The papers are sorted in 6 categories: Overview and Highlight papers; Instrument, software and techniques; Galactic sources; Extragalactic sources; Multi-wavelength and joint campaigns; Fundamental physics, dark matter and cosmic rays.

We report on the observation of $\gamma$-rays above 25\,GeV from the Crab pulsar (PSR B0532+21) using the MAGIC I telescope. Two data sets from observations during the winter period 2007/2008 and 2008/2009 are used. In order to discuss the spectral shape from 100\,MeV to 100\,GeV, one year of public {\it Fermi} Large Area Telescope ({\it Fermi}-LAT) data are also analyzed to complement the MAGIC data. The extrapolation of the exponential cutoff spectrum determined with the Fermi-LAT data is inconsistent with MAGIC measurements, which requires a modification of the standard pulsar emission models. In the energy region between 25 and 100\,GeV, the emission in the P1 phase (from -0.06 to 0.04, location of the main pulse) and the P2 phase (from 0.32 to 0.43, location of the interpulse) can be described by power laws with spectral indices of $-3.1 \pm 1.0_{stat} \pm 0.3_{syst}$ and $-3.5 \pm 0.5_{stat} \pm 0.3_{syst}$, respectively. Assuming an asymmetric Lorentzian for the pulse shape, the peak positions of the main pulse and the interpulse are estimated to be at phases $-0.009 \pm 0.007$ and $0.393 \pm 0.009$, while the full widths at half maximum are $0.025 \pm 0.008$ and $0.053 \pm 0.015$, respectively.

We present the results of 2.5 years of AGILE observations of PSR B1509-58 and of the same interval of Fermi observations. The modulation significance of AGILE light curve above 30 MeV is at a 5 sigma confidence level and the light curve shows a broad asymmetric first peak reaching its maximum $0.39 \pm 0.02$ cycles after the radio peak plus a second peak at$0.94 \pm 0.03$. The gamma-ray spectral energy distribution of pulsed flux is well described by a power-law (photon index $\alpha=1.87\pm0.09$) with a remarkable cutoff at below 100 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce AGILE light curves and spectra match the polarization measurements.

The observation of Gamma Ray Bursts (GRBs) in the gamma-ray band has been advanced by the AGILE and Fermi satellites after the era of the Compton Gamma-Ray Observatory. AGILE and Fermi are showing that the GeV-bright GRBs share a set of common features, particularly the high fluence from the keV up to the GeV energy bands, the high value of the minimum Lorentz factor, the presence of an extended emission of gamma-rays, often delayed with respect to lower energies, and finally the possible presence of multiple spectral components. GRB 100724B, localised in a joint effort by Fermi and the InterPlanetary Newtork, is the brightest burst detected in gamma-rays so far by AGILE. Characteristic features of GRB 100724B are the simultaneous emissions at MeV and GeV, without delayed onset nor time lag as shown by the analysis of the cross correlation function, and the significant spectral evolution in hard X-rays over the event duration. In this paper we show the analysis of the AGILE data of GRB 100724B and we discuss its features in the context of the bursts observed so far in gamma-rays and the recently proposed models.

We report on the extremely intense and fast gamma-ray are above 100 MeV detected by AGILE from the Crab Nebula in mid-April 2011. This event is the fourth of a sequence of reported major gamma-ray flares produced by the Crab Nebula in the period 2007/mid-2011. These events are attributed to strong radiative and plasma instabilities in the inner Crab Nebula, and their properties are crucial for theoretical studies of fast and efficient particle acceleration up to 10^15 eV. Here we study the very rapid flux and spectral evolution of the event that reached on April 16, 2011 the record-high peak flux of F = (26 +/- 5) x 10^-6 ph cm^-2 s^-1 with a risetime timescale that we determine to be in the range 6-10 hrs. The peak flaring gamma-ray spectrum reaches a distinct maximum near 500 MeV with no substantial emission above 1 GeV. The very rapid risetime and overall evolution of the Crab Nebula are strongly constrain the acceleration mechanisms and challenge MHD models. We briefly discuss the theoretical implications of our observations.

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

The acceleration of particles up to GeV or higher energies in microquasars has been the subject of considerable theoretical and observational efforts in the past few years. Sco X-1 is a microquasar from which evidence of highly energetic particles in the jet has been found when it is in the so-called Horizontal Branch (HB), a state when the radio and hard X-ray fluxes are higher and a powerful relativistic jet is present. Here we present the first very high energy gamma-ray observations of Sco X-1 obtained with the MAGIC telescopes. An analysis of the whole dataset does not yield a significant signal, with 95% CL flux upper limits above 300 GeV at the level of 2.4x10^{-12} ph/cm^2/s. Simultaneous RXTE observations were conducted to search for TeV emission during particular X-ray states of the source. A selection of the gamma-ray data obtained during the HB based on the X-ray colors did not yield a signal either, with an upper limit of 3.4x10^{-12} ph/cm^2/s. These upper limits place a constraint on the maximum TeV luminosity to non-thermal X-ray luminosity of L_{VHE}/L_{ntX}<0.02, that can be related to a maximum TeV luminosity to jet power ratio of L_{VHE}/L_{j}<10^{-3}. Our upper limits indicate that the underlying high-energy emission physics in Sco X-1 must be inherently different from that of the hitherto detected gamma-ray binaries.

We present the results of 2.5 years AGILE observations of PSR B1509-58. The modulation significance of the light-curve above 30 MeV is at a sigma confidence level and the light-curve is similar to those found earlier up to 30 MeV by \comp: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of pulsed flux is well described by a power-law (photon index alpha=1.87 +/-0.09) with a remarkable cutoff at E_c=81+/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce AGILE light-curves and spectra match the polarization measurements.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

Very High Energy (VHE) gamma-ray emission from the flat spectrum radio quasar (FSRQ) PKS 1222+21 (4C 21.35, z=0.432) was detected with the MAGIC Cherenkov telescopes during a short observation (~0.5 hr) performed on 2010 June 17. The MAGIC detection coincides with high energy MeV/GeV gamma-ray activity measured by the Large Area Telescope (LAT) on board the Fermi satellite. The VHE spectrum measured by MAGIC extends from about 70 GeV up to at least 400 GeV and can be well described by a power law dN/dE \propto E^-Gamma with a photon index Gamma= 3.75+/-0.27stat +/-0.2syst. The averaged integral flux above 100 GeV is (4.56+/-0.46)x10^(-10) cm^-2 s^-1 (~1 Crab Nebula flux). The VHE flux measured by MAGIC varies significantly within the 30 min exposure implying a flux doubling time of about 10 min. The VHE and MeV/GeV spectra, corrected for the absorption by the extragalactic background light (EBL), can be described by a single power law with photon index 2.72+/-0.34 between 3 GeV and 400 GeV, and is consistent with emission belonging to a single component in the jet. The absence of a spectral cutoff constrains the gamma-ray emission region outside the Broad Line Region, which would otherwise absorb the VHE gamma-rays. Together with the detected fast variability, this challenges present emission models from jets in FSRQ. Moreover, the combined Fermi/LAT and MAGIC spectral data yield constraints on the density of the Extragalactic Background Light in the UV-optical to near-infrared range that are compatible with recent models.

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

The well known Crab Nebula is at the center of the SN1054 supernova remnant. It consists of a rotationally-powered pulsar interacting with a surrounding nebula through a relativistic particle wind. The emissions originating from the pulsar and nebula have been considered to be essentially stable. Here we report the detection of strong gamma-ray (100 MeV-10 GeV) flares observed by the AGILE satellite in September, 2010 and October, 2007. In both cases, the unpulsed flux increased by a factor of 3 compared to the non-flaring flux. The flare luminosity and short timescale favor an origin near the pulsar, and we discuss Chandra Observatory X-ray and HST optical follow-up observations of the nebula. Our observations challenge standard models of nebular emission and require power-law acceleration by shock-driven plasma wave turbulence within a ~1-day timescale.

The AGILE gamma-ray satellite provides large sky exposure levels ($\geq 10^9$ cm$^2$ s per year on the Galactic Plane) with sensitivity peaking at $E\sim$100 MeV where the bulk of pulsar energy output is typically released. Its $\sim$1 $\mu$s absolute time tagging capability makes it perfectly suited for the study of gamma-ray pulsars. AGILE collected a large number of gamma-ray photons from EGRET pulsars ($\geq$40,000 pulsed counts for Vela) in two years of observations unveiling new interesting features at sub-millisecond level in the pulsars' high-energy light-curves, gamma-ray emission from pulsar glitches and Pulsar Wind Nebulae. AGILE detected about 20 nearby and energetic pulsars with good confidence through timing and/or spatial analysis. Among the newcomers we find pulsars with very high rotational energy losses, such as the remarkable PSR B1509--58 with a magnetic field in excess of 10$^{13}$ Gauss, and PSR J2229+6114 providing a reliable identification for the previously unidentified EGRET source 3EG 2227+6122. Moreover, the powerful millisecond pulsar B1821--24, in the globular cluster M28, is detected during a fraction of the observations.

The neutron star spin down imposes a balance between the energy and angular momentum (E/L) losses of a pulsar. This translates into constraints on the region of emission. The E/L balance may be cogent in discriminating among the models of gamma-ray production. Hence, models which require the release of the entire luminosity at the polar caps should be excluded. Also models where the radiative zone is well inside the speed of light radius have some difficulties. It is argued that a local unbalance of E/L should generate a global instability of the magnetosphere, possibly quenching the emission at the polar caps.

Identification of gamma-ray-emitting Galactic sources is a long-standing problem in astrophysics. One such source, 1AGL J2022+4032, coincident with the interior of the radio shell of the supernova remnant Gamma Cygni (SNR G78.2+2.1) in the Cygnus Region, has recently been identified by Fermi as a gamma-ray pulsar, LAT PSR J2021+4026. We present long-term observations of 1AGL J2022+4032 with the AGILE gamma-ray telescope, measuring its flux and light curve. We compare the light curve of 1AGL J2022+4032 with that of 1AGL J2021+3652 (PSR J2021+3651), showing that the flux variability of 1AGL J2022+4032 appears to be greater than the level predicted from statistical and systematic effects and producing detailed simulations to estimate the probability of the apparent observed variability. We evaluate the possibility that the gamma-ray emission may be due to the superposition of two or more point sources, some of which may be variable, considering a number of possible counterparts. We consider the possibility of a nearby X-ray quiet microquasar contributing to the flux of 1AGL J2022+4032 to be more likely than the hypotheses of a background blazar or intrinsic gamma-ray variabilty of LAT PSR J2021+4026.

Terrestrial Gamma-Ray Flashes (TGFs) are very short bursts of high energy photons and electrons originating in Earth's atmosphere. We present here a localization study of TGFs carried out at gamma-ray energies above 20 MeV based on an innovative event selection method. We use the AGILE satellite Silicon Tracker data that for the first time have been correlated with TGFs detected by the AGILE Mini-Calorimeter. We detect 8 TGFs with gamma-ray photons of energies above 20 MeV localized by the AGILE gamma-ray imager with an accuracy of 5-10 degrees at 50 MeV. Remarkably, all TGF-associated gamma rays are compatible with a terrestrial production site closer to the sub-satellite point than 400 km. Considering that our gamma rays reach the AGILE satellite at 540 km altitude with limited scattering or attenuation, our measurements provide the first precise direct localization of TGFs from space.

We present the results of new AGILE observations of PSR B1509-58 performed over a period of $\sim$2.5 years following the detection obtained with preliminary data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier up to 30 MeV by COMPTEL: a broad asymmetric first peak reaching its maximum $0.39 \pm 0.02$ cycles after the radio peak plus a second peak at $0.94 \pm 0.03$. The $\gamma$-ray spectral energy distribution of pulsed flux is well described by a power-law (photon index $\alpha=1.87\pm0.09$) with a remarkable cutoff at $E_c=81\pm 20$ MeV, representing the softest spectrum observed among $\gamma$-ray pulsars so far. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce AGILE lightcurves and spectra match the polarization measurements.

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

We present the gamma-ray data of the extraordinary flaring activity above 100 MeV from the flat spectrum radio quasar 3C 454.3 detected by AGILE during the month of December 2009. 3C 454.3, that has been among the most active blazars of the FSRQ type since 2007, was detected in the gamma-ray range with a progressively rising flux since November 10, 2009. The gamma-ray flux reached a value comparable with that of the Vela pulsar on December 2, 2009. Remarkably, between December 2 and 3, 2009 the source more than doubled its gamma-ray emission and became the brightest gamma-ray source in the sky with a peak flux of F_{\gamma,p} = (2000 \pm 400) x 10^-8 ph cm^-2 s^-1 for a 1-day integration above 100 MeV. The gamma-ray intensity decreased in the following days with the source flux remaining at large values near F \simeq (1000 \pm 200) x 10^-8 ph cm^-2 s^-1 for more than a week. This exceptional gamma-ray flare dissipated among the largest ever detected intrinsic radiated power in gamma-rays above 100 MeV (L_{\gamma, source, peak} \simeq 3 x 10^46 erg s^-1, for a relativistic Doppler factor of {\delta} \simeq 30). The total isotropic irradiated energy of the month-long episode in the range 100 MeV - 3 GeV is E_{\gamma,iso} \simeq 10^56 erg. We report the intensity and spectral evolution of the gamma-ray emission across the flaring episode. We briefly discuss the important theoretical implications of our detection.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

Supernova remnants (SNRs) are believed to be the main sources of Galactic cosmic rays. Molecular clouds associated with SNRs can produce gamma-ray emission through the interaction of accelerated particles with the concentrated gas. The middle aged SNR W28, for its associated system of dense molecular clouds, provides an excellent opportunity to test this hypothesis. We present the AGILE/GRID observations of SNR W28, and compare them with observations at other wavelengths (TeV and 12CO J=1-->0 molecular line emission). The gamma-ray flux detected by AGILE from the dominant source associated with W28 is (14 +- 5) 10^-8 ph cm^-2 s^-1 for E > 400 MeV. This source is positionally well correlated with the TeV emission observed by the HESS telescope. The local variations of the GeV to TeV flux ratio suggest a difference between the CR spectra of the north-west and south molecular cloud complexes. A model based on a hadronic-induced interaction and diffusion with two molecular clouds at different distances from the W28 shell can explain both the morphological and spectral features observed by AGILE in the MeV-GeV energy range and by the HESS telescope in the TeV energy range. The combined set of AGILE and H.E.S.S. data strongly support a hadronic model for the gamma-ray production in W28.

We present the observation of Cyg X-1 in hard spectral state performed during the AGILE Science Verification Phase and Observing Cycle 1 in hard X-rays (with SuperAGILE) and gamma rays (with the GRID) and lasting for about 160 days with a live time of $\sim 6$ Ms. We investigate the variability of Cyg X-1 in hard X-rays at different timescales, from $\sim 300$ s up to one day, and we apply different tools of timing analysis, such as the autocorrelation function, the first order structure function and the Lomb-Scargle periodogram, to our data (from SuperAGILE) and to the simultaneous data in soft X-rays (from RXTE/ASM). We conclude our investigation with a search for emission in the energy range above 100 MeV with the maximum likelihood technique. In the hard X-ray band the flux of Cyg X-1 shows its typical erratic fluctuations at all timescales with variations of about a factor of two that do not affect significantly the shape of the energy spectrum. From the first order structure function we find that the X-ray emission of Cyg X-1 is characterized by \textit{antipersistence}, indication of a negative feedback mechanism at work. In the gamma ray data a statistically significant point-like source at the position of Cyg X-1 is not found and the upper limit on the flux is $\mathrm{5 \times 10^{-8} \; ph \; cm^{-2} \; s^{-1}}$, over the whole observation (160 days). Finally we compare our upper limit in gamma rays with the expectation of various models of the Cyg X-1 emission, of both hadronic and leptonic origin, in the GeV -- TeV band. The time history of Cyg X-1 in the hard X-ray band over 13 months (not continuous) is shown. Different tools of analysis do not provide fully converging results of the characteristic timescales in the system, suggesting that the timescales found in the structure function are not intrinsic to the physics of the source.