Source author record

S. Covino

S. Covino 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

238works
8topics
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

238 published item(s)

preprint2022arXiv

A misfired outburst in the neutron star X-ray binary Centaurus X-4

We report on a long-term optical monitoring of the neutron star X-ray binary Centaurus X-4 performed during the last 13.5 years. This source has been in quiescence since its outburst in 1979. Our monitoring reveals the overall evolution of the accretion disc; we detect short-duration flares, likely originating also in the disc, superimposed with a small-amplitude (< 0.1 mag) ellipsoidal modulation from the companion star due to geometrical effects. A long-term (~2300 days) downward trend, followed by a shorter (~1000 days) upward one, is observed in the disc light curve. Such a rise in the optical has been observed for other X-ray binaries preceding outbursts, as predicted by the disc instability model. For Cen X-4, the rise of the optical flux proceeded for ~3 years, and culminated in a flux increase at all wavelengths (optical-UV-X-rays) at the end of 2020. This increase faded after ~2 weeks, without giving rise to a full outburst. We suggest that the propagation of an inside-out heating front was ignited due to a partial ionization of hydrogen in the inner disc. The propagation might have stalled soon after the ignition due to the increasing surface density in the disc that the front encountered while propagating outwards. The stall was likely eased by the low level irradiation of the outer regions of the large accretion disc, as shown by the slope of the optical/X-ray correlation, suggesting that irradiation does not play a strong role in the optical, compared to other sources of emission.

preprint2022arXiv

CUBES, the Cassegrain U-Band Efficient Spectrograph

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028.

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.

preprint2022arXiv

Short timescale imaging polarimetry of geostationary satellite Thor-6: the nature of micro-glints

Large constellations of orbiting communication satellites will become an important source of noise for present and future astronomical observatories. Mitigation measures rely on high quality predictive models of the position and expected brightness of these objects. Optical linear imaging polarimetry holds promise as a quantitative tool to improve our understanding of the physics of reflection of sunlight off satellite components and through which models of expected brightness can be improved. We present the first simultaneous short-timescale linear polarimetry and optical photometry observations of a geostationary satellite, using the new MOPTOP imaging polarimeter on the 2m Liverpool Telescope. Our target, telecommunication satellite Thor-6, shows prominent short timescale glint-like features in the lightcurve, some as short as seconds. Our polarimetric observations overlap with several of these micro-glints, and have the cadence required to resolve them. We find that the polarisation lightcurve is remarkably smooth, the short time scale glints are not seen to produce strong polarimetric features in our observation. We show how short timescale polarimetry can further constrain the properties of the components responsible for these micro-glints.

preprint2022arXiv

The peculiar short-duration GRB 200826A and its supernova

Gamma-ray bursts (GRBs) are classified as long and short events. Long GRBs (LGRBs) are associated with the end states of very massive stars, while short GRBs (SGRBs) are linked to the merger of compact objects. GRB 200826A was a peculiar event, because by definition it was a SGRB, with a rest-frame duration of ~ 0.5 s. However, this event was energetic and soft, which is consistent with LGRBs. The relatively low redshift (z = 0.7486) motivated a comprehensive, multi-wavelength follow-up campaign to characterize its host, search for a possible associated supernova (SN), and thus understand the origin of this burst. To this aim we obtained a combination of deep near-infrared (NIR) and optical imaging together with spectroscopy. Our analysis reveals an optical and NIR bump in the light curve whose luminosity and evolution is in agreement with several LGRB-SNe. Analysis of the prompt GRB shows that this event follows the $E_{\rm p,i}-E_{\rm iso}$ relation found for LGRBs. The host galaxy is a low-mass star-forming galaxy, typical for LGRBs, but with one of the highest star-formation rates (SFR), especially with respect to its mass ($\log M_\ast/M_\odot = 8.6$, SFR $\sim 4.0 \,M_\odot$/yr). We conclude that GRB 200826A is a typical collapsar event in the low tail of the duration distribution of LGRBs. These findings support theoretical predictions that events produced by collapsars can be as short as 0.5 s in the host frame and further confirm that duration alone is not an efficient discriminator for the progenitor class of a GRB.

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

eROSITA detection of flares from the Be/X-ray binary A0538-66

In 2018, XMM-Newton observed the awakening in X-rays of the Be/X-ray binary (Be/XRB) A0538-66. It showed bright and fast flares close to periastron with properties that had never been observed in other Be/XRBs before. We report the results from the observations of A0538-66 collected during the first all-sky survey of eROSITA, an X-ray telescope (0.2-10 keV) on board the Spektrum-Roentgen-Gamma (SRG) satellite. eROSITA caught two flares within one orbital cycle at orbital phases $ϕ= 0.29$ and $ϕ= 0.93$ (where $ϕ=0$ corresponds to periastron), with peak luminosities of $\sim 2-4 \times 10^{36}$ erg/s (0.2-10 keV) and durations of $42 \leq Δt_{\rm fl} \leq 5.7\times 10^4$ s. The flare observed at $ϕ\approx 0.29$ shows that the neutron star can accrete considerably far from periastron, although it is expected to be outside of the circumstellar disk, thus providing important new information about the plasma environment surrounding the binary system. We also report the results from the photometric monitoring of A0538-66 carried out with the REM, OGLE, and MACHO telescopes from January 1993 until March 2020. We found that the two sharp peaks that characterize the orbital modulation in the optical occur asymmetrically in the orbit, relative to the position of the donor star.

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 for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology.

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 comparison between short GRB afterglows and AT2017gfo: shedding light on kilonovae properties

Multi-messenger astronomy received a great boost following the discovery of kilonova AT2017gfo, the optical counterpart of the gravitational wave source GW170817 associated with the short gamma-ray burst GRB 170817A. AT2017gfo was the first kilonova that could be extensively monitored in time both photometrically and spectroscopically. Previously, only few candidates have been observed against the glare of short GRB afterglows. In this work, we aim to search the fingerprints of AT2017gfo-like kilonova emissions in the optical/NIR light curves of 39 short GRBs with known redshift. For the first time, our results allow us to study separately the range of luminosity of the blue and red components of AT2017gfo-like kilonovae in short GRBs. In particular, the red component is similar in luminosity to AT2017gfo, while the blue kilonova can be more than 10 times brighter. Finally, we find further evidence to support all the claimed kilonova detections and we exclude an AT2017gfo-like kilonova in GRBs 050509B and 061201.

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

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

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

Lyman continuum leakage in faint star-forming galaxies at redshift z=3-3.5 probed by gamma-ray bursts

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$α$ emission. From its Ly$α$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$α$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged]

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

Search for the optical counterpart of the GW170814 gravitationalwave event with the VLT Survey Telescope

We report on the search for the optical counterpart of the gravitational event GW170814, which was carried out with the VLT Survey Telescope (VST) by the GRAvitational Wave Inaf TeAm (GRAWITA). Observations started 17.5 hours after the LIGO/Virgo alert and we covered an area of 99 deg$^2$ that encloses $\sim 77\%$ and $\sim 59\%$ of the initial and refined localization probability regions, respectively. A total of six epochs were secured over nearly two months. The survey reached an average limiting magnitude of 22 AB mag in the $r-$band. After assuming the model described in Perna et al. 2019 that derives as possible optical counterpart of a BBH event a transient source declining in about one day, we have computed a survey efficiency of about $5\%$. This paper describes the VST observational strategy and the results obtained by our analysis pipelines developed to search for optical transients in multi-epoch images. We report the catalogue of the candidates with possible identifications based on light-curve fitting. We have identified two dozens of SNe, nine AGNs, one QSO. Nineteen transients characterized by a single detection were not classified. We have restricted our analysis only to the candidates that fall into the refined localization map. None out of 39 left candidates could be positively associated with GW170814. This result implies that the possible emission of optical radiation from a BBH merger had to be fainter than r $\sim$ 22 ($L_{optical}$ $\sim$ $1.4 \times 10^{42}$ erg/s) on a time interval ranging from a few hours up to two months after the GW event.

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

preprint2019arXiv

GRB 190114C in the nuclear region of an interacting galaxy -- A detailed host analysis using ALMA, HST and VLT

GRB 190114C is the first GRB for which the detection of very-high energy emission up to the TeV range has been reported. It is still unclear whether environmental properties might have contributed to the production of these very high-energy photons, or if it is solely related to the released GRB emission. The relatively low redshift of the GRB (z=0.425) allows us to study the host galaxy of this event in detail, and to potentially identify idiosyncrasies that could point to progenitor characteristics or environmental properties responsible for such a unique event. We use ultraviolet, optical, infrared and submillimetre imaging and spectroscopy obtained with HST, VLT and ALMA to obtain an extensive dataset on which the analysis of the host galaxy is based. The host system is composed of a close pair of interacting galaxies (Delta v = 50 km s^-1), both of which are well-detected by ALMA in CO(3-2). The GRB occurred within the nuclear region (~170 pc from the centre) of the less massive but more star-forming galaxy of the pair. The host is more massive (log(M/M_odot)=9.3) than average GRB hosts at that redshift and the location of the GRB is rather unique. The enhanced star-formation rate was probably triggered by tidal interactions between the two galaxies. Our ALMA observations indicate that both host galaxy and companion have a high molecular gas fraction, as has been observed before in interacting galaxy pairs. The location of the GRB within the core of an interacting galaxy with an extinguished line-of-sight is indicative of a denser environment than typically observed for GRBs and could have been crucial for the generation of the very-high-energy photons that were observed.

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

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

Different twins in the millisecond pulsar recycling scenario: optical polarimetry of PSR J1023+0038 and XSS J12270-4859

We present the first optical polarimetric study of the two transitional pulsars PSR J1023+0038 and XSS J12270-4859. This work is focused on the search for intrinsical linear polarisation (LP) in the optical emission from the two systems. We carried out multiband optical and NIR photo-polarimetry of the two systems using the ESO NTT at La Silla (Chile), equipped with the EFOSC2 and the SOFI instruments. XSS J12270-4859 was observed during its radio-pulsar state; we did not detect LP in all bands, with 3 sigma upper limits of, e.g., 1.4% in the R-band. We built the NIR-optical averaged spectral energy distribution (SED) of the system, that could be well described by an irradiated black body with radius $R_{*} = 0.33\pm0.03\,R_{\odot}$ and albedo $η=0.32\pm0.05$, without the need of further components (thus excluding the visible presence of an extended accretion disc and/or of relativistic jets). The case was different for PSR J1023+0038, that was in its accretion phase during our campaign. We measured a LP of $1.09\pm0.27\%$ and $0.90\pm 0.17\%$ in the V and R bands, respectively. The phase-resolved polarimetric curve of the source in the R-band reveals a hint of a sinusoidal modulation at the source 4.75 hr orbital period, peaked at the same orbital phase as the light curve. The measured LP of PSR J1023+0038 could in principle be interpreted as scattering with free electrons (that can be found in the accretion disc of the system or even in the hot corona that surrounds the disc itself) or to synchrotron emission from a relativistic particles jet or outflow. However, the NIR-optical SED of the system built starting from our dataset did not suggest the presence of a jet. We conclude that the optical LP observed for PSR J1023+0038 is possibly due to Thomson scattering with electrons in the disc, as also suggested from the possible modulation of the R-band LP at the system orbital period.

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

Microlensing due to both gravitation and refraction as a further probe of universe evolution

Microlensings events are predicted for the light coming from cosmological sources. In addition to the microlensing due to gravitation lensing, microlensing produced also by refraction of light due to either ionized, or not, gas clouds can be considered. A detailed prediction is here given assuming that the ray of light coming from the distant source traverses a gas cloud with a King's density profile for various possible environments. We conclude that the additional deviation due to relativistic refraction is in most cases negligible compared to the gravitational deviation. Deviation due to refraction can anyway become an interesting analysis tool for future facility with great resolving power and the effects can be singled out with dedicated surveys.

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

Optical and near-infrared photometric monitoring of the transient X-ray binary A0538-66 with REM

The transient Be/X-ray binary A0538-66 shows peculiar X-ray and optical variability. Despite numerous studies, the intrinsic properties underlying its anomalous behaviour remain poorly understood. Since 2014 September we are conducting the first quasi-simultaneous optical and near-infrared photometric monitoring of A0538-66 in seven filters with the Rapid Eye Mount (REM) telescope, aiming to understand the properties of this binary system. We found that the REM lightcurves show fast flares lasting one or two days that repeat almost regularly every ~16.6 days, the orbital period of the neutron star. If the optical flares are powered by X-ray outbursts through photon reprocessing, the REM lightcurves indicate that A0538-66 is still active in X-rays: bright X-ray flares (L_x > 1E37 erg/s) could be observable during the periastron passages. The REM lightcurves show a long-term variability that is especially pronounced in the g band and decreases with increasing wavelength, until it no longer appears in the near-infrared lightcurves. In addition, A0538-66 is fainter with respect to previous optical observations most likely due to the higher absorption of the stellar radiation of a denser circumstellar disc. On the basis of the current models, we interpret these observational results with a circumstellar disc around the Be star observed nearly edge-on during a partial depletion phase. The REM lightcurves also show short-term variability on timescales of ~1 day possibly indicative of perturbations in the density distribution of the circumstellar disc caused by the tidal interaction with the neutron star.

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

Short GRBs at the dawn of the gravitational wave era

We derive the luminosity function and redshift distribution of short Gamma Ray Bursts (SGRBs) using (i) all the available observer-frame constraints (i.e. peak flux, fluence, peak energy and duration distributions) of the large population of Fermi SGRBs and (ii) the rest-frame properties of a complete sample of Swift SGRBs. We show that a steep $ϕ(L)\propto L^{-a}$ with a>2.0 is excluded if the full set of constraints is considered. We implement a Monte Carlo Markov Chain method to derive the $ϕ(L)$ and $ψ(z)$ functions assuming intrinsic Ep-Liso and Ep-Eiso correlations or independent distributions of intrinsic peak energy, luminosity and duration. To make our results independent from assumptions on the progenitor (NS-NS binary mergers or other channels) and from uncertainties on the star formation history, we assume a parametric form for the redshift distribution of SGRBs. We find that a relatively flat luminosity function with slope ~0.5 below a characteristic break luminosity ~3$\times10^{52}$ erg/s and a redshift distribution of SGRBs peaking at z~1.5-2 satisfy all our constraints. These results hold also if no Ep-Liso and Ep-Eiso correlations are assumed. We estimate that, within ~200 Mpc (i.e. the design aLIGO range for the detection of GW produced by NS-NS merger events), 0.007-0.03 SGRBs yr$^{-1}$ should be detectable as gamma-ray events. Assuming current estimates of NS-NS merger rates and that all NS-NS mergers lead to a SGRB event, we derive a conservative estimate of the average opening angle of SGRBs: $θ_{jet}$~3-6 deg. Our luminosity function implies an average luminosity L~1.5$\times 10^{52}$ erg/s, nearly two orders of magnitude higher than previous findings, which greatly enhances the chance of observing SGRB "orphan" afterglows. Efforts should go in the direction of finding and identifying such orphan afterglows as counterparts of GW events.

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

Accessing the population of high redshift Gamma Ray Bursts

Gamma Ray Bursts (GRBs) are a powerful probe of the high redshift Universe. We present a tool to estimate the detection rate of high-z GRBs by a generic detector with defined energy band and sensitivity. We base this on a population model that reproduces the observed properties of GRBs detected by Swift, Fermi and CGRO in the hard X-ray and gamma-ray bands. We provide the expected cumulative distributions of the flux and fluence of simulated GRBs in different energy bands. We show that scintillator detectors, operating at relatively high energies (e.g. tens of keV to the MeV), can detect only the most luminous GRBs at high redshifts due to the link between the peak spectral energy and the luminosity (Ep-Liso) of GRBs. We show that the best strategy for catching the largest number of high-z bursts is to go softer (e.g. in the soft X-ray band) but with a very high sensitivity. For instance, an imaging soft X-ray detector operating in the 0.2-5 keV energy band reaching a sensitivity, corresponding to a fluence of ~10^-8 erg cm^-2, is expected to detect ~40 GRBs yr^-1 sr^-1 at z>5 (~3 GRBs yr^-1 sr^-1 at z>10). Once high-z GRBs are detected the principal issue is to secure their redshift. To this aim we estimate their NIR afterglow flux at relatively early times and evaluate the effectiveness of following them up and construct usable samples of events with any forthcoming GRB mission dedicated to explore the high-z Universe.

preprint2015arXiv

Are LGRBs biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of LGRBs. I: Stellar mass at z<1

LGRBs are associated with massive stars and are therefore linked to star formation. The conditions necessary to produce LGRBs can affect the relation between the LGRB rate and star formation. By using the power of a complete LGRB sample, our aim is to understand whether such a bias exists and, if it does, what is its origin. In this first paper, we build the SED of the z<1 host galaxies of the BAT6 LGRB sample, and determine their stellar masses from SED fitting. We compare the resulting stellar mass distribution (i) with star-forming galaxies observed in deep surveys (UltraVISTA); (ii) with semi-analitical models of the z<1 star forming galaxy population and (iii) with numerical simulations of LGRB hosts having different metallicity thresholds for the progenitor star environment. We find that at z<1 LGRBs tend to avoid massive galaxies and are powerful in selecting faint low-mass star-forming galaxies. The stellar mass distribution of the hosts is not consistent with that of the UltraVISTA star-forming galaxies weighted for their SFR. This implies that, at least at z<1, LGRBs are not unbiased tracers of star formation. To make the two distributions consistent, a much steeper faint-end of the mass function would be required, or a very shallow SFR-Mass relation for the low mass galaxy population. GRB host galaxy simulations indicates that, to reproduce the stellar mass distribution, a metallicity threshold of the order of Z_th=0.3-0.5Z_sun is necessary. Models without a metallicity threshold or with an extreme threshold of Z_th = 0.1Z_sun are excluded at z<1. The use of the BAT6 complete sample makes this result not affected by possible biases which could have influenced past results based on incomplete samples. The preference for low metallicities (Z<~0.5Z_sun) can be a consequence of the particular conditions needed for the progenitor star to produce a GRB. (Abridged)

preprint2015arXiv

Awakening of The High Redshift Blazar CGRaBS J0809+5341

CGRaBS J0809+5341, a high redshift blazar at z = 2.144, underwent a giant optical outburst on 2014 April 19 when it brightened by $\sim$5 mag and reached an unfiltered apparent magnitude of 15.7 mag. This implies an absolute magnitude of -30.5 mag, making it one of the brightest quasars in the Universe. This optical flaring triggered us to carry out observations during the decaying part of the flare covering a wide energy range using the {\it Nuclear Spectroscopic Telescope Array}, {\it Swift}, and ground based optical facilities. For the first time, the source is detected in $γ$-rays by the Large Area Telescope onboard the {\it Fermi Gamma-Ray Space Telescope}. A high optical polarization of $\sim$10\% is also observed. Using the Sloan Digital Sky Survey spectrum, accretion disk luminosity and black hole mass are estimated as $1.5 \times 10^{45}$ erg s$^{-1}$ and $10^{8.4}~M_{\odot}$ respectively. Using a single zone leptonic emission model, we reproduce the spectral energy distribution of the source during the flaring activity. This analysis suggests that the emission region is probably located outside the broad line region, and the jet becomes radiatively efficient. We also show that the overall properties of CGRaBS J0809+5341 seems not to be in agreement with the general properties observed in high redshift blazars up to now.

preprint2015arXiv

Electromagnetic follow-up of gravitational wave candidates: perspectives in INAF

The electromagnetic (EM) emission associated with a gravitational wave (GW) signal is one of the main goal of future astronomy. Merger of neutron stars and/or black holes and core-collapse of massive stars are expected to cause rapid transient electromagnetic signals. The EM follow-up of GW signals will have to deal with large position uncertainties. The gravitational sky localization is expected to be tens to hundreds of square degrees. Wide-field cameras and rapid follow-up observations will be crucial to characterize the EM candidates for the first EM counterpart identification. We present some of the activities that we are currently carrying on to optimize the response of the INAF network of facilities to expected GW triggers. The INAF network will represent an efficient operational framework capable of fast reaction on large error box triggers and direct identification and characterization of the candidates.

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

GRB host galaxies with VLT/X-Shooter: properties at 0.8 < z < 1.3

Long gamma-ray bursts (LGRBs) are associated with the death of massive stars. Their host galaxies therefore represent a unique class of objects tracing star formation across the observable Universe. Indeed, recently accumulated evidence shows that GRB hosts do not differ substantially from general population of galaxies at high (z > 2) redshifts. However, it has been long recognised that the properties of z < 1.5 hosts, compared to general star-forming population, are unusual. To better understand the reasons for the supposed difference in LGRB hosts properties at z < 1.5, we obtained VLT/X- Shooter spectra of six hosts lying in the redshift range of 0.8 < z < 1.3. Some of these hosts have been observed before, yet we still lack well constrained information on their characteristics such as metallicity, dust extinction and star formation rate. We search for emission lines in the VLT/X-Shooter spectra of the hosts and measure their fluxes. We perform a detailed analysis, estimating host average extinction, star-formation rates, metallicities and electron densities where possible. Measured quantities of our hosts are compared to a larger sample of previously observed GRB hosts at z < 2. Star-formation rates and metallicities are measured for all the hosts analyzed in this paper and metallicities are well determined for 4 hosts. The mass-metallicity relation, the fundamental metallicity relation and SFRs derived from our hosts occupy similar parameter space as other host galaxies investigated so-far at the same redshift. We therefore conclude that GRB hosts in our sample support the found discrepancy between the properties of low-redshift GRB hosts and the general population of star- forming galaxies.

preprint2015arXiv

GRB hosts through cosmic time - VLT/X-Shooter emission-line spectroscopy of 96 GRB-selected galaxies at 0.1 < z < 3.6

[Abridged] We present data and initial results from VLT/X-Shooter emission-line spectroscopy of 96 GRB-selected galaxies at 0.1<z<3.6, the largest sample of GRB host spectroscopy available to date. Most of our GRBs were detected by Swift and 76% are at 0.5<z<2.5 with a median z~1.6. Based on Balmer and/or forbidden lines of oxygen, nitrogen, and neon, we measure systemic redshifts, star formation rates (SFRs), visual attenuations, oxygen abundances (12+log(O/H)), and emission-line widths. We find a strong change of the typical physical properties of GRB hosts with redshift. The median SFR, for example, increases from ~0.6 M_sun/yr at z~0.6 up to ~15 M_sun/yr at z~2. A higher ratio of [OIII]/[OII] at higher redshifts leads to an increasing distance of GRB-selected galaxies to the locus of local galaxies in the BPT diagram. Oxygen abundances of the galaxies are distributed between 12+log(O/H)=7.9 and 12+log(O/H)=9.0 with a median of 12+log(O/H)~8.5. The fraction of GRB-selected galaxies with super-solar metallicities is around 20% at z<1 in the adopted metallicity scale. This is significantly less than the fraction of star formation in similar galaxies, illustrating that GRBs are scarce in high-metallicity environments. At z~3, sensitivity limits us to probing only the most luminous GRB hosts for which we derive metallicities of Z ~< 0.5 Z_sun. Together with a high incidence of galaxies with similar metallicity in our sample at z~1.5, this indicates that the metallicity dependence observed at low redshift will not be dominant at z~3.

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

NuSTAR and multifrequency study of the two high-redshift blazars S5 0836+710 and PKS 2149-306

The most powerful blazars are the flat spectrum radio quasars whose emission is dominated by a Compton component peaking between a few hundred keV and a few hundred MeV. We selected two bright blazars, PKS 2149-306 at redshift z=2.345 and S5 0836+710 at z=2.172, in order to observe them in the hard X-ray band with the NuSTAR satellite. In this band the Compton component is rapidly rising almost up to the peak of the emission. Simultaneous soft-X-rays and UV-optical observations were performed with the Swift satellite, while near-infrared (NIR) data were obtained with the REM telescope. To study their variability, we repeated these observations for both sources on a timescale of a few months. While no fast variability was detected during a single observation, both sources were found to be variable in the X-ray band, up to 50%, between the two observations, with larger variability at higher energies. No variability was detected in the optical/NIR band. These data together with Fermi-LAT, WISE and other literature data are then used to study the overall spectral energy distributions (SEDs) of these blazars. Although the jet non-thermal emission dominates the SED, it leaves the UV band unhidden, allowing us to detect the thermal emission of the disc and to estimate the mass of the black hole. The non-thermal emission is well reproduced by a one-zone leptonic model. The non-thermal radiative processes are synchrotron, self-Compton and external Compton using seed photons from both the broad-line region (BLR) and the torus. We find that our data are better reproduced if we assume that the location of the dissipation region of the jet, R_diss, is in-between the torus, (at R_torus), and the BLR (R_torus>R_diss>R_BLR). The observed variability is explained by changing a minimum number of model parameters by a very small amount.

preprint2015arXiv

Short Timescale Photometric and Polarimetric Behavior of two BL Lacertae Type Objects

Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as synchrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects are still unexplored. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6m Telescopio Nazionale Galileo. Several hours of almost continuous observations were obtained for both sources. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Essentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our relatively short monitoring, different regimes can be singled out. (abridged)

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

Spectrophotometric analysis of GRB afterglow extinction curves with X-shooter

In this work we use gamma-ray burst (GRB) afterglow spectra observed with the VLT/X-shooter spectrograph to measure rest-frame extinction in GRB lines-of-sight by modeling the broadband near-infrared (NIR) to X-ray afterglow spectral energy distributions (SEDs). Our sample consists of nine Swift GRBs, eight of them belonging to the long-duration and one to the short-duration class. Dust is modeled using the average extinction curves of the Milky Way and the two Magellanic Clouds. We derive the rest-frame extinction of the entire sample, which fall in the range $0 \lesssim {\it A}_{\rm V} \lesssim 1.2$. Moreover, the SMC extinction curve is the preferred extinction curve template for the majority of our sample, a result which is in agreement with those commonly observed in GRB lines-of-sights. In one analysed case (GRB 120119A), the common extinction curve templates fail to reproduce the observed extinction. To illustrate the advantage of using the high-quality X-shooter afterglow SEDs over the photometric SEDs, we repeat the modeling using the broadband SEDs with the NIR-to-UV photometric measurements instead of the spectra. The main result is that the spectroscopic data, thanks to a combination of excellent resolution and coverage of the blue part of the SED, are more successful in constraining the extinction curves and therefore the dust properties in GRB hosts with respect to photometric measurements. In all cases but one the extinction curve of one template is preferred over the others. We show that the modeled values of the extinction and the spectral slope, obtained through spectroscopic and photometric SED analysis, can differ significantly for individual events. Finally we stress that, regardless of the resolution of the optical-to-NIR data, the SED modeling gives reliable results only when the fit is performed on a SED covering a broader spectral region.

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 high-redshift gamma-ray burst GRB140515A

High-redshift gamma-ray bursts have several advantages for the study of the distant universe, providing unique information about the structure and properties of the galaxies in which they exploded. Spectroscopic identification with large ground-based telescopes has improved our knowledge of the class of such distant events. We present the multi-wavelength analysis of the high-$z$ Swift gamma-ray burst GRB140515A ($z = 6.327$). The best estimate of the neutral hydrogen fraction of the intergalactic medium (IGM) towards the burst is $x_{HI} \leq 0.002$. The spectral absorption lines detected for this event are the weakest lines ever observed in gamma-ray burst afterglows, suggesting that GRB140515A exploded in a very low density environment. Its circum-burst medium is characterised by an average extinction (A$_{\rm V} \sim 0.1$) that seems to be typical of $z \ge 6$ events. The observed multi-band light curves are explained either with a very flat injected spectrum ($p = 1.7$) or with a multi-component emission ($p = 2.1$). In the second case a long-lasting central engine activity is needed in order to explain the late time X-ray emission. The possible origin of GRB140515A from a Pop III (or from a Pop II stars with local environment enriched by Pop III) massive star is unlikely.

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

The rate and luminosity function of long Gamma Ray Bursts

We derive, adopting a direct method, the luminosity function and the formation rate of long Gamma Ray Bursts through a complete, flux-limited, sample of Swift bursts which has a high level of completeness in redshift z (~82%). We parametrise the redshift evolution of the GRB luminosity as L = L_0(1+ z)^k and we derive k = 2.5, consistently with recent estimates. The de-evolved luminosity function of GRBs can be represented by a broken power law with slopes a = -1.32 +- 0.21 and b = -1.84 +- 0.24 below and above, respectively, a characteristic break luminosity L_0,b = 10^51.45+-0.15 erg/s. Under the hypothesis of luminosity evolution we find that the GRB formation rate increases with redshift up to z~2, where it peaks, and then decreases in agreement with the shape of the cosmic star formation rate. We test the direct method through numerical simulations and we show that if it is applied to incomplete (both in redshift and/or flux) GRB samples it can misleadingly result in an excess of the GRB formation rate at low redshifts.

preprint2015arXiv

There is a short gamma-ray burst prompt phase at the beginning of each long one

We compare the prompt intrinsic spectral properties of a sample of short Gamma--ray Burst (GRB) with the first 0.3 seconds (rest frame) of long GRBs observed by Fermi/GBM. We find that short GRBs and the first part of long GRBs lie on the same E_p--E_iso correlation, that is parallel to the relation for the time averaged spectra of long GRBs. Moreover, they are indistinguishable in the E_p--L_iso plane. This suggests that the emission mechanism is the same for short and for the beginning of long events, and both short and long GRBs are very similar phenomena, occurring on different timescales. If the central engine of a long GRB would stop after ~0.3 * (1+z) seconds the resulting event would be spectrally indistinguishable from a short GRB.

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.

preprint2015arXiv

VLT/X-shooter spectroscopy of the afterglow of the Swift GRB 130606A: Chemical abundances and reionisation at $z\sim6$

The reionisation of the Universe is thought to have ended around z~6, as inferred from spectroscopy of distant bright background sources, such as quasars (QSO) and gamma-ray burst (GRB) afterglows. Furthermore, spectroscopy of a GRB afterglow provides insight in its host galaxy, which is often too dim and distant to study otherwise. We present the high S/N VLT/X-shooter spectrum of GRB130606A at z=5.913. We aim to measure the degree of ionisation of the IGM between 5.02<z<5.84 and to study the chemical abundance pattern and dust content of its host galaxy. We measured the flux decrement due to absorption at Ly$α$, $β$ and $γ$ wavelength regions. The hydrogen and metal absorption lines formed in the host galaxy were fitted with Voigt profiles to obtain column densities. Our measurements of the Ly$α$-forest optical depth are consistent with previous measurements of QSOs, but have a much smaller uncertainty. The analysis of the red damping wing yields a neutral fraction $x_{HI}<0.05$ (3$σ$). We obtain column density measurements of several elements. The ionisation corrections due to the GRB is estimated to be negligible (<0.03 dex), but larger corrections may apply due to the pre-existing radiation field (up to 0.4 dex based on sub-DLA studies). Our measurements confirm that the Universe is already predominantly ionised over the redshift range probed in this work, but was slightly more neutral at z>5.6. GRBs are useful probes of the ionisation state of the IGM in the early Universe, but because of internal scatter we need a larger statistical sample to draw robust conclusions. The high [Si/Fe] in the host can be due to dust depletion, alpha-element enhancement, or a combination of both. The very high value of [Al/Fe]=2.40+/-0.78 might connected to the stellar population history. We estimate the host metallicity to be -1.7<[M/H]<-0.9 (2%-13% of solar). (trunc.)

preprint2014arXiv

A complete sample of bright Swift short Gamma-Ray Bursts

We present a carefully selected sample of short gamma-ray bursts (SGRBs) observed by the Swift satellite up to June 2013. Inspired by the criteria we used to build a similar sample of bright long GRBs (the BAT6 sample), we selected SGRBs with favorable observing conditions for the redshift determination on ground, ending up with a sample of 36 events, almost half of which with a redshift measure. The redshift completeness increases up to about 70% (with an average redshift value of z = 0.85) by restricting to those events that are bright in the 15-150 keV Swift Burst Alert Telescope energy band. Such flux-limited sample minimizes any redshift-related selection effects, and can provide a robust base for the study of the energetics, redshift distribution and environment of the Swift bright population of SGRBs. For all the events of the sample we derived the prompt and afterglow emission in both the observer and (when possible) rest frame and tested the consistency with the correlations valid for long GRBs. The redshift and intrinsic X-ray absorbing column density distributions we obtain are consistent with the scenario of SGRBs originated by the coalescence of compact objects in primordial binaries, with a possible minor contribution (~10%-25%) of binaries formed by dynamical capture (or experiencing large natal kicks). This sample is expected to significantly increase with further years of Swift activity.

preprint2014arXiv

A magnetar powering the ordinary monster GRB 130427A?

We present the analysis of the extraordinarily bright Gamma-Ray Burst (GRB) 130427A under the hypothesis that the GRB central engine is an accretion-powered magnetar. In this framework, initially proposed to explain GRBs with precursor activity, the prompt emission is produced by accretion of matter onto a newly-born magnetar, and the observed power is related to the accretion rate. The emission is eventually halted if the centrifugal forces are able to pause accretion. We show that the X-ray and optical afterglow is well explained as the forward shock emission with a jet break plus a contribution from the spin-down of the magnetar. Our modelling does not require any contribution from the reverse shock, that may still influence the afterglow light curve at radio and mm frequencies, or in the optical at early times. We derive the magnetic field ($B\sim 10^{16}$ G) and the spin period ($P\sim 20$ ms) of the magnetar and obtain an independent estimate of the minimum luminosity for accretion. This minimum luminosity results well below the prompt emission luminosity of GRB 130427A, providing a strong consistency check for the scenario where the entire prompt emission is the result of continuous accretion onto the magnetar. This is in agreement with the relatively long spin period of the magnetar. GRB 130427A was a well monitored GRB showing a very standard behavior and, thus, is a well-suited benchmark to show that an accretion-powered magnetar gives a unique view of the properties of long GRBs.

preprint2014arXiv

A quiescent galaxy at the position of the long GRB 050219A

Long-duration gamma-ray bursts (LGRBs) are produced by the collapse of very massive stars. Due to the short lifetime of their progenitors, LGRBs pinpoint star-forming galaxies. We present here a multi-band search for the host galaxy of the long dark GRB 050219A within the enhanced Swift/XRT error circle. We used spectroscopic observations acquired with VLT/X-shooter to determine the redshift and star-formation rate of the putative host galaxy. We compared the results with the optical/IR spectral energy distribution obtained with different facilities. Surprisingly, the host galaxy is a old and quiescent early-type galaxy at z = 0.211 characterised by an unprecedentedly low specific star-formation rate. It is the first LGRB host to be also an early-type post-starburst galaxy. This is further evidence that GRBs can explode in all kind of galaxies, with the only requirement being an episode of star-formation.

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

Circular polarization in the optical afterglow of GRB 121024A

Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets.

preprint2014arXiv

Comparing the spectral lag of short and long gamma-ray bursts and its relation with the luminosity

We investigated the rest frame spectral lags of two complete samples of bright long (50) and short (6) gamma-ray bursts (GRB) detected by Swift. We analysed the Swift/BAT data through a discrete cross-correlation function (CCF) fitted with an asymmetric Gaussian function to estimate the lag and the associated uncertainty. We find that half of the long GRBs have a positive lag and half a lag consistent with zero. All short GRBs have lags consistent with zero. The distributions of the spectral lags for short and long GRBs have different average values. Limited by the small number of short GRBs, we cannot exclude at more than 2 sigma significance level that the two distributions of lags are drawn from the same parent population. If we consider the entire sample of long GRBs, we do not find evidence for a lag-luminosity correlation, rather the lag-luminosity plane appears filled on the left hand side, thus suggesting that the lag-luminosity correlation could be a boundary. Short GRBs are consistent with the long ones in the lag-luminosity plane.

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

Diversity of GRB energetics vs. SN homogeneity: supernova 2013cq associated with the gamma-ray burst 130427A

Long-duration gamma-ray bursts (GRBs) have been found to be associated with broad-lined type-Ic supernovae (SNe), but only a handful of cases have been studied in detail. Prompted by the discovery of the exceptionally bright, nearby GRB130427A (redshift z=0.3399), we aim at characterising the properties of its associated SN2013cq. This is the first opportunity to test directly the progenitors of high-luminosity GRBs. We monitored the field of the Swift long duration GRB130427A using the 3.6-m TNG and the 8.2-m VLT during the time interval between 3.6 and 51.6 days after the burst. Photometric and spectroscopic observations revealed the presence of the type Ic SN2013cq. Spectroscopic analysis suggests that SN2013cq resembles two previous GRB-SNe, SN1998bw and SN2010bh associated with GRB980425 and XRF100316D, respectively. The bolometric light curve of SN2013cq, which is significantly affected by the host galaxy contribution, is systematically more luminous than that of SN2010bh ($\sim$ 2 mag at peak), but is consistent with SN1998bw. The comparison with the light curve model of another GRB-connected SN2003dh, indicates that SN2013cq is consistent with the model when brightened by 20%. This suggests a synthesised radioactive $^{56}$Ni mass of $\sim 0.4 M_\odot$. GRB130427A/SN2013cq is the first case of low-z GRB-SN connection where the GRB energetics are extreme ($E_{\rm γ, iso} \sim 10^{54}$ erg). We show that the maximum luminosities attained by SNe associated with GRBs span a very narrow range, but those associated with XRFs are significantly less luminous. On the other hand the isotropic energies of the accompanying GRBs span 6 orders of magnitude (10$^{48}$ erg $< E_{\rm γ, iso} <$ 10$^{54}$ erg), although this range is reduced when corrected for jet collimation. The GRB total radiated energy is in fact a small fraction of the SN energy budget.

preprint2014arXiv

Engulfing a radio pulsar: the case of PSR J1023+0038

The binary millisecond radio pulsar PSR J1023+0038 has been recently the subject of multiwavelength monitoring campaigns which revealed that an accretion disc has formed around the neutron star (since 2013 June). We present here the results of X-ray and UV observations carried out by the Swift satellite between 2013 October and 2014 May, and of optical and NIR observations performed with the REM telescope, the Liverpool Telescope, the 2.1-m telescope at the San Pedro Mártir Observatory and the 1.52-m telescope at the Loiano observing station. The X-ray spectrum is well described by an absorbed power law, which is softer than the previous quiescent epoch (up to 2013 June). The strong correlation between the X-ray and the UV emissions indicates that the same mechanism should be responsible for part of the emission in these bands. Optical and infrared photometric observations show that the companion star is strongly irradiated. Double-peaked emission lines in the optical spectra provide compelling evidence for the presence of an outer accretion disc too. The spectral energy distribution from IR to X-rays is well modelled if the contributions from the companion, the disc and the intra-binary shock emission are all considered. Our extensive data set can be interpreted in terms of an engulfed radio pulsar: the radio pulsar is still active, but undetectable in the radio band due to a large amount of ionized material surrounding the compact object. X-rays and gamma-rays are produced in an intra-binary shock front between the relativistic pulsar wind and matter from the companion and an outer accretion disc. The intense spin-down power irradiates the disc and the companion star, accounting for the UV and optical emissions.

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

GRB 120422A/SN 2012bz: Bridging the Gap between Low- And High-Luminosity GRBs

At low redshift, a handful of gamma-ray bursts (GRBs) have been discovered with peak luminosities ($L_{\rm iso} < 10^{48.5}~\rm{erg\,s}^{-1}$) substantially lower than the average of the more distant ones ($L_{\rm iso} > 10^{49.5}~\rm{erg\,s}^{-1}$). The properties of several low-luminosity (low-$L$) GRBs indicate that they can be due to shock break-out, as opposed to the emission from ultrarelativistic jets. Owing to this, it is highly debated how both populations are connected, and whether there is a continuum between them. The burst at redshift $z=0.283$ from 2012 April 22 is one of the very few examples of intermediate-$L$ GRBs with a $γ$-ray luminosity of $L\sim10^{48.9}~\rm{erg\,s}^{-1}$ that have been detected up to now. Together with the robust detection of its accompanying supernova SN 2012bz, it has the potential to answer important questions on the origin of low- and high-$L$ GRBs and the GRB-SN connection. We carried out a spectroscopy campaign using medium- and low-resolution spectrographs at 6--10-m class telescopes, covering the time span of 37.3 days, and a multi-wavelength imaging campaign from radio to X-ray energies over a duration of $\sim270$ days. Furthermore, we used a tuneable filter centred at H$α$ to map star formation in the host galaxy and the surrounding galaxies. We used these data to extract and model the properties of different radiation components and incorporate spectral-energy-distribution fitting techniques to extract the properties of the host galaxy. Modelling the light curve and spectral energy distribution from the radio to the X-rays revealed the blast-wave to expand with an initial Lorentz factor of $Γ_0\sim60$, low for a high-$L$ GRB, and that the afterglow had an exceptional low peak luminosity-density of $\lesssim2\times10^{30}~\rm{erg\,s}^{-1}\,\rm{Hz}^{-1}$ in the sub-mm. [Abridged]

preprint2014arXiv

GRB 120711A: an intense INTEGRAL burst with long-lasting soft gamma-ray emission and a powerful optical flash

A long and intense gamma-ray burst (GRB) was detected by INTEGRAL on July 11 2012 with a duration of ~115s and fluence of 2.8x10^-4 erg cm^-2 in the 20 keV-8 MeV energy range. GRB 120711A was at z~1.405 and produced soft gamma-ray emission (>20 keV) for at least ~10 ks after the trigger. The GRB was observed by several ground-based telescopes that detected a powerful optical flash peaking at an R-band brightness of ~11.5 mag at ~126 s after the trigger. We present a comprehensive temporal and spectral analysis of the long-lasting soft gamma-ray emission detected in the 20-200 keV band with INTEGRAL, the Fermi/LAT post-GRB detection above 100 MeV, the soft X-ray afterglow from XMM-Newton, Chandra, and Swift and the optical/NIR detections from Watcher, Skynet, GROND, and REM. We modelled the long-lasting soft gamma-ray emission using the standard afterglow scenario, which indicates a forward shock origin. The combination of data extending from the NIR to GeV energies suggest that the emission is produced by a broken power-law spectrum consistent with synchrotron radiation. The afterglow is well modelled using a stratified wind-like environment with a density profile k~1.2, suggesting a massive star progenitor (i.e. Wolf-Rayet). The analysis of the reverse and forward shock emission reveals an initial Lorentz factor of ~120-340, a jet half-opening angle of ~2deg-5deg, and a baryon load of ~10^-5-10^-6 Msun consistent with the expectations of the fireball model when the emission is highly relativistic. Long-lasting soft gamma-ray emission from other INTEGRAL GRBs with high peak fluxes, such as GRB 041219A, was not detected, suggesting that a combination of high Lorentz factor, emission above 100 MeV, and possibly a powerful reverse shock are required. Similar long-lasting soft gamma-ray emission has recently been observed from the nearby and extremely bright Fermi/LAT burst GRB 130427A.

preprint2014arXiv

GRB 140206A: the most distant polarized Gamma-Ray Burst

The nature of the prompt gamma-ray emission of Gamma-Ray Bursts (GRBs) is still far from being completely elucidated. The measure of linear polarization is a powerful tool that can be used to put further constraints on the content and magnetization of the GRB relativistic outflows, as well as on the radiation processes at work. To date only a handful of polarization measurements are available for the prompt emission of GRBs. Here we present the analysis of the prompt emission of GRB 140206A, obtained with INTEGRAL/IBIS, Swift/BAT, and Fermi/GBM. Using INTEGRAL/IBIS as a Compton polarimeter we were able to constrain the linear polarization level of the second peak of this GRB as being larger than 28% at 90% c.l. We also present the GRB afterglow optical spectroscopy obtained at the Telescopio Nazionale Galileo (TNG), which allowed us the measure the distance of this GRB, z=2.739. This distance value together with the polarization measure obtained with IBIS, allowed us to derive the deepest and most reliable limit to date (xi <1x10-16) on the possibility of Lorentz Invariance Violation, measured through the vacuum birefringence effect on a cosmological source.

preprint2014arXiv

GRB orphan afterglows in present and future radio transient surveys

Orphan Afterglows (OA) are slow transients produced by Gamma Ray Bursts seen off-axis that become visible on timescales of days/years at optical/NIR and radio frequencies, when the prompt emission at high energies (X and gamma rays) has already ceased. Given the typically estimated jet opening angle of GRBs theta_jet ~ 3 deg, for each burst pointing to the Earth there should be a factor ~ 700 more GRBs pointing in other directions. Despite this, no secure OAs have been detected so far. Through a population synthesis code we study the emission properties of the population of OA at radio frequencies. OAs reach their emission peak on year-timescales and they last for a comparable amount of time. The typical peak fluxes (which depend on the observing frequency) are of few micro Jy in the radio band with only a few OA reaching the mJy level. These values are consistent with the upper limits on the radio flux of SN Ib/c observed at late times. We find that the OA radio number count distribution has a typical slope -1.7 at high fluxes and a flatter (-0.4) slope at low fluxes with a break at a frequency-dependent flux. Our predictions of the OA rates are consistent with the (upper) limits of recent radio surveys and archive searches for radio transients. Future radio surveys like VAST/ASKAP at 1.4 GHz should detect ~ 3x10^-3 OA deg^-2 yr-1, MeerKAT and EVLA at 8.4 GHz should see ~ 3x10^-1 OA deg-2 yr-1. The SKA, reaching the micro Jy flux limit, could see up to ~ 0.2-1.5 OA deg^-2 yr^-1. These rates also depend on the duration of the OA above a certain flux limit and we discuss this effect with respect to the survey cadence.

preprint2014arXiv

How to switch on and off a Gamma-ray burst through a magnetar

One of the most elusive features of Gamma Ray Bursts (GRBs) is the sporadic emission prior to the main prompt event observed in at least $\sim 15\%$ of cases. These precursors have spectral and temporal properties similar to the main prompt emission, and smaller, but comparable, energetics. They are separated from the main event by a quiescent time that may be extremely long and, in some cases, more than one precursor has been observed in the same burst. Precursors are still a puzzle: despite many attempts none of the proposed models can account for all the observed features. Based on the complete sample of bright long GRBs observed by Swift (BAT6), we propose a new scenario for which precursors are explained by assuming that the central GRB engine is a newly born magnetar. In this model the precursor and the prompt emission arise from accretion of matter onto the surface of the magnetar. The accretion process can be halted by the centrifugal drag exerted by the rotating magnetosphere onto the in-falling matter, allowing for multiple precursors and very long quiescent times.

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

New light on gamma-ray burst host galaxies with Herschel

Until recently, dust emission has been detected in very few host galaxies of gamma-ray bursts (GRBHs). With Herschel, we have now observed 17 GRBHs up to redshift z~3 and detected seven of them at infrared (IR) wavelengths. This relatively high detection rate (41%) may be due to the composition of our sample which at a median redshift of 1.1 is dominated by the hosts of dark GRBs. Although the numbers are small, statistics suggest that dark GRBs are more likely to be detected in the IR than their optically-bright counterparts. Combining our IR data with optical, near-infrared, and radio data from our own datasets and from the literature, we have constructed spectral energy distributions (SEDs) which span up to 6 orders of magnitude in wavelength. By fitting the SEDs, we have obtained stellar masses, dust masses, star-formation rate (SFR), and extinctions for our sample galaxies. We find that GRBHs are galaxies that tend to have a high specfic SFR (sSFR), and like other star-forming galaxies, their ratios of dust-to-stellar mass are well correlated with sSFR. We incorporate our Herschel sample into a larger compilation of GRBHs, and compare this combined sample to SFR-weighted median stellar masses of the widest, deepest galaxy survey to date. This is done in order to establish whether or not GRBs can be used as an unbiased tracer of cosmic comoving SFR density (SFRD) in the universe. In contrast with previous results, this comparison shows that GRBHs are medium-sized galaxies with relatively high sSFRs; stellar masses and sSFRs of GRBHs as a function of redshift are similar to what is expected for star-forming galaxy populations at similar redshifts. We conclude that there is no strong evidence that GRBs are biased tracers of SFRD; thus they should be able to reliably probe the SFRD to early epochs.

preprint2014arXiv

Optical and infrared polarimetry of the transient LMXB Cen X-4 in quiescence

We present the first optical and infrared polarimetric study of the low mass transient X-ray binary Cen X-4 during its quiescent phase. This work is aimed to search for an intrinsic linear polarisation component in the system emitted radiation that might be due, e.g., to synchrotron emission from a compact jet, or to Thomson scattering with free electrons in an accretion disc. Multiband (BVRI) optical polarimetric observations were obtained during two nights in 2008 at the ESO La Silla 3.6 m telescope (EFOSC2) in polarimetric mode. These observations cover about the 30% of the 15.1 hours orbital period. J-band observations were obtained in 2007 with the NICS (TNG) instrument at La Palma, for a totality of 1 hour observation. We obtained 3-sigma upper limits to the polarisation degree in all the optical bands, with the most constraining one being in the I-band (P<0.5%). No phase-correlated variability has been noticed in all the filters. The J-band observations provided a 6% upper limit on the polarisation level. The constraining upper limits to the polarisation in the optical allowed us to evaluate the contribution of the possible emission of a relativistic particles jet to the total system radiation to be less then the 10%. This is in agreement with the observation of a spectral energy distribution typical of a single black body of a K-spectral type main sequence star irradiated from the compact object. Due to the low S/N ratio it was not possible to investigate the possible dependency of the polarisation degree from the wavelength, that could be suggestive of polarisation induced by Thomson scattering of radiation with free electrons in the outer part of the accretion disc. Observations with higher S/N ratio are required to examine in depth this hypothesis, searching for significant phase-correlated variability.

preprint2014arXiv

Optical and X-ray Rest-frame Light Curves of the BAT6 sample

We present the rest-frame light curves in the optical and X-ray bands of an unbiased and complete sample of Swift long Gamma-Ray Bursts (GRBs), namely the BAT6 sample. The unbiased BAT6 sample (consisting of 58 events) has the highest level of completeness in redshift ($\sim$ 95%), allowing us to compute the rest-frame X-ray and optical light curves for 55 and 47 objects, respectively. We compute the X-ray and optical luminosities accounting for any possible source of absorption (Galactic and intrinsic) that could affect the observed fluxes in these two bands. We compare the behaviour observed in the X-ray and in the optical bands to assess the relative contribution of the emission during the prompt and afterglow phases. We unarguably demonstrate that the GRBs rest-frame optical luminosity distribution is not bimodal, being rather clustered around the mean value Log(L$_{R}$) = 29.9 $\pm$ 0.8 when estimated at a rest frame time of 12 hr. This is in contrast with what found in previous works and confirms that the GRB population has an intrinsic unimodal luminosity distribution. For more than 70% of the events the rest-frame light curves in the X-ray and optical bands have a different evolution, indicating distinct emitting regions and/or mechanisms. The X-ray light curves normalised to the GRB isotropic energy (E$_{\rm iso}$), provide evidence for X-ray emission still powered by the prompt emission until late times ($\sim$ hours after the burst event). On the other hand, the same test performed for the E$_{\rm iso}$-normalised optical light curves shows that the optical emission is a better proxy of the afterglow emission from early to late times.

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

Polarimetric and spectroscopic optical observations of the ultra-compact X-ray binary 4U 0614+091

Aims: We present a polarimetric and spectroscopic study of the persistent ultra compact X-ray binary 4U 0614+091 aimed at searching for the emission of a relativistic particle jet and at unveiling the orbital period of the system. Methods: We obtained r-band polarimetric observations with the Telescopio Nazionale Galileo (TNG) equipped with the PAOLO polarimeter and with the Nordic Optical Telescope (NOT) equipped with the ALFOSC instrument, covering ~ 2 hours and ~ 0.5 hours observations, respectively. We carried out low resolution spectroscopy of the system using the ESO Very Large Telescope equipped with FORS1 for ~ 1.5 hours (16 spectra covering the range 430-800 nm). Results: The polarimetric analysis performed starting from the TNG dataset revealed a polarisation degree in the r-band of 3 % +/- 1 %. From the NOT dataset, due to the lower S/N ratio, we could obtain only a 3 sigma upper limit of 3.4 %. From the joining of a spectroscopic and photometric analysis, through the study of the equivalent width variations of the CII 724 nm line and the r-band light curve, we could find a hint of a ~ 45 min periodicity. Conclusions: A polarisation degree P of ~ 3 % in the r-band is consistent with the emission of a relativistic particle jet, which is supposed to emit intrinsically linearly polarised synchrotron radiation. Since no variations of P with time have been detected, and the accretion disc of the system does not contain ionised hydrogen, scattering by free electrons in the accretion disc has been rejected. The period of ~ 45 min obtained through the analysis of the system light curve and of the equivalent width variations of the selected spectral line is probably linked to the presence of a hot spot or a superhump in the accretion disc, and lead to an orbital period > 1 hour for the binary system.

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.

preprint2014arXiv

Spatially-resolved dust properties of the GRB 980425 host galaxy

Gamma-ray bursts (GRBs) have been proposed as a tool for studying star formation in the Universe, so it is crucial to investigate whether their host galaxies and immediate environments are in any way special compared with other star-forming galaxies. Here we present spatially resolved maps of dust emission of the host galaxy of the closest known GRB 980425 at z=0.0085 using our new high-resolution observations from Herschel, Atacama Pathfinder Experiment (APEX), Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). We modelled the spectral energy distributions of the host and of the star-forming region displaying the Wolf-Rayet signatures in the spectrum (WR region), located 800 pc from the GRB position. The host is characterised by low dust content and a high fraction of UV-visible star formation, similar to other dwarf galaxies. These galaxies are abundant in the local universe, so it is not surprising to find a GRB in one of them, assuming the correspondence between the GRB rate and star formation. The WR region contributes substantially to the host emission at the far-infrared, millimetre, and radio wavelengths and we propose that this is a consequence of its high gas density. If dense environments are also found close to the positions of other GRBs, then the ISM density should also be considered, along with metallicity, an important factor influencing whether a given stellar population can produce a GRB.

preprint2014arXiv

The metallicity and dust content of a redshift 5 gamma-ray burst host galaxy

Observations of the afterglows of long gamma-ray bursts (GRBs) allow the study of star-forming galaxies across most of cosmic history. Here we present observations of GRB 111008A from which we can measure metallicity, chemical abundance patterns, dust-to-metals ratio and extinction of the GRB host galaxy at z=5.0. The host absorption system is a damped Lyman-alpha absorber (DLA) with a very large neutral hydrogen column density of log N(HI)/cm^(-2) = 22.30 +/- 0.06, and a metallicity of [S/H]= -1.70 +/- 0.10. It is the highest redshift GRB with such a precise metallicity measurement. The presence of fine-structure lines confirms the z=5.0 system as the GRB host galaxy, and makes this the highest redshift where Fe II fine-structure lines have been detected. The afterglow is mildly reddened with A_V = 0.11 +/- 0.04 mag, and the host galaxy has a dust-to-metals ratio which is consistent with being equal to or lower than typical values in the Local Group.

preprint2014arXiv

The mysterious optical afterglow spectrum of GRB140506A at z=0.889

Context. Gamma-ray burst (GRBs) afterglows probe sightlines to star-forming regions in distant star-forming galaxies. Here we present a study of the peculiar afterglow spectrum of the z = 0.889 Swift GRB 140506A. Aims. Our aim is to understand the origin of the very unusual properties of the absorption along the line-of-sight. Methods. We analyse spectroscopic observations obtained with the X-shooter spectrograph mounted on the ESO/VLT at two epochs 8.8 h and 33 h after the burst as well as imaging from the GROND instrument. We also present imaging and spectroscopy of the host galaxy obtained with the Magellan telescope. Results. The underlying afterglow appears to be a typical afterglow of a long-duration GRB. However, the material along the line-of- sight has imprinted very unusual features on the spectrum. Firstly, there is a very broad and strong flux drop below 8000 AA (4000 AA in the rest frame), which seems to be variable between the two spectroscopic epochs. We can reproduce the flux-drops both as a giant 2175 AA extinction bump and as an effect of multiple scattering on dust grains in a dense environment. Secondly, we detect absorption lines from excited H i and He i. We also detect molecular absorption from CH+ . Conclusions. We interpret the unusual properties of these spectra as reflecting the presence of three distinct regions along the line-of-sight: the excited He i absorption originates from an H ii-region, whereas the Balmer absorption must originate from an associated photodissociation region. The strong metal line and molecular absorption and the dust extinction must originate from a third, cooler region along the line-of-sight. The presence of (at least) three separate regions is reflected in the fact that the different absorption components have different velocities relative to the systemic redshift of the host galaxy.

preprint2014arXiv

VLT/X-shooter spectroscopy of the GRB 120327A afterglow

We present a study of the environment of the Swift long gamma-ray burst GRB 120327A at z ~2.8 through optical spectroscopy of its afterglow. We analyzed medium-resolution, multi-epoch spectroscopic observations (~7000 - 12000, corresponding to ~ 15 - 23 km/s, S/N = 15- 30 and wavelength range 3000-25000AA) of the optical afterglow of GRB 120327A, taken with X-shooter at the VLT 2.13 and 27.65 hr after the GRB trigger. The first epoch spectrum shows that the ISM in the GRB host galaxy at z = 2.8145 is extremely rich in absorption features, with three components contributing to the line profiles. The hydrogen column density associated with GRB 120327A has log NH / cm^(-2) = 22.01 +/- 0.09, and the metallicity of the host galaxy is in the range [X/H] = -1.3 to -1.1. In addition to the ground state lines, we detect absorption features associated with excited states of CII, OI, SiII, FeII, and NiII, which we used to derive information on the distance between the host absorbing gas and the site of the GRB explosion. The variability of the FeI\lambda2396 excited line between the two epochs proves that these features are excited by the GRB UV flux. Moreover, the distance of component I is found to be dI=200+100-60 pc, while component II is located closer to the GRB, at dII=100+40-30 pc. These values are among the lowest found in GRBs. Component III does not show excited transitions, so it should be located farther away from the GRB. The presence of H2 molecules is firmly established, with a molecular fraction in the range f=4 X 10^(-7) - 10^(-4). This particularly low value can be attributed to the small dust content. This represents the third positive detection of molecules in a GRB environment.

preprint2013arXiv

A Complete Sample of Long Bright Swift GRBs

Starting from the Swift sample we define a complete sub-sample of 58 bright long Gamma Ray Bursts (GRB), 55 of them (95%) with a redshift determination, in order to characterize their properties. Our sample (BAT6) allows us to study the properties of the long GRB population and their evolution with cosmic time. We focus in particular on the GRB luminosity function, on the spectral-energy correlations of their prompt emission, on the nature of dark bursts, on possible correlations between the prompt and the X-ray afterglow properties, and on the dust extinction.

preprint2013arXiv

A new population of ultra-long duration gamma-ray bursts

We present comprehensive multiwavelength observations of three gamma-ray bursts (GRBs) with durations of several thousand seconds. We demonstrate that these events are extragalactic transients; in particular we resolve the long-standing conundrum of the distance of GRB 101225A (the "Christmas-day burst"), finding it to have a redshift z=0.847, and showing that two apparently similar events (GRB 111209A and GRB 121027A) lie at z=0.677 and z=1.773 respectively. The systems show extremely unusual X-ray and optical lightcurves, very different from classical GRBs, with long lasting highly variable X-ray emission and optical light curves that exhibit little correlation with the behaviour seen in the X-ray. Their host galaxies are faint, compact, and highly star forming dwarf galaxies, typical of "blue compact galaxies". We propose that these bursts are the prototypes of a hitherto largely unrecognized population of ultra-long GRBs, that while observationally difficult to detect may be astrophysically relatively common. The long durations may naturally be explained by the engine driven explosions of stars of much larger radii than normally considered for GRB progenitors which are thought to have compact Wolf-Rayet progenitor stars. However, we cannot unambiguously identify supernova signatures within their light curves or spectra. We also consider the alternative possibility that they arise from the tidal disruption of stars by supermassive black holes.

preprint2013arXiv

Discovery of the broad-lined Type Ic SN 2013cq associated with the very energetic GRB 130427A

Long-duration gamma-ray bursts (GRBs) at z < 1 are in most cases found to be accompanied by bright, broad-lined Type Ic supernovae (SNe Ic-BL). The highest-energy GRBs are mostly located at higher redshifts, where the associated SNe are hard to detect observationally. Here we present early and late observations of the optical counterpart of the very energetic GRB 130427A. Despite its moderate redshift z = 0.3399+/-0.0002, GRB 130427A is at the high end of the GRB energy distribution, with an isotropic-equivalent energy release of Eiso ~ 9.6x10^53 erg, more than an order of magnitude more energetic than other GRBs with spectroscopically confirmed SNe. In our dense photometric monitoring, we detect excess flux in the host-subtracted r-band light curve, consistent with what expected from an emerging SN, ~0.2 mag fainter than the prototypical SN 1998bw. A spectrum obtained around the time of the SN peak (16.7 days after the GRB) reveals broad undulations typical of SNe Ic-BL, confirming the presence of a SN, designated SN 2013cq. The spectral shape and early peak time are similar to those of the high expansion velocity SN 2010bh associated with GRB 100316D. Our findings demonstrate that high-energy long-duration GRBs, commonly detected at high redshift, can also be associated with SNe Ic-BL, pointing to a common progenitor mechanism.

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

Dust extinction for an unbiased sample of GRB afterglows

In this paper we compute rest-frame extinctions for the afterglows of a sample of gamma-ray bursts complete in redshift. The selection criteria of the sample are based on observational high-energy parameters of the prompt emission and therefore our sample should not be biased against dusty sight-lines. It is therefore expected that our inferences hold for the general population of gamma-ray bursts. Our main result is that the optical/near-infrared extinction of gamma-ray burst afterglows in our sample does not follow a single distribution. 87% of the events are absorbed by less than 2 mag, and 50% suffer from less than 0.3-0.4 mag extinction. The remaining 13% of the afterglows are highly absorbed. The true percentage of gamma-ray burst afterglows showing high absorption could be even higher since a fair fraction of the events without reliable redshift measurement are probably part of this class. These events may be due to highly dusty molecular clouds/star forming regions associated with the gamma-ray burst progenitor or along the afterglow line of sight, and/or to massive dusty host galaxies. No clear evolution in the dust extinction properties is evident within the redshift range of our sample, although the largest extinctions are at z~1.5-2, close to the expected peak of the star formation rate. Those events classified as dark are characterized, on average, by a higher extinction than typical events in the sample. A correlation between optical/near-infrared extinction and hydrogen-equivalent column density based on X-ray studies is shown although the observed NH appears to be well in excess compared to those observed in the Local Group. Dust extinction does not seem to correlate with GRB energetics or luminosity.

preprint2013arXiv

Light from the Cosmic Frontier: Gamma-Ray Bursts

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

preprint2013arXiv

Long-term monitoring of PKS 0537-441 with Fermi-LAT and multiwavelength observations

We report on multiwavelength observations of the blazar PKS 0537-441 (z = 0.896) obtained from microwaves through gamma rays by SMA, REM, ATOM, Swift and Fermi during 2008 August-2010 April. Strong variability has been observed in gamma rays, with two major flaring episodes (2009 July and 2010 March) and a harder-when-brighter behaviour, quite common for FSRQs and low-synchrotron-peaked BL Lacs, in 2010 March. In the same way the SED of the source cannot be modelled by a simple synchrotron self-Compton model, as opposed to many BL Lacs, but the addition of an external Compton component of seed photons from a dust torus is needed. The 230 GHz light curve showed an increase simultaneous with the gamma-ray one, indicating co-spatiality of the mm and gamma-ray emission region likely at large distance from the central engine. The low, average, and high activity SED of the source could be fit changing only the electron distribution parameters, but two breaks in the electron distribution are necessary. The ensuing extra spectral break, located at NIR-optical frequencies, together with that in gamma rays seem to indicate a common origin, most likely due to an intrinsic feature in the underlying electron distribution. An overall correlation between the gamma-ray band with the R-band and K-band has been observed with no significant time lag. On the other hand, when inspecting the light curves on short time scales some differences are evident. In particular, flaring activity has been detected in NIR and optical bands with no evident gamma-ray counterparts in 2009 September and November. Moderate variability has been observed in X-rays with no correlation between flux and photon index. An increase of the detected X-ray flux with no counter part at the other wavelengths has been observed in 2008 October, suggesting once more a complex correlation between the emission at different energy bands.

preprint2013arXiv

Molecular Hydrogen in the Damped Lyman-alpha System towards GRB 120815A at z=2.36

We present the discovery of molecular hydrogen (H_2), including the presence of vibrationally-excited H_2^* in the optical spectrum of the afterglow of GRB 120815A at z=2.36 obtained with X-shooter at the VLT. Simultaneous photometric broad-band data from GROND and X-ray observations by Swift/XRT place further constraints on the amount and nature of dust along the sightline. The galactic environment of GRB 120815A is characterized by a strong DLA with log(N(H)/cm^-2) = 21.95 +/- 0.10, prominent H_2 absorption in the Lyman-Werner bands (log(N(H_2)/cm^-2) = 20.53 +/- 0.13) and thus a molecular gas fraction log f(H_2)=-1.14 +/- 0.15. The distance d between the absorbing neutral gas and GRB 120815A is constrained via photo-excitation modeling of fine-structure and meta-stable transitions of FeII and NiII to d = 0.5 +/- 0.1 kpc. The DLA metallicity ([Zn/H] = -1.15 +/- 0.12), visual extinction (A_V < 0.15 mag) and dust depletion ([Zn/Fe] = 1.01 +/- 0.10) are intermediate between the values of well-studied, H_2-deficient GRB-DLAs observed at high spectral resolution, and the approximately solar metallicity, highly-obscured and H_2-rich GRB 080607 sightline. With respect to N(H), metallicity, as well as dust-extinction and depletion, GRB 120815A is fairly representative of the average properties of GRB-DLAs. This demonstrates that molecular hydrogen is present in at least a fraction of the more typical GRB-DLAs, and H_2 and H_2^* are probably more wide-spread among GRB-selected systems than the few examples of previous detections would suggest.

preprint2013arXiv

On the relation between dynamics and geometry in GRBs

We estimate the initial bulk Lorentz factors Gamma_0 for GRBs that show the onset of the afterglow in their optical light curves. We find that Gamma_0 is strongly correlated with both the isotropic equivalent luminosity L_iso and energy E_iso and, with a larger scatter, also with the rest frame peak energy E_peak. These new correlations allow us to interpret the spectral energy correlations E_peak-L_iso (-E_iso) as a sequence of Gamma_0 factors. By accounting for the beaming effects, we find that the comoving frame properties of GRBs result clustered around typical values (e.g. L'_iso~5x10^48 erg/s). Moreover, it is theoretically predicted that there should be a link between the jet dynamics (Gamma_0) and its geometry (theta_jet). Through a population synthesis code we reconstruct the Gamma_0 and theta_jet distributions and search for a possible link between them. We find that Gamma_0 and theta_jet in GRBs should have log-normal distributions and they should be anti correlated (i.e. theta_jet^2*Gamma_0=const).

preprint2013arXiv

PAOLO: a Polarimeter Add-On for the LRS Optics at a Nasmyth focus of the TNG

We describe a new polarimetric facility available at the Istituto Nazionale di AstroFisica / Telescopio Nazionale Galileo at La Palma, Canary islands. This facility, PAOLO (Polarimetric Add-On for the LRS Optics), is located at a Nasmyth focus of an alt-az telescope and requires a specific modeling in order to remove the time- and pointing position-dependent instrumental polarization. We also describe the opto-mechanical structure of the instrument and its calibration and present early examples of applications.

preprint2013arXiv

Radio afterglows of a complete sample of bright Swift GRBs: predictions from present days to the SKA era

Radio observations of Gamma Ray Bursts afterglows are fundamental in providing insights into their physics and environment, and in constraining the true energetics of these sources. Nonetheless, radio observations of GRB afterglows are presently sparse in the time/frequency domain. Starting from a complete sample of 58 bright Swift long bursts (BAT6), we constructed a homogeneous sub-sample of 38 radio detections/upper limits which preserves all the properties of the parent sample. One half of the bursts have detections between 1 and 5 days after the explosion with typical fluxes F>100 muJy at 8.4 GHz. Through a Population SYnthesis Code coupled with the standard afterglow Hydrodynamical Emission model (PSYCHE) we reproduce the radio flux distribution of the radio sub-sample. Based on these results we study the detectability in the time/frequency domain of the entire long GRB population by present and future radio facilities. We find that the GRBs that typically trigger Swift can be detected at 8.4 GHz by JVLA within few days with modest exposures even at high redshifts. The final SKA can potentially observe the whole GRB population provided that there will be a dedicated GRB gamma-ray detector more sensitive than Swift. For a sizable fraction (50%) of these bursts, SKA will allow us to perform radio-calorimetry, after the trans-relativistic transition (occurring ~100 d), providing an estimate of the true (collimation corrected) energetics of GRBs.

preprint2013arXiv

Spectroscopy of the short-hard GRB 130603B: The host galaxy and environment of a compact object merger

Short duration gamma-ray bursts (SGRBs) are thought to be related to the violent merger of compact objects, such as neutron stars or black holes, which makes them promising sources of gravitational waves. The detection of a 'kilonova'-like signature associated to the Swift-detected GRB 130603B has suggested that this event is the result of a compact object merger. Our knowledge on SGRB has been, until now, mostly based on the absence of supernova signatures and the analysis of the host galaxies to which they cannot always be securely associated. Further progress has been significantly hampered by the faintness and rapid fading of their optical counterparts (afterglows), which has so far precluded spectroscopy of such events. Afterglow spectroscopy is the key tool to firmly determine the distance at which the burst was produced, crucial to understand its physics, and study its local environment. Here we present the first spectra of a prototypical SGRB afterglow in which both absorption and emission features are clearly detected. Together with multiwavelength photometry we study the host and environment of GRB 130603B. From these spectra we determine the redshift of the burst to be z = 0.3565+/-0.0002, measure rich dynamics both in absorption and emission, and a substantial line of sight extinction of A_V = 0.86+/-0.15 mag. The GRB was located at the edge of a disrupted arm of a moderately star forming galaxy with near-solar metallicity. Unlike for most long GRBs (LGRBs), N_HX / A_V is consistent with the Galactic ratio, indicating that the explosion site differs from those found in LGRBs. The merger is not associated with the most star-forming region of the galaxy; however, it did occur in a dense region, implying a rapid merger or a low natal kick velocity for the compact object binary.

preprint2013arXiv

Swift X-ray and ultraviolet observations of the shortest orbital period double-degenerate system RX J0806.3+1527 (HM Cnc)

The system RX J0806.3+1527 (HM Cnc) is a pulsating X-ray source with 100 per cent modulation on a period of 321.5 s (5.4 min). This period reflects the orbital motion of a close binary consisting of two interacting white dwarfs. Here we present a series of simultaneous X-ray (0.2-10 keV) and near-ultraviolet (2600 angstrom and 1928 angstrom) observations that were carried out with the Swift satellite. In the near-ultraviolet, the counterpart of RX J0806.3+1527 was detected at flux densities consistent with a blackbody with temperature 27E+3 K. We found that the emission at 2600 angstrom is modulated at the 321.5-s period with the peak ahead of the X-ray one by 0.28 cycles and is coincident within 0.05 cycles with the optical. This phase-shift measurement confirms that the X-ray hot spot (located on the primary white dwarf) is at about 80-100 degrees from the direction that connects the two white dwarfs. Albeit at lower significance, the 321.5-s signature is present also in the 1928-angstrom data; at this wavelength, however, the pulse peak is better aligned with that observed at X-rays. We use the constraints on the source luminosity and the geometry of the emitting regions to discuss the merits and limits of the main models for RX J0806.3+1527.

preprint2013arXiv

The obscured hyper-energetic GRB 120624B hosted by a luminous compact galaxy at z = 2.20

Gamma-ray bursts are the most luminous explosions that we can witness in the Universe. Studying the most extreme cases of these phenomena allows us to constrain the limits for the progenitor models. In this Letter, we study the prompt emission, afterglow, and host galaxy of GRB 120624B, one of the brightest GRBs detected by Fermi, to derive the energetics of the event and characterise the host galaxy in which it was produced. Following the high-energy detection we conducted a multi-wavelength follow-up campaign, including near-infrared imaging from HAWKI/VLT, optical from OSIRIS/GTC, X-ray observations from the Chandra X-ray Observatory and at sub-millimetre/millimetre wavelengths from SMA. Optical/nIR spectroscopy was performed with X-shooter/VLT. We detect the X-ray and nIR afterglow of the burst and determine a redshift of z = 2.1974 +/- 0.0002 through the identification of emission lines of [OII], [OIII] and H-alpha from the host galaxy of the GRB. This implies an energy release of Eiso = (3.0+/-0.2)x10^54 erg, amongst the most luminous ever detected. The observations of the afterglow indicate high obscuration with AV > 1.5. The host galaxy is compact, with R1/2 < 1.6 kpc, but luminous, at L ~ 1.5 L* and has a star formation rate of 91 +/- 6 Msol/yr as derived from H-alpha. As other highly obscured GRBs, GRB 120624B is hosted by a luminous galaxy, which we also proof to be compact, with a very intense star formation. It is one of the most luminous host galaxies associated with a GRB, showing that the host galaxies of long GRBs are not always blue dwarf galaxies, as previously thought.

preprint2013arXiv

The polarized Gamma-Ray Burst GRB 061122

We report on the polarization measure, obtained with IBIS on board INTEGRAL, of the prompt emission of GRB 061122. Over an 8 s interval containing the brightest part of the Gamma-Ray Burst we put a lower limit on its polarization fraction of 60% at 68% c.l. and of 33% at 90% c.l. on the 250-800 keV energy range. We performed late time optical and near infra-red imaging observations of the GRB field using the Telescopio Nazionale Galileo (TNG), and the Canada-France-Hawaii Telescope (CFHT). Our multi-band (ugrizYJHK) photometry allowed us to identify the host galaxy of GRB 061122 and to build its SED. Using a photometric redshift code we fitted these data, and derived the basic properties of the galaxy, including its type and redshift, that we could constrain to the interval [0.57, 2.10] at a 90% c.l., with a best fit value of z=1.33. The polarization measurement in different energy bands, together with the distance determination, allowed us to put the most stringent limit (xi < 3.4x10-16) to date to a possible Lorentz Invariance Violation based on the vacuum birefringence effect, predicted by some quantum-gravity theories.

preprint2013arXiv

X-Shooter spectroscopy of young stellar objects III. Photospheric and chromospheric properties of Class III objects

We analyzed X-Shooter/VLT spectra of 24 ClassIII sources from three nearby star-forming regions (sigmaOrionis, LupusIII, and TWHya). We determined the effective temperature, surface gravity, rotational velocity, and radial velocity by comparing the observed spectra with synthetic BT-Settl model spectra. We investigated in detail the emission lines emerging from the stellar chromospheres and combined these data with archival X-ray data to allow for a comparison between chromospheric and coronal emissions. Both X-ray and Halpha luminosity as measured in terms of the bolometric luminosity are independent of the effective temperature for early-M stars but decline toward the end of the spectral M sequence. For the saturated early-M stars the average emission level is almost one dex higher for X-rays than for Halpha: log(L_x/L_bol) = -2.85 +- 0.36 vs. log(L_Halpha/L_bol) = -3.72 +- 0.21. When all chromospheric emission lines (including the Balmer series up to H11, CaII HK, the CaII infrared triplet, and several HeI lines) are summed up the coronal flux still dominates that of the chromosphere, typically by a factor 2-5. Flux-flux relations between activity diagnostics that probe different atmospheric layers (from the lower chromosphere to the corona) separate our sample of active pre-main sequence stars from the bulk of field M dwarfs studied in the literature. Flux ratios between individual optical emission lines show a smooth dependence on the effective temperature. The Balmer decrements can roughly be reproduced by an NLTE radiative transfer model devised for another young star of similar age. Future, more complete chromospheric model grids can be tested against this data set.

preprint2013arXiv

X-Shooter spectroscopy of young stellar objects: II. Impact of chromospheric emission on accretion rate estimates

Context. The lack of knowledge of photospheric parameters and the level of chromospheric activity in young low-mass pre-main sequence stars introduces uncertainties when measuring mass accretion rates in accreting (Class II) Young Stellar Objects. A detailed investigation of the effect of chromospheric emission on the estimates of mass accretion rate in young low-mass stars is still missing. This can be undertaken using samples of young diskless (Class III) K and M-type stars. Aims. Our goal is to measure the chromospheric activity of Class III pre main sequence stars to determine its effect on the estimates of accretion luminosity (Lacc) and mass accretion rate (Macc) in young stellar objects with disks. Methods. Using VLT/X-Shooter spectra we have analyzed a sample of 24 non-accreting young stellar objects of spectral type between K5 and M9.5. We identify the main emission lines normally used as tracers of accretion in Class II objects, and we determine their fluxes in order to estimate the contribution of the chromospheric activity to the line luminosity. Results. We have used the relationships between line luminosity and accretion luminosity derived in the literature for Class II objects to evaluate the impact of chromospheric activity on the accretion rate measurements. We find that the typical chromospheric activity would bias the derived accretion luminosity by Lacc,noise< 10-3Lsun, with a strong dependence with the Teff of the objects. The noise on Macc depends on stellar mass and age, and the typical values of log(Macc,noise) range between -9.2 to -11.6Msun/yr. Conclusions. Values of Lacc< 10-3Lsun obtained in accreting low-mass pre main sequence stars through line luminosity should be treated with caution as the line emission may be dominated by the contribution of chromospheric activity.

preprint2012arXiv

A complete sample of bright Swift Gamma-Ray Bursts: X-ray afterglow luminosity and its correlation with the prompt emission

We investigate wheter there is any correlation between the X-ray afterglow luminosity and the prompt emission properties of a carefully selected sub-sample of bright Swift long Gamma-Ray Bursts (GRBs) nearly complete in redshift (~90%). Being free of selection effects (except flux limit), this sample provides the possibility to compare the rest frame physical properties of GRB prompt and afterglow emission in an unbiased way. The afterglow X-ray luminosities are computed at four different rest frame times (5 min, 1 hr, 11 hr and 24 hr after trigger) and compared with the prompt emission isotropic energy E_iso, the isotropic peak luminosity L_iso and the rest frame peak energy E_peak. We find that the rest frame afterglow X-ray luminosity do correlate with these prompt emission quantities, but the significance of each correlation decreases over time. This result is in agreement with the idea that the GRB X-ray light curve can be described as the result of a combination of different components whose relative contribution and weight change with time, with the prompt and afterglow emission dominating at early and late time, respectively. In particular, we found evidence that the plateau and the shallow decay phase often observed in GRB X-ray light curves are powered by activity from the central engine. The existence of the L_X-E_iso correlation at late times (t_rf > 11 hr) suggests a similar radiative efficiency among different bursts with on average about 6% of the total kinetic energy powering the prompt emission.

preprint2012arXiv

A complete sample of bright Swift Long Gamma Ray Bursts: testing the spectral-energy correlations

We use a nearly complete sample of Gamma Ray Bursts (GRBs) detected by the Swift satellite to study the correlations between the spectral peak energy Ep of the prompt emission, the isotropic energetics Eiso and the isotropic luminosity Liso. This GRB sample is characterized by a high level of completeness in redshift (90%). This allows us to probe in an unbiased way the issue related to the physical origin of these correlations against selection effects. We find that one burst, GRB 061021, is an outlier to the Ep-Eiso correlation. Despite this case, we find strong Ep-Eiso and Ep-Liso correlations for the bursts of the complete sample. Their slopes, normalisations and dispersions are consistent with those found with the whole sample of bursts with measured redshift and Ep. This means that the biases present in the total sample commonly used to study these correlations do not affect their properties. Finally, we also find no evolution with redshift of the Ep-Eiso and Ep-Liso correlations.

preprint2012arXiv

A complete sample of bright Swift Long Gamma-Ray Bursts: Sample presentation, Luminosity Function and evolution

We present a carefully selected sub-sample of Swift Long Gamma-ray Bursts (GRBs), that is complete in redshift. The sample is constructed by considering only bursts with favorable observing conditions for ground-based follow-up searches, that are bright in the 15-150 keV Swift/BAT band, i.e. with 1-s peak photon fluxes in excess to 2.6 ph s^-1 cm^-2. The sample is composed by 58 bursts, 52 of them with redshift for a completeness level of 90%, while another two have a redshift constraint, reaching a completeness level of 95%. For only three bursts we have no constraint on the redshift. The high level of redshift completeness allows us for the first time to constrain the GRB luminosity function and its evolution with cosmic times in a unbiased way. We find that strong evolution in luminosity (d_l=2.3\pm 0.6) or in density (d_d=1.7\pm 0.5) is required in order to account for the observations. The derived redshift distribution in the two scenarios are consistent with each other, in spite of their different intrinsic redshift distribution. This calls for other indicators to distinguish among different evolution models. Complete samples are at the base of any population studies. In future works we will use this unique sample of Swift bright GRBs to study the properties of the population of long GRBs.

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

Detailed optical and near-infrared polarimetry, spectroscopy and broadband photometry of the afterglow of GRB 091018: Polarisation evolution

[Abridged] A number of phenomena have been observed in GRB afterglows that defy explanation by simple versions of the standard fireball model, leading to a variety of new models. Polarimetry can be a major independent diagnostic of afterglow physics, probing the magnetic field properties and internal structure of the GRB jets. In this paper we present the first high quality multi-night polarimetric light curve of a Swift GRB afterglow, aimed at providing a well calibrated dataset of a typical afterglow to serve as a benchmark system for modelling afterglow polarisation behaviour. In particular, our dataset of the afterglow of GRB 091018 (at redshift z=0.971) comprises optical linear polarimetry (R band, 0.13 - 2.3 days after burst); circular polarimetry (R band) and near-infrared linear polarimetry (Ks band). We add to that high quality optical and near-infrared broadband light curves and spectral energy distributions as well as afterglow spectroscopy. The linear polarisation varies between 0 and 3%, with both long and short time scale variability visible. We find an achromatic break in the afterglow light curve, which corresponds to features in the polarimetric curve. We find that the data can be reproduced by jet break models only if an additional polarised component of unknown nature is present in the polarimetric curve. We probe the ordered magnetic field component in the afterglow through our deep circular polarimetry, finding P_circ < 0.15% (2 sigma), the deepest limit yet for a GRB afterglow, suggesting ordered fields are weak, if at all present. Our simultaneous R and Ks band polarimetry shows that dust induced polarisation in the host galaxy is likely negligible.

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

GRB 100219A with X-shooter - abundances in a galaxy at z = 4.7

Abundances of galaxies at redshifts z > 4 are difficult to obtain from damped Ly α (DLA) systems in the sightlines of quasars (QSOs) due to the Ly α forest blanketing and the low number of high-redshift quasars detected. Gamma-ray bursts (GRBs) with their higher luminosity are well suited to study galaxies out to the formation of the first stars at z > 10. Its large wavelength coverage makes the X-shooter spectrograph an excellent tool to study the interstellar medium (ISM) of high redshift galaxies, in particular if the redshift is not known beforehand. Here we determine the properties of a GRB host at z = 4.66723 from a number of resonant low- and high ionization and fine-structure absorption lines. This is one of the highest redshifts where a detailed analysis with medium-resolution data has been possible. We detect one intervening system at z = 2.18. The velocity components of the absorption lines are fitted with Voigt-profiles and we determine a metallicity of [M/H] = -1.0 \pm 0.1 using S. The absorption lines show a complicated kinematic structure which could point to a merger in progress. Si II* together with the restframe UV energy release determined from GROND data gives us the distance of 0.3 to 1 kpc of the absorbing material from the GRB. We measure a low extinction of AV = 0.24 \pm 0.06 mag using X-ray spectral information and the flux calibrated X-shooter spectrum. GRB-DLAs have a shallower evolution of metallicity with redshift than QSO absorbers and no evolution in HI column density or ionization fraction. GRB hosts at high redshift might continue the trend towards lower metallicities in the LZ-relation with redshift, but the sample is still too small to draw a definite conclusion. While the detection of GRBs at z > 4 with current satellites is still difficult, they are very important for our understanding of the early epochs of star- and galaxy-formation.

preprint2012arXiv

GRBs have preferred jet opening angles and bulk Lorentz factors

We recently found that Gamma Ray Burst energies and luminosities, in their comoving frame, are remarkably similar. This, coupled with the clustering of energetics once corrected for the collimation factor, suggests the possibility that all bursts, in their comoving frame, have the same peak energy E'peak (of the order of a few keV) and the same energetics of the prompt emission E'gamma (of the order of 2e48 erg). The large diversity of bursts energies is then due to the different bulk Lorentz factor Gamma and jet aperture angle theta_jet. We investigated, through a population synthesis code, what are the distributions of Gamma and theta_jet compatible with the observations. Both quantities must have preferred values, with log-normal best fitting distributions and <Gamma0> ~ 275 and <θ_jet> ~ 8.7 degree. Moreover, the peak values of the Gamma and theta_jet distributions must be related - theta_jet^2.5 Gamma =const: the narrower the jet angle, the larger the bulk Lorentz factor. We predict that ~6% of the bursts that point to us should not show any jet break in their afterglow light curve since they have sin(theta_jet)<1/Gamma. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and ~2.5% of local hypernovae, i.e. SNIb/c with broad absorption lines.

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

On the Environment of Short Gamma-ray Bursts

In this paper we present a sample of 10 short gamma-ray bursts (GRBs) with a robust redshift determination, discovered by the Swift satellite up to January 2011. We measure their X-ray absorbing column densities and collect data on the host galaxy offsets. We find evidence for intrinsic absorption and no correlation between the intrinsic absorbing column density and the projected offset of the GRB from its host galaxy center. We find that the properties in the gamma regime (T90, fluence and 1-s peak photon flux) of short GRBs with "bright" and "faint" X-ray afterglow likely disfavour different prompt emission mechanisms. Host galaxy offset and GRB duration (T90) do not correlate. Instead, there is a hint of anti-correlation between the effective radius normalised host galaxy offset and T90. Finally, we examine the properties of short GRBs with short-lived and long-lived X-ray afterglows, finding that some short GRBs with short-lived X-ray afterglows have their optical afterglow detected. In light of this, the X-ray afterglow duration does not seem to be an unique indicator of a specific progenitor and/or environment for short GRBs.

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 Dark Bursts population in a complete sample of bright Swift Long Gamma-Ray Bursts

We study the properties of the population of optically dark events present in a carefully selected complete sample of bright Swift long gamma-ray bursts. The high level of completeness in redshift of our sample (52 objects out of 58) allow us to establish the existence of a genuine dark population and we are able to estimate the maximum fraction of dark burst events (~30%) expected for the whole class of long gamma-ray burst. The redshift distribution of this population of dark bursts is similar to the one of the whole sample. Interestingly, the rest-frame X-ray luminosity (and the de-absorbed X-ray flux) of the sub-class of dark bursts is slightly higher than the average luminosity of the non-dark events. At the same time the prompt properties do not differ and the optical flux of dark events is at the lower tail of the optical flux distribution, corrected for Galactic absorption. All these properties suggest that dark bursts events generate in much denser environments with respect to normal bright events. We can therefore exclude the high-z and the low-density scenarios and conclude that the major cause of the origin of optically dark events is the dust extinction.

preprint2012arXiv

The impact of selection biases on the Ep-Liso correlation of Gamma Ray Bursts

We study the possible effects of selection biases on the Ep-Liso correlation caused by the unavoidable presence of flux-limits in the existing samples of Gamma Ray Bursts (GRBs). We consider a well defined complete sample of Swift GRBs and perform Monte Carlo simulations of the GRB population under different assumptions for their luminosity functions. If we assume that there is no correlation between the peak energy Ep and the isotropic luminosity Liso, we are unable to reproduce it as due to the flux limit threshold of the Swift complete sample. We can reject the null hypothesis that there is no intrinsic correlation between Ep and Liso at more than 2.7 sigma level of confidence. This result is robust against the assumptions of our simulations and it is confirmed if we consider, instead of Swift, the trigger threshold of the Batse instrument. Therefore, there must be a physical relation between these two quantities. Our simulations seem to exclude, at a lower confidence level of 1.6 sigma, the possibility that the observed Ep-Liso correlation among different bursts is caused by a boundary, i.e. such that for any given Ep, we see only the largest Liso, which has a flux above the threshold of the current instruments.

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.

preprint2012arXiv

The origin of the early time optical emission of Swift GRB 080310

We present broadband multi-wavelength observations of GRB 080310 at redshift z = 2.43. This burst was bright and long-lived, and unusual in having extensive optical and near IR follow-up during the prompt phase. Using these data we attempt to simultaneously model the gamma-ray, X-ray, optical and IR emission using a series of prompt pulses and an afterglow component. Initial attempts to extrapolate the high energy model directly to lower energies for each pulse reveal that a spectral break is required between the optical regime and 0.3 keV to avoid over predicting the optical flux. We demonstrate that afterglow emission alone is insufficient to describe all morphology seen in the optical and IR data. Allowing the prompt component to dominate the early-time optical and IR and permitting each pulse to have an independent low energy spectral indices we produce an alternative scenario which better describes the optical light curve. This, however, does not describe the spectral shape of GRB 080310 at early times. The fit statistics for the prompt and afterglow dominated models are nearly identical making it difficult to favour either. However one enduring result is that both models require a low energy spectral index consistent with self absorption for at least some of the pulses identified in the high energy emission model.

preprint2012arXiv

The X-ray absorbing column density of a complete sample of bright Swift Gamma-Ray Bursts

A complete sample of bright Swift Gamma-ray Bursts (GRBs) has been recently selected by Salvaterra et al. (2011). The sample has a high level of completeness in redshift (91%). We derive here the intrinsic absorbing X-ray column densities of these GRBs making use of the Swift X-ray Telescope data. This distribution has a mean value of log(NH/cm-2)=21.7+-0.5. This value is consistent with the distribution of the column densities derived from the total sample of GRBs with redshift. We find a mild increase of the intrinsic column density with redshift. This can be interpreted as due to the contribution of intervening systems along the line of sight. Making use of the spectral index connecting optical and X-ray fluxes at 11 hr (beta_OX), we investigate the relation of the intrinsic column density and the GRB `darkness'. We find that there is a very tight correlation between dark GRBs and high X-ray column densities. This clearly indicates that the dark GRBs are formed in a metal-rich environment where dust must be present.

preprint2011arXiv

A detailed spectral study of GRB 041219A and its host galaxy

GRB 041219A is one of the longest and brightest gamma-ray bursts (GRBs) ever observed. It was discovered by the INTEGRAL satellite, and thanks to a precursor happening about 300 s before the bulk of the burst, ground based telescopes were able to catch the rarely-observed prompt emission in the optical and in the near infrared bands. Here we present the detailed analysis of its prompt gamma-ray emission, as observed with IBIS on board INTEGRAL, and of the available X-ray afterglow data collected by XRT on board Swift. We then present the late-time multi-band near infrared imaging data, collected at the TNG, and the CFHT, that allowed us to identify the host galaxy of the GRB as an under-luminous, irregular galaxy of about 5x10^9 M_Sun at best fit redshift of z=0.31 -0.26 +0.54. We model the broad-band prompt optical to gamma-ray emission of GRB 041219A within the internal shock model. We were able to reproduce the spectra and light curve invoking the synchrotron emission of relativistic electrons accelerated by a series of propagating shock waves inside a relativistic outflow. On the other hand, it is less easy to simultaneously reproduce the temporal and spectral properties of the infrared data.

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

Constraints on Lorentz Invariance Violation using INTEGRAL/IBIS observations of GRB041219A

One of the experimental tests of Lorentz invariance violation is to measure the helicity dependence of the propagation velocity of photons originating in distant cosmological obejcts. Using a recent determination of the distance of the Gamma-Ray Burst GRB 041219A, for which a high degree of polarization is observed in the prompt emission, we are able to improve by 4 orders of magnitude the existing constraint on Lorentz invariance violation, arising from the phenomenon of vacuum birefringence.

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

Diversity of multiwavelength emission bumps in the GRB 100219A afterglow

Context. Multi-wavelength observations of gamma-ray burst (GRB) afterglows provide important information about the activity of their central engines and their environments. In particular, the short timescale variability, such as bumps and/or rebrightening features visible in the multi-wavelength light curves, is still poorly understood. Aims. We analyze the multi-wavelength observations of the GRB100219A afterglow at redshift 4.7. In particular, we attempt to identify the physical origin of the late achromatic flares/bumps detected in the X-ray and optical bands. Methods. We present ground-based optical photometric data and Swift X-ray observations on GRB100219A. We analyzed the temporal behavior of the X-ray and optical light curves, as well as the X-ray spectra. Results. The early flares in the X-ray and optical light curves peak simultaneously at about 1000 s after the burst trigger, while late achromatic bumps in the X-ray and optical bands appear at about 20000 s after the burst trigger. These are uncommon features in the afterglow phenomenology. Considering the temporal and spectral properties, we argue that both optical and X-ray emissions come from the same mechanism. The late flares/bumps may be produced by late internal shocks from long-lasting activity of the central engine. An off-axis origin for a structured jet model is also discussed to interpret the bump shapes. The early optical bump can be interpreted as the afterglow onset, while the early X-ray flare could be caused by the internal activity. GRB 100219A exploded in a dense environment as revealed by the strong attenuation of X-ray emission and the optical-to-X-ray spectral energy distribution.

preprint2011arXiv

Gamma Ray Bursts in the comoving frame

We estimate the bulk Lorentz factor Gamma_0 of 31 GRBs using the measured peak time of their afterglow light curves. We consider two possible scenarios for the estimate of Gamma_0: the case of a homogeneous circumburst medium or a wind density profile. The values of Gamma_0 are broadly distributed between few tens and several hundreds with average values ~138 and ~66 for the homogeneous and wind density profile, respectively. We find that the isotropic energy and luminosity correlate in a similar way with Gamma_0, i.e. Eiso Gamma_0^2 and Liso Gamma_0^2, while the peak energy Epeak Gamma_0. These correlations are less scattered in the wind density profile than in the homogeneous case. We then study the energetics, luminosities and spectral properties of our bursts in their comoving frame. The distribution of Liso' is very narrow with a dispersion of less than a decade in the wind case, clustering around Liso'=5x10^48 erg/s. Peak photon energies cluster around Epeak'=6 keV. The newly found correlations involving Gamma_0 offer a general interpretation scheme for the spectral-energy correlations of GRBs. The Epeak-Eiso and Epeak-Liso correlations are due to the different Gamma_0 factors and the collimation-corrected correlation, Epeak-Egamma (obtained by correcting the isotropic quantities for the jet opening angle theta_j), can be explained if theta_j^2*Gamma_0=constant. Assuming the Epeak-Egamma correlation as valid, we find a typical value of theta_j*Gamma_0 ~ 6-20, in agreement with the predictions of magnetically accelerated jet models.

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

GRB 081029: Understanding Multiple Afterglow Components

We present an analysis of the unusual optical light curve of the gamma-ray burst GRB~081029, which occurred at a redshift of z = 3.8479$. We combine X-ray and optical observations from the Swift X-Ray Telescope and the Swift UltraViolet/Optical Telescope with optical and infrared data obtained using the REM and ROTSE telescopes to construct a detailed data set extending from 86 s to approximately 100,000 s after the BAT trigger. Our data also cover a wide energy range, from 10 keV to 0.77 eV (1.24 Angstrom to 16,000 Angstrom). The X-ray afterglow shows a shallow initial decay followed by a rapid decay starting at about 18,000s. The optical and infrared afterglow, however, shows an uncharacteristic rise at about 5000 s that does not correspond to any feature in the X-ray light curve. Our data are not consistent with synchrotron radiation from a single-component jet interacting with an external medium. We do, however, find that the observed light curve can be explained using multi-component model for the jet.

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 discovery of VHE Emission from the FSRQ PKS 1222+21

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

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

Pre-ALMA observations of GRBs in the mm/submm range

GRBs generate an afterglow emission that can be detected from radio to X-rays during days, or even weeks after the initial explosion. The peak of this emission crosses the mm/submm range during the first hours to days, making their study in this range crucial for constraining the models. Observations have been limited until now due to the low sensitivity of the observatories in this range. We present observations of 10 GRB afterglows obtained from APEX and SMA, as well as the first detection of a GRB with ALMA, and put them into context with all the observations that have been published until now in the spectral range that will be covered by ALMA. The catalogue of mm/submm observations collected here is the largest to date and is composed of 102 GRBs, of which 88 had afterglow observations, whereas the rest are host galaxy searches. With our programmes, we contributed with data of 11 GRBs and the discovery of 2 submm counterparts. In total, the full sample, including data from the literature, has 22 afterglow detections with redshift ranging from 0.168 to 8.2. GRBs have been detected in mm/submm wavelengths with peak luminosities spanning 2.5 orders of magnitude, the most luminous reaching 10^33erg s^-1 Hz^-1. We observe a correlation between the X-ray brightness at 0.5 days and the mm/submm peak brightness. Finally we give a rough estimate of the distribution of peak flux densities of GRB afterglows, based on the current mm/submm sample. Observations in the mm/submm bands have been shown to be crucial for our understanding of the physics of GRBs, but have until now been limited by the sensitivity of the observatories. With the start of the operations at ALMA, the sensitivity will be increased by more than an order of magnitude. Our estimates predict that, once completed, ALMA will detect up to 98% of the afterglows if observed during the passage of the peak synchrotron emission.

preprint2011arXiv

Probing the ambient medium of GRB 090618 with XMM-Newton observations

Long Gamma-ray Bursts (GRBs) signal the death of massive stars. The afterglow emission can be used to probe the progenitor ambient through a detailed study of the absorption pattern imprinted by the circumburst material as well as the host galaxy interstellar medium on the continuum spectrum. This has been done at optical wavelengths with impressive results. Similar studies can in principle be carried out in the X-ray band, allowing us to shed light on the material metallicity, composition and distance of the absorber. We start exploiting this route through high resolution spectroscopy XMM-Newton observations of GRB 090618. We find a high metallicity absorbing medium (Z> 0.2 Zsun) with possible enhancements of S and Ne with respect to the other elements (improving the fit at a level of >3.4 sigma). Including the metallicity effects on the X-ray column density determination, the X-ray and optical evaluations of the absorption are in agreement for a Small Magellanic Cloud extinction curve.

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

Spectroscopic evidence for SN 2010ma associated with GRB 101219B

We report on the spectroscopic detection of supernova SN 2010ma associated with the long gamma-ray burst GRB 101219B. We observed the optical counterpart of the GRB on three nights with the X-shooter spectrograph at the VLT. From weak absorption lines, we measure a redshift of z = 0.55. The first epoch UV--near-infrared afterglow spectrum, taken 11.6 hr after the burst, is well fit by a power law consistent with the slope of the X-ray spectrum. The second and third epoch spectra (obtained 16.4 and 36.7 days after the burst), however, display clear bumps closely resembling those of the broad-lined type-Ic SN 1998bw if placed at z=0.55. Apart from demonstrating that spectroscopic SN signatures can be observed for GRBs at these large distances, our discovery makes a step forward in establishing a general connection between GRBs and SNe. In fact, unlike most previous unambiguous GRB-associated SNe, GRB 101219B has a large gamma-ray energy (E_iso = 4.2 x 10^{51} erg), a bright afterglow, and obeys the "Amati" relation, thus being fully consistent with the cosmological population of GRBs.

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

UVES and FORS2 spectroscopy of GRB081008

We study the GRB 081008 environment with simultaneous high- and low-resolution spectroscopy using UVES and FORS2 data acquired ~5 hr after the Swift trigger. The interstellar medium (ISM) of the host galaxy at z=1.9683 is constituted by at least three components which contribute to the line profiles. Component I is the redmost one, and is 20 km/s and 78 km/s redward component II and III, respectively. We detect several ground state and excited absorption features in components I and II. These features have been used to compute the distances between the GRB and the absorbers. Component I is found to be 52 +/-6 pc away from the GRB, while component II presents few excited transitions and its distance is 200+60-80 pc. Component III only features a few, low ionization and saturated lines suggesting that it is even farther from the GRB. The hydrogen column density associated to GRB 081008 is log N_H/cm^-2 = 21.11 +/- 0.10, and the metallicity of the host galaxy is in the range [X/H] = -1.29 to -0.52. In particular, we found [Fe/H]=-1.19 +/- 0.11 and [Zn/H]=-0.52+/- 0.11 with respect to solar values. This discrepancy can be explained by the presence of dust in the GRB ISM, given the opposite refractory properties of iron and zinc. By deriving the depletion pattern for GRB 081008, we find the optical extinction in the visual band to be A_V ~0.19 mag. The Curve of Growth analysis applied to the FORS2 spectra brings column densities consistent at the 3 sigma level to that evaluated from the UVES data using the line fitting procedure. This reflects the low saturation of the detected GRB 081008 absorption features.

preprint2011arXiv

VLT/UVES and FORS2 spectroscopy of the GRB 081008 afterglow

We study the gamma-ray burst GRB 081008 environment by analysing the spectra of its optical afterglow. UVES high resolution spectroscopy of GRB 081008 was secured ~5 hr after the Swift trigger. Our dataset comprises also three FORS2 nearly simultaneous spectra of the same source. The availability of nearly simultaneous high and low resolution spectra for a GRB afterglow is an extremely rare event. The GRB-Damped Lyman Alpha system at z = 1.9683 shows that the interstellar medium (ISM) of the host galaxy is constituted by at least three components which contribute to the line profiles. Component I is the redmost one, and is 20 km/s and 78 km/s redward component II and III, respectively. We detect several ground state and excited absorption features in components I and II. These features have been used to compute the distances between the GRB and the absorbers. Component I is found to be 52 +/- 6 pc away from the GRB, while component II presents few excited transitions and its distance is 200+60-80} pc. The hydrogen column density associated to GRB 081008 is log NH/cm^-2 = 21.11 +/- 0.10, and the metallicity of the host galaxy is in the range [X/H] = -1.29 to -0.52. In particular, we found [Fe/H]=-1.19 +/- 0.11 and [Zn/H]=-0.52 +/- 0.11 with respect to solar values. This discrepancy can be explained by the presence of dust in the GRB ISM, given the opposite refractory properties of iron and zinc. By deriving the depletion pattern for GRB 081008, we find the optical extinction in the visual band to be AV ~ 0.19 mag. The Curve of Growth analysis applied to the FORS2 spectra brings column densities consistent at the 3 sigma level to that evaluated from the UVES data using the line fitting procedure. This reflects the low saturation of the detected GRB\,081008 absorption features.

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

Discovery of the nearby long, soft GRB 100316D with an associated supernova

We report the Swift discovery of nearby long, soft gamma-ray burst GRB 100316D, and the subsequent unveiling of its low redshift host galaxy and associated supernova. We derive the redshift of the event to be z = 0.0591 +/- 0.0001 and provide accurate astrometry for the GRB-SN. We study the extremely unusual prompt emission with time-resolved gamma-ray to X-ray spectroscopy, and find that the spectrum is best modelled with a thermal component in addition to a synchrotron emission component with a low peak energy. The X-ray light curve has a remarkably shallow decay out to at least 800 s. The host is a bright, blue galaxy with a highly disturbed morphology and we use Gemini South, VLT and HST observations to measure some of the basic host galaxy properties. We compare and contrast the X-ray emission and host galaxy of GRB 100316D to a subsample of GRB-SNe. GRB 100316D is unlike the majority of GRB-SNe in its X-ray evolution, but resembles rather GRB 060218, and we find that these two events have remarkably similar high energy prompt emission properties. Comparison of the host galaxies of GRB-SNe demonstrates, however, that there is a great diversity in the environments in which GRB-SNe can be found. GRB 100316D is an important addition to the currently sparse sample of spectroscopically confirmed GRB-SNe, from which a better understanding of long GRB progenitors and the GRB--SN connection can be gleaned.

preprint2010arXiv

Evidence for Supernova-Synthesised Dust from the Rising Afterglow of GRB 071025 at z~5

We present observations and analysis of the broadband afterglow of Swift GRB 071025. Using optical and infrared (RIYJHK) photometry, we derive a photometric redshift of 4.4 < z < 5.2; at this redshift our simultaneous multicolour observations begin at ~30 s after the GRB trigger in the host frame and during the initial rising phase of the afterglow. We associate the light curve peak at 580 s in the observer frame with the formation of the forward shock, giving an estimate of the initial Lorentz factor Gamma_0 ~ 200. The red spectral energy distribution (even in regions not affected by the Lyman-alpha break) provides secure evidence of a large dust column. However, the inferred extinction curve shows a prominent flat component between 2000-3000 Angstroms in the rest-frame, inconsistent with any locally observed template but well-fit by models of dust formed by supernovae. Time-dependent fits to the extinction profile reveal no evidence of dust destruction and limit the decrease in the extinction column to Delta A_3000 < 0.54 mag after t = 50 s in the rest frame. Our observations provide evidence of a transition in dust properties at z~5, in agreement with studies of high-z quasars, and suggest that SN-formed dust continues to dominate the opacity of typical galaxies at this redshift.

preprint2010arXiv

GRB 050502B optical afterglow: a jet break at high redshift

Aims: GRB 050502B is well known for the very bright flare displayed in its X-ray light curve. Despite extensive studies, however, the optical light curve has never been discussed and its redshift is unconstrained. Possible correlations between optical and X-ray data are analysed. Methods: Photometric data from TNG in the R and I bands were used to compare the optical afterglow with the X-ray light curve. The HyperZ package and a late time VLT host observation were used to derive redshift estimates. Results: The I-band afterglow decay followed a power-law of index α = 2.1 $\pm$ 0.6, after a late break at ~ $1.3 \times 10^5$ s. The R - I color is remarkably red and the broadband spectral index β_OX = $0.9 \pm 0.1$ is consistent with the X-ray spectral slope β_X. Although a photometric redshift of z > 4 is the most conservative result to consider, a photometric redshift of z = 5.2 $\pm$ 0.3 is suggested with no extinction in the host, based on which an isotropic energy E_γ,iso = $(3.8 \pm 0.7) \times 10^52$ erg and a jet opening angle θ ~ $3.7^{\circ}$ are subsequently derived. Conclusions: The combined X-ray and optical data suggest an achromatic break, which we interpret as a jet break. The post jet break slope obeys roughly the closure relation for the jet slow cooling model. Because of the afterglow's very red color, in order for the redshift to be low (z < 1), extinction, if present in the host, must be significantly high. Since the optical-to-X-ray index is consistent with the X-ray spectrum, and there is no XRT evidence for excess N_H, GRB 050502B was likely at high redshift.

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

LOFT - a Large Observatory For x-ray Timing

The high time resolution observations of the X-ray sky hold the key to a number of diagnostics of fundamental physics, some of which are unaccessible to other types of investigations, such as those based on imaging and spectroscopy. Revealing strong gravitational field effects, measuring the mass and spin of black holes and the equation of state of ultradense matter are among the goals of such observations. At present prospects for future, non-focused X-ray timing experiments following the exciting age of RXTE/PCA are uncertain. Technological limitations are unavoidably faced in the conception and development of experiments with effective area of several square meters, as needed in order to meet the scientific requirements. We are developing large-area monolithic Silicon Drift Detectors offering high time and energy resolution at room temperature, which require modest resources and operation complexity (e.g., read-out) per unit area. Based on the properties of the detector and read-out electronics that we measured in the lab, we developed a realistic concept for a very large effective area mission devoted to X-ray timing in the 2-30 keV energy range. We show that effective areas in the range of 10-15 square meters are within reach, by using a conventional spacecraft platform and launcher of the small-medium class.

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

Optical and infrared photometry of the blazar PKS0537-441

We present a large collection of photometric data on the Blazar PKS 0537-441 in the VRIJHK bands taken in 2004-2009. At least three flare-like episodes with months duration, and >3 mag amplitude are apparent. The spectral energy distribution is consistent with a power law, and no indication of a thermal component is found. We searched for short time scale variability, and an interesting event was identified in the J band, with a duration of ~25 minutes.

preprint2010arXiv

Optical Observations of GRB Afterglows

GRB afterglows are among the best examples of astrophysical sources requiring a true multiwavelength observational approach. Radiation processes and main physical ingredients can only be disentangled with knowledge of their spectral and temporal properties through the largest possible band, i.e. from the X-ray down to radio. We now briefly review some of the most relevant observational findings recently obtained through optical observations.

preprint2010arXiv

Photometry and spectroscopy of GRB 060526: A detailed study of the afterglow and host galaxy of a z=3.2 gamma-ray burst

Aims: With this paper we want to investigate the highly variable afterglow light curve and environment of gamma-ray burst (GRB) 060526 at $z=3.221$. Methods: We present one of the largest photometric datasets ever obtained for a GRB afterglow, consisting of multi-color photometric data from the ultraviolet to the near infrared. The data set contains 412 data points in total to which we add additional data from the literature. Furthermore, we present low-resolution high signal-to-noise spectra of the afterglow. The afterglow light curve is modeled with both an analytical model using broken power law fits and with a broad-band numerical model which includes energy injections. The absorption lines detected in the spectra are used to derive column densities using a multi-ion single-component curve-of-growth analysis from which we derive the metallicity of the host of GRB 060526. Results: The temporal behaviour of the afterglow follows a double broken power law with breaks at $t=0.090\pm0.005$ and $t=2.401\pm0.061$ days. It shows deviations from the smooth set of power laws that can be modeled by additional energy injections from the central engine, although some significant microvariability remains. The broadband spectral-energy distribution of the afterglow shows no significant extinction along the line of sight. The metallicity derived from \ion{S}{II} and \ion{Fe}{II} of [S/H] = --0.57 $\pm$0.25 and [Fe/H] = --1.09$\pm$0.24 is relatively high for a galaxy at that redshift but comparable to the metallicity of other GRB hosts at similar redshifts. At the position of the afterglow, no host is detected to F775W(AB) = 28.5 mag with the HST, implying an absolute magnitude of the host M(1500 Å)$>$--18.3 mag which is fainter than most long-duration hosts, although the GRB may be associated with a faint galaxy at a distance of 11 kpc.

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

Searching for differences in Swift's intermediate GRBs

Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.

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

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 afterglow and host galaxy of GRB 090205: evidence for a Ly-alpha emitter at z=4.65

Gamma-ray bursts have been proved to be detectable up to distances much larger than any other astrophysical object, providing the most effective way, complementary to ordinary surveys, to study the high redshift universe. To this end, we present here the results of an observational campaign devoted to the study of the high-z GRB 090205. We carried out optical/NIR spectroscopy and imaging of GRB 090205 with the ESO-VLT starting from hours after the event up to several days later to detect the host galaxy. We compared the results obtained from our optical/NIR observations with the available Swift high-energy data of this burst. Our observational campaign led to the detection of the optical afterglow and host galaxy of GRB 090205 and to the first measure of its redshift, z=4.65. Similar to other, recent high-z GRBs, GRB 090205 has a short duration in the rest-frame with T_{90,rf}=1.6 s, which suggests the possibility that it might belong to the short GRBs class. The X-ray afterglow of GRB 090205 shows a complex and interesting behaviour with a possible rebrightening at 500-1000s from the trigger time and late flaring activity. Photometric observations of the GRB 090205 host galaxy argue in favor of a starburst galaxy with a stellar population younger than ~ 150 Myr. Moreover, the metallicity of Z > 0.27 Z_Sun derived from the GRB afterglow spectrum is among the highest derived from GRB afterglow measurement at high-z, suggesting that the burst occurred in a rather enriched envirorment. Finally, a detailed analysis of the afterglow spectrum shows the existence of a line corresponding to Lyman-alpha emission at the redshift of the burst. GRB 090205 is thus hosted in a typical Lyman-alpha emitter (LAE) at z=4.65. This makes the GRB 090205 host the farthest GRB host galaxy, spectroscopically confirmed, detected to date.

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 Sooner: a Large Robotic Telescope

The approach of Observational Astronomy is mainly aimed at the construction of larger aperture telescopes, more sensitive detectors and broader wavelength coverage. Certainly fruitful, this approach turns out to be not completely fulfilling the needs when phenomena related to the formation of black holes (BH), neutron stars (NS) and relativistic stars in general are concerned. Recently, mainly through the Vela, Beppo-SAX and Swift satellites, we reached a reasonable knowledge of the most violent events in the Universe and of some of the processes we believe are leading to the formation of black holes (BH). We plan to open a new window of opportunity to study the variegated physics of very fast astronomical transients, particularly the one related to extreme compact objects. The innovative approach is based on three cornerstones: 1) the design (the conceptual design has been already completed) of a 3m robotic telescope and related focal plane instrumentation characterized by the unique features: "No telescope points faster"; 2) simultaneous multi-wavelengths observations (photometry, spectroscopy o\& polarimetry); 3) high time resolution observations. The conceptual design of the telescope and related instrumentation is optimized to address the following topics: High frequency a-periodic variability, Polarization, High z GRBs, Short GRBs, GRB-Supernovae association, Multi-wavelengths simultaneous photometry and rapid low dispersion spectroscopy. This experiment will turn the "exception" (like the optical observations of GRB 080319B) to "routine".

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.

preprint2010arXiv

VLT Phase Resolved Optical Spectroscopy of the Ultra-Compact Binary HM Cnc

A 321.5 s modulation was discovered in 1999 in the X-ray light curve of HM Cnc. In 2001 and 2002, optical photometric and spectroscopic observations revealed that HM Cnc is a very blue object with no intrinsic absorptions but broad (FWHM 1500 km s^-1) low equivalent width emission lines (EW 1-6A), which were first identified with the HeII Pickering series. The combination of X-ray and optical observations pictures HM Cnc as a double degenerate binary hosting two white dwarfs, and possibly being the shortest orbital period binary discovered so far. The present work is aimed at studying the orbital motion of the two components by following the variations of the shape, centroid and intensity of the emission lines through the orbit. In February 2007, we carried out the first phase resolved optical spectroscopic study with the VLT/FORS2 in the High Time Resolution (HIT) mode, yielding five phase bins in the 321 s modulation. Despite the low SNR, the data show that the intensity of the three most prominent emission lines, already detected in 2001, varies with the phase. These lines are detected at phases 0.2-0.6 where the optical emission peaks, and marginally detected or not detected at all elsewhere. Moreover, the FWHM of the emission lines in the phase resolved spectra is smaller, by almost a factor 2, than that in the the phase-averaged 2001 spectrum. Our results are consistent with both the pulsed optical component and emission lines originating in the same region which we identify with the irradiated surface of the secondary. Moreover, regardless of the origin of the un-pulsed optical continuum, we note that the EWs of the emission lines might be up to -15 / -25A, larger than thought before; these values are more similar to those detected in cataclysmic variables. All the findings further confirm that the 321s modulation observed in HM Cnc is the orbital period of the system, the shortest known to date.

preprint2010arXiv

VLT/X-shooter spectroscopy of the GRB 090926A afterglow

The aim of this paper is to study the environment and intervening absorbers of the gamma-ray burst GRB 090926A through analysis of optical spectra of its afterglow. We analyze medium resolution spectroscopic observations (R=10 000, corresponding to 30 km/s, S/N=15 - 30 and wavelength range 3000-25000) of the optical afterglow of GRB 090926A, taken with X-shooter at the VLT ~ 22 hr after the GRB trigger. The spectrum shows that the ISM in the GRB host galaxy at z = 2.1071 is rich in absorption features, with two components contributing to the line profiles. In addition to the ground state lines, we detect C II, O I, Si II, Fe II and Ni II excited absorption features. No host galaxy emission lines, molecular absorption features nor diffuse interstellar bands are detected in the spectrum. The Hydrogen column density associated to GRB 090926A is log N_H/cm^{-2} = 21.60 +/- 0.07, and the metallicity of the host galaxy is in the range [X/H] =3.2X10^{-3}-1.2X10^{-2} with respect to the solar values, i.e., among the lowest values ever observed for a GRB host galaxy. A comparison with galactic chemical evolution models has suggested that the host of GRB090926A is likely to be a dwarf irregular galaxy. We put an upper limit to the Hydrogen molecular fraction of the host galaxy ISM, which is f < 7X10^{-7}. We derive information on the distance between the host absorbing gas and the site of the GRB explosion. The distance of component I is found to be 2.40 +/- 0.15 kpc, while component II is located far away from the GRB, possibly at ~ 5 kpc. These values are compatible with that found for other 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

Near-infrared follow-up to the May 2008 activation of SGR 1627-41

On 28 May 2008, the Swift satellite detected the first reactivation of SGR 1627-41 since its discovery in 1998. Following this event we began an observing campaign in near infrared wavelengths to search for a possible counterpart inside the error circle of this SGR, which is expected to show flaring activity simultaneous to the high energy flares or at least some variability as compared to the quiescent state. For the follow-up we used the 0.6m REM robotic telescope at La Silla Observatory, which allowed a fast response within 24 hours and, through director discretionary time, the 8.2m Very Large Telescope at Paranal Observatory. There, we observed with NACO to produce high angular resolution imaging with the aid of adaptive optics. These observations represent the fastest near infrared observations after an activation of this SGR and the deepest and highest spatial resolution observations of the Chandra error circle. 5 sources are detected in the immediate vicinity of the most precise X-ray localisation of this source. For 4 of them we do not detect variability, although the X-ray counterpart experimented a significant decay during our observation period. The 5th source is only detected in one epoch, where we have the best image quality, so no variability constrains can be imposed and remains as the only plausible counterpart. We can impose a limit of Ks > 21.6 magnitudes to any other counterpart candidate one week after the onset of the activity. Our adaptive optics imaging, with a resolution of 0.2" provides a reference frame for subsequent studies of future periods of activity.

preprint2009arXiv

Non variability of intervening absorbers observed in the UVES spectra of the "naked-eye" GRB080319

The aim of this paper is to investigate the properties of the intervening absorbers lying along the line of sight of Gamma-Ray Burst (GRB) 080319B through the analysis of its optical absorption features. To this purpose, we analyze a multi-epoch, high resolution spectroscopic observations (R=40000, corresponding to 7.5 km/s) of the optical afterglow of GRB080319B (z=0.937), taken with UVES at the VLT. Thanks to the rapid response mode (RRM), we observed the afterglow just 8m:30s after the GRB onset when the magnitude was R ~ 12. This allowed us to obtain the best signal-to-noise, high resolution spectrum of a GRB afterglow ever (S/N per resolution element ~ 50). Two further RRM and target of opportunity observations were obtained starting 1.0 and 2.4 hours after the event, respectively. Four MgII absorption systems lying along the line of sight to the afterglow have been detected in the redshift range 0.5 < z < 0.8, most of them showing a complex structure featuring several components. Absorptions due to FeII, MgI and MnII are also present; they appear in four, two and one intervening absorbers, respectively. One out of four systems show a MgII2796 rest frame equivalent width larger than 1A. This confirms the excess of strong MgII absorbers compared to quasars, with dn/dz = 0.9, ~ 4 times larger than the one observed along quasar lines of sight. In addition, the analysis of multi-epoch, high-resolution spectra allowed us to exclude a significant variability in the column density of the single components of each absorber. Combining this result with estimates of the size of the emitting region, we can reject the hypothesis that the difference between GRB and QSO MgII absorbers is due to a different size of the emitting regions.

preprint2009arXiv

The optical counterparts of Accreting Millisecond X-Ray Pulsars during quiescence

Eight Accreting Millisecond X-ray Pulsars (AMXPs) are known to date. Optical and NIR observations carried out during quiescence give a unique opportunity to constrain the nature of the donor star and to investigate the origin of the observed quiescent luminosity at long wavelengths. Using data obtained with the ESO-Very Large Telescope, we performed a deep optical and NIR photometric study of the fields of XTE J1814-338 and of the ultracompact systems XTE J0929-314 and XTE J1807-294 during quiescence in order to look for the presence of a variable counterpart. If suitable candidates were found, we also carried out optical spectroscopy. We present here the first multi-band (VR) detection of the optical counterpart of XTE J1814-338 in quiescence together with its optical spectrum. The optical light curve shows variability in both bands consistent with a sinusoidal modulation at the known 4.3 hr orbital period and presents a puzzling decrease of the V-band flux around superior conjunction that may be interpreted as a partial eclipse. The marginal detection of the very faint counterpart of XTE J0929-314 and deep upper limits for the optical/NIR counterpart of XTE J1807-294 are also reported. We also briefly discuss the results reported in the literature for the optical/NIR counterpart of XTE J1751-305. Our findings are consistent with AMXPs being systems containing an old, weakly magnetized neutron star, reactivated as a millisecond radio pulsar during quiescence which irradiates the low-mass companion star. The absence of type I X-ray bursts and of hydrogen and helium lines in outburst spectra of ultracompact (P_orb < 1 hr) AMXPs suggests that the companion stars are likely evolved dwarf stars.

preprint2009arXiv

The X-ray absorbing column densities of Swift Gamma-ray bursts

Long gamma-ray bursts (GRBs) are associated with the explosion of massive stars in star forming regions. A large fraction of GRBs show intrinsic absorption as detected in optical spectra but absorption signatures are also detectable in afterglow X-ray spectra. We present here a comprehensive analysis the full sample of 93 GRBs with known redshift promptly observed by Swift XRT up to June 2009. The distribution of X-ray column densities clearly shows that GRBs are heavily absorbed indicating that they indeed occur in dense environments. Furthermore, there is a lack of heavily absorbed GRBs at low redshift (z<1-2) that might therefore be candidates for the missing `dark' GRB population. However, there is no statistically significant correlation between the amount of X-ray absorption and the `darkness' of a GRB. Finally, we compare the hydrogen column densities derived in the optical with those derived from X-ray absorption. The two distributions are different, with the optical column densities being lower than the X-ray ones, which is even more apparent when correcting for metallicity effects. The most likely explanation is photoionization of hydrogen in the circumburst material caused by the radiation field of the burst.

preprint2009arXiv

The X-ray afterglow of GRB 081109A: clue to the wind bubble structure

We present the prompt BAT and afterglow XRT data of Swift-discovered GRB081109A up to ~ 5\times 10^5 sec after the trigger, and the early ground-based optical follow-ups. The temporal and spectral indices of the X-ray afterglow emission change remarkably. We interpret this as the GRB jet first traversing the freely expanding supersonic stellar wind of the progenitor with density varying as $ρ\propto r^{-2}$. Then after approximately 300 sec the jet traverses into a region of apparent constant density similar to that expected in the stalled-wind region of a stellar wind bubble or the interstellar medium (ISM). The optical afterglow data are generally consistent with such a scenario. Our best numerical model has a wind density parameter {$A_{*} \sim 0.02$, a density of the stalled wind $n\sim 0.12 {\rm cm}^{-3}$, and a transition radius $ \sim 4.5 \times 10^{17}$ cm}. Such a transition radius is smaller than that predicted by numerical simulations of the stellar wind bubbles and may be due to a rapidly evolving wind of the progenitor close to the time of its core-collapse.

preprint2009arXiv

UVES/VLT high resolution absorption spectroscopy of the GRB080330 afterglow: a study of the GRB host galaxy and intervening absorbers

We study the Gamma Ray Burst (GRB) environment and intervening absorbers by analyzing the optical absorption features produced by gas surrounding the GRB or along its line of sight. We analyzed high resolution spectroscopic observations (R=40000, S/N=3 - 6) of the optical afterglow of GRB080330, taken with UVES at the VLT ~ 1.5 hours after the GRB trigger. The spectrum illustrates the complexity of the ISM of the GRB host galaxy at z = 1.51 which has at least four components in the main absorption system. We detect strong FeII, SiII, and NiII excited absorption lines associated with the bluemost component only. In addition to the host galaxy, at least two more absorbers lying along the line of sight to the afterglow have been detected in the redshift range 0.8 < z < 1.1, each exhibiting MgII absorption. For the bluemost component in the host galaxy, we derive information about its distance from the site of the GRB explosion. We do so by assuming that the excited absorption lines are produced by indirect UV pumping, and compare the data with a time dependent photo-excitation code. The distance of this component is found to be 280+40-50 pc, which is lower than found for other GRBs (1 - 6 kpc). We identify two additional MgII absorbers, one of them with a rest frame equivalent width larger than 1A. The distance between the GRB and the absorber measured in this paper confirms that the power of the GRB radiation can influence the conditions of the interstellar medium up to a distance of at least several hundred pc. For the intervening absorbers, we confirm the trend that on average one strong intervening system is found per afterglow, as has been noted in studies exhibiting an excess of strong MgII absorbers along GRB sightlines compared to quasars.

preprint2008arXiv

Observations of X-ray Flares from Gamma Ray Bursts

Swift-XRT observations of the X-ray emission from gamma ray bursts (GRBs) and during the GRB afterglow have led to many new results during the past two years. One of these exciting results is that approximately 1/3-1/2 of GRBs contain detectable X-ray flares. The mean fluence of the X-ray flares is ~10 times less than that of the initial prompt emission, but in some cases the flare is as energetic as the prompt emission itself. The flares display fast rises and decays, and they sometimes occur at very late times relative to the prompt emission (sometimes as late as 10^5 s after T_0) with very high peak fluxes relative to the underlying afterglow decay that has clearly begun prior to some flares. The temporal and spectral properties of the flares are found to favor models in which flares arise due to the same GRB internal engine processes that spawned the prompt GRB emission. Therefore, both long and short GRB internal engine models must be capable of producing high fluences in the X-ray band at very late times.

preprint2008arXiv

Probing the complex environments of GRB host galaxies and intervening systems: high resolution spectroscopy of GRB050922C

The aim of this paper is to investigate the environment of gamma ray bursts (GRBs) and the interstellar matter of their host galaxies. We use to this purpose high resolution spectroscopic observations of the afterglow of GRB050922C, obtained with UVES/VLT 3.5 hours after the GRB event. We found that, as for most high resolution spectra of GRBs, the spectrum of the afterglow of GRB050922C is complex. At least seven components contribute to the main absorption system at z=2.1992. The detection of lines of neutral elements like MgI and the detection of fine-structure levels of the ions FeII, SiII and CII allows us to separate components in the GRB ISM along the line of sight. Moreover, in addition to the main system, we have analyzed the five intervening systems between z = 2.077 and z = 1.5664 identified along the GRB line of sight. GRB afterglow spectra are very complex, but full of information. This can be used to disentangle the contribution of the different parts of the GRB host galaxy and to study their properties. Our metallicity estimates agree with the scenario of GRBs exploding in low metallicity galaxies

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

Attributes of flares in Gamma Ray Bursts: sample I

We discuss some of the preliminary results and findings derived from the analysis of a first sample of flares detected by the XRT on board Swift. The analysis shows that the morphology of flares is the one we expect from the collision of ultra-relativistic shells as it happens during the internal shock model proposed by Rees and Meszaros in 1994. Furthermore the Delta(t)/t ratio and the decay-time to rise-time ratio have mean values that are in good agreement with the values observed in the prompt emission pulses that are believed to originate from internal shocks. The conclusion is that the flare analysis favors the internal shock as due to shells that have been ejected by the central engine after the prompt emission. The central engine seems to remain active and capable of generating large amounts of energy also at later times. More data are needed to establish whether or not some of the flares could be due to shells that have been emitted with small Lorentz factor at the time of the prompt emission and generate later time flares due to the catch up of these shells.

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

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

Long-term monitoring of the X-ray afterglow of GRB 050408 with Swift/XRT

We present observations of the X-ray afterglow of GRB 050408, a gamma-ray burst discovered by HETE-II. Swift began observing the field 42 min after the burst, performing follow-up over a period of 38 d (thus spanning three decades in time).The X-ray light curve showed a steepening with time, similar to many other afterglows. However, the steepening was unusually smooth, over the duration of the XRT observation, with no clear break time. The early decay was too flat to be described in terms of standard models. We therefore explore alternative explanations, such as the presence of a structured afterglow or of long-lasting energy injection into the fireball from the central GRB engine. The lack of a sharp break puts constraints on these two models. In the former case, it may indicate that the angular energy profile of the jet was not a simple power law, while in the second model it implies that injection did not stop abruptly. The late decay may be due either to a standard afterglow (that is, with no energy injection), or to a jetted outflow still being refreshed. A significant amount of absorption was present in the X-ray spectrum, corresponding to a rest-frame Hydrogen column density NH = 1.2 (-0.3,+0.4)*10^22 cm^-2, indicative of a dense environment.

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.

preprint2006arXiv

X-ray flare in XRF 050406: evidence for prolonged engine activity

We present observations of XRF 050406, the first burst detected by Swift showing a flare in its X-ray light curve. During this flare, which peaks at t_peak ~210s after the BAT trigger, a flux variation of (delta F)/F~6 in a very short time (delta t)/t_peak<<1 was observed. Its measured fluence in the 0.2-10 keV band was ~1.4x10^-8 erg cm^-2, which corresponds to 1-15% of the prompt fluence. We present indications of spectral variations during the flare. We argue that the producing mechanism is late internal shocks, which implies that the central engine is still active at 210s, though with a reduced power with respect to the prompt emission. The X-ray light curve flattens to a very shallow slope with decay index of ~0.5 after ~4400s, which also supports continued central engine activity at late times. This burst is classified as an X-ray flash, with a relatively low fluence (~10^-7 erg cm^-2 in the 15-350 keV band, E_iso~10^51 erg), a soft spectrum (photon index 2.65), no significant flux above ~50 keV and a peak energy E_p<15 keV. XRF 050406 is one of the first examples of a well-studied X-ray light curve of an XRF. We show that the main afterglow characteristics are qualitatively similar to those of normal GRBs. In particular, X-ray flares superimposed on a power-law light curve have now been seen in both XRFs and GRBs. This indicates that a similar mechanism may be at work for both kinds of events.

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

Open issues in gamma-ray bursts: polarimetry and dark GRBs

We review some open problems in the physics of afterglows, namely their polarization properties and the existence of dark/faint bursts. Polarization studies yield precious insights in the physical structure and dynamical evolution of GRB jets, revealing their magnetization properties and their energy profile. Polarimetric observations of GRB 020813 already allowed to exclude a homogeneous jet for this event. We then present observations of faint/dark bursts, showing that some of them may be obscured by dust, while others are possibly just intrinsically dim.

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.

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.

preprint1995arXiv

SPECTROPHOTOMETRIC INDICES AND METAL CONTENT OF GALACTIC GLOBULAR CLUSTERS

Spectrophotometric indices for $18$ Galactic globular clusters, obtained from CCD observations and careful reductions, were used to determine reliable calibrations on metallicity $[Fe/H]$. The indices were measured in the bandpasses adopted by Burnstein et al. (\cite{BFGK84}). Adding other observations of Burnstein et al. (\cite{BFGK84}) we obtained our results from an homogeneous sample of indices for $26$ globular clusters. Relations with indices defined by other Authors and with metallicity photometric indices or parameters were also computed. In each case the relations are quite satisfactory. Observational data were compared with synthetic indices derived from Buzzoni's (\cite{B89}) models and detailed discussions were performed for $Mg_2$, $Fe_{52}$, and $H_β$. The observational points seem to be systematically shifted with respect to the fiducial lines traced by the models. The scenario confirms that a certain degree of oxygen enhancement would be necessary to obtain a better agreement between observed data and theoretical predictions. This enhancement, however, removes some of the disagreement, but not all of it. The dependence of the observed $Fe_{52}$ and $H_β$ indices on the metal content for different HB morphologies was considered. Finally, some results were also discussed from a statistical point of view. A principal component analysis was applied to the index sample to study the number of independent parameters necessary to reproduce the observations. The whole index set is completely consistent with a one-parameter family.