Source author record

L. A. Antonelli

L. A. Antonelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

169works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

169 published item(s)

preprint2022arXiv

Accretion geometry of the neutron star low mass X-ray binary Cyg X-2 from X-ray polarization measurements

We report spectro-polarimetric results of an observational campaign of the bright neutron star low-mass X-ray binary Cyg X-2 simultaneously observed by IXPE, NICER and INTEGRAL. Consistently with previous results, the broad-band spectrum is characterized by a lower-energy component, attributed to the accretion disc with $kT_{\rm in} \approx$ 1 keV, plus unsaturated Comptonization in thermal plasma with temperature $kT_{\rm e} = 3$ keV and optical depth $τ\approx 4$, assuming a slab geometry. We measure the polarization degree in the 2-8 keV band $P=1.8 \pm 0.3$ per cent and polarization angle $ϕ= 140^{\circ} \pm 4^{\circ}$, consistent with the previous X-ray polarimetric measurements by OSO-8 as well as with the direction of the radio jet which was earlier observed from the source. While polarization of the disc spectral component is poorly constrained with the IXPE data, the Comptonized emission has a polarization degree $P =4.0 \pm 0.7$ per cent and a polarization angle aligned with the radio jet. Our results strongly favour a spreading layer at the neutron star surface as the main source of the polarization signal. However, we cannot exclude a significant contribution from reflection off the accretion disc, as indicated by the presence of the iron fluorescence line.

preprint2022arXiv

ASTRI Mini-Array Core Science at the Observatorio del Teide

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily. The ASTRI Mini-Array will surpass the current Cherenkov telescope array differential sensitivity above a few tera-electronvolt (TeV), extending the energy band well above hundreds of TeV. This will allow us to explore a new window of the electromagnetic spectrum, by convolving the sensitivity performance with excellent angular and energy resolution figures. In this paper we describe the Core Science that we will address during the first four years of operation, providing examples of the breakthrough results that we will obtain when dealing with current open questions, such as the acceleration of cosmic rays, cosmology and fundamental physics and the new window, for the TeV energy band, of the time-domain astrophysics.

preprint2022arXiv

Extragalactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observational feasibility of extragalactic sources and on astrophysical processes that best complement and expand the ASTRI Mini-Array core science, presenting the most relevant examples that are at reach of detection over long-term time scales and whose observation can provide breakthrough achievements in the very-high energy extragalactic science. Such examples cover a wide range of $γ$-ray emitters, including the study of AGN low states in the multi-TeV energy range, the possible detection of Seyfert galaxies with long exposures and the searches of dark matter lines above 10 TeV. Simulations of the presented objects show that the instrument performance will be competitive at multi-TeV energies with respect to current arrays of Cherenkov telescopes.

preprint2022arXiv

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

preprint2021arXiv

Combined searches for dark matter in dwarf spheroidal galaxies observed with the MAGIC telescopes, including new data from Coma Berenices and Draco

Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particles (WIMPs) in dSphs. Results on the observations of Segue 1 and Ursa Major II dSphs have already been published and include some of the most stringent upper limits (ULs) on the velocity-averaged cross-section $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP annihilation from observations of dSphs. In this work, we report on the analyses of 52.1 h of data of Draco dSph and 49.5 h of Coma Berenices dSph observed with the MAGIC telescopes in 2018 and in 2019 respectively. No hint of a signal has been detected from either of these targets and new constraints on the $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP candidates have been derived. In order to improve the sensitivity of the search and reduce the effect of the systematic uncertainties due to the $J$-factor estimates, we have combined the data of all dSphs observed with the MAGIC telescopes. Using 354.3 h of dSphs good quality data, 95 % CL ULs on $\langle σ_{\mathrm{ann}} v \rangle$ have been obtained for 9 annihilation channels. For most of the channels, these results reach values of the order of $10^{-24} $cm$^3$/s at ${\sim}1$ TeV and are the most stringent limits obtained with the MAGIC telescopes so far.

preprint2021arXiv

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.

preprint2021arXiv

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

preprint2021arXiv

The supernova of the MAGIC GRB190114C

We observed GRB190114C (redshift z = 0.4245), the first GRB ever detected at TeV energies, at optical and near-infrared wavelengths with several ground-based telescopes and the Hubble Space Telescope, with the primary goal of studying its underlying supernova, SN2019jrj. The monitoring spanned the time interval between 1.3 and 370 days after the burst, in the observer frame. We find that the afterglow emission can be modelled with a forward shock propagating in a uniform medium modified by time-variable extinction along the line of sight. A jet break could be present after 7 rest-frame days, and accordingly the maximum luminosity of the underlying SN ranges between that of stripped-envelope corecollapse supernovae (SNe) of intermediate luminosity, and that of the luminous GRB-associated SN2013dx. The observed spectral absorption lines of SN2019jrj are not as broad as in classical GRB-SNe, and are rather more similar to those of less-luminous core-collapse SNe. Taking the broad-lined stripped-envelope core-collapse SN2004aw as an analogue, we tentatively derive the basic physical properties of SN2019jrj. We discuss the possibility that a fraction of the TeV emission of this source might have had a hadronic origin and estimate the expected high-energy neutrino detection level with IceCube.

preprint2020arXiv

A search for dark matter in Triangulum II with the MAGIC telescopes

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

preprint2020arXiv

AGILE Observations of Two Repeating Fast Radio Bursts with Small Intrinsic Dispersion Measures

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 $γ$-ray emission for time intervals up to 100 days, obtaining 3$\,σ$ upper limits (ULs) for the average isotropic luminosity above 50 MeV, $L_{γ,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_{γ,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.

preprint2020arXiv

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

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

preprint2020arXiv

An X-Ray Burst from a Magnetar Enlightening the Mechanism of Fast Radio Bursts

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.

preprint2020arXiv

Bounds on Lorentz invariance violation from MAGIC observation of GRB 190114C

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

preprint2020arXiv

Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016

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

preprint2020arXiv

First detection of the Crab Nebula at TeV energies with a Cherenkov telescope in a dual-mirror Schwarzschild-Couder configuration: the ASTRI-Horn telescope

We report on the first detection of very high-energy (VHE) gamma-ray emission from the Crab Nebula by a Cherenkov telescope in dual-mirror Schwarzschild-Couder (SC) configuration. The result has been achieved by means of the 4 m size ASTRI-Horn telescope, operated on Mt. Etna (Italy) and developed in the context of the Cherenkov Telescope Array Observatory preparatory phase. The dual-mirror SC design is aplanatic and characterized by a small plate scale, allowing us to implement large field of view cameras with small-size pixel sensors and a high compactness. The curved focal plane of the ASTRI camera is covered by silicon photo-multipliers (SiPMs), managed by an unconventional front-end electronics based on a customized peak-sensing detector mode. The system includes internal and external calibration systems, hardware and software for control and acquisition, and the complete data archiving and processing chain. The observations of the Crab Nebula were carried out in December 2018, during the telescope verification phase, for a total observation time (after data selection) of 24.4 h, equally divided into on- and off-axis source exposure. The camera system was still under commissioning and its functionality was not yet completely exploited. Furthermore, due to recent eruptions of the Etna Volcano, the mirror reflection efficiency was reduced. Nevertheless, the observations led to the detection of the source with a statistical significance of 5.4 sigma above an energy threshold of ~3 TeV. This result provides an important step towards the use of dual-mirror systems in Cherenkov gamma-ray astronomy. A pathfinder mini-array based on nine large field-of-view ASTRI-like telescopes is under implementation.

preprint2020arXiv

Gamma-Ray and X-Ray Observations of the Periodic-Repeater FRB 180916 During Active Phases

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 $γ$-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 $γ$-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 $τ_d$ to satisfy the relation $R_{m,6}^3 B_{m,16}^2 τ_{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 $τ_{d,8}$ in units of $10^8$ s.

preprint2020arXiv

MAGIC observations of the nearby short gamma-ray burst GRB 160821B

The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only $\sim$20 events detected in the GeV band, and none in the TeV band. GRB~160821B is one of the nearest short GRBs known at $z=0.162$. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 seconds after the burst trigger. Evidence of a gamma-ray signal is found above $\sim$0.5 TeV at a significance of $\sim3\,σ$ during observations that lasted until 4 hours after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.

preprint2020arXiv

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

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

preprint2020arXiv

Monitoring of the radio galaxy M87 during a low emission state from 2012 to 2015 with MAGIC

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

preprint2020arXiv

New hard-TeV extreme blazars detected with the MAGIC telescopes

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

preprint2020arXiv

Observation of inverse Compton emission from a long $γ$-ray burst

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

preprint2020arXiv

Statistics of VHE gamma-Rays in Temporal Association with Radio Giant Pulses from the Crab Pulsar

The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows we find upper limits at 95% confidence level on an increase in VHE gamma-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges respectively. This puts upper limits on the flux increase during a radio GP of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.}

preprint2020arXiv

Studying the nature of the unidentified gamma-ray source HESS J1841-055 with the MAGIC telescopes

We investigate the physical nature and origin of the gamma-ray emission from the extended source HESS J1841-055 observed at TeV and GeV energies. We observed HESS J1841-055 at TeV energies for a total effective time of 43 hours with the MAGIC telescopes, in 2012 and 2013. Additionally, we analysed the GeV counterpart making use of about 10 years of Fermi-LAT data. Using both Fermi-LAT and MAGIC, we study both the spectral and energy-dependent morphology of the source for almost four decades of energy. The origin of the gamma-ray emission from this region is investigated using multi-waveband information on sources present in this region, suggested to be associated with this unidentified gamma-ray source. We find that the extended emission at GeV-TeV energies is best described by more than one source model. We also perform the first energy-dependent analysis of the HESS J1841-055 region at GeV-TeV. We find that the emission at lower energies comes from a diffuse or extended component, while the major contribution of gamma rays above 1 TeV arises from the southern part of the source. Moreover, we find that a significant curvature is present in the combined observed spectrum of MAGIC and Fermi-LAT. The first multi-wavelength spectral energy distribution of this unidentified source shows that the emission at GeV-TeV energies can be well explained with both leptonic and hadronic models. For the leptonic scenario, bremsstrahlung is the dominant emission compared to inverse Compton. On the other hand, for the hadronic model, gamma-ray resulting from the decay of neutral pions ($π^0$) can explain the observed spectrum. The presence of dense molecular clouds overlapping with HESS J1841-055 makes both bremsstrahlung and $π^0$-decay processes the dominant emission mechanisms for the source.

preprint2020arXiv

Teraelectronvolt emission from the $γ$-ray burst GRB 190114C

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

preprint2020arXiv

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

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

preprint2020arXiv

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

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

preprint2020arXiv

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

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

preprint2019arXiv

A Fast Very High Energy $γ$-ray Flare from BL Lacertae during a Period of Multiwavelength activity in June 2015

The mechanisms producing fast variability of the $γ$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $γ$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $γ$-ray bands. In this paper we present the MAGIC VHE $γ$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $γ$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $γ$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $γ$-rays is very similar to two other observed VHE $γ$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $γ$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $γ$-ray flare, but all have some tension with the multiwavelength observations.

preprint2019arXiv

Constraints on gamma-ray and neutrino emission from NGC 1068 with the MAGIC telescopes

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

preprint2018arXiv

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$σ$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $σ$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

preprint2018arXiv

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

preprint2018arXiv

The blazar TXS 0506+056 associated with a high-energy neutrino: insights into extragalactic jets and cosmic ray acceleration

A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3σ$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $γγ$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton up-scattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV.

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

preprint2017arXiv

A cut-off in the TeV gamma-ray spectrum of the SNR Cassiopeia A

It is widely believed that the bulk of the Galactic cosmic rays are accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between December 2014 and October 2016 we observed Cas A with the MAGIC telescopes, accumulating 158 hours of good-quality data. We derived the spectrum of the source from 100 GeV to 10 TeV. We also analysed $\sim$8 years of $Fermi$-LAT to obtain the spectral shape between 60 MeV and 500 GeV. The spectra measured by the LAT and MAGIC telescopes are compatible within the errors and show a clear turn off (4.6 $σ$) at the highest energies, which can be described with an exponential cut-off at $E_c = 3.5\left(^{+1.6}_{-1.0}\right)_{\textit{stat}} \left(^{+0.8}_{-0.9}\right)_{\textit{sys}}$ TeV. The gamma-ray emission from 60 MeV to 10 TeV can be attributed to a population of high-energy protons with spectral index $\sim$2.2 and energy cut-off at $\sim$10 TeV. This result indicates that Cas A is not contributing to the high energy ($\sim$PeV) cosmic-ray sea in a significant manner at the present moment. A one-zone leptonic model fails to reproduce by itself the multi-wavelength spectral energy distribution. Besides, if a non-negligible fraction of the flux seen by MAGIC is produced by leptons, the radiation should be emitted in a region with a low magnetic field (B$\lessapprox$100$μ$G) like in the reverse shock.)

preprint2017arXiv

Constraining Lorentz invariance violation using the Crab Pulsar emission observed up to TeV energies by MAGIC

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the MAGIC collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent~CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, found to worsen the statistical limits by about 36--42\%. Our constraints would have resulted much more competitive if the intrinsic pulse shape of the pulsar between 200GeV and 400GeV was understood in sufficient detail and allowed inclusion of events well below 400GeV.

preprint2017arXiv

First Multi-wavelength Campaign on the Gamma-ray-loud Active Galaxy IC 310

The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars.

preprint2017arXiv

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab Nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for several nights during the period June 18-27, for more than 10 hours. One hour of observation was conducted simultaneously to a giant 22 GHz radio flare and a hint of signal at GeV energies seen by Fermi-LAT. The MAGIC observations did not show significant emission in any of the analysed time intervals. The derived flux upper limit, in the energy range 200--1250 GeV, is 4.8$\times 10^{-12}$ ph cm$^{-2}$ s$^{-1}$. We estimate the gamma-ray opacity during the flaring period, which along with our non-detection, points to an inefficient acceleration in the V404\,Cyg jets if VHE emitter is located further than $1\times 10^{10}$ cm from the compact object.

preprint2017arXiv

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

preprint2017arXiv

Performance of the MAGIC telescopes under moonlight

MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor ~1.7 (Reduced HV settings) with respect to standard settings (Nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15-30% and between 60 and 80% for Nominal HV, Reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from ~18%, under dark nights only, to up to ~40% in total with only moderate performance degradation.

preprint2017arXiv

Prospects for CTA observations of the young SNR RX J1713.7-3946

We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.

preprint2017arXiv

Search for very-high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $Γ=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

preprint2017arXiv

Very-high-energy gamma-ray observations of the Type Ia Supernova SN 2014J with the MAGIC telescopes

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-Γ}$, with a spectral index of $Γ=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

preprint2016arXiv

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

preprint2016arXiv

AGILE Observations of the Gravitational Wave Event GW150914

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.

preprint2016arXiv

Deep observation of the NGC 1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse gamma-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E>100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6$^\circ$ from the centre, have been detected as point-like VHE gamma-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is well described by a power law with a photon index of $3.6\pm0.2_{stat}\pm0.2_{syst}$ between 90 GeV and 1.2 TeV. We do not detect any diffuse gamma-ray emission from the cluster and set stringent constraints on its CR population. In order to bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes $α$. For $α=2.2$, the CR-to-thermal pressure within the cluster virial radius is constrained to be below 1-2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to <20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, $B_0$, that depend on the CR distribution. For $α=2.2$, $B_0\gtrsim5-8 μ$G, which is below the 25 $μ$G inferred from Faraday rotation measurements, whereas, for $α\lesssim2.1$, the hadronic interpretation of the diffuse radio emission is in contrast with our gamma-ray flux upper limits independently of the magnetic field strength.

preprint2016arXiv

Detection of very high energy gamma-ray emission from the gravitationally-lensed blazar QSO B0218+357 with the MAGIC telescopes

Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC, are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray sources detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broad band emission can be modeled in the framework of a two zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.

preprint2016arXiv

Insights into the emission of the blazar 1ES 1011+496 through unprecedented broadband observations during 2011 and 2012

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $γ$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Metsähovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.

preprint2016arXiv

Investigating the peculiar emission from the new VHE gamma-ray source H1722+119

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $γ$-ray band with a statistical significance of 5.9 $σ$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $γ$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

preprint2016arXiv

Limits to dark matter annihilation cross-section from a combined analysis of MAGIC and Fermi-LAT observations of dwarf satellite galaxies

We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV - the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

MAGIC observations of the February 2014 flare of 1ES 1011+496 and ensuing constraint of the EBL density

In February-March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z=0.212) in flaring state at very-high energy (VHE, E>100GeV). The flux reached a level more than 10 times higher than any previously recorded flaring state of the source. We present the description of the characteristics of the flare presenting the light curve and the spectral parameters of the night-wise spectra and the average spectrum of the whole period. From these data we aim at detecting the imprint of the Extragalactic Background Light (EBL) in the VHE spectrum of the source, in order to constrain its intensity in the optical band. For this we implement the method developed by the H.E.S.S. collaboration in which the intrinsic energy spectrum of the source is modeled with a simple function, and the EBL-induced optical depth is calculated using a template EBL model. The likelihood of the observed spectrum is then maximized, including a normalization factor for the EBL opacity among the free parameters. From the data collected differential energy spectra was produced for all nights of the observed period. Evaluating the changes in the fit parameters we conclude that the spectral shape for most of the nights were compatible, regardless of the flux level, which enabled us to produce an average spectrum from which the EBL imprint could be constrained. The likelihood ratio test shows that the model with an EBL density 1.07(-0.20,+0.24)_{stat+sys}, relative to the one in the tested EBL template (Dominguez et al.2011), is preferred at the 4.6 sigma level to the no-EBL hypothesis, with the assumption that the intrinsic source spectrum can be modeled as a log-parabola. This would translate into a constraint of the EBL density in the wavelength range [0.24 um,4.25 um], with a peak value at 1.4 um of F=12.27_{-2.29}^{+2.75} nW m^{-2} sr^{-1}, including systematics.

preprint2016arXiv

Observations of Sagittarius A* during the pericenter passage of the G2 object with MAGIC

Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; $\geq$100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 hours of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 year observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio Arc.

preprint2016arXiv

Search for VHE gamma-ray emission from Geminga pulsar and nebula with the MAGIC telescopes

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2015arXiv

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsähovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

preprint2015arXiv

MAGIC observations of MWC 656, the only known Be/BH system

Context: MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims: Our aim is to evaluate if the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods. We have observed MWC 656 with the MAGIC telescopes for $\sim$23 hours during two observation periods: between May and June 2012 and June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results: We have not detected the MWC 656 binary system at TeV energies with the MAGIC Telescopes in either of the two campaigns carried out. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of $\sim$5% of the Crab Nebula flux. The results obtained from the MAGIC observations do not support persistent emission of very high energy gamma rays from this system at a level of 2.4% the Crab flux.

preprint2015arXiv

SN 2013dx associated with GRB 130702A: a detailed photometric and spectroscopic monitoring and a study of the environment

Long duration gamma-ray bursts (GRBs) and broad-line, type Ic supernovae (SNe) are strongly connected. We aim at characterizing SN 2013dx, associated with GRB\,130702A, through sensitive and extensive ground-based observational campaigns in the optical-IR band. We monitored the field of the Swift GRB 130702A (redshift z = 0.145) using the 8.2-m VLT, the 3.6-m TNG and the 0.6-m REM telescopes during the time interval between 4 and 40 days after the burst. Photometric and spectroscopic observations revealed the presence of the associated Type Ic SN 2013dx. Our multi-band photometry allowed the construction of a bolometric light curve.} The bolometric light curve of SN 2013dx resembles that of 2003dh (associated with GRB 030329), but is ~10% faster and ~25% dimmer. From this we infer a synthesized 56Ni mass of ~0.2 solar masses. The multi-epoch optical spectroscopy shows that the SN 2013dx behavior is best matched by SN 1998bw, among the other well-known low-redshift SNe associated with GRBs and XRFs, and by SN 2010ah, an energetic Type Ic SN not associated with any GRB. The photospheric velocity of the ejected material declines from ~2.7X10^4 km/s at 8 rest frame days from the explosion, to ~3.5X10^3 km/s at 40 days. These values are extremely close to those of SN1998bw and 2010ah. We deduce for SN 2013dx a kinetic energy of ~35X10^51 erg, and an ejected mass of ~7 solar masses. This suggests that the progenitor of SN2013dx had a mass of ~25 solar masses, i.e., 15-20% less massive than that of SN 1998bw. Finally, we performed a study of the SN 2013dx environment, through spectroscopy of the closeby galaxies. 9 out of the 14 inspected galaxies lie within 0.03 in redshift from z=0.145, indicating that the host of GRB 130702A/SN 2013dx belongs to a group of galaxies, an unprecedented finding for a GRB-associated SN and, to our knowledge, for long GRBs in general.

preprint2015arXiv

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

preprint2015arXiv

The major upgrade of the MAGIC telescopes, Part I: The hardware improvements and the commissioning of the system

The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers initiated by gamma rays in the energy regime from around 50 GeV to more than 50 TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011-2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, dead time of the system, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system.

preprint2015arXiv

The major upgrade of the MAGIC telescopes, Part II: A performance study using observations of the Crab Nebula

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ~50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66+/-0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the sigma of a 2-dimensional Gaussian distribution, at those energies is < 0.07 degree, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11-18% in flux normalization and +/-0.15 for the energy spectrum power-law slope.

preprint2015arXiv

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($δ$, $B$ and $R$), in addition to the parameters related to the particle population.

preprint2015arXiv

Very-high-energy γ-ray observations of novae and dwarf novae with the MAGIC telescopes

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV γ-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the γ-ray spectrum at TeV energies. Aims. We aim to explore the very-high-energy domain to search for γ-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, triggered by optical and subsequent GeV γ-ray detec- tions. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp<~0.15 Le.

preprint2014arXiv

Black hole lightning due to particle acceleration at subhorizon scales

Supermassive black holes with masses of millions to billions of solar masses are commonly found in the centers of galaxies. Astronomers seek to image jet formation using radio interferometry, but still suffer from insufficient angular resolution. An alternative method to resolve small structures is to measure the time variability of their emission. Here, we report on gamma-ray observations of the radio galaxy IC 310 obtained with the MAGIC telescopes revealing variability with doubling time scales faster than 4.8 min. Causality constrains the size of the emission region to be smaller than 20\% of the gravitational radius of its central black hole. We suggest that the emission is associated with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet.

preprint2014arXiv

Contemporaneous observations of the radio galaxy NGC 1275 from radio to very high energy gamma-rays

The radio galaxy NGC 1275, recently identified as a very high energy (VHE, >100 GeV) gamma-ray emitter by MAGIC, is one of the few non-blazar AGN detected in the VHE regime. In order to better understand the origin of the gamma-ray emission and locate it within the galaxy, we studied contemporaneous multi-frequency observations of NGC 1275 and modeled the overall spectral energy distribution (SED). We analyzed unpublished MAGIC observations carried out between Oct. 2009 and Feb. 2010, and the previously published ones taken between Aug. 2010 and Feb. 2011. We studied the multi-band variability and correlations by analyzing data of Fermi-LAT (0.1-100 GeV), as well as Chandra (X-ray), KVA (optical) and MOJAVE (radio) data taken during the same period. Using customized Monte Carlo simulations corresponding to early MAGIC stereo data, we detect NGC 1275 also in the earlier campaign. The flux level and energy spectra are similar to the results of the second campaign. The monthly light curve >100 GeV shows a hint of variability at the 3.6 sigma level. In the Fermi-LAT band, both flux and spectral shape variabilities are reported. The optical light curve is variable and shows a clear correlation with the gamma-ray flux >100 MeV. In radio, 3 compact components are resolved in the innermost part of the jet. One of them shows a similar trend as the LAT and KVA light curves. The 0.1-650 GeV spectra measured simultaneously with MAGIC and Fermi-LAT can be well fitted either by a log-parabola or by a power-law with a sub-exponential cutoff for both campaigns. A single-zone synchrotron-self-Compton model, with an electron spectrum following a power-law with an exponential cutoff, can explain the broadband SED and the multi-band behavior of the source. However, this model suggests an untypical low bulk-Lorentz factor or a velocity alignment closer to the line of sight than the pc-scale radio jet.

preprint2014arXiv

Detection of bridge emission above 50 GeV from the Crab pulsar with the MAGIC telescopes

The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma-rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of {\it Fermi}-LAT data were also analyzed. The known two pulses per period were detected with a significance of $8.0 σ$ and $12.6 σ$. In addition, significant emission was found between the two pulses with $6.2 σ$. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.

preprint2014arXiv

Discovery of TeV gamma-ray emission from the pulsar wind nebula 3C 58 by MAGIC

The pulsar wind nebula (PWN) 3C 58 is one of the historical very-high-energy (VHE; E>100 GeV) gamma-ray source candidates. It is energized by one of the highest spin-down power pulsars known (5% of Crab pulsar) and it has been compared to the Crab Nebula due to their morphological similarities. This object was previously observed by imaging atmospheric Cherenkov telescopes (Whipple, VERITAS and MAGIC), although not detected, with an upper limit of 2.4% Crab Unit (C.U.) at VHE. It was detected by Fermi-LAT with a spectrum extending beyond 100 GeV. We analyzed 81 hours of 3C 58 data taken with the MAGIC telescopes and we detected VHE gamma-ray emission with a significance of 5.7 sigma and an integral flux of 0.65% C.U. above 1 TeV. The differential energy spectrum between 400 GeV and 10 TeV is well described by a power-law function dϕ/dE=f_0(E/1TeV)^{-Gamma} with f_0=(2.0\pm0.4_{stat}\pm0.6_{sys})\times10^{-13}cm^{-2}s^{-1}TeV^{-1} and Gamma=2.4\pm0.2_{stat}\pm0.2_{sys}. The skymap is compatible with an unresolved source. We report the first significant detection of PWN 3C 58 at TeV energies. According to our results 3C 58 is the least luminous VHE gamma-ray PWN ever detected at VHE and the one with the lowest flux at VHE to date. We compare our results with the expectations of time-dependent models in which electrons up-scatter photon fields. The best representation favors a distance to the PWN of 2 kpc and Far Infrared (FIR) comparable to CMB photon fields. If we consider an unexpectedly high FIR density, the data can also be reproduced by models assuming a 3.2 kpc distance. A low magnetic field, far from equipartition, is required to explain the VHE data. Hadronic contribution from the hosting supernova remnant (SNR) requires unrealistic energy budget given the density of the medium, disfavoring cosmic ray acceleration in the SNR as origin of the VHE gamma-ray emission.

preprint2014arXiv

Discovery of very high energy gamma-ray emission from the blazar 1ES 0033+595 by the MAGIC telescopes

The number of known very high energy (VHE) blazars is $\sim\,50$, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the $γ$-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy $γ$-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES\,0033+595 with a statistical significance of 5.5\,$σ$. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parameterized with a power law with an integral flux above 150 GeV of $(7.1\pm1.3)\times 10^{-12} {\mathrm{ph\,cm^{-2}\,s^{-1}}}$ and a photon index of ($3.8\pm0.7$). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED we used simultaneous optical R-band data from the KVA telescope, archival X-ray data by \textit{Swift} as well as \textit{INTEGRAL}, and simultaneous high energy (HE, $300$\,MeV~--~$10$\,GeV) $γ$-ray data from the \textit{Fermi} LAT observatory. Using the empirical approach of Prandini et al. (2010) and the \textit{Fermi}-LAT and MAGIC spectra for this object, we estimate the redshift of this source to be $0.34\pm0.08\pm0.05$. This is a relevant result because this source is possibly one of the ten most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.

preprint2014arXiv

First broadband characterization and redshift determination of the VHE blazar MAGIC J2001+439

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 σ (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

preprint2014arXiv

MAGIC gamma-ray and multifrequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

preprint2014arXiv

MAGIC long-term study of the distant TeV blazar PKS 1424+240 in a multiwavelength context

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

preprint2014arXiv

MAGIC reveals a complex morphology within the unidentified gamma-ray source HESS J1857+026

HESS J1857+026 is an extended TeV gamma-ray source that was discovered by H.E.S.S. as part of its Galactic plane survey. Given its broadband spectral energy distribution and its spatial coincidence with the young energetic pulsar PSR J1856+0245, the source has been put forward as a pulsar wind nebula (PWN) candidate. MAGIC has performed follow-up observations aimed at mapping the source down to energies approaching 100 GeV in order to better understand its complex morphology. HESS J1857+026 was observed by MAGIC in 2010, yielding 29 hours of good quality stereoscopic data that allowed us to map the source region in two separate ranges of energy. We detected very-high-energy gamma-ray emission from HESS J1857+026 with a significance of $12 σ$ above $150$ GeV. The differential energy spectrum between $100$ GeV and $13$ TeV is well described by a power law function $dN/dE = N_0(E/1\textrm{TeV})^{-Γ}$ with $N_0 = (5.37 \pm0.44_{stat} \pm1.5_{sys}) \times 10^{-12} (\textrm{TeV}^{-1} \textrm{cm}^{-2}$ $\textrm{ s}^{-1})$ and $Γ= 2.16\pm0.07_{stat} \pm0.15_{sys}$, which bridges the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S.. In addition, we present a detailed analysis of the energy-dependent morphology of this region. We couple these results with archival multi-wavelength data and outline evidence in favor of a two-source scenario, whereby one source is associated with a PWN, while the other could be linked with a molecular cloud complex containing an HII region and a possible gas cavity.

preprint2014arXiv

MAGIC search for VHE $γ$-ray emission from AE Aquarii in a multiwavelength context

It has been claimed that the nova-like cataclysmic variable (CV) AE Aquarii (AE Aqr) is a very-high-energy (VHE, $E>$100 GeV) source both on observational and theoretical grounds. We aim to search for VHE gamma-ray emission from AE Aqr during different states of the source at several wavelengths to confirm or rule out previous claims of detection of gamma-ray emission from this object. We report on observations of AE Aqr performed by MAGIC. The source was observed during 12 hours as part of a multiwavelength campaign carried out between May and June 2012 covering the optical, X-ray, and gamma-ray ranges. Besides MAGIC, the other facilities involved were the KVA, Skinakas, and Vidojevica telescopes in the optical and Swift in X-rays. We calculated integral upper limits coincident with different states of the source in the optical. We computed upper limits to the pulsed emission limiting the signal region to 30% of the phaseogram and we also searched for pulsed emission at different frequencies applying the Rayleigh test. AE Aqr was not detected at VHE energies during the multiwavelength campaign. We establish integral upper limits at the 95\% confidence level for the steady emission assuming the differential flux proportional to a power-law function dϕ/dE \propto E^{-Gamma}, with a Crab-like photon spectral index of Gamma=2.6. The upper limit above 200 GeV is 6.4\times10^{-12} cm^{-2}s^{-1} and above 1 TeV is 7.4\times10^{-13} cm^{-2}s^{-1}. We obtained an upper limit for the pulsed emission of 2.6\times10^{-12} cm^{-2}s^{-1} for energies above 200 GeV. Applying the Rayleigh test for pulsed emission at different frequencies we did not find any significant signal. Our results indicate that AE Aqr is not a VHE gamma-ray emitter at the level of emission previously claimed. We have established the most constraining upper limits for the VHE gamma-ray emission of AE Aqr.

preprint2014arXiv

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

preprint2014arXiv

Multiwavelength observations of Mrk 501 in 2008

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

preprint2014arXiv

Optimized dark matter searches in deep observations of Segue 1 with MAGIC

We present the results of stereoscopic observations of the satellite galaxy Segue 1 with the MAGIC Telescopes, carried out between 2011 and 2013. With almost 160 hours of good-quality data, this is the deepest observational campaign on any dwarf galaxy performed so far in the very high energy range of the electromagnetic spectrum. We search this large data sample for signals of dark matter particles in the mass range between 100 GeV and 20 TeV. For this we use the full likelihood analysis method, which provides optimal sensitivity to characteristic gamma-ray spectral features, like those expected from dark matter annihilation or decay. In particular, we focus our search on gamma-rays produced from different final state Standard Model particles, annihilation with internal bremsstrahlung, monochromatic lines and box-shaped signals. Our results represent the most stringent constraints to the annihilation cross-section or decay lifetime obtained from observations of satellite galaxies, for masses above few hundred GeV. In particular, our strongest limit (95% confidence level) corresponds to a ~500 GeV dark matter particle annihilating into tau+tau-, and is of order <sigma v> ~ 1.2x10^{-24} cm^3 s^{-1} - a factor ~40 above the <sigma v> thermal value.

preprint2014arXiv

Rapid and multi-band variability of the TeV-bright active nucleus of the galaxy IC 310

Context. The radio galaxy IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E>100GeV) with the MAGIC telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. Aims. In order to understand the nature of IC 310 and the origin of the VHE emission we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE gamma-ray regime. Methods. The light curve of IC 310 above 300GeV has been measured with the MAGIC telescopes from Oct. 2009 to Feb. 2010. Fermi-LAT data (2008-2011) in the 10-500GeV energy range were also analyzed. In X-ray, archival observations from 2003 to 2007 with XMM, Chandra, and Swift-XRT in the 0.5-10keV band were studied. Results. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present. The photon index between 120GeV and 8TeV remains at the value $Γ\sim2.0$ during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1TeV changes by a factor of $\sim7$. Fermi-LAT detected only eight gamma-ray events in the energy range 10GeV-500GeV in three years of observation. The measured photon index of $Γ=1.3\pm0.5$ in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in flux and photon index. The latter varied from $1.76\pm0.07$ to $2.55\pm0.07$. Conclusion. The rapid variability measured confirms the blazar-like behavior of IC 310. The TeV emission seems to originate from scales of less than 80 Schwarzschild radii within the compact core of its FRI radio jet with orientation angle 10deg-38deg. The SED resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.

preprint2014arXiv

Search for Very-High-Energy Gamma Rays from the z = 0.896 Quasar 4C +55.17 with the MAGIC telescopes

The bright gamma-ray quasar 4C +55.17 is a distant source ($z = 0.896$) with a hard spectrum at GeV energies as observed by the Large Area Telescope (LAT) on board the {\it Fermi} satellite. This source is identified as a good source candidate for very-high-energy (VHE; $> 30$ GeV) gamma rays. In general VHE gamma rays from distant sources provide an unique opportunity to study the extragalactic background light (EBL) and underlying astrophysics. The flux intensity of this source in the VHE range is investigated. Then, constraints on the EBL are derived from the attenuation of gamma-ray photons coming from the distant blazar. We searched for a gamma-ray signal from this object using the 35-hour observations taken by the MAGIC telescopes between November 2010 and January 2011. No significant VHE gamma-ray signal was detected. We computed the upper limits of the integrated gamma-ray flux at $95\%$ confidence level of $9.4 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ and $2.5 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ above $100$ GeV and $200$ GeV, respectively. The differential upper limits in four energy bins in the range from $80$ GeV to $500$ GeV are also derived. The upper limits are consistent with the attenuation predicted by low-flux EBL models on the assumption of a simple power-law spectrum extrapolated from LAT data.

preprint2013arXiv

An updated list of AGILE bright gamma-ray sources and their variability in pointing mode

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°303, Cyg X-3, and 1AGLR J2021+4030.

preprint2013arXiv

Discovery of very high energy gamma-ray emission from the blazar 1ES 1727+502 with the MAGIC Telescopes

Motivated by the Costamante & Ghisellini (2002) predictions we investigated if the blazar 1ES 1727+502 (z=0.055) is emitting very high energy (VHE, E>100 GeV) gamma rays. We observed the BL Lac object 1ES 1727+502 in stereoscopic mode with the two MAGIC telescopes during 14 nights between May 6th and June 10th 2011, for a total effective observing time of 12.6 hours. For the study of the multiwavelength spectral energy distribution (SED) we use simultaneous optical R-band data from the KVA telescope, archival UV/optical and X-ray observations by instruments UVOT and XRT on board of the Swift satellite and high energy (HE, 0.1 GeV - 100 GeV) gamma-ray data from the Fermi-LAT instrument. We detect, for the first time, VHE gamma-ray emission from 1ES 1727+502 at a statistical significance of 5.5 sigma. The integral flux above 150 GeV is estimated to be (2.1\pm0.4)% of the Crab Nebula flux and the de-absorbed VHE spectrum has a photon index of (2.7\pm0.5). No significant short-term variability was found in any of the wavebands presented here. We model the SED using a one-zone synchrotron self-Compton model obtaining parameters typical for this class of sources.

preprint2013arXiv

Expected performance of the ASTRI-SST-2M telescope prototype

ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) is an Italian flagship project pursued by INAF (Istituto Nazionale di Astrofisica) strictly linked to the development of the Cherenkov Telescope Array, CTA. Primary goal of the ASTRI program is the design and production of an end-to-end prototype of a Small Size Telescope for the CTA sub-array devoted to the highest gamma-ray energy region. The prototype, named ASTRI SST-2M, will be tested on field in Italy during 2014. This telescope will be the first Cherenkov telescope adopting the double reflection layout in a Schwarzschild-Couder configuration with a tessellated primary mirror and a monolithic secondary mirror. The collected light will be focused on a compact and light-weight camera based on silicon photo-multipliers covering a 9.6 deg full field of view. Detailed Monte Carlo simulations have been performed to estimate the performance of the planned telescope. The results regarding its energy threshold, sensitivity and angular resolution are shown and discussed.

preprint2013arXiv

The ASTRI Mini-Array Science Case

ASTRI is a Flagship Project financed by the Italian Ministry of Education, University and Research, and led by INAF, the Italian National Institute of Astrophysics. Within this framework, INAF is currently developing an end-to-end prototype of a Small Size Telescope in a dual-mirror configuration (SST-2M) for the Cherenkov Telescope Array (CTA), scheduled to start data acquisition in 2014. Although the ASTRI SST-2M prototype is mainly a technological demonstrator, it will perform scientific observations of the Crab Nebula, Mrk 421 and Mrk 501 at E>1 TeV. A remarkable improvement in terms of performance could come from the operation, in 2016, of a SST-2M mini-array, composed of a few SST-2M telescopes to be placed at final CTA Southern Site. The SST mini-array will be able to study in great detail relatively bright sources (a few x 10E-12 erg/cm2/s at 10 TeV) with angular resolution of a few arcmin and energy resolution of about 10-15%. Thanks to the stereo approach, it will be possible to verify the wide field of view (FoV) performance through the detections of very high-energy showers with core located at a distance up to 500 m, to compare the mini-array performance with the Monte Carlo expectations by means of deep observations of selected targets, and to perform the first CTA science at the beginning of the mini-array operations. Prominent sources such as extreme blazars, nearby well-known BL Lac objects and radio-galaxies, galactic pulsar wind nebulae, supernovae remnants, micro-quasars, and the Galactic Center can be observed in a previously unexplored energy range, in order to investigate the electron acceleration and cooling, relativistic and non relativistic shocks, the search for cosmic-ray (CR) Pevatrons, the study of the CR propagation, and the impact of the extragalactic background light on the spectra of the sources.

preprint2013arXiv

The ASTRI Mini-Array Software System

ASTRI ("Astrofisica con Specchi a Tecnologia Replicante Italiana") is a Flagship Project financed by the Italian Ministry of Education, University and Research, and led by INAF, the Italian National Institute of Astrophysics. Main goals of the ASTRI project are the realization of an end-to-end prototype of a Small Size Telescope (SST) for the Cherenkov Telescope Array (CTA) in a dual-mirror configuration (SST-2M) and, subsequently, of a mini-array composed of a few SST-2M telescopes to be placed at the final CTA Southern Site. Here we present the main features of the Mini-Array Software System (MASS) that has a central role in the success of the ASTRI Project and will also serve as a prototype for the CTA software system. The MASS will provide a set of tools to prepare an observing proposal, to perform the observations specified therein (monitoring and controlling all the hardware components of each telescope), to analyze the acquired data online and to store/retrieve all the data products to/from the archive.}

preprint2013arXiv

The ASTRI Project: a mini-array of dual-mirror small Cherenkov telescopes for CTA

ASTRI is a flagship project of the Italian Ministry of Education, University and Research, which aims to develop an end-to-end prototype of the CTA small-size telescope. The proposed design is characterized by a dual-mirror Schwarzschild-Couder configuration and a camera based on Silicon photo-multipliers, two challenging but innovative technological solutions which will be adopted for the first time on a Cherenkov telescope. Here we describe the current status of the project, the expected performance and the possibility to realize a mini-array composed by a few small-size telescopes, which shall be placed at the final CTA Southern Site.

preprint2013arXiv

Towards the ASTRI mini-array

The Cherenkov Telescope Array (CTA) will consist of an array of three types of telescopes covering a wide energy range, from tens of GeV up to more than 100 TeV. The high energy section (> 3 TeV) will be covered by the Small Size Telescopes (SST). ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) is a flagship project of the Italian Ministry of Research and Education led by INAF, aiming at the design and construction of a prototype of the Dual Mirror SST. In a second phase the ASTRI project foresees the installation of the first elements of the SST array at the CTA southern site, a mini-array of 5-7 telescopes. The optimization of the layout of this mini-array embedded in the SST array of the CTA Observatory has been the object of an intense simulation effort. In this work we present the expected mini-array performance in terms of energy threshold, angular and energy resolution and sensitivity.

preprint2012arXiv

Constraining Cosmic Rays and Magnetic Fields in the Perseus Galaxy Cluster with TeV observations by the MAGIC telescopes

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.

preprint2012arXiv

Detection of very-high energy γ-ray emission from NGC 1275 by the MAGIC telescopes

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 Γ=-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.

preprint2012arXiv

Discovery of VHE gamma-ray emission from the BL Lac object B3 2247+381 with the MAGIC telescopes

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) $γ$-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 $γ$-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 $γ$-rays at a statistical significance of 5.6 $σ$. 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.

preprint2012arXiv

Discovery of VHE gamma-rays from the blazar 1ES 1215+303 with the MAGIC Telescopes and simultaneous multi-wavelength observations

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ähovi 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.

preprint2012arXiv

Galactic Sources Science With Agile: The Case Of The Carina Region

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.

preprint2012arXiv

High zenith angle observations of PKS 2155-304 with the MAGIC-I telescope

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.

preprint2012arXiv

MAGIC observations of the giant radio galaxy M87 in a low-emission state between 2005 and 2007

We present the results of a long M87 monitoring campaign in very high energy $γ$-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 $γ$-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σ$. We present the spectrum between 100\,GeV and 2\,TeV, which is consistent with a simple power law with a photon index $Γ=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.

preprint2012arXiv

Morphological and spectral properties of the W51 region measured with the MAGIC telescopes

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.

preprint2012arXiv

Mrk 421 active state in 2008: the MAGIC view, simultaneous multi-wavelength observations and SSC model constrained

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.

preprint2012arXiv

Observations of the magnetars 4U 0142+61 and 1E 2259+586 with the MAGIC telescopes

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.

preprint2012arXiv

PG 1553+113: five years of observations with MAGIC

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 Γ= 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.

preprint2012arXiv

Phase-resolved energy spectra of the Crab Pulsar in the range of 50-400 GeV measured with the MAGIC Telescopes

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.

preprint2012arXiv

The 2010 very high energy gamma-ray flare & 10 years of multi-wavelength observations of M 87

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.

preprint2012arXiv

The Optical SN 2012bz Associated with the Long GRB 120422A

The association of Type Ic SNe with long-duration GRBs is well established. We endeavor, through accurate ground-based observational campaigns, to characterize these SNe at increasingly high redshifts. We obtained a series of optical photometric and spectroscopic observations of the Type Ic SN2012bz associated with the Swift long-duration GRB120422A (z=0.283) using the 3.6-m TNG and the 8.2-m VLT telescopes. The peak times of the light curves of SN2012bz in various optical filters differ, with the B-band and i'-band light curves reaching maximum at ~9 and ~23 rest-frame days, respectively. The bolometric light curve has been derived from individual bands photometric measurements, but no correction for the unknown contribution in the near-infrared (probably around 10-15%) has been applied. Therefore, the present light curve should be considered as a lower limit to the actual UV-optical-IR bolometric light curve. This pseudo-bolometric curve reaches its maximum (Mbol = -18.56 +/- 0.06) at 13 +/- 1 rest-frame days; it is similar in shape and luminosity to the bolometric light curves of the SNe associated with z<0.2 GRBs and more luminous than those of SNe associated with XRFs. A comparison with the model generated for the bolometric light curve of SN2003dh suggests that SN2012bz produced only about 15% less 56Ni than SN2003dh, about 0.35 Msol. Similarly the VLT spectra of SN2012bz, after correction for Galactic extinction and for the contribution of the host galaxy, suggest comparable explosion parameters with those observed in SN2003dh (EK~3.5 x 10^52 erg, Mej~7 Msol) and a similar progenitor mass (~25-40 Msol). GRB120422A is consistent with the Epeak-Eiso and the EX,iso-Egamma,iso-E_peak relations. GRB120422A/SN2012bz shows the GRB-SN connection at the highest redshift so far accurately monitored both photometrically and spectroscopically.

preprint2011arXiv

A minor body falling onto a neutron star as an explanation for the unusual gamma-ray burst GRB 101225A

The tidal disruption of a solar mass star around a supermassive black hole has been extensively studied analytically and numerically. In these events the star develops into an elongated banana-shaped structure. After completing an eccentric orbit, the bound debris fall onto the black hole, forming an accretion disk and emitting radiation. The same process may occur on planetary scales, if a minor body orbits too close to its star. In the Solar System, comets fall directly onto our Sun or onto planets. If the star is a compact object, the minor body can become tidally disrupted. Indeed, one of the first mechanisms invoked to produce strong gamma-ray emission involved accretion of comets onto neutron stars in our Galaxy. Here we report that the peculiarities of the `Christmas' gamma-ray burst (GRB 101225A) can be explained by a tidal disruption event of a minor body around a Galactic isolated neutron star. This result would indicate either that minor bodies can be captured by compact stellar remnants more frequently than it occurs in our Solar System or that minor body formation is relatively easy around millisecond radio pulsars. A peculiar Supernova associated to a GRB may provide an alternative explanation.

preprint2011arXiv

A search for Very High Energy gamma-ray emission from Scorpius X-1 with the MAGIC telescopes

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.

preprint2011arXiv

AGILE detection of extreme gamma-ray activity from the blazar PKS 1510-089 during March 2009. Multifrequency analysis

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.

preprint2011arXiv

Detection of the gamma-ray binary LS I +61 303 in a low flux state at Very High Energy gamma-rays with the MAGIC Telescopes in 2009

We present very high energy (VHE, E > 100 GeV) γ-ray observations of the γ-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3σ γ-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.

preprint2011arXiv

Discovery of the Onset of Rapid Accretion by a Dormant Massive Black Hole

Massive black holes are believed to reside at the centres of most galaxies. They can be- come detectable by accretion of matter, either continuously from a large gas reservoir or impulsively from the tidal disruption of a passing star, and conversion of the gravitational energy of the infalling matter to light. Continuous accretion drives Active Galactic Nuclei (AGN), which are known to be variable but have never been observed to turn on or off. Tidal disruption of stars by dormant massive black holes has been inferred indirectly but the on- set of a tidal disruption event has never been observed. Here we report the first discovery of the onset of a relativistic accretion-powered jet in the new extragalactic transient, Swift J164449.3+573451. The behaviour of this new source differs from both theoretical models of tidal disruption events and observations of the jet-dominated AGN known as blazars. These differences may stem from transient effects associated with the onset of a powerful jet. Such an event in the massive black hole at the centre of our Milky Way galaxy could strongly ionize the upper atmosphere of the Earth, if beamed towards us.

preprint2011arXiv

Gamma-ray excess from a stacked sample of high- and intermediate-frequency peaked blazars observed with the MAGIC telescope

Between 2004 and 2009 a sample of 28 X-ray selected high- and intermediate-frequency peaked blazars with a X-ray flux larger than 2 uJy at 1 keV in the redshift range from 0.018 to 0.361 was observed with the MAGIC telescope at energies above 100 GeV. Seven among them were detected and the results of these observations are discussed elsewhere. Here we concentrate on the remaining 21 blazars which were not detected during this observation campaign and present the 3 sigma (99.7 %) confidence upper limits on their flux. The individual flux upper limits lie between 1.6 % and 13.6 % of the integral flux from the Crab Nebula. Applying a stacking method to the sample of non-detections with a total of 394.1 hours exposure time, we find evidence for an excess with a cumulative significance of 4.9 standard deviations. It is not dominated by individual objects or flares, but increases linearly with the observation time as for a constant source with an integral flux level of ~1.5 % of that observed from the Crab Nebula above 150 GeV.

preprint2011arXiv

GRB091127/SN2009nz and the VLT/X-shooter spectroscopy of its host galaxy: probing the faint end of the mass-metallicity relation

We perform a detailed study of the gamma-ray burst GRB091127/SN2009nz host galaxy at z=0.490 using the VLT/X-shooter spectrograph in slit and integral-field unit (IFU). From the analysis of the optical and X-ray afterglow data obtained from ground-based telescopes and Swift-XRT we confirm the presence of a bump associated with SN2009nz and find evidence of a possible jet break in the afterglow lightcurve. The X-shooter afterglow spectra reveal several emission lines from the underlying host, from which we derive its integrated properties. These are in agreement with those of previously studied GRB-SN hosts and, more generally, with those of the long GRB host population. We use the Hubble Space Telescope and ground based images of the host to determine its stellar mass (M_star). Our results extend to lower M_star values the M-Z plot derived for the sample of long GRB hosts at 0.3<z<1.0 adding new information to probe the faint end of the M-Z relation and the shift of the LGRB host M-Z relation from that found from emission line galaxy surveys. Thanks to the IFU spectroscopy we can build the 2D velocity, velocity dispersion and star formation rate (SFR) maps. They show that the host galaxy has a perturbed rotation kinematics with evidence of a SFR enhancement consistent with the afterglow position.

preprint2011arXiv

MAGIC contributions to the 32nd International Cosmic Ray Conference

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.

preprint2011arXiv

MAGIC observations and multiwavelength properties of the quasar 3C279 in 2007 and 2009

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.

preprint2011arXiv

MAGIC Upper Limits for two Milagro-detected, Bright Fermi Sources in the Region of SNR G65.1+0.6

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

preprint2011arXiv

Observations of the Crab pulsar between 25 GeV and 100 GeV with the MAGIC I telescope

We report on the observation of $γ$-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.

preprint2011arXiv

Performance of the MAGIC stereo system obtained with Crab Nebula data

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.

preprint2011arXiv

Searches for Dark Matter annihilation signatures in the Segue 1 satellite galaxy with the MAGIC-I telescope

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.

preprint2011arXiv

The Afterglows of Swift-era Gamma-Ray Bursts II.: Type I GRB versus Type II GRB Optical Afterglows

We use a large sample of GRB afterglow and prompt-emission data (adding further GRB afterglow observations in this work) to compare the optical afterglows (or the lack thereof) of Type I GRBs with those of Type II GRBs. In comparison to the afterglows of Type II GRBs, we find that those of Type I GRBs have a lower average luminosity and show an intrinsic spread of luminosities at least as wide. From late and deep upper limits on the optical transients, we establish limits on the maximum optical luminosity of any associated supernova, confirming older works and adding new results. We use deep upper limits on Type I GRB optical afterglows to constrain the parameter space of possible mini-SN emission associated with a compact-object merger. Using the prompt emission data, we search for correlations between the parameters of the prompt emission and the late optical afterglow luminosities. We find tentative correlations between the bolometric isotropic energy release and the optical afterglow luminosity at a fixed time after trigger (positive), and between the host offset and the luminosity (negative), but no significant correlation between the isotropic energy release and the duration of the GRBs. We also discuss three anomalous GRBs, GRB 060505, GRB 060614, and GRB 060121, in the light of their optical afterglow luminosities. (Abridged)

preprint2011arXiv

The AGILE observations of the hard and bright GRB 100724B

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.

preprint2011arXiv

The flaring blazars of the first 1.5 years of the AGILE mission

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.

preprint2011arXiv

The observation of Gamma Ray Bursts and Terrestrial Gamma-ray Flashes with AGILE

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.

preprint2010arXiv

A year-long AGILE observation of Cygnus X-1 in hard spectral state

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.

preprint2010arXiv

Challenging GRB models through the broadband dataset of GRB060908

Context: Multiwavelength observations of gamma-ray burst prompt and afterglow emission are a key tool to disentangle the various possible emission processes and scenarios proposed to interpret the complex gamma-ray burst phenomenology. Aims: We collected a large dataset on GRB060908 in order to carry out a comprehensive analysis of the prompt emission as well as the early and late afterglow. Methods: Data from Swift-BAT, -XRT and -UVOT together with data from a number of different ground-based optical/NIR and millimeter telescopes allowed us to follow the afterglow evolution from about a minute from the high-energy event down to the host galaxy limit. We discuss the physical parameters required to model these emissions. Results: The prompt emission of GRB060908 was characterized by two main periods of activity, spaced by a few seconds of low intensity, with a tight correlation between activity and spectral hardness. Observations of the afterglow began less than one minute after the high-energy event, when it was already in a decaying phase, and it was characterized by a rather flat optical/NIR spectrum which can be interpreted as due to a hard energy-distribution of the emitting electrons. On the other hand, the X-ray spectrum of the afterglow could be fit by a rather soft electron distribution. Conclusions: GRB060908 is a good example of a gamma-ray burst with a rich multi-wavelength set of observations. The availability of this dataset, built thanks to the joint efforts of many different teams, allowed us to carry out stringent tests for various interpretative scenarios showing that a satisfactorily modeling of this event is challenging. In the future, similar efforts will enable us to obtain optical/NIR coverage comparable in quality and quantity to the X-ray data for more events, therefore opening new avenues to progress gamma-ray burst research.

preprint2010arXiv

Detection of very high energy gamma-ray emission from the Perseus cluster head-tail galaxy IC 310 by the MAGIC telescopes

We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy.

preprint2010arXiv

Direct Evidence for Hadronic Cosmic-Ray Acceleration in the Supernova Renmant IC 443

The Supernova Remnant (SNR) IC 443 is an intermediate-age remnant well known for its radio, optical, X-ray and gamma-ray energy emissions. In this Letter we study the gamma-ray emission above 100 MeV from IC 443 as obtained by the AGILE satellite. A distinct pattern of diffuse emission in the energy range 100 MeV-3 GeV is detected across the SNR with its prominent maximum (source "A") localized in the Northeastern shell with a flux F = (47 \pm 10) 10^{-8} photons cm^{-2} s^{-1} above 100 MeV. This location is the site of the strongest shock interaction between the SNR blast wave and the dense circumstellar medium. Source "A" is not coincident with the TeV source located 0.4 degree away and associated with a dense molecular cloud complex in the SNR central region. From our observations, and from the lack of detectable diffuse TeV emission from its Northeastern rim, we demonstrate that electrons cannot be the main emitters of gamma-rays in the range 0.1-10 GeV at the site of the strongest SNR shock. The intensity, spectral characteristics, and location of the most prominent gamma-ray emission together with the absence of co-spatial detectable TeV emission are consistent only with a hadronic model of cosmic-ray acceleration in the SNR. A high-density molecular cloud (cloud "E") provides a remarkable "target" for nucleonic interactions of accelerated hadrons: our results show enhanced gamma-ray production near the molecular cloud/shocked shell interaction site. IC 443 provides the first unambiguous evidence of cosmic-ray acceleration by SNRs.

preprint2010arXiv

Episodic Transient Gamma-Ray Emission from the Microquasar Cygnus X-1

Cygnus X-1 is the archetypal black hole (BH) binary system in our Galaxy. We report the main results of an extensive search for transient gamma-ray emission from Cygnus X-1 carried out in the energy range 100 MeV - 3 GeV by the AGILE satellite, during the period 2007 July - 2009 October. The total exposure time is about 300 days, during which the source was in the "hard" X-ray spectral state. We divided the observing intervals in 2 or 4 week periods, and searched for transient and persistent emission. We report an episode of significant transient gamma-ray emission detected on 2009, October 16 in a position compatible with Cygnus X-1 optical position. This episode, occurred during a hard spectral state of Cygnus X-1, shows that a 1-2 day time variable emission above 100 MeV can be produced during hard spectral states, having important theoretical implications for current Comptonization models for Cygnus X-1 and other microquasars. Except for this one short timescale episode, no significant gamma-ray emission was detected by AGILE. By integrating all available data we obtain a 2$σ$ upper limit for the total integrated flux of $F_{γ,U.L.} = 3 \times 10^{-8} \rm ph cm^{-2} s^{-1}$ in the energy range 100 MeV - 3 GeV. We then clearly establish the existence of a spectral cutoff in the energy range 1-100 MeV that applies to the typical hard state outside the flaring period and that confirms the historically known spectral cutoff above 1 MeV.

preprint2010arXiv

Gamma-Ray Localization of Terrestrial Gamma-Ray Flashes

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.

preprint2010arXiv

GRB 090313: X-shooter's first shot at a GRB

Context. X-shooter is the first second-generation instrument to become operative at the ESO Very Large Telescope (VLT). It is a broad-band medium-resolution spectrograph designed with gamma-ray burst (GRB) afterglow spectroscopy as one of its main science drivers. Aims. During the first commissioning night on sky with the instrument fully assembled, X-shooter observed the afterglow of GRB 090313 as a demonstration of the instrument's capabilities. Methods. GRB 090313 was observed almost two days after the burst onset, when the object had already faded to R~21.6. Furthermore, the 90% illuminated Moon was just 30 degrees away from the field. In spite of the adverse conditions, we obtained a spectrum that, for the first time in GRB research, covers simultaneously the range from 5700 to 23000 Angstroms. Results. The spectrum shows multiple absorption features at a redshift of 3.3736, the redshift of the GRB. These features are composed of 3 components with different ionisation levels and velocities. Some of the features have never been observed before in a GRB at such a high redshift. Furthermore, we detect two intervening systems at redshifts of 1.8005 and 1.9597. Conclusions. These results demonstrate the potential of X-shooter in the GRB field, as it was capable of observing a GRB down to a magnitude limit that would include 72% of long GRB afterglows 2 hours after the burst onset. Coupled with the rapid response mode available at VLT, allowing reaction times of just a few minutes, X-shooter constitutes an important leap forward on medium resolution spectroscopic studies of GRBs, their host galaxies and intervening systems, probing the early history of the Universe.

preprint2010arXiv

Magic constraints on Gamma-ray emission from Cygnus X-3

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

preprint2010arXiv

MAGIC Gamma-ray Telescope Observation of the Perseus Cluster of Galaxies: implications for cosmic rays, dark matter, and NGC1275

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

preprint2010arXiv

MAGIC TeV Gamma-Ray Observations of Markarian 421 during Multiwavelength Campaigns in 2006

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.

preprint2010arXiv

Observations of the Blazar 3C 66A with the MAGIC Telescopes in Stereoscopic Mode

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A with the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play an essential role for the separation between 3C 66A and the nearby radio galaxy 3C 66B, which is at a distance of only $6^\prime$. The derived integral flux above $100\eh{GeV}$ is 8.3% of Crab Nebula flux and the energy spectrum is reproduced by a power law of photon index $3.64 \pm 0.39_{\rm stat} \pm 0.25_{\rm sys}$. Within errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected for absorption by the extragalactic background light, we only find small differences between the four models that we applied, and constrain the redshift of the blazar to $z < 0.68$.

preprint2010arXiv

Search for an extended VHE gamma-ray emission from Mrk 421 and Mrk 501 with the MAGIC Telescope

Context: Part of the very high energy $γ$-ray radiation coming from extragalactic sources is absorbed through the pair production process on the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright VHE sources. Aims: We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. Methods: If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival directions for a point source. Results: We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We discuss possible constraints on the extragalactic magnetic fields strength. We obtain upper limits on the extended emission around the Mrk~421 (Mrk~501) on the level of < 5% (< 4%) of the Crab Nebula flux above the energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times $10^{-15}$~G.

preprint2010arXiv

Simultaneous multi-frequency observation of the unknown redshift blazar PG 1553+113 in March-April 2008

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

preprint2010arXiv

Spectral Energy Distribution of Markarian 501: Quiescent State vs. Extreme Outburst

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including γ-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

preprint2010arXiv

Study of the gamma-ray source 1AGL J2022+4032 in the Cygnus Region

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.

preprint2010arXiv

Testing GRB models with the afterglow of GRB 090102

We present the observations of the afterglow of gamma-ray burst GRB 090102. Optical data taken by the TAROT, REM, GROND, together with publicly available data from Palomar, IAC and NOT telescopes, and X-ray data taken by the XRT instrument on board the Swift spacecraft were used. This event features an unusual light curve. In X-rays, it presents a constant decrease with no hint of temporal break from 0.005 to 6 days after the burst. In the optical, the light curve presents a flattening after 1 ks. Before this break, the optical light curve is steeper than that of the X-ray. In the optical, no further break is observed up to 10 days after the burst. We failed to explain these observations in light of the standard fireball model. Several other models, including the cannonball model were investigated. The explanation of the broad band data by any model requires some fine tuning when taking into account both optical and X-ray bands.

preprint2010arXiv

The Afterglows of Swift-era Gamma-Ray Bursts. I. Comparing pre-Swift and Swift era Long/Soft (Type II) GRB Optical Afterglows

We have gathered optical photometry data from the literature on a large sample of Swift-era gamma-ray burst (GRB) afterglows including GRBs up to September 2009, for a total of 76 GRBs, and present an additional three pre-Swift GRBs not included in an earlier sample. Furthermore, we publish 840 additional new photometry data points on a total of 42 GRB afterglows, including large data sets for GRBs 050319, 050408, 050802, 050820A, 050922C, 060418, 080413A and 080810. We analyzed the light curves of all GRBs in the sample and derived spectral energy distributions for the sample with the best data quality, allowing us to estimate the host galaxy extinction. We transformed the afterglow light curves into an extinction-corrected z=1 system and compared their luminosities with a sample of pre-Swift afterglows. The results of a former study, which showed that GRB afterglows clustered and exhibited a bimodal distribution in luminosity space, is weakened by the larger sample. We found that the luminosity distribution of the two afterglow samples (Swift-era and pre-Swift) are very similar, and that a subsample for which we were not able to estimate the extinction, which is fainter than the main sample, can be explained by assuming a moderate amount of line-of-sight host extinction. We derived bolometric isotropic energies for all GRBs in our sample, and found only a tentative correlation between the prompt energy release and the optical afterglow luminosity at one day after the GRB in the z=1 system. A comparative study of the optical luminosities of GRB afterglows with echelle spectra (which show a high number of foreground absorbing systems) and those without reveals no indication that the former are statistically significantly more luminous. (abridged)

preprint2010arXiv

The extraordinary gamma-ray flare of the blazar 3C 454.3

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_{γ,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_{γ, source, peak} \simeq 3 x 10^46 erg s^-1, for a relativistic Doppler factor of δ \simeq 30). The total isotropic irradiated energy of the month-long episode in the range 100 MeV - 3 GeV is E_{γ,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.

preprint2010arXiv

The two-component jet of GRB 080413B

The quick and precise localization of GRBs by the Swift telescope allows the early evolution of the afterglow light curve to be captured by ground-based telescopes. With GROND measurements we can investigate the optical/near-infrared light curve of the afterglow of gamma-ray burst 080413B in the context of late rebrightening. Multi-wavelength follow-up observations were performed on the afterglow of GRB 080413B. X-ray emission was detected by the X-ray telescope onboard the Swift satellite and obtained from the public archive. Optical and near-infrared photometry was performed with the seven-channel imager GROND mounted at the MPG/ESO 2.2 m telescope and additionally with the REM telescope, both in La Silla, Chile. The light curve model was constructed using the obtained broad-band data. The broad-band light curve of the afterglow of GRB 080413B is well fitted with an on-axis two-component jet model. The narrow ultra-relativistic jet is responsible for the initial decay, while the rise of the moderately relativistic wider jet near its deceleration time is the cause of the rebrightening of the light curve. The later evolution of the optical/NIR light curve is then dominated by the wide component, the signature of which is almost negligible in the X-ray wavelengths. These components have opening angles of theta(narrow) ~1.7 degrees and theta(wide) ~9 degrees, and Lorentz factors of Gamma(narrow) >188 and Gamma(wide) ~18.5. We calculated the beaming-corrected energy release to be 7.9 x 10^48 erg.

preprint2009arXiv

AGILE detection of a rapid gamma-ray flare from the blazar PKS 1510-089 during the GASP-WEBT monitoring

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

preprint2009arXiv

Agile Detection of Delayed Gamma-Ray Emission from the Short Gamma-Ray Burst GRB 090510

Short gamma-ray bursts (GRBs), typically lasting less than 2 s, are a special class of GRBs of great interest. We report the detection by the AGILE satellite of the short GRB 090510 which shows two clearly distinct emission phases: a prompt phase lasting ~ 200 msec and a second phase lasting tens of seconds. The prompt phase is relatively intense in the 0.3-10 MeV range with a spectrum characterized by a large peak/cutoff energy near 3 MeV, in this phase, no significant high-energy gamma-ray emission is detected. At the end of the prompt phase, intense gamma-ray emission above 30 MeV is detected showing a power-law time decay of the flux of the type t^-1.3 and a broad-band spectrum remarkably different from that of the prompt phase. It extends from sub-MeV to hundreds of MeV energies with a photon index alpha ~ 1.5. GRB 090510 provides the first case of a short GRB with delayed gamma-ray emission. We present the timing and spectral data of GRB 090510 and briefly discuss its remarkable properties within the current models of gamma-ray emission of short GRBs.

preprint2009arXiv

GRB 090423 at a redshift of z~8.1

Gamma-ray bursts (GRBSs) are produced by rare types of massive stellar explosions. Their rapidly fading afterglows are often bright enough at optical wavelengths, that they are detectable up to cosmological distances. Hirtheto, the highest known redshift for a GRB was z=6.7, for GRB 080913, and for a galaxy was z=6.96. Here we report observations of GRB 090423 and the near-infrared spectroscopic measurement of its redshift z=8.1^{+0.1}_{-0.3}. This burst happened when the Universe was only ~4% of its current age. Its properties are similar to those of GRBs observed at low/intermediate redshifts, suggesting that the mechanisms and progenitors that gave rise to this burst about 600 million years after the Big Bang are not markedly different from those producing GRBs ~10 billion years later.

preprint2009arXiv

GRB090426: the farthest short gamma-ray burst?

Aims. With an observed and rest-frame duration of < 2s and < 0.5s, respectively, GRB090426 could be classified as a short GRB. The prompt detection, both from space and ground-based telescopes, of a bright optical counterpart to this GRB offered a unique opportunity to complete a detailed study. Methods. Based on an extensive ground-based observational campaign, we obtained the spectrum of the optical afterglow of GRB090426, measuring its redshift and obtaining information about the medium in which the event took place. We completed follow-up observation of the afterglow optical light curve down to the brightness level of the host galaxy that we firmly identified and studied. We also retrieved and analyzed all the available high-energy data of this event, and compared the results with our findings in the optical. This represents one of the most detailed studies of a short-duration event presented so far. Results. The time properties qualify GRB090426 as a short burst. In this case, its redshift of z = 2.61 would be the highest yet found for a GRB of this class. On the other hand, the spectral and energy properties are more similar to those of long bursts. LBT late-time deep imaging identifies a star-forming galaxy at a redshift consistent with that of the GRB. The afterglow lies within the light of its host and shows evidence of local absorption.

preprint2009arXiv

Multiwavelength observations of 3C 454.3 II. The AGILE 2007 December campaign

We report on the second AGILE multiwavelength campaign of the blazar 3C 454.3 during the first half of December 2007. This campaign involved AGILE, Spitzer, Swift,Suzaku,the WEBT consortium,the REM and MITSuME telescopes,offering a broad band coverage that allowed for a simultaneous sampling of the synchrotron and inverse Compton (IC) emissions.The 2-week AGILE monitoring was accompanied by radio to optical monitoring by WEBT and REM and by sparse observations in mid-Infrared and soft/hard X-ray energy bands performed by means of Target of Opportunity observations by Spitzer, Swift and Suzaku, respectively.The source was detected with an average flux of~250x10^{-8}ph cm^-2s^-1 above 100 MeV,typical of its flaring states.The simultaneous optical and gamma-ray monitoring allowed us to study the time-lag associated with the variability in the two energy bands, resulting in a possible ~1-day delay of the gamma-ray emission with respect to the optical one. From the simultaneous optical and gamma-ray fast flare detected on December 12, we can constrain the delay between the gamma-ray and optical emissions within 12 hours. Moreover, we obtain three Spectral Energy Distributions (SEDs) with simultaneous data for 2007 December 5, 13, 15, characterized by the widest multifrequency coverage. We found that a model with an external Compton on seed photons by a standard disk and reprocessed by the Broad Line Regions does not describe in a satisfactory way the SEDs of 2007 December 5, 13 and 15. An additional contribution, possibly from the hot corona with T=10^6 K surrounding the jet, is required to account simultaneously for the softness of the synchrotron and the hardness of the inverse Compton emissions during those epochs.

preprint2009arXiv

Temporal properties of GX 301-2 over a year-long observation with SuperAGILE

We present the long-term monitoring of the High Mass X-ray Binary GX 301-2 performed with the SuperAGILE instrument on-board the AGILE mission. The source was monitored in the 20-60 keV energy band during the first year of the mission from 2007 July 17 to 2008 August 31, covering about one whole orbital period and three more pre-periastron passages for a total net observation time of about 3.7 Ms. The SuperAGILE dataset represents one of the most continuous and complete monitoring at hard X-ray energies of the 41.5 day long binary period available to date. The source behavior was characterized at all orbital phases in terms of hard X-ray flux, spectral hardness, spin period history, pulsed fraction and pulse shape profile. We also complemented the SuperAGILE observations with the soft X-ray data of the RossiXTE/ASM. Our analysis shows a clear orbital modulation of the spectral hardness, with peaks in correspondence with the pre-periastron flare and near phase 0.25. The hardness peaks we found could be related with the wind-plus-stream accretion model proposed in order to explain the orbital light curve modulation of GX 301-2. Timing analysis of the pulsar spin period shows that the secular trend of the about 680 s pulse period is consistent with the previous observations, although there is evidence of a slight decrease in the spin-down rate. The analysis of the hard X-ray pulsed emission also showed a variable pulse shape profile as a function of the orbital phase, with substructures detected near the passage at the periastron, and a clear modulation of the pulsed fraction, which appears in turn strongly anti-correlated with the source intensity.

preprint2008arXiv

Gamma-ray burst detection with the AGILE mini-calorimeter

The Mini-Calorimeter (MCAL) instrument on-board the AGILE satellite is a non-imaging gamma-ray scintillation detector sensitive in the 300keV-100MeV energy range with a total on-axis geometrical area of 1400cm^2. Gamma-Ray Bursts (GRBs) are one of the main scientific targets of the AGILE mission and the MCAL design as an independent self-triggering detector makes it a valuable all-sky monitor for GRBs. Furthermore MCAL is one of the very few operative instruments with microsecond timing capabilities in the MeV range. In this paper the results of GRB detections with MCAL after one year of operation in space are presented and discussed. A flexible trigger logic implemented in the AGILE payload data-handling unit allows the on-board detection of GRBs. For triggered events, energy and timing information are sent to telemetry on a photon-by-photon basis, so that energy and time binning are limited by counting statistics only. When the trigger logic is not active, GRBs can be detected offline in ratemeter data, although with worse energy and time resolution. Between the end of June 2007 and June 2008 MCAL detected 51 GRBs, with a detection rate of about 1 GRB/week, plus several other events at a few milliseconds timescales. Since February 2008 the on-board trigger logic has been fully active. Comparison of MCAL detected events and data provided by other space instruments confirms the sensitivity and effective area estimations. MCAL also joined the 3rd Inter-Planetary Network, to contribute to GRB localization by means of triangulation.

preprint2008arXiv

High energy variability of 3C 273 during the AGILE multiwavelength campaign of December 2007 - January 2008

We report the results of a 3-weeks multi-wavelength campaign on the flat spectrum radio quasar 3C 273 carried out with the AGILE gamma-ray mission, covering the 30 MeV -50 GeV and 18-60 keV, the REM observatory (covering the near-IR and optical), Swift (near-UV/Optical, 0.2-10 keV and 15-50 keV), INTEGRAL (3 - 200 keV) and Rossi XTE (2-12 keV). This is the first observational campaign including gamma-ray data, after the last EGRET observations, more than 8 years ago. Our study was carried out using simultaneous light curves of the source flux from all the involved instruments, in the different energy ranges, in search for correlated variability. Then a time-resolved SED was used for a detailed physical modelling of the emission mechanisms. The source was detected in gamma-rays only in the second week of our campaign. We found indication of a possible anti-correlation between the emission at gamma-rays and at soft and hard X-rays, supported by the complete set of instruments. Instead, optical data do not show short term variability. Only in two EGRET observations (in 1993 and 1997) 3C 273 showed intra-observation variability in gamma-rays. The SED with almost simultaneous data, partially covers the regions of the synchrotron emission, the big blue bump, and the inverse-Compton. The soft X-ray emission is consistent with combined SSC and EC mechanisms, while hard X and gamma-ray emissions are compatible with EC from thermal photons of the disk. The time evolution of the SED is well interpreted and modelled in terms of an acceleration episode of the electrons population, leading to a shift in the IC peak towards higher energies.

preprint2007arXiv

A study of the prompt and afterglow emission of the Short GRB 061201

Our knowledge of the intrinsic properties of short duration Gamma-Ray Bursts has relied, so far, only upon a few cases for which the estimate of the distance and an extended, multiwavelength monitoring of the afterglow have been obtained. We carried out multiwavelength observations of the short GRB 061201 aimed at estimating its distance and studying its properties. We performed a spectral and timing analysis of the prompt and afterglow emission and discuss the results in the context of the standard fireball model. A clear temporal break was observed in the X-ray light curve about 40 minutes after the burst trigger. We find that the spectral and timing behaviour of the X-ray afterglow is consistent with a jet origin of the observed break, although the optical data can not definitively confirm this and other scenarios are possible. No underlying host galaxy down to R~26 mag was found after fading of the optical afterglow. Thus, no secure redshift could be measured for this burst. The nearest galaxy is at z=0.111 and shows evidence of star formation activity. We discuss the association of GRB 061201 with this galaxy and with the ACO S 995 galaxy cluster, from which the source is at an angular distance of 17'' and 8.5', respectively. We also test the association with a possible undetected, positionally consistent galaxy at z~1. In all these cases, in the jet interpretation, we find a jet opening angle of 1-2 degrees.

preprint2007arXiv

Swift GRBs: the early afterglow spectral energy distribution

We present the first results of a program to systematically study the optical-to-X-ray spectral energy distribution (SED) of Swift GRB afterglows with known redshift. The goal is to study the properties of the GRB explosion and of the intervening absorbing material. In this report we present the preliminary analysis on 23 afterglows. Thanks to Swift, we could build the SED at early times after the GRB (minutes to hours). We derived the Hydrogen column densities and the spectral slopes from the X-ray spectrum. We then constrained the visual extinction by requiring that the combined optical/X-ray SED is due to synchrotron, namely either a single power law or a broken power law with a slope change by 0.5. We confirm a low dust-to-metal ratio, smaller than in the SMC, even from the analysis of data taken significantly earlier than previously possible. Our analysis does not support the existence of ``grey'' dust. We also find that the synchrotron spectrum works remarkably well to explain afterglow SEDs. We clearly see, however, that during the X-ray steep decay phases and the flares, the X-ray radiation cannot be due only to afterglow emission.

preprint2007arXiv

The afterglow onset for GRB060418 and GRB060607A

Gamma-ray burst are thought to be produced by highly relativistic outflows. Although upper and lower limits for the outflow initial Lorentz factor $Γ_0$ are available, observational efforts to derive a direct determination of $Γ_0$ have so far failed or provided ambiguous results. As a matter of fact, the shape of the early-time afterglow light curve is strongly sensitive on $Γ_0$ which determines the time of the afterglow peak, i.e. when the outflow and the shocked circumburst material share a comparable amount of energy. We now comment early-time observations of the near-infrared afterglows of GRB 060418 and GRB 060607A performed by the REM robotic telescope. For both events, the afterglow peak was singled out and allowed us to determine the initial fireball Lorentz, $Γ_0\sim 400$.

preprint2007arXiv

UVES/VLT high resolution spectroscopy of GRB 050730 afterglow: probing the features of the GRB environment

We analyze high resolution spectroscopic observations of the optical afterglow of GRB050730, obtained with UVES@VLT about hours after the GRB trigger. The spectrum shows that the ISM of the GRB host galaxy at z = 3.967 is complex, with at least five components contributing to the main absorption system. We detect strong CII*, SiII*, OI* and FeII* fine structure absorption lines associated to the second and third component. For the first three components we derive information on the relative distance from the site of the GRB explosion. Component 1, which has the highest redshift, does not present any fine structure nor low ionization lines; it only shows very high ionization features, such as CIV and OVI, suggesting that this component is very close to the GRB site. From the analysis of low and high ionization lines and fine structure lines, we find evidences that the distance of component 2 from the site of the GRB explosion is 10-100 times smaller than that of component 3. We evaluated the mean metallicity of the z=3.967 system obtaining values about 0.01 of the solar metallicity or less. However, this should not be taken as representative of the circumburst medium, since the main contribution to the hydrogen column density comes from the outer regions of the galaxy while that of the other elements presumably comes from the ISM closer to the GRB site. Furthermore, difficulties in evaluating dust depletion correction can modify significantly these values. The mean [C/Fe] ratio agrees well with that expected by single star-formation event models. Interestingly the [C/Fe] of component 2 is smaller than that of component 3, in agreement with GRB dust destruction scenarios, if component 2 is closer than component 3 to the GRB site.

preprint2006arXiv

GRB 050223: A dark GRB in a dusty starburst galaxy

Aims: We aim at detecting and determining the properties of the host galaxy of the dark GRB 050223. Methods: We use VLT optical/NIR images coupled to Swift X-ray positioning, and optical spectra of the host galaxy to measure its properties. Results: We find a single galaxy within the Swift error box of GRB 050223. It is located at z = 0.584 and its luminosity is L ~ 0.4 L*. Emission lines in the galaxy spectrum imply an intrinsic SFR > 7 Msun/yr, and a large extinction A_V > 2 mag within it. We also detect absorption lines, which reveal an underlying stellar population with an age between 40 Myr and 1.5 Gyr. Conclusions: The identification of a host galaxy with atypical properties using only the X-ray transient suggests that a bias may be present in the former sample of host galaxies. Dust obscuration together with intrinsic faintness are the most probable causes for the darkness of this burst.

preprint2006arXiv

Optical emission from GRB 050709: a short/hard GRB in a star forming galaxy

We present optical observations of the short/hard gamma-ray burst GRB 050709, the first such event with an identified optical counterpart. The object is coincident with a weak X-ray source and is located inside a galaxy at redshift z = 0.1606 +- 0.0002. Multiband photometry allowed us to study the broad-band spectral energy distribution. Late-time monitoring places strong limits on any supernova simultaneous with the GRB. The host galaxy is not of early type. Spectra show that the dominant stellar population is relatively young (~ 1 Gyr), and that ongoing star formation is present at a level of 2-3 (L/L*) Msun/yr. This is at least 2 orders of magnitude larger than that observed in the elliptical hosts of the short GRB 050509B and GRB 050724. This shows that at least some short GRBs originate in a young population. Short/hard GRB models based on the merger of a binary degenerate system are compatible with the host galaxy characteristics, although there is still the possibility of a connection between young stars and at least a fraction of such events.

preprint2006arXiv

The multiwavelength afterglow of GRB 050721: a puzzling rebrightening seen in the optical but not in the X-ray

Context. GRB050721 was detected by Swift and promptly followed-up, in the X-ray by Swift itself and, in the optical band, by the VLT operated, for the first time, in rapid response mode. A multiwavelength monitoring campaign was performed in order to study its afterglow behavior. Aims. We present the analysis of the early and late afterglow emission in both the X-ray and optical bands, as observed by Swift, a robotic telescope, and the VLT. We compare early observations with late afterglow observations obtained with Swift and the VLT in different bands in order to constrain the density of the medium in which the fireball is expanding. Methods. We have analyzed both the X-ray and the optical light curves and compared the spectral energy distribution of the afterglow at two different epochs. Results. We observed an intense rebrightening in the optical band at about one day after the burst which was not seen in the X-ray band. This is the first observation of a GRB afterglow in which a rebrightening is observed in the optical but not in the X-ray band. The lack of detection in X-ray of such a strong rebrightening at lower energies can be described with a variable external density profile. In such a scenario, the combined X-ray and optical observations allow us to derive the matter density at 10^15 cm from the burst. This is about a factor of 10 higher than in the inner region.

preprint2005arXiv

An origin for short g-ray bursts unassociated with current star formation

Two short (<2 s) g-ray bursts (GRBs) have recently been localized and fading afterglow counterparts detected. The combination of these two results left unclear the nature of the host galaxies of the bursts, because one was a star-forming dwarf, while the other was probably an elliptical galaxy. Here we report the X-ray localization of a short burst (GRB 050724) with unusual g-ray and X-ray properties. The X-ray afterglow lies off the centre of an elliptical galaxy at a redshift of z=0.258, coincident with the position determined by ground-based optical and radio observations. The low level of star formation typical for elliptical galaxies makes it unlikely that the burst originated n a supernova explosion. A supernova origin was also ruled out for GRB 050709, even though that burst took place in a galaxy with current star formation. The isotropic energy for the short bursts is 2-3 orders of magnitude lower han that for the long bursts. Our results therefore suggest that an alternative source of bursts -- the of binary systems of neutron stars or a neutron star-black hole pair -- are the progenitors of short bursts.

preprint2005arXiv

Evidence for intrinsic absorption in the Swift X-ray afterglows

Gamma-ray burst (GRB) progenitors are observationally linked to the death of massive stars. X-ray studies of the GRB afterglows can deepen our knowledge of the ionization status and metal abundances of the matter in the GRB environment. Moreover, the presence of local matter can be inferred through its fingerprints in the X-ray spectrum, i.e. the presence of absorption higher than the Galactic value. A few studies based on BeppoSAX and XMM-Newton found evidence of higher than Galactic values for the column density in a number of GRB afterglows. Here we report on a systematic analysis of 17 GRBs observed by Swift up to April 15, 2005. We observed a large number of GRBs with an excess of column density. Our sample, together with previous determinations of the intrinsic column densities for GRBs with known redshift, provides evidence for a distribution of absorption consistent with that predicted for randomly occurring GRB within molecular clouds.

preprint2005arXiv

GRB 050904 at redshift 6.3: observations of the oldest cosmic explosion after the Big Bang

We present optical and near-infrared observations of the afterglow of the gamma-ray burst GRB 050904. We derive a photometric redshift z = 6.3, estimated from the presence of the Lyman break falling between the I and J filters. This is by far the most distant GRB known to date. Its isotropic-equivalent energy is 3.4x10^53 erg in the rest-frame 110-1100 keV energy band. Despite the high redshift, both the prompt and the afterglow emission are not peculiar with respect to other GRBs. We find a break in the J-band light curve at t_b = 2.6 +- 1.0 d (observer frame). If we assume this is the jet break, we derive a beaming-corrected energy E_gamma = (4-12)x10^51 erg. This limit shows that GRB 050904 is consistent with the Amati and Ghirlanda relations. This detection is consistent with the expected number of GRBs at z > 6 and shows that GRBs are a powerful tool to study the star formation history up to very high redshift.

preprint2005arXiv

GRB environment properties through X and optical afterglow observations

We present the spectral analysis of 14 gamma-ray bursts (GRB) X-ray afterglows in order to investigate the properties of interstellar matter (ISM) along the line of sight of GRB. We carried out a simultaneous analysis of the NIR-optical and X-band for those afterglows with an optical counterpart too, in order to evaluate and strongly constrain the absorption effect on the spectral energy distribution due to dust extinction from GRB environment. We evaluated the equivalent hydrogen column density Nh from X-ray spectroscopy and rest frame visual extinction Av by assuming different type of ISM composition and dust grain size distribution. From our analysis we obtained a distribution of the GRB rest frame consistent with the one expected if GRB were embedded in a Galactic-like molecular cloud. Moreover, values of the visual extinction estimated from the simultaneous analysis of NIR-to-X band favour an environment where small dust grain are destroyed by the interaction with the X-ray and UV photons from GRB.

preprint2005arXiv

GRB990123: multiwavelength afterglow study

We report on the BeppoSAX and multiwavelength data analysis of the afterglow of Gamma-Ray Burst (GRB) 990123. Mainly due to its exceptional brightness, this is the only source for which the Wide Field Cameras have allowed an early detection of the X-ray afterglow between about 20 and 60 min after the GRB trigger. Besides, again for the first time, high-energy emission from the afterglow was detected up to 60 keV. The backwards extrapolation of the afterglow decay smoothly reconnects with the late GRB emission, thus indicating that both are consistent with being produced by the same phenomenon. For the X-ray afterglow we found a power-law decay with index alpha_X = 1.46 +/- 0.04; the spectrum has a power-law shape with photon index Gamma about 2. An extensive set of multiwavelength observations on the afterglow, collected from the literature and made during the BeppoSAX pointing, allowed constructing a Spectral Flux Distribution. We performed an analysis of the GRB afterglow emission in the context of the "fireball" model.

preprint2005arXiv

Swift observations of GRB050128: the early X-ray afterglow

Swift discovered GRB050128 with the Burst Alert Telescope and promptly pointed its narrow field instruments to monitor the afterglow. X-ray observations started 108 s after the trigger time. The early decay of the afterglow is relatively flat with a temporal decay modeled with a power law with index ~ -0.3. A steepening occurs at later times (~ 1500 s) with a power law index of ~ -1.3. During this transition, the observed X-ray spectrum does not change. We interpret this behaviour as either an early jet break or evidence for a transition from the fast cooling regime to the slow cooling regime in a wind environment.

preprint2005arXiv

The GRB of 1999 January 23: prompt emission and broad-band afterglow modeling

We report on BeppoSAX simultaneous X- and gamma-ray observations of the bright gamma-ray burst (GRB) 990123. We present the broad-band spectrum of the prompt emission, including optical, X- and gamma-rays, confirming the suggestion that the emission mechanisms at low and high frequencies must have different physical origins. We discuss the X-ray afterglow observed by the Narrow Field Instruments (NFIs) on board BeppoSAX and its hard X-ray emission up to 60 keV several hours after the burst, in the framework of the standard fireball model. The effects of including an important contribution of Inverse Compton scattering or modifying the hydrodynamics are studied.

preprint2004arXiv

Probing a GRB progenitor at a redshift of z=2: a comprehensive observing campaign of the afterglow of GRB 030226

We report results from a comprehensive follow-up observing campaign of the afterglow of GRB 030226, including VLT spectroscopy, VLT polarimetry, and Chandra X-ray observations. In addition, we present BOOTES-1 wide-field observations at the time of the occurrence of the burst. First observations at ESO started 0.2 days after the event when the GRB afterglow was at a magnitude of R~19 and continued until the afterglow had faded below the detection threshold (R>26). No underlying host galaxy was found. The optical light curve shows a break around 0.8 days after the burst, which is achromatic within the observational errors, supporting the view that it was due to a jetted explosion. Close to the break time the degree of linear polarization of the afterglow light was less than 1.1%, which favors a uniform jet model rather than a structured one. VLT spectra show two absorption line systems at redshifts z=1.962+/-0.001 and at z=1.986+/-0.001, placing the lower limit for the redshift of the GRB close to 2. We emphasize that the kinematics and the composition of the absorbing clouds responsible for these line systems is very similar to those observed in the afterglow of GRB 021004. This corroborates the picture in which at least some GRBs are physically related to the explosion of a Wolf-Rayet star.

preprint2003arXiv

Evidence for Supernova Signatures in the Spectrum of the Late-time Bump of the Optical Afterglow of GRB 021211

We present photometric and spectroscopic observations of the gamma-ray burst GRB 021211 obtained during the late stages of its afterglow. The light curve shows a rebrightening occurring ~25 days after the GRB. The analysis of a VLT spectrum obtained during the bump (27 days after the GRB) reveals a suggestive resemblance with the spectrum of the prototypical type-Ic SN 1994I, obtained about ~10 days past maximum light. Particularly we have measured a strong, broad absorption feature at 3770 A, which we have identified with Ca II blueshifted by ~14400 km/s, thus indicating that a supernova (SN) component is indeed powering the `bump' in the afterglow decay. Assuming SN 1994I as a template, the spectroscopic and photometric data together indicate that the SN and GRB explosions were at most separated by a few days. Our results suggest that GRBs might be associated also to standard type-Ic supernovae.

preprint2003arXiv

The dark burst 010214 with BeppoSAX: possible variable absorption and jet X-ray emission

We report on the prompt and afterglow emission observations of the dark burst GRB010214 with BeppoSAX. The prompt emission shows possible evidence of variable absorption from N_H = 3.0^{+5.1}_{-2.0} x 10^23 cm-2 in the first 6 s of the event to a value consistent with the Galactic column density (N_H^(G) = 2.66 x 10^20 cm-2) in the GRB direction. An X-ray afterglow emission in the 2-10 keV energy band was detected with BeppoSAX, but an analogue search at lower wavelengths (optical, IR and radio) was unsuccessful. The X-ray afterglow spectrum is consistent with a power-law with Galactic absorption. The light curve shows a complex decay, if the tail of the prompt emission is assumed as the onset of the afterglow: if the origin of the afterglow is coincident with the GRB onset, a bump before ~ 3 x 10^4 s is inferred, while if the afterglow is assumed to start later, a steepening of the power-law light curve at t ~ 3 x 10^4 s is deduced. We discuss these results in the light of the current models of afterglows and the possible origin of the GRB darkness. Finally, we tentatively derive an estimate of the burst redshift.

preprint2002arXiv

The bright Gamma-Ray Burst of February 10, 2000: a case study of an optically dark GRB

The gamma-ray burst GRB000210 had the highest gamma-ray peak flux of any event localized by BeppoSAX as yet but it did not have a detected optical afterglow. It is therefore one of the events recently classified as dark GRBs or GHOST (GRB Hiding Optical Source Transient), whose origin is still unclear. Chandra observations allowed us to localize this GRB within ~1" and a radio transient was detected with the VLA. We identify the likely (P=0.01) host galaxy of this burst at z=0.846. The X-ray spectrum of the afterglow shows intrinsic absorption N_H=5x10**21 cm-2. The amount of dust needed to absorb the optical flux of this object is consistent with the above HI column density, given a dust-to-gas ratio similar to that of our Galaxy. We do not find evidence for a partially ionized absorber expected if the absorption takes place in a Giant Molecular Cloud. We therefore conclude that either the gas is local to the GRB, but is condensed in small-scale high-density (n>~10**9 cm-3) clouds, or that the GRB is located in a dusty, gas-rich region of the galaxy. Finally, if GRB000210 lies at z>5, its X-ray absorbing medium would have to be substantially different from that observed in GRBs with optical afterglows.

preprint2001arXiv

GRB 990704: the most X-ray rich BeppoSAX gamma-ray burst

We present the X- and $γ$-ray detection of GRB 990704 and the discovery and study of its X-ray afterglow, 1SAX J1219.5-0350. Two pointed BeppoSAX observations with the narrow field instruments were performed on this source, separated in time by one week. The decay of the X-ray flux within the first observation appears unusually slow, being best-fit by a power law with negative index 0.83$\pm$0.16. Such a slow decay is consistent with the non-detection in our second observation, but its back-extrapolation to the time of the GRB largely underestimates the detected GRB X-ray prompt emission. In addition, the GRB prompt event shows, among the BeppoSAX-WFC detected sample, unprecedentedly high ratios of X- and gamma-ray peak fluxes (F$_{2-10 keV}$/F$_{40-700 keV}$$\sim$0.6, and F$_{2-26 keV}$/F$_{40-700 keV}$$\sim$1.6) and fluences (S$_{2-10 keV}$/S$_{40-700 keV}$$\sim$1.5 and S$_{2-26 keV}$/S$_{40-700 keV}$$\sim$2.8), making it, among the BeppoSAX arcminute-localized GRBs, the closest to the recently discovered class of Fast X-ray Transients.

preprint2000arXiv

Discovery of a redshifted Iron K-line in the X-ray afterglow of GRB 000214

We report the detection (3 sigma significance level) of a strong iron emission line in the X-ray spectrum of the afterglow of GRB 000214 (``Valentine's Day Burst'') observed by BeppoSAX. An emission line feature was observed with a centroid energy of 4.7+/-0.2 keV which, if interpreted as K-alpha emission from hydrogen-like iron, corresponds to a redshift of z=0.47. The observed intensity (EW ~2 keV) and duration (tens of hours) of the line give information on the distance, from the burst region, of the emitting material (R>3x10^15 cm) and its mass (M>1.4 solar mass). These results are not easily reconciled with the binary merger and hypernova models for gamma ray bursts.

preprint1999arXiv

The X-ray, Optical and Infrared Counterpart to GRB 980703

We report on X-ray, optical and infrared follow-up observations of GRB 980703. We detect a previously unknown X-ray source in the GRB error box; assuming a power law decline we find for its decay index alpha < -0.91 (3-sigma). We invoke host galaxy extinction to match the observed spectral slope with the slope expected from `fireball' models. We find no evidence for a spectral break in the infrared to X-ray spectral range on 1998 July 4.4, and determine a lower limit of the cooling break frequency: nu_c > 1.3 x 10^{17} Hz. For this epoch we obtain an extinction of A_V = 1.50 +/- 0.11. From the X-ray data we estimate the optical extinction to be A_V = 20.2 +12.3 -7.3, inconsistent with the former value. Our optical spectra confirm the redshift of z = 0.966 found by Djorgovski et al. (1998). We compare the afterglow of GRB 980703 with that of GRB 970508 and find that the fraction of the energy in the magnetic field, epsilon_B < 6x10^{-5}, is much lower in the case of GRB 980703, which is a consequence of the high frequency of the cooling break.