Source author record

F. Tavecchio

F. Tavecchio 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

185works
5topics
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

185 published item(s)

preprint2025arXiv

How do recollimation-induced instabilities shape the propagation of hydrodynamic relativistic jets?

Recollimation is a phenomenon of particular importance in the dynamic evolution of jets and in the emission of high-energy radiation. Additionally, the full comprehension of this phenomenon provides insights into fundamental properties of jets in the vicinity of the Active Galactic Nucleus (AGN). Three-dimensional (magneto-)hydrodynamic simulations revealed that the jet conditions at recollimation favor the growth of strong instabilities, challenging the traditional view-supported from two-dimensional simulations-of confined jets undergoing a series of recollimation and reflection shocks. To investigate the stability of relativistic jets in AGNs at recollimation sites, we perform a set of long duration three-dimensional relativistic hydrodynamic simulations with the state-of-the-art PLUTO code, to focus on the development of hydrodynamical instabilities. We explore the non-linear growth of the instabilities and their effects on the physical jet properties as a function of the initial jet parameters: jet Lorentz factor, temperature, opening angle and jet-environment density-contrast. The parameter space is designed to describe low-power, weakly magnetized jets at small distances from the core (around the parsec scale). All collimating jets we simulated develop instabilities. Recollimation instabilities decelerate the jet, heat it, entrain external material, and move the recollimation point to shorter distances from the core. This is true for both conical and cylindrical jets. The instabilities, that are first triggered by the centrifugal instability, appear to be less disruptive in the case of narrower, denser, more relativistic, and warmer jets. These results provide valuable insights into the complex processes governing AGN jets and could be used to model the properties of low-power, weakly magnetized jetted AGNs.

preprint2022arXiv

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

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

preprint2022arXiv

ASTRI Mini-Array Core Science at the Observatorio del Teide

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

preprint2022arXiv

Direct observation of an extended X-ray jet at $z$=6.1

We report on the direct observation of an extended X-ray jet in the $z$=6.1 radio-loud Active Galactic Nucleus PSO J030947.49+271757.31 from a deep Chandra X-ray observation (128 ksec). This detection represents the most distant kpc off-nuclear emission resolved in the X-rays to date. The angular distance of the emission is $\sim$4" (corresponding to $\sim$20 kpc at $z$=6.1), along the same direction of the jet observed at parsec scales in previous VLBA high-resolution radio observations. Moreover, the 0.5-7.0 keV isophotes coincide with the extended radio emission as imaged by the VLA Sky Survey at 3 GHz. The rest-frame 2-10 keV luminosity of the extended component is L$_{2-10keV}$=5.9$\times$10$^{44}$ erg s$^{-1}$, about 8% of the core: this makes it one of the most luminous jets resolved in the X-rays so far. Through Spectral Energy Distribution modelling we find that this emission can be explained by the Inverse Compton interaction with the photons of the Cosmic Microwave Background assuming that the jet's physical parameters are similar to those in the local Universe. At the same time, we find that the radiation produced by a putative population of high-energetic electrons through the synchrotron process observed at low redshift is quenched at high redshift, hence becoming negligible.

preprint2022arXiv

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

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

preprint2022arXiv

Extreme blazars: the result of unstable recollimated jets?

Extreme BL Lacs (EHBL) form a subclass of blazars which challenge standard emission scenarios. In a recent study it has been argued that their peculiar properties can be explained if emitting electrons are accelerated in a series of oblique shocks induced by the recollimation of the relativistic jet. However, new 3D simulations of recollimated, weakly magnetized jets reveal that, in correspondence with the first recollimation shock, the flow develops a rapidly growing instability, becomes highly turbulent and decelerates, effectively hampering the formation of the multiple shock structure routinely observed in 2D simulations. Building on these new findings, we propose here a revised scenario for EHBL, in which the emission is produced by electrons accelerated at the recollimation shock and subsequently further energized through stochastic acceleration in the turbulent downstream flow. We apply a simple version of this scenario to the prototypical EHBL 1ES 0229-200, showing that the SED can be satisfactorily reproduced with standard values of the main physical parameters.

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

Testing particle acceleration models for BL LAC jets with the Imaging X-ray Polarimetry Explorer

To determine if the recently launched Imaging X-ray Polarimetry Explorer (IXPE) can follow the polarization variations induced by different particle acceleration mechanisms in blazars jets we simulate observations of a high synchrotron peak (HSP) blazar variable in polarization degree and angle according to the theoretical predictions of Tavecchio et al., 2020 (shock acceleration model), Bodo et al., 2021 (magnetic reconnection model) and Marscher \& Jorstad, 2021 (Turbulent Extreme Multi-Zone Model). We use the Monte Carlo tool ixpeobssim to create realistic IXPE data products. We create simulated light-curves of polarization degree and angle by time-slicing the simulated data into arbitrary short time bins. We use a $χ^2$ test to assess the performance of the observations in detecting the time variability of the polarization properties. A series of 10 ks long observations permits IXPE to follow the time variability of the polarization degree in the case of the shock acceleration model. In the case of the magnetic reconnection model, the nominal injected model provides the best fit of the simulated IXPE data for time bins of 5-10 ks, depending on the tested flux level. For the TEMZ model, shorter time slices of 0.5 ks are needed for obtaining a formally good fit of the simulated IXPE data with the injected model. On the other hand, we find that a fit with a constant model provides a $χ^2$ lower than the fit with the nominal injected model when using time slices of 20 ks, 60/70 ks, and 5 ks for the case of the shock acceleration, magnetic reconnection and TEMZ model, respectively. We conclude that, provided that the statistics of the observation allows slicing of the data in adequately short time bins, IXPE observations of an HSP blazar at a typical flux level can detect the time variability predicted by popular models for particle acceleration in jets.

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

Conceptual Design of BabyIAXO, the intermediate stage towards the International Axion Observatory

This article describes BabyIAXO, an intermediate experimental stage of the International Axion Observatory (IAXO), proposed to be sited at DESY. IAXO is a large-scale axion helioscope that will look for axions and axion-like particles (ALPs), produced in the Sun, with unprecedented sensitivity. BabyIAXO is conceived to test all IAXO subsystems (magnet, optics and detectors) at a relevant scale for the final system and thus serve as prototype for IAXO, but at the same time as a fully-fledged helioscope with relevant physics reach itself, and with potential for discovery. The BabyIAXO magnet will feature two 10 m long, 70 cm diameter bores, and will host two detection lines (optics and detector) of dimensions similar to the final ones foreseen for IAXO. BabyIAXO will detect or reject solar axions or ALPs with axion-photon couplings down to $g_{aγ} \sim 1.5 \times 10^{-11}$ GeV$^{-1}$, and masses up to $m_a\sim 0.25$ eV. BabyIAXO will offer additional opportunities for axion research in view of IAXO, like the development of precision x-ray detectors to identify particular spectral features in the solar axion spectrum, and the implementation of radiofrequency-cavity-based axion dark matter setups.

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

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

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

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

Constraining the shear acceleration model for the X-ray emission of large-scale extragalactic jets

The nature of the intense X-ray emission from powerful extragalactic jets at large ($>1$ kpc) scale is still debated. The scenario that invokes the inverse Compton scattering of the CMB by electrons is challenged by the lack of gamma-ray emission in the GeV band. An alternative assumes synchrotron emission by a distinct population of ultra-high energy electrons. Here we present a concrete attempt to apply this scenario, exploring the specific model in which the ultra-high energy electrons are accelerated in a shear layer surrounding the jet. We limit the study to non-relativistic flows and particle acceleration is treated by a Fokker-Planck equation. The observed relation between low energy (radio, optical) and X-ray emission prompts us to assume that the required population of pre-accelerated particles is provided by a shock responsible for the acceleration of the electrons emitting at low frequencies. We apply the model to the emission of the principal knots of the jets of PKS 0637-752 and PKS 1136-135, two of the best studied objects. For the set of fiducial parameters adopted, the condition that the jet power does not exceeds a limiting value of $10^{48}$ erg s$^{-1}$ constrains the magnetic field above $10$ $μ$G and indicates moderate beaming ($δ\simeq 2$) for PKS 0637-752. For both sources, the requirement that acceleration of the electrons proceeds faster than radiative cooling can be met if the magnetic turbulence in the shear layer follows a Kolmogorov spectrum, $I(k)\propto k^{-q}$ with $q=5/3$, but cannot satisfied in the Bohm-like case ($q=1$).

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

On the distribution of fluxes of gamma-ray blazars: hints for a stochastic process?

We examine a model for the observed temporal variability of powerful blazars in the $γ$-ray band in which the dynamics is described in terms of a stochastic differential equation, including the contribution of a deterministic drift and a stochastic term. The form of the equation is motivated by the current astrophysical framework, accepting that jets are powered through the extraction of the rotational energy of the central supermassive black hole mediated by magnetic fields supported by a so-called \emph{magnetically arrested} accretion disk. We apply the model to the $γ$-ray light curves of several bright blazars and we infer the parameters suitable to describe them. In particular, we examine the differential distribution of fluxes ($dN/dF_γ$) and we show that the predicted probability density function for the assumed stochastic equation naturally reproduces the observed power law shape at large fluxes $dN/dF_γ \propto F_γ^{-α}$ with $α>2$.

preprint2020arXiv

Probing shock acceleration in BL Lac jets through X-ray polarimetry: the time-dependent view

Polarimetric measurements, especially if extended at high energy, are expected to provide important insights into the mechanisms underlying the acceleration of relativistic particles in jets. In a previous work we have shown that the polarization of the synchrotron X-ray emission produced by highly energetic electrons accelerated by a mildly relativistic shock carries essential imprints of the geometry and the structure of the magnetic fields in the downstream region. Here we present the extension of our analysis to the non-stationary case, especially suitable to model the highly variable emission of high-energy emitting BL Lacs. We anticipate a large ($Π\approx 40\%$), almost time-independent degree of polarization in the hard/medium X-ray band, a prediction soon testable with the upcoming mission {\it IXPE}. The situation in other bands, in particular in the optical, is more complex. A monotonic decrease of the optical degree of polarization is observed during the development of a flare. At later stages $Π$ reaches zero and then it starts to increase, recovering large values at late times. The instant at which $Π=0$ is marked by a rotation of the polarization angle by $90$ degrees. However, at optical frequencies it is likely that more than one region contributes to the observed emission, potentially making it difficult to detect the predicted behavior.

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

Anisotropic electron populations in BL Lac jets: consequences for the observed emission

We investigate the impact on the properties of high-energy emitting BL Lac objects of a population of electrons with an anisotropic momentum distribution. We adopt a simple phenomenological description of the anisotropy, in which the most energetic electrons have a small pitch angle and the least energetic electrons are isotropic, as proposed by Sobacchi \& Lyubarsky (2019). We explore (i) a simple model that assumes a phenomenological shape for the electron energy distribution, and (ii) a self-consistent scheme in which the electrons follow a distribution which is the result of the balance between injection and radiative losses (we include the effects of the anisotropy on the synchrotron cooling rate). Considering the BL Lac object Mkn 421 as representative of the entire class, we show that in both cases the emission can be satisfactorily reproduced under equipartition between the magnetic field and the relativistic electrons. This is in better agreement with the idea that jets are launched as Poynting dominated flows with respect to the standard isotropic scenario, which requires both a low magnetization and a low radiative efficiency to reproduce the observed emission. The hard spectrum predicted for the inverse Compton continuum at TeV energies could be used as a potential test of the anisotropic model.

preprint2019arXiv

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

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

preprint2018arXiv

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

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

preprint2018arXiv

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

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

preprint2018arXiv

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

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

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

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

preprint2017arXiv

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

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

preprint2017arXiv

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

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

preprint2017arXiv

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

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

preprint2017arXiv

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

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

preprint2017arXiv

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

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

preprint2017arXiv

Performance of the MAGIC telescopes under moonlight

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

preprint2017arXiv

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

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

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

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

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2015arXiv

Blazar candidates beyond redshift 4 observed by Swift

We have selected SDSS J222032.50+002537.5 and SDSS J142048.01+120545.9 as best blazar candidates out of a complete sample of extremely radio-loud quasars at z>4, with highly massive black holes. We observed them and a third serendipitous candidate with similar features (PMN J2134-0419) in the X-rays with the Swift/XRT telescope, to confirm their blazar nature. We observed strong and hard X-ray fluxes (i.e. $α_X<0.6$, where $F(ν)\proptoν^{-α_X}$ in the 0.3-10keV observed energy range, ~1-40keV rest frame) in all three cases. This allowed us to classify our candidates as real blazars, being characterized by large Lorentz factors (~13) and very small viewing angles (~3deg). All three sources have black hole masses exceeding 10^9Msun and their classification provides intriguing constraints on supermassive black hole formation and evolution models. We confirm our earlier suggestion that there are different formation epochs of extremely massive black holes hosted in jetted (z~4) and non-jetted systems (z~2.5).

preprint2015arXiv

CMB quenching of high-redshift radio-loud AGNs

The very existence of more than a dozen of high-redshift (z>4) blazars indicates that a much larger population of misaligned powerful jetted AGN was already in place when the Universe was <1.5 Gyr old. Such parent population proved to be very elusive, and escaped direct detection in radio surveys so far. High redshift blazars themselves seem to be failing in producing extended radio-lobes, raising questions about the connection between such class and the vaster population of radio-galaxies. We show that the interaction of the jet electrons with the intense cosmic microwave background (CMB) radiation explains the lack of extended radio emission in high redshift blazars and in their parent population, helping to explain the apparently missing misaligned counterparts of high redshift blazars. On the other hand, the emission from the more compact and more magnetised hot spots are less affected by the enhanced CMB energy density. By modelling the spectral energy distribution of blazar lobes and hot spots we find that most of them should be detectable by low frequency deep radio observations, e.g., by LOw-Frequency ARray for radio astronomy (LOFAR) and by relatively deep X-ray observations with good angular resolution, e.g., by the Chandra satellite. At high redshifts, the emission of a misaligned relativistic jet, being de-beamed, is missed by current large sky area surveys. The isotropic flux produced in the hot spots can be below ~1 mJy and the isotropic lobe radio emission is quenched by the CMB cooling. Consequently, even sources with very powerful jets can go undetected in current radio surveys, and misclassified as radio-quiet AGNs.

preprint2015arXiv

Extreme BL Lacs: probes for cosmology and UHECR candidates

High-energy observations of extreme BL Lac objects, such as 1ES0229+200 or 1ES 0347-121, recently focused interest both for blazar and jet physics and for the implication on the extragalactic background light and intergalactic magnetic field estimate. Moreover, their enigmatic properties have been interpreted in a scenario in which their primary high- energy output is through a beam of high-energy hadrons. However, despite their possible important role in all these topics, the number of these extreme highly peaked BL Lac objects (EHBL) is still rather small. Aiming at increase their number, we selected a group of EHBL candidates considering those undetected (or only barely detected) by the LAT onboard Fermi and characterized by a high X-ray versus radio flux ratio. We assembled the multi-wavelength spectral energy distribution of the resulting 9 sources, using available archival data of Swift, GALEX, and Fermi satellites, confirming their nature. Through a simple one-zone synchrotron self-Compton model we estimate the expected very high energy flux, finding that in the majority of cases it is within the reach of present generation of Cherenkov arrays or of the forthcoming CTA.

preprint2015arXiv

Fermi/LAT broad emission line blazars

We study the broad emission line blazars detected in the gamma-ray band by the Large Area Telescope onboard the Fermi satellite and with the optical spectrum studied by Shaw et al. (2012, 2013). The observed broad line strength provides a measure of the ionizing luminosity of the accretion disk, while the gamma-luminosity is a proxy for the bolometric non-thermal beamed jet emission. The resulting sample, composed by 217 blazars, is the best suited to study the connection between accretion and jet properties. We compare the broad emission line properties of these blazars with those of radio-quiet and radio-loud quasars present in the Sloan Digital Sky Survey, to asses differences and similarities of the disc luminosity and the virial black hole mass. For most sources, we could derive the black hole mass by reproducing the IR-optical-UV data with a standard accretion disc spectrum, and we compared the black hole masses derived with the two methods. The distributions of the masses estimated in the two ways agree satisfactorily. We then apply a simple, one-zone, leptonic model to all the 217 objects of our sample. The knowledge of the black hole mass and disc luminosity helps to constrain the jet parameters. On average they are similar to what found by previous studies of smaller samples of sources.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

No axions from the Sun

Preliminary evidence of solar axions in XMM-Newton observations has quite recently been claimed by Fraser et al. as an interpretation of their detection of a seasonally-modulated excess of the X-ray background. Within such an interpretation, these authors also estimate the axion mass to be $m_a \simeq 2.3 \cdot 10^{- 6}$ eV. Since an axion with this mass behaves as a cold dark matter particle, according to the proposed interpretation the considered detection directly concerns cold dark matter as well. So, the suggested interpretation would lead to a revolutionary discovery if confirmed. Unfortunately, we have identified three distinct problems in this interpretation of the observed result of Fraser et al. which ultimately imply that the detected signal - while extremely interesting in itself - cannot have any relation with hypothetical axions produced by the Sun. Thus, a physically consistent interpretation of the observed seasonally-modulated X-ray excess still remains an exciting challenge.

preprint2015arXiv

On the detectability of Lorentz invariance violation through anomalous multi-TeV $γ$-ray spectra of blazars

Cosmic opacity for very high-energy gamma rays ($E>10$ TeV) due to the interaction with the extragalactic background light can be strongly reduced because of possible Lorentz-violating terms in the particle dispersion relations expected, e.g., in several versions of quantum gravity theories. We discuss the possibility to use very high energy observations of blazars to detect anomalies of the cosmic opacity induced by LIV, considering in particular the possibility to use -- besides the bright and close-by BL Lac Mkn 501 -- extreme BL Lac objects. We derive the modified expression for the optical depth of $γ$ rays considering also the redshift dependence and we apply it to derive the expected high-energy spectrum above 10 TeV of Mkn 501 in high and low state and the extreme BL Lac 1ES 0229+200. We find that, besides the nearby and well studied BL Lac Mkn 501 -- especially in high state --, suitable targets are extreme BL Lac objects, characterized by quite hard TeV intrinsic spectra likely extending at the energies relevant to detect LIV features.

preprint2015arXiv

On the magnetization of BL Lac jets

The current paradigm foresees that relativistic jets are launched as magnetically dominated flows, whose magnetic power is progressively converted to kinetic power of of the matter of the jet, until equipartition is reached. Therefore, at the end of the acceleration phase, the jet should still carry a substantial fraction ($\approx$ half) of its power in the form of a Poynting flux. It has been also argued that, in these conditions, the best candidate particle acceleration mechanism is efficient reconnection of magnetic field lines, for which it is predicted that magnetic field and accelerated relativistic electron energy densities are in equipartition.Through the modeling of the jet non--thermal emission, we explore if equipartition is indeed possible in BL Lac objects, i.e. low-power blazars with weak or absent broad emission lines. We find that one-zone models (for which only one region is involved in the production of the radiation we observe) the particle energy density is largely dominating (by 1-2 orders of magnitude) over the magnetic one. As a consequence, the jet kinetic power largely exceeds the magnetic power. Instead, if the jet is structured (i.e. made by a fast spine surrounded by a slower layer), the amplification of the IC emission due to the radiative interplay between the two components allows us to reproduce the emission in equipartition conditions.

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

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

Unveiling the nature of the gamma-ray emitting AGN PKS 0521-36

PKS 0521-36 is an Active Galactic Nucleus (AGN) with uncertain classification. We investigate the properties of this source from radio to gamma rays. The broad emission lines in the optical and UV bands and steep radio spectrum indicate a possible classification as an intermediate object between broad-line radio galaxies (BLRG) and steep spectrum radio quasars (SSRQ). On pc-scales PKS 0521-36 shows a knotty structure similar to misaligned AGN. The core dominance and the gamma-ray properties are similar to those estimated for other SSRQ and BLRG detected in gamma rays, suggesting an intermediate viewing angle with respect to the observer. In this context the flaring activity detected from this source by Fermi-LAT between 2010 June and 2012 February is very intriguing. We discuss the gamma-ray emission of this source in the framework of the structured jet scenario, comparing the spectral energy distribution (SED) of the flaring state in 2010 June with that of a low state. We present three alternative models corresponding to three different choices of the viewing angles theta_v = 6, 15, and 20 deg. We obtain a good fit for the the first two cases, but the SED obtained with theta_v = 15 deg if observed at a small angle does not resemble that of a typical blazar since the synchrotron emission should dominate by a large factor (about 100) the inverse Compton component. This suggests that a viewing angle between 6 and 15 deg is preferred, with the rapid variability observed during gamma-ray flares favouring a smaller angle. However, we cannot rule out that PKS 0521-36 is the misaligned counterpart of a synchrotron-dominated blazar.

preprint2015arXiv

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

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

preprint2014arXiv

An active state of the BL Lac Object Markarian 421 detected by INTEGRAL in April 2013

Multiwavelength variability of blazars offers indirect insight into their powerful engines and on the mechanisms through which energy is propagated from the centre down the jet. The BL Lac object Mkn 421 is a TeV emitter, a bright blazar at all wavelengths, and therefore an excellent target for variability studies. Mkn 421 was observed by INTEGRAL and Fermi-LAT in an active state on 16-21 April 2013. Well sampled optical, soft, and hard X-ray light curves show the presence of two flares. The average flux in the 20-100 keV range is 9.1e-11 erg/s/cm2 (~4.5 mCrab) and the nuclear average apparent magnitude, corrected for Galactic extinction, is V ~12.2. In the time-resolved X-ray spectra (3.5-60 keV), which are described by broken power laws and, marginally better, by log-parabolic laws, we see a hardening that correlates with flux increase, as expected in refreshed energy injections in a population of electrons that later cool via synchrotron radiation. The hardness ratios between the JEM-X fluxes in two different bands and between the JEM-X and IBIS/ISGRI fluxes confirm this trend. During the observation, the variability level increases monotonically from the optical to the hard X-rays, while the large LAT errors do not allow a significant assessment of the MeV-GeV variability. The cross-correlation analysis during the onset of the most prominent flare suggests a monotonically increasing delay of the lower frequency emission with respect to that at higher frequency, with a maximum time-lag of about 70 minutes, that is however not well constrained. The spectral energy distributions from the optical to the TeV domain are satisfactorily described by homogeneous models of blazar emission based on synchrotron radiation and synchrotron self-Compton scattering, except in the state corresponding to the LAT softest spectrum and highest flux.

preprint2014arXiv

Black hole lightning due to particle acceleration at subhorizon scales

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Exploring the blazar zone in High Energy flares of FSRQs

The gamma-ray emission offers a powerful diagnostic tool to probe jets and their surroundings in flat spectrum radio quasars (FSRQ). In particular, sources emitting at high energies (>10 GeV) give us the strongest constraints. This motivates us to start a systematic study of flares with bright emission above 10 GeV, examining archival data of Fermi-LAT gamma-ray telescope. At the same time, we began to trigger Target of Opportunity observations to the Swift observatory at the occurrence of high-energy flares, obtaining a wide coverage of the spectral energy distributions for several FSRQs during flares. Among the others we investigate the SED of a peculiar flare of 3C 454.3, showing a remarkable hard gamma-ray spectrum, quite different from the brightest flares of this source, and a bright flare of CTA 102. We modeled the SED in the framework of the one--zone leptonic model, using also archival optical spectroscopic data to derive the luminosity of the broad lines and thus estimate the disk luminosity, from which the structural parameters of the FSRQ nucleus can be inferred. The model allowed us to evaluate the magnetic field intensity in the blazar zone, and to locate the emitting region of gamma rays in the particular case in which gamma-ray spectra show neither absorption from the BLR, nor the Klein-Nishina curvature expected in leptonic models assuming the BLR as source of seed photons for the External Compton. For FSRQs bright above 10 GeV, we where able to identify short periods lasting less than 1 day characterized by high rate of high energy gamma rays, and hard gamma-ray spectra. We discussed the observed spectra and variability timescales in terms of injection and cooling of energetic particles, arguing that these flares could be triggered by magnetic reconnections events or turbulence in the flow.

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

High-energy cosmic neutrinos from spine-sheath BL Lac jets

We recently proposed that structured (spine-sheath) jets associated to BL Lac objects offer a suitable environment for the production of the extragalactic high-energy ($E>100$ TeV) neutrino recently revealed by IceCube. Our previous analysis was limited to low-power BL Lac objects. We extend our preliminary study to the entire BL Lac population. We assume that the power of cosmic rays as well as the radiative luminosity of the sheath depend linearly on the the jet power. In turn, we assume that the latter is well traced by the $γ$-ray luminosity. We exploit the BL Lac $γ$-ray luminosity function and its cosmic evolution as recently inferred from Fermi-LAT data to derive the expected neutrino cumulative intensity from the entire BL Lac population. When considering only the low-power BL Lacs, a large cosmic ray power for each source is required to account for the neutrino flux. Instead, if BL Lacs of all powers produce neutrinos, the power demand decreases, and the required cosmic ray power becomes of the same order of the radiative jet power. We also discuss the prospects for the direct association of IceCube events with BL Lacs, providing an estimate of the expected counts for the most promising sources.

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

On the spine-layer scenario for the very high-energy emission of NGC 1275

We discuss the $γ$-ray emission of the radiogalaxy NGC 1275 (the central galaxy of the Perseus Cluster), detected by Fermi-LAT and MAGIC, in the framework of the "spine-layer" scenario, in which the jet is assumed to be characterized by a velocity structure, with a fast spine surrounded by a slower layer. The existence of such a structure in the parsec scale jet of NGC 1275 has been recently proved through VLBI observations. We discuss the constraints that the observed spectral energy distribution imposes to the parameters and we present three alternative models, corresponding to three different choices of the angles between the jet and the line of sight ($θ_{\rm v}=6^{\circ}, 18^{\circ}$ and 25$^{\circ}$). While for the the case with $θ_{\rm v}=6^{\circ}$ we obtain an excellent fit, we consider this solution unlikely, since such small angles seems to be excluded by radio observations of the large-scale jet. For $θ_{\rm v}=25^{\circ}$ the required large intrinsic luminosity of the soft (IR--optical) component of the spine determines a large optical depth for $γ$-rays through the pair production scattering $γγ\rightarrow e^+ e^-$, implying a narrow cut--off at $\sim50$~GeV. We conclude that intermediate angles are required. In this case the low frequency and the high--energy emissions are produced by two separate regions and, in principle, a full variety of correlations is expected. The correlation observed between the optical and the $γ$-ray flux, close to linearity, is likely linked to variations of the emissivity of the spine.

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

Photon-ALP conversions inside AGN

An intriguing possibility to partially circumvent extragalactic background light (EBL) absorption in very-high-energy (VHE) observations of blazars is that photons convert into axion-like particles (ALPs) $γ\to a$ inside or close to a blazar and reconvert into photons $a \to γ$ in the Milky Way magnetic field. This idea has been put forward in 2008 and has attracted a considerable interest. However, while the probability for the back-conversion $a \to γ$ has been computed in detail (using the maps of the Galatic magnetic field), regretfully no realistic estimate of the probability for the conversion $γ\to a$ inside a blazar has been performed, in spite of the fact that the present-day knowledge allows this task to be accomplished in a reliable fashion. We present a detailed calculation that fills this gap, considering both types of blazars, namely BL Lac objects (BL Lacs) and flat spectrum radio quasars (FSRQ) with their specific structural and environmental properties. We also include the host elliptical galaxy into account. Our somewhat surprising results show that the conversion probability in BL Lacs is strongly dependent on the source parameters -- like the position of the emission region along the jet and the strength of the magnetic field therein -- making it effectivelly unpredictable. On the other hand, the lobes at the termination of FSRQ jets lead to an effective "equipartition" between photons and ALPs due to its chaotic nature, thereby allowing us to make a clear-cut prediction. These results are quite important in view of the planned VHE detectors like the CTA, HAWK and HiSCORE.

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

SDSS J114657.79+403708.6: the third most distant blazar at z=5.0

The radio-loud quasar SDSS J114657.79+403708.6 at a redshift z=5.0 is one of the most distant radio-loud objects. The IR-optical luminosity and spectrum suggest that its black hole has a very large mass: M=(5+-1)x 1e9 Msun. The radio-loudness (ratio of the radio to optical flux) of the source is large (around 100), suggesting that the source is viewed at small angles from the jet axis, and could be a blazar. The X-ray observations fully confirm this hypothesis, due to the high level and hardness of the flux. This makes SDSS J114657.79+403708.6 the third most distant blazar known, after Q0906+693 (z=5.47) and B2 1023+25 (z=5.3). Among those, SDSS J114657.79+403708.6 has the largest black hole mass, setting interesting constraints on the mass function of heavy (larger than one billion solar masses) black holes at high redshifts.

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

Structured jets in BL Lac objects: efficient PeV neutrino factories?

The origin of the high-energy neutrinos (0.1$-$1 PeV range) detected by IceCube remains a mystery. In this work we explore the possibility that efficient neutrino production can occur in structured jets of BL Lac objects, characterized by a fast inner spine surrounded by a slower layer. This scenario has been widely discussed in the framework of the high-energy emission models for BL Lacs and radiogalaxies. One of the relevant consequences of a velocity structure is the enhancement of the inverse Compton emission caused by the radiative coupling of the two zones. We show that a similar boosting could occur for the neutrino output of the spine through the photo-meson reaction of high-energy protons scattering off the amplified soft target photon field of the layer. Assuming the local density and the cosmological evolution of $γ$-ray BL Lac derived from {\it Fermi}-LAT data, we calculate the expected diffuse neutrino intensity, that can match the IceCube data for a reasonable choice of the parameters.

preprint2014arXiv

The power of relativistic jets is larger than the luminosity of their accretion disks

Theoretical models for the production of relativistic jets from active galactic nuclei predict that jet power arises from the spin and mass of the central black hole, as well as the magnetic field near the event horizon. The physical mechanism mechanism underlying the contribution from the magnetic field is the torque exerted on the rotating black hole by the field amplified by the accreting material. If the squared magnetic field is proportional to the accretion rate, then there will be a correlation between jet power and accretion luminosity. There is evidence for such a correlation, but inadequate knowledge of the accretion luminosity of the limited and inhomogeneous used samples prevented a firm conclusion. Here we report an analysis of archival observations of a sample of blazars (quasars whose jets point towards Earth) that overcomes previous limitations. We find a clear correlation between jet power as measured through the gamma-ray luminosity, and accretion luminosity as measured by the broad emission lines, with the jet power dominating over the disk luminosity, in agreement with numerical simulations. This implies that the magnetic field threading the black hole horizon reaches the maximum value sustainable by the accreting matter.

preprint2013arXiv

Basic properties of Narrow-Line Seyfert 1 Galaxies with Relativistic Jets

We present the preliminary results of a survey performed with Swift to observe a sample of radio-loud Narrow-Line Seyfert 1 Galaxies (RLNLS1s). Optical-to-X-ray data from Swift are complemented with gamma-ray observations from Fermi/LAT and radio measurements available in the literature. The comparison with a sample of bright Fermi blazars indicates that RLNLS1s seem to be the low-power tail of the distribution.

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

Fermi LAT detection of extraordinary variability in the gamma-ray emission of the blazar PKS 1510-089

We have reanalyzed the giant outburst of the blazar PKS 1510-089 (z=0.36) that occurred on 2011 October-November. The gamma-ray flux in the 0.1-100 GeV energy range exceeded the value of 10^-5 ph/cm2/s for several days. The peak flux was reached on 2011 October 19, with a value of 4.4 x 10^-5 ph/cm2/s, which in turn corresponds to a luminosity of 2 x 10^49 erg/s. A very short timescale variability was measured. Particularly on 2011 October 18, the flux-doubling time was as short as 20 minutes. This is the shortest variability ever detected in the MeV-GeV energy band. We compared our analysis with two other outbursts observed in 2009 March and 2012 February-March, when the blazar was also detected by HESS and MAGIC to infer information about the emission at hundreds of GeV.

preprint2013arXiv

Radio-loud AGNs at high redshifts and the cosmic microwave background

We discuss how the interaction between the electrons in a relativistic jet and the Cosmic Microwave Background (CMB) affects the observable properties of radio-loud AGN at early epochs. At high z the magnetic energy density in the radio lobes of powerful radio-loud quasars can be exceeded by the energy density of the CMB (because of its (1+z)^4 dependance). In this case, relativistic electrons cool preferentially by scattering off CMB photons, rather than by synchrotron. Thus, sources sharing the same intrinsic properties have different extended radio and X-ray luminosities when located at different z: more distant sources are less luminous in radio and more luminous in X-rays than their closer counterparts. Instead, in compact regions where the local magnetic field still exceeds the CMB in terms of energy density, synchrotron radiation would be unaffected by the presence of the CMB. Such regions include the compact inner jet and the so-called hot spots in the radio lobes. The decrease in radio luminosity is larger in misaligned sources, whose radio flux is dominated by the extended isotropic component. These sources can fail detection in current flux limited radio surveys, and therefore they are possibly under-represented in the associated samples. As the cooling time is longer for lower energy electrons, the radio luminosity deficit due to the CMB photons is less important at low radio frequencies. Therefore objects not detected so far in current surveys at a few GHz could be picked up by low frequency deep surveys, such as LOFAR and SKA. Until then, we can estimate the number of high redshift radio-loud AGNs through the census of their aligned proxies, i.e., blazars. Indeed, their observed radio emission arises in the inner and strongly magnetized compact core of the relativistic jet, and not affected by inverse Compton scattering off CMB photons.

preprint2013arXiv

The ASTRI Mini-Array Science Case

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

preprint2013arXiv

The far emission region of the gamma-ray blazar PKS B1424-418

We report multi-wavelength (IR-O, UV, X-ray, gamma-ray) data of the flat spectrum radio quasar (FSRQ) PKS B1424-418 (z=1.52) gathered during an active phase in 2013 April . As for a few other cases of FSRQs reported in literature, the detection by the Large Area Telescope onboard Fermi of gamma rays with energy above 10 GeV indicates that the emission likely occurs beyond the highly opaque (tau ~10) broad line region. This conclusion is strengthened by a model fit to the spectral energy distribution, which allows us to locate the emission region even beyond the distance generally assumed for the dusty torus. The consequent large size (~1 pc) inferred for the emission region cannot account for the observed daily-scale variability of the gamma-ray flux. We discuss the possibility that short-term variability results from fast magnetic reconnection events, as proposed in recent works.

preprint2013arXiv

The red blazar PMN J2345-1555 becomes blue

The Flat Spectrum Radio Quasar PMN J2345-1555 is a bright gamma-ray source, that recently underwent a flaring episode in the IR, UV and gamma-ray bands. The flux changed quasi simultaneously at different frequencies, suggesting that it was produced by a single population of emitting particles, hence by a single and well localized region of the jet. While the overall Spectral Energy Distribution (SED) before the flare was typical of powerful blazars (namely two broad humps peaking in the far IR and below 100 MeV bands, respectively), during the flare the peaks moved to the optical-UV and to energies larger than 1 GeV, to resemble low power BL Lac objects, even if the observed bolometric luminosity increased by more than one order of magnitude. We interpret this behavior as due to a change of the location of the emission region in the jet, from within the broad line region, to just outside. The corresponding decrease of the radiation energy density as seen in the comoving frame of the jet allowed the relativistic electrons to be accelerated to higher energies, and thus produce a "bluer" SED.

preprint2012arXiv

"Flat" broad line region and gamma-ray absorption in blazars

We study the impact of the geometry of the broad line region (BLR) on the expected absorption, through the γγ\rightarrow e^{\pm} process, of gamma rays produced in the relativistic jet of flat spectrum radio quasars (FSRQ). We consider "flat" (or "disky") BLR models, and use BLR spectra calculated with the photoionization code CLOUDY, already used to investigate the emission and the absorption of high-energy photons in FSRQ. We characterize the energy-dependent optical depth of the process, τ(E), for different accretion disk luminosities, aperture angles of the BLR (α, as measured from the equatorial plane), and initial injection eighths of the high-energy photons, $R_{\rm o}$. We study in particular how the change of these parameters influences the spectral break at GeV energies, predicted if the emission occurs within the BLR. We found a well defined relation between the break energy and the post-break slope, both uniquely determined by α. We finally find that even a rather disk-like BLR (α\sim 25^\circ) corresponds to important absorption (τ>1) of photons above few tens of GeV produced within the BLR. We therefore conclude that the VHE emission detected from FSRQs occurs beyond the BLR.

preprint2012arXiv

Blue Fermi Flat Spectrum Radio Quasars

Many blazars detected by the Fermi satellite, observed spectroscopically in the optical, are line-less, and have been classified as BL Lac objects. Optical-UV photometry of nearly one hundred of them allowed to determine the redshift for a handful of objects and redshift upper limits for the great majority. A few of these are candidates to be "blue quasars", namely flat spectrum radio quasars whose broad emission lines are hidden by an overwhelming synchrotron emission peaking in the UV. This implies that the emitting electrons have high energies. In turn, this requires relatively weak radiative cooling, a condition that can be met if the main radiative dissipation of the jet power occurs outside the broad line region. We confirm this hypothesis by studying and modelling the spectral energy distributions of the 4 "blue quasars" recently discovered. Furthermore, we discuss the distribution of Fermi blazars in the gamma-ray spectral index -- gamma-ray luminosity plane, and argue that "blue quasars" objects are a minority within the blazar populations.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Evidence for an axion-like particle from PKS 1222+216?

The surprising discovery by MAGIC of an intense, rapidly varying emission in the energy range 70 - 400 GeV from the flat spectrum radio quasar PKS 1222+216 represents a challenge for all interpretative scenarios. Indeed, in order to avoid absorption of γrays in the dense ultraviolet radiation field of the broad line region (BLR), one is forced to invoke some unconventional astrophysical picture, like for instance the existence of a very compact (r\sim 10^{14} cm) emitting blob at a large distance (R \sim10^{18} cm) from the jet base. We offer the investigation of a scenario based on the standard blazar model for PKS 1222+216 where γrays are produced close to the central engine, but we add the new assumption that inside the source photons can oscillate into axion-like particles (ALPs), which are a generic prediction of several extensions of the Standard Model of elementary particle interactions. As a result, a considerable fraction of very-high-energy photons can escape absorption from the BLR through the mechanism of photon-ALP oscillations much in the same way as they largely avoid absorption from extragalactic background light when propagating over cosmic distances in the presence of large-scale magnetic fields in the nG range. In addition we show that the above MAGIC observations and the simultaneous Fermi/LAT observations in the energy range 0.3 - 3 GeV can both be explained by a standard spectral energy distribution for experimentally allowed values of the model parameters. In particular, we need a very light ALP just like in the case of photon-ALP oscillations in cosmic space. Moreover, we find it quite tantalizing that the most favorable value of the photon-ALP coupling happens to be the same in both situations. Although our ALPs cannot contribute to the cold dark matter, they are a viable candidate for the quintessential dark energy. [abridged]

preprint2012arXiv

High redshift Fermi blazars observed by GROND and Swift

We observed 5 gamma-ray loud blazars at redshift greater than 2 with the X-Ray Telescope (XRT) and the UltraViolet and Optical Telescope (UVOT) onboard the Swift satellite, and the Gamma-Ray burst Optical Near-Infrared Detector (GROND) instrument. These observations were quasi simultaneous, usually within a few hours. For 4 of these blazars the near-IR to UV data show the presence of an accretion disc, and we could reliably estimate its accretion rate and black hole mass. One of them, PKS 1348+007, was found in an extraordinarily high IR-optical state, almost two orders of magnitude brighter than at the epoch of the Sloan Digital Sky Survey observations. For all the 5 quasars the physical parameters of the jet emitting zone, derived by applying a one-zone emission model, are similar to that found for the bulk of other gamma-ray loud quasars. With our observations we have X-ray data for the full sample of blazars at z>2 present in the Fermi 2-yrs (2LAC) catalog. This allows to have a rather complete view of the spectral energy distribution of all high-redshift Fermi blazars, and to draw some conclusions about their properties, and especially about the relation between the accretion rate and the jet power.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

PG 1553+113: five years of observations with MAGIC

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

preprint2012arXiv

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

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

preprint2012arXiv

Radio-to-gamma-ray Monitoring of the Narrow-Line Seyfert 1 Galaxy PMN J0948+0022 from 2008 to 2011

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

preprint2012arXiv

SDSS J102623.61+254259.5: the second most distant blazar at z=5.3

The radio-loud quasar SDSS J102623.61+254259.5, at a redshift z=5.3, is one of the most distant radio-loud objects. Since its radio flux exceeds 100 mJy at a few GHz, it is also one of the most powerful radio-loud sources. We propose that this source is a blazar, i.e. we are seeing its jet at a small viewing angle. This claim is based on the spectral energy distribution of this source, and especially on its strong and hard X-ray spectrum, as seen by Swift, very typical of powerful blazars. Observations by the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) and by theWide-field Infrared Survey Explorer (WISE) allow to establish the thermal nature of the emission in the near IR-optical band. Assuming that this is produced by a standard accretion disk, we derive that it emits a luminosity of L_d \simeq 9 \times 10^46 erg s^{-1} and that the black hole has a mass between 2 and 5 billion solar masses. This poses interesting constraints on the mass function of heavy (> 10^9 M_sun) black holes at high redshifts.

preprint2012arXiv

Short time scale variability at gamma rays in FSRQs and implications on the current models

We studied the rapid variability at GeV gamma rays of the flat-spectrum radio quasar PKS 1222+216, which was recently found by the MAGIC Cerenkov telescope to display very short variability (minutes time scale) at hundreds of GeV. We analyzed the time period between 2010 April 29 and June 20, when the source generated a few gamma-ray flares with flux in the MeV-GeV band in excess of 10^-5 ph cm^-2 s^-1 on daily basis. We set tight upper limits on the observed doubling time scale (about 1 hour on 2010 April 30), the smallest measured to date at MeV-GeV energies, which can constrain the size of the gamma-ray emitting region. We also studied the spectra measured during two flares (2010 April 30 and June 17-18). The combination of spectral and variability studies obtained in the present work favors the hypothesis that gamma rays are generally produced inside the broad-line region (BLR), but sometimes the dissipation can occur at larger distances, nearby the infrared torus.

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.

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

Constraining blazars distances with combined GeV and TeV data

Recently, a new method to constrain the distances of blazars with unknown redshift using combined observations in the GeV and TeV regimes has been developed. The underlying assumption is that the Very High Energy (VHE, E > 100 GeV) spectrum corrected for the absorption of TeV photons by the Extragalactic Background Light (EBL) via photon-photon interaction should still be softer than the gamma-ray spectrum observed by Fermi/LAT. The constraints found are related to the real redshifts by a simple linear relation, that has been used to infer the unknown or uncertain distance of blazars. The sample is revised with the up-to-date spectra in both TeV and GeV bands and the method applied to the unknown distance blazar PKS 1424+240 detected at VHE.

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

Extreme TeV blazars and the intergalactic magnetic field

We study the four BL Lac objects (RGB J0152+017, 1ES 0229+200, 1ES 0347-121 and PKS 0548-322) detected in the TeV band but not present in the 1FGL catalogue of the Fermi/Large Area Telescope. We analize the 24 months of LAT data deriving gamma-ray fluxes or upper limits that we use to assemble their spectral energy distributions (SED). We model the SEDs with a standard one-zone leptonic model, also including the contribution of the reprocessed radiation in the multi GeV band, emitted by the pairs produced through the conversion of the primary TeV emission by interaction with the cosmic optical-IR background. For simplicity, in the calculation of this component we adopt an analytical approach including some simplifying assumptions, in particular i) the blazar high energy emission is considered on average stable over times of the order of 10^7 years and ii) the observer is exactly on-axis. We compare the physical parameters derived by the emission model with those of other high-energy emitting BL Lacs, confirming that TeV BL Lacs with a rather small GeV flux are characterized by extremely low values of the magnetic field and large values of the electron energies. The comparison between the flux in the GeV band and that expected from the reprocessed TeV emission allows us to confirm and strengthen the lower limit of B >10^{-15} G for the intergalactic magnetic field using a theoretically motivated spectrum for the primary high-energy photons.

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

Gamma-ray variability of radio-loud narrow-line Seyfert 1 galaxies

The recent detection of gamma-ray emission from four radio-loud narrow-line Seyfert 1 galaxies suggests that the engine driving the AGN activity of these objects share some similarities with that of blazars, namely the presence of a gamma-ray emitting, variable, jet of plasma closely aligned to the line of sight. In this work we analyze the gamma-ray light curves of the four radio-loud narrow-line Seyfert 1 galaxies for which high-energy gamma-ray emission has been discovered by Fermi/LAT, in order to study their variability. We find significant flux variability in all the sources. This allows us to exclude a starburst origin of the gamma-ray photons and confirms the presence of a relativistic jet. Furthermore we estimate the minimum e-folding variability timescale (3 - 30 days) and infer an upper limit for the size of the emitting region (0.2 - 2 pc, assuming a relativistic Doppler factor delta=10 and a jet aperture of theta=0.1 rad).

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

On the origin of the gamma-ray emission from the flaring blazar PKS 1222+216

The flat spectrum radio quasar PKS 1222+216 was detected in the very high energy gamma-ray band by MAGIC during a highly active gamma-ray phase following an alert by the LAT onboard Fermi. Its relatively hard spectrum without a cut off, together with the observed variability on timescale of ~10 min challenges standard emission models. If the emission originates in a portion of the relativistic jet located inside the BLR, severe absorption of gamma rays above few tens of GeV is expected due to the pair production process. These observations imply the existence of a very compact (R_b ~5 x 10^{14} cm) and very fast blob located far beyond the BLR radius, responsible for the rapidly varying high energy flux. However the long term (days-weeks) coherent evolution of the GeV flux recorded by LAT indicates that there could be also the substantial contribution from another larger emission region. We model the spectral energy distribution of PKS 1222+216 during the epoch of the MAGIC detection assuming three different scenarios,: (1) a one-zone model considering only the emission from a compact blob outside the BLR; (2) a two-zone model considering the compact blob plus an emitting region filling the whole jet cross-section located outside the BLR and (3) a two zone model with the jet emitting region inside the BLR. We find that in all cases the high-energy emission from the compact blob is dominated by the inverse Compton scattering of the IR thermal radiation of the dusty torus. Both regions are matter dominated, with a negligible Poynting flux. These results disfavor models in which the compact blob is the result of reconnection events inside the jet or "needles" of high-energy electrons. Instead, the observational framework and our radiative models could be compatible with scenarios in which the jet is re-collimated and focussed at large distances from the central BH.

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

Search for the shortest variability at gamma rays in flat-spectrum radio quasars

We report about the search for short-term variability in the high-energy gamma-ray energy band of three flat-spectrum radio quasars (3C 454.3, 3C 273, PKS B1222+216), whose flux at E > 100 MeV exceeded the value of 10^-5 ph cm^-2 s^-1 for at least one day. Although, the statistics was not yet sufficient to effectively measure the characteristic time scale, it allowed us to set tight upper limits on the observed doubling time scale (< 2-3 hours) -- the smallest measured to date at MeV energies --, which can constrain the size of the gamma-ray emitting region. The results obtained in the present work favor the hypothesis that gamma rays are generated inside the broad-line region.

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

Study of the variability of Blazars gamma-ray emission

The gamma-ray emission of blazar jets shows a pronounced variability and this feature provides limits to the size and to the speed of the emitting region. We study the gamma-ray variability of bright blazars using data from the first 18 months of activity of the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. From the daily light-curves of the blazars characterized by a remarkable activity, we firstly determine the minimum variability time-scale, giving an upper limit for the size of the emitting region of the sources, assumed to be spheroidal blobs in relativistic motion. These regions must be smaller than ~10^-3 parsec. Another interesting time-scale is the duration of the outbursts. We conclude that they cannot correspond to radiation produced by a single blob moving relativistically along the jet, but they are either the signature of emission from a standing shock extracting energy from a modulated jet, or the superposition of a number of flares occurring on a shorter time-scale. We also derive lower limits on the bulk Lorentz factor needed to make the emitting region transparent for gamma-rays interacting through photon-photon collisions.

preprint2011arXiv

The first gamma-ray outburst of a Narrow-Line Seyfert 1 Galaxy: the case of PMN J0948+0022 in July 2010

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.5846) performed in 2010 July-September and triggered by a high-energy gamma-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. The peak flux in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^-6 ph cm^-2 s^-1, corresponding to an observed luminosity of 10^48 erg s^-1. Although the source was too close to the Sun position to organize a densely sampled follow-up, it was possible to gather some multiwavelength data that confirmed the state of high activity across the sampled electromagnetic spectrum. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

preprint2011arXiv

The July 2010 outburst of the NLS1 PMN J0948+0022

We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

preprint2011arXiv

The transition between BL Lac objects and Flat Spectrum Radio Quasars

We study the BL Lac objects detected in the one year all sky survey of the Fermi satellite, with a energy spectral slope alpha_gamma in the [0.1-100 GeV] band greater than 1.2. In the alpha_gamma vs gamma-ray luminosity plane, these BL Lacs occupy the region populated by Flat Spectrum Radio Quasars (FSRQs). Studying the properties of their spectral energy distributions (SED) and of their emitting lines, we find that several of these BL Lacs have a SED similar to FSRQs and that they do have broad lines of large equivalent width, and should be reclassified as FSRQs even adopting the current phenomenological definition (i.e. equivalent width EW of the emitting line greater than 5 A). In other cases, even if the EW width is small, the emitting lines can be as luminous as in quasars, and again their SED is similar to the SED of FSRQs. Sources classified as BL Lacs with a SED appearing as intermediate between BL Lacs and FSRQs also have relatively weak broad emission lines and small EW, and can be considered as transition sources. These properties are confirmed also by model fitting, that allows to derive the relevant intrinsic jet parameters and the jet power. This study leads us to propose a physical distinction between the two classes of blazars, based on the luminosity of the broad line region measured in Eddington units. The dividing line is of the order of L_BLR/L_Edd ~ 5e-4, in good agreement with the idea that the presence of strong emitting lines is related to a transition in the accretion regime, becoming radiatively inefficient below a disk luminosity of the order of one per cent of the Eddington one.

preprint2010arXiv

Chasing the heaviest black holes of jetted Active Galactic Nuclei

We investigate the physical properties of the 10 blazars at redshift greater than 2 detected in the 3-years all sky survey performed by the Burst Alert Telescope (BAT) onboard the Swift satellite. We find that the jets of these blazars are among the most powerful known. Furthermore, the mass of their central black hole, inferred from the optical-UV bump, exceeds a few billions of solar masses, with accretion luminosities being a large fraction of the Eddington one. We compare their properties with those of the brightest blazars of the 3-months survey performed by the Large Area Telescope (LAT) onboard the Fermi satellite. We find that the BAT blazars have more powerful jets, more luminous accretion disks and larger black hole masses than LAT blazars. These findings can be simply understood on the basis of the blazar sequence, that suggests that the most powerful blazars have a spectral energy distribution with a high energy peak at MeV (or even sub-MeV) energies. This implies that the most extreme blazars can be found more efficiently in hard X-rays, rather than in the high energy gamma-ray band. We then discuss the implications of our findings for future missions, such as the New Hard X-ray Mission (NHXM) and especially the Energetic X-ray Imaging Survey Telescope (EXIST) mission which, during its planned 2 years all sky survey, is expected to detect thousands of blazars, with a few of them at z greater than 6.

preprint2010arXiv

Compton rockets and the minimum power of relativistic jets

The power of a relativistic jet depends on the number of leptons and protons carried by the jet itself. We have reasons to believe that powerful gamma-ray flat spectrum radio sources emit most of their radiation where radiative cooling is severe. This helps to find the minimum number of emitting leptons needed to explain the radiation we see. The number of protons is more uncertain. If there is one proton per electron, they dominate the jet power, but they could be unimportant if the emission is due to electron-positron pairs. In this case the total jet power could be much smaller. However, if the gamma-ray flux is due to inverse Compton scattering with seed photons produced outside the jet, the radiation is anisotropic also in the comoving frame, making the jet to recoil. This Compton rocket effect is strong for light, electron-positron jets, and negligible for heavy, proton dominated jets. No significant deceleration, required by fast superluminal motion, requires a minimum number of protons per lepton, and thus a minimum jet power. We apply these ideas to the blazar 3C 454.3, to find a robust lower limit to its total jet power: if the viewing angle theta_v ~ 1/Gamma the jet power is larger than the accretion luminosity L_d for any bulk Lorentz factor Gamma. For theta_v =0, instead, the minimum jet power can be smaller than L_d for Gamma<25. No more than ~10 pairs per proton are allowed.

preprint2010arXiv

Constraining the location of the emitting region in Fermi blazars through rapid gamma-ray variability

We consider the 1.5 years Fermi Large Area Telescope light curves (E > 100 MeV) of the flat spectrum radio quasars 3C 454.3 and PKS 1510-089, which show high activity in this period of time. We characterise the duty cycle of the source by comparing the time spent by the sources at different flux levels. We consider in detail the light curves covering periods of extreme flux. The large number of high-energy photons collected by LAT in these events allows us to find evidence of variability on timescales of few hours. We discuss the implications of significant variability on such short timescales, that challenge the scenario recently advanced in which the bulk of the gamma-ray luminosity is produced in regions of the jet at large distances (tens of parsec) from the black hole.

preprint2010arXiv

Correlation of Fermi-LAT sources with the AT20GHz radio survey

We cross correlate the Fermi 11 months survey catalogue (1FGL) with the 20 GHz Australia Telescope Compact Array radio survey catalogue (AT20G) composed by 5890 sources at declination <0 deg. Among the 738 Fermi sources distributed in the southern sky we find 230 highly probable candidate counterparts in the AT20G survey. Of these, 222 are already classified as blazars (176 of known type and 46 of unknown optical class) in the Fermi 1-year LAT AGN Catalogue (1LAC) and 8 are new associations. By studying the gamma-ray and radio properties of these associations we find a strong correlation between the gamma-ray flux (above 100 MeV) and the 20 GHz flux density. This correlation is more than 3 sigma statistically significant both for the population of BL Lacs and of FSRQ considered separately. We also find that the radio counterparts associated to the Fermi sources have on average flat radio spectra between 5 and 20 GHz and that Fermi gamma-ray sources are not preferentially associated with "ultra inverted spectrum" radio sources. For 2 of the 8 new associations we build the broad band spectral energy distribution combining Fermi, Swift and radio observations. One of these two sources is identified with the high redshift FSRQ Swift J1656.3-3302 (z=2.4) and we classify the other source as a candidate new FSRQ. We also study the brightest radio source of the 46 associations without an optical classification and classify it as a new BL Lac candidate "twin" of the prototypical BL Lac OJ 287 if its redshift is somewhat larger, z~0.4.

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

Does the gamma-ray flux of the blazar 3C 454.3 vary on sub-hour timescales?

In the early days of April 2010, the blazar 3C 454.3 (z=0.859) underwent a strong gamma-ray outburst, reaching fluxes (E > 100 MeV) in excess of 10^-5 ph cm^-2 s^-1. The Fermi Gamma ray Space Telescope performed a 200 ks long pointed observation starting from 5 April 2010 19:38 UTC. This allowed us to try probing the variability of the gamma-ray emission on timescales of hours or less. We found the variability on a few hours timescale. On sub-hour timescale we found no evidence of significant variability, although the present statistics is not yet conclusive and further observations are needed.

preprint2010arXiv

High redshift Fermi blazars

With the release of the first year Fermi catalogue, the number of blazars detected above 100 MeV lying at high redshift has been largely increased. There are 28 blazars at z>2 in the "clean" sample. All of them are Flat Spectrum Radio Quasars (FSRQs). We study and model their overall spectral energy distribution in order to find the physical parameters of the jet emitting region, and for all of them we estimate their black hole masses and accretion rates. We then compare the jet with the accretion disk properties, setting these sources in the broader context of all the other bright gamma-ray or hard X-ray blazars. We confirm that the jet power correlates with the accretion luminosity. We find that the high energy emission peak shifts to smaller frequencies as the observed luminosity increases, according to the blazar sequence, making the hard X-ray band the most suitable for searching the most luminous and distant blazars.

preprint2010arXiv

INTEGRAL observations of the GeV blazar PKS1502+106 and of the hard X-ray bright Seyfert Galaxy Mkn841

Extragalactic nuclear activity is best explored with observations at high energies, where the most extreme flux and spectral variations are expected to occur, witnessing changes in the accretion flow or in the kinematics of the plasma. In active galactic nuclei of blazar type, these variations are the most dramatic. By following blazar outbursts from their onset and by correlating the observed variations at many different wavelengths we can reconstruct the behavior of the plasma and map out the development of the flare within the jet. The advent of the Fermi satellite has allowed the start of a systematic and intensive monitoring program of blazars. Blazar outbursts are very effectively detected by the LAT instrument in the MeV-GeV domain, and these can be promptly followed up with other facilities. Based on a Fermi LAT detection of a high MeV-GeV state, we have observed the blazar PKS 1502+106 with the INTEGRAL satellite between 9 and 11 August 2008. Simultaneous Swift observations have been also accomplished, as well as optical follow-up at the Nordic Optical Telescope. The IBIS instrument onboard INTEGRAL detected a source at a position inconsistent with the optical coordinates of PKS 1502+106, but consistent with those of the Seyfert 1 galaxy Mkn 841, located at 6.8 arcmin south-west from the blazar, which is therefore responsible for all the hard X-ray flux detected by IBIS. At the location of the blazar, IBIS sets an upper limit of ~10^{-11} erg/s/cm2 on the 15-60 keV flux, that turns out to be consistent with a model of inverse Compton scattering accounting for the soft X-ray and gamma-ray spectra measured by Swift XRT and Fermi LAT, respectively. The gamma-ray spectrum during the outburst indicates substantial variability of the characteristic energy of the inverse Compton component in this blazar.(abridged).

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

Relativistic jets in Narrow-Line Seyfert 1

Narrow-Line Seyfert 1 (NLS1) class of active galactic nuclei (AGNs) is generally radio-quiet, but a small percent of them are radio-loud. The recent discovery by Fermi/LAT of high-energy gamma-ray emission from 4 NLS1s proved the existence of relativistic jets in these systems. It is therefore important to study this new class of gamma-ray emitting AGNs. Here we report preliminary results about the observations of the July 2010 gamma-ray outburst of PMN J0948+0022, when the source flux exceeded for the first time 10^-6 ph cm^-2 s^-1 (E > 100 MeV).

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

Suzaku Observations of Luminous Quasars: Revealing the Nature of High-Energy Blazar Emission in Quiescent States

We present the results from the Suzaku X-ray observations of five flat-spectrum radio quasars (FSRQs), namely PKS0208-512, Q0827+243, PKS1127-145, PKS1510-089 and 3C 454.3. All these sources were additionally monitored simultaneously or quasi-simultaneously by the Fermi satellite in gamma-rays and the Swift UVOT in the UV and optical bands, respectively. We constructed their broad-band spectra covering the frequency range from 10^14 Hz up to 10^25 Hz, and those reveal the nature of high-energy emission of luminous blazars in their low-activity states. The analyzed X-ray spectra are well fitted by a power-law model with photoelectric absorption. In the case of PKS0208-512, PKS1127-145, and 3C 454.3, the X-ray continuum showed indication of hard-ening at low-energies. Moreover, when compared with the previous X-ray observations, we see a significantly increasing contribution of low-energy photons to the total X-ray fluxes when the sources are getting fainter. The same behavior can be noted in the Suzaku data alone. A likely explanation involves a variable, flat-spectrum component produced via inverse-Compton (IC) emission, plus an additional, possibly steady soft X-ray component prominent when the source gets fainter. This soft X-ray excess is represented either by a steep powerlaw (photon indices Gamma ~ 3 - 5) or a blackbody-type emission with temperatures kT ~ 0.1-0.2 keV. We model the broad-band spectra spectra of the five observed FSRQs using synchrotron self-Compton (SSC) and/or external-Compton radiation (ECR) models. Our modeling suggests that the difference between the low and high-activity states in luminous blazars is due to the different total kinetic power of the jet, most likely related to varying bulk Lorentz factor of the outflow within the blazar emission zone.

preprint2010arXiv

The gamma-ray brightest days of the blazar 3C 454.3

In the first week of December 2009, the blazar 3C 454.3 became the brightest high energy source in the sky. Its photon flux reached and surpassed the level of 1e-5 ph/cm2/s above 100 MeV. The Swift satellite observed the source several times during the period of high gamma-ray flux, and we can construct really simultaneous spectral energy distributions (SED) before, during, and after the luminosity peak. Our main findings are: i) the optical, X-ray and gamma-ray fluxes correlate; ii) the gamma-ray flux varies quadratically (or even more) with the optical flux; iii) a simple one-zone synchrotron inverse Compton model can account for all the considered SED; iv) in this framework the gamma-ray vs optical flux correlation can be explained if the magnetic field is slightly fainter when the overall jet luminosity is stronger; v) the power that the jet spent to produce the peak gamma-ray luminosity is of the same order, or larger, than the accretion disk luminosity. During the flare, the total jet power surely surpassed the accretion power.

preprint2010arXiv

The intergalactic magnetic field constrained by Fermi/LAT observations of the TeV blazar 1ES 0229+200

TeV photons from blazars at relatively large distances, interacting with the optical-IR cosmic background, are efficiently converted into electron-positron pairs. The produced pairs are extremely relativistic (Lorentz factors of the order of 1e6 1e7 and promptly loose their energy through inverse Compton scatterings with the photons of the microwave cosmic background, producing emission in the GeV band. The spectrum and the flux level of this reprocessed emission is critically dependent on the intensity of the intergalactic magnetic field, B, that can deflect the pairs diluting the intrinsic emission over a large solid angle. We derive a simple relation for the reprocessed spectrum expected from a steady source. We apply this treatment to the blazar 1ES 0229+200, whose intrinsic very hard TeV spectrum is expected to be approximately steady. Comparing the predicted reprocessed emission with the upper limits measured by the Fermi/Large Area Telescope, we constrain the value of the intergalactic magnetic field to be larger than $B \simeq 5\times 10^{-15}$ Gauss, depending on the model of extragalactic background light.

preprint2009arXiv

Canonical high power blazars

The jets of powerful blazars propagate within regions relatively dense of radiation produced externally to the jet. This radiation is a key ingredient to understand the origin of the high energy emission of blazars, from the X-ray to the gamma-ray energy band. These external radiation fields control the amount of the inverse Compton radiation with respect to the synchrotron flux. Therefore the predicted spectral energy distribution (SED) will depend on where the jet dissipates part of its energy to produce the observed radiation. We investigate in detail how the SED changes as a function of the location of the jet dissipation region, by assuming rather "standard" (i.e. "canonical") prescriptions for the accretion disk and its X-ray corona, the profile of the jet magnetic field and the external radiation. The magnetic energy density of a "canonical" jet almost never dominates the radiative cooling of the emitting electrons, and consequently the inverse Compton flux almost always dominates the bolometric output. This is more so for large black hole masses. Dissipation taking place beyond the broad line region is particularly interesting, since it accounts in a simple way for the largest inverse Compton to synchrotron flux ratios accompanied by an extremely hard X-ray spectrum. Furthermore it makes the high power blazars at high redshift useful tools to study the optical to UV cosmic backgrounds.

preprint2009arXiv

Fermi/LAT discovery of gamma-ray emission from a relativistic jet in the narrow-line quasar PMN J0948+0022

We report the discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope of high-energy gamma-ray emission from the peculiar quasar PMN J0948+0022 (z=0.5846). The optical spectrum of this object exhibits rather narrow Hbeta (FWHM(Hbeta) ~ 1500 km s^-1), weak forbidden lines and is therefore classified as a narrow-line type I quasar. This class of objects is thought to have relatively small black hole mass and to accrete at high Eddington ratio. The radio loudness and variability of the compact radio core indicates the presence of a relativistic jet. Quasi simultaneous radio-optical-X-ray and gamma-ray observations are presented. Both radio and gamma-ray emission (observed over 5-months) are strongly variable. The simultaneous optical and X-ray data from Swift show a blue continuum attributed to the accretion disk and a hard X-ray spectrum attributed to the jet. The resulting broad band spectral energy distribution (SED) and, in particular, the gamma-ray spectrum measured by Fermi are similar to those of more powerful FSRQ. A comparison of the radio and gamma-ray characteristics of PMN J0948+0022 with the other blazars detected by LAT shows that this source has a relatively low radio and gamma-ray power, with respect to other FSRQ. The physical parameters obtained from modelling the SED also fall at the low power end of the FSRQ parameter region discussed in Celotti & Ghisellini (2008). We suggest that the similarity of the SED of PMN J0948+0022 to that of more massive and more powerful quasars can be understood in a scenario in which the SED properties depend on the Eddington ratio rather than on the absolute power.

preprint2009arXiv

General physical properties of bright Fermi blazars

We studied all blazars of known redshift detected by the Fermi satellite during its first three months survey. For the majority of them, pointed Swift observations ensures a good multiwavelength coverage, enabling us to to reliably construct their spectral energy distributions (SED). We model the SEDs using a one-zone leptonic model and study the distributions of the derived interesting physical parameters as a function of the observed gamma-ray luminosity. We confirm previous findings concerning the relation of the physical parameters with source luminosity which are at the origin of the blazar sequence. The SEDs allow to estimate the luminosity of the accretion disk for the majority of broad emitting line blazars, while for the line-less BL Lac objects in the sample upper limits can be derived. We find a positive correlation between the jet power and the luminosity of the accretion disk in broad line blazars. In these objects we argue that the jet must be proton-dominated, and that the total jet power is of the same order of (or slightly larger than) the disk luminosity. We discuss two alternative scenarios to explain this result.

preprint2009arXiv

Jet and accretion power in the most powerful Fermi blazars

Among the blazars detected by the Fermi satellite, we have selected the 23 blazars that in the three months of survey had an average gamma-ray luminosity above 1e48 erg/s. For 17 out of the 23 sources we found and analysed X-ray and optical-UV data taken by the Swift satellite. With these data, implemented by archival and not simultaneous data, we construct the spectral energy distributions, and interpreted them with a simple one-zone, leptonic, synchrotron and inverse Compton model. When possible, we also compare different high energy states of single sources, like 0528+134 and 3C 454.3, for which multiple good sets of multi-wavelength data are available. In our powerful blazars the high energy emission always dominates the electromagnetic output, and the relatively low level of the synchrotron radiation often does not hide the accretion disk emission. We can then constrain the black hole mass and the disk luminosity. Both are large (i.e. masses equal or greater than 1e9 solar masses and disk luminosities above 0.1 Eddington). By modelling the non-thermal continuum we derive the power that the jet carries in the form of bulk motion of particles and fields. On average, the jet power is found to be slightly larger than the disk luminosity, and proportional to the mass accretion rate.

preprint2009arXiv

Radio-Loud Narrow-Line Seyfert 1 as a New Class of Gamma-Ray AGN

We report the discovery with Fermi/LAT of gamma-ray emission from three radio-loud narrow-line Seyfert 1 galaxies: PKS 1502+036 (z=0.409), 1H 0323+342 (z=0.061) and PKS 2004-447 (z=0.24). In addition to PMN J0948+0022 (z=0.585), the first source of this type to be detected in gamma rays, they may form an emerging new class of gamma-ray active galactic nuclei (AGN). These findings can have strong implications on our knowledge about relativistic jets and the unified model of AGN.

preprint2009arXiv

TeV BL Lac objects at the dawn of the Fermi era

We reconsider the emission properties of the BL Lac objects emitting in the high-energy gamma-ray band exploiting the new information in the MeV-GeV band obtained by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope in its first three months of operation. To this aim we construct the spectral energy distribution of all the BL Lacs revealed by LAT and of the known TeV BL Lacs not detected by LAT, also including data from the Swift satellite, and model them with a simple one-zone leptonic model. The analysis shows that the BL Lacs detected by LAT (being or not already detected in the TeV band) share similar physical parameters. While some of the TeV BL Lacs not revealed by LAT have spectral energy distributions and physical parameters very similar to the LAT BL Lacs, a group of objects displays peculiar properties (larger electron energies and smaller magnetic fields) suggesting different physical conditions in the emission region. Finally, we discuss possible criteria to effectively select good new candidates for the Cherenkov telescopes among the LAT sources, presenting a list of predicted fluxes in the very high-energy band calculated including the effect of the absorption by the extragalactic background light.

preprint2009arXiv

The blazar S5 0014+813: a real or apparent monster?

A strong hard X-ray luminosity from a blazar flags the presence of a very powerful jet. If the jet power is in turn related to the mass accretion rate, the most luminous hard X-ray blazars should pinpoint the largest accretion rates, and therefore the largest black hole masses. These ideas are confirmed by the Swift satellite observations of the blazar S5 0014+813, at the redshift z=3.366. Swift detected this source with all its three instruments, from the optical to the hard X-rays. Through the construction of its spectral energy distribution we are confident that its optical-UV emission is thermal in origin. Associating it to the emission of a standard optically thick geometrically thin accretion disk, we find a black hole mass of 40 billion solar masses, radiating at 40% the Eddington value. The derived mass is among the largest ever found. Super-Eddington slim disks or thick disks with the presence of a collimating funnel can in principle reduce the black hole mass estimate, but tends to produce spectra bluer than observed.

preprint2009arXiv

The Fermi blazars divide

Flat Spectrum Radio Quasars (FSRQs) and BL Lac objects detected in the first three months of the Fermi survey neatly separate in the gamma-ray spectral index vs gamma-ray luminosity plane. BL Lac objects are less luminous and have harder spectra than broad line blazars. We suggest that this division has its origin in the different accretion regimes of the two classes of objects. Using the gamma-ray luminosity as a proxy for the observed bolometric one we show that the boundary between the two subclasses of blazars can be associated with the threshold between the regimes of optically thick accretion disks and of radiatively inefficient accretion flows, which lies at an accretion rate of the order of 0.01 the Eddington rate. The spectral separation in hard (BL Lacs) and soft (FSRQs) objects can then result from the different radiative cooling suffered by the relativistic electrons in jets propagating in different ambients. We argue that the bulk of the most luminous blazars alread detected by Fermi should be characterised by large black hole masses, around 10^9 solar masses, and predict that lowering the gamma-ray flux threshold the region of the alpha_gamma-L_gamma plane corresponding to steep spectral indices and lower luminosities will be progressively populated by FSRQs with lower mass black holes, while the region of hard spectra and large luminosities will remain forbidden.

preprint2009arXiv

The hard TeV spectrum of 1ES 0229+200: new clues from Swift

The BL Lac object 1ES 0229+200 (z=0.14) has been detected by HESS during observations taking place in 2005-2006. The TeV spectrum, when corrected for the absorption of gamma-ray photons through the interaction with the extragalactic background light, is extremely hard, even if the most conservative level for the background is considered. The case of 1ES 0229+200 is very similar to that of 1ES 1101-232, for which a possible explanation, in the framework of the standard one-zone synchrotron-self Compton model, is that the high-energy emission is synchrotron-self Compton radiation of electrons distributed as a power law with a large value of the minimum energy. In this scenario the hard TeV spectrum is accompanied by a very hard synchrotron continuum below the soft X-ray band. We will show that recent Swift observations of 1ES 0229+200 in the critical UV-X-ray band strongly support this model, showing the presence of the expected spectral break and hard continuum between the UV and the X-ray bands.

preprint2009arXiv

The jet of the BL Lac object PKS 0521 -365 in the near-IR : MAD adaptive optics observations

BL Lac objects are low--power active nuclei exhibiting a variety of peculiar properties that are caused by the presence of a relativistic jet and orientation effects. We present here adaptive optics near-IR images at high spatial resolution of the nearby BL Lac object PKS 0521-365, which is known to display a prominent jet both at radio and optical frequencies. The observations were obtained in Ks--band using the ESO multi-conjugated adaptive optics demonstrator at the Very Large Telescope. This allowed us to obtain images with 0.1 arcsec effective resolution. We performed a detailed analysis of the jet and its related features from the near-IR images, and combined them with images previously obtained with HST in the R band and by a re-analysis of VLA radio maps. We find a remarkable similarity in the structure of the jet at radio, near-IR, and optical wavelengths. The broad--band emission of the jet knots is dominated by synchrotron radiation, while the nucleus also exhibits a significant inverse Compton component. We discovered the near-IR counterpart of the radio hotspot and found that the near-IR flux is consistent with being a synchrotron emission from radio to X-ray. The bright red object (red-tip), detached but well aligned with the jet, is well resolved in the near-IR and has a linear light profile. Since it has no radio counterpart, we propose that it is a background galaxy not associated with the jet. The new adaptive optics near-IR images and previous observations at other frequencies allow us to study the complex environment around the remarkable BL Lac object PKS 0521-365. These data exemplify the capabilities of multi conjugate adaptive optics observations of extragalactic extended sources.

preprint2007arXiv

An Infrared Study of the Large-scale Jet in Quasar PKS 1136-135

We present Spitzer IRAC imaging of the large-scale jet in the quasar PKS 1136-135 at wavelengths of 3.6 and 5.8 um, combined with previous VLA, HST, and Chandra observations. We clearly detect infrared emission from the jet, resulting in the most detailed multifrequency data among the jets in lobe-dominated quasars. The spectral energy distributions of the jet knots have significant variations along the jet, like the archetypal jet in 3C 273. The infrared measurements with IRAC are consistent with the previous idea that the jet has two spectral components, namely (1) the low-energy synchrotron spectrum extending from radio to infrared, and (2) the high-energy component responsible for the X-ray flux. The optical fluxes may be a mixture of the two components. We consider three radiation models for the high-energy component: inverse Compton scattering off CMB photons by radio-emitting electrons in a highly relativistic jet, synchrotron radiation by a second distinct electron population, and synchrotron radiation by ultra high energy protons. Each hypothesis leads to important insights into and constraints on particle acceleration in the jet, as well as the basic physical properties of the jet such as bulk velocity, transporting power, and particle contents.

preprint2007arXiv

Blazar microvariability at hard X-rays

Blazars are known to display strong and erratic variability at almost all the wavelengths of electromagnetic spectrum. Presently, variability studies at high-energies (hard X-rays, gamma-rays) are hampered by low sensitivity of the instruments. Nevertheless, the latest generation of satellites (INTEGRAL, Swift) have given suggestions not yet fully confirmed of variability on intraday timescales. Some specific cases recently observed are presented and physical implications are discussed (e.g. NRAO 530). The contribution that SIMBOL-X could give in this research topic is also outlined.

preprint2007arXiv

Chandra and HST observations of gamma-ray blazars: comparing jet emission at small and large scales

We present new Chandra and HST data for four gamma-ray blazars selected on the basis of radio morphology with the aim of revealing X-ray and optical emission from their jets at large scales. All the sources have been detected. Spectral Energy Distributions of the large scale jets are obtained as well as new X-ray spectra for the blazar cores. Modeling for each object the core (sub-pc scale) and large-scale (>100 kpc) jet SEDs, we derive the properties of the same jet at the two scales. The comparison of speeds and powers at different scales supports a simple scenario for the dynamics and propagation of high power relativistic jets.

preprint2007arXiv

Wind circumburst density profile: a linear E_p-E_gamma correlation

Ghirlanda et al. (2004) derived the collimation-corrected energy E_gamma for a sample of 15 bursts under the assumption of a homogeneous circumburst density profile. They found a correlation (the so-called Ghirlanda correlation) between E_gamma and the rest frame peak energy of the nuF_nu prompt spectrum (E_p). Nava et al. (2006) showed that, assuming a circumburst density distribuited with a r^-2 wind profile, the Ghirlanda correlation remains tight and becomes linear. This implies that: i) it remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors, ii) it entails that different bursts have the same number of relevant photons. We have updated these findings including recently detected bursts (21 in total), stressing the two important implications.

preprint2006arXiv

The jet/disk connection in AGN: Chandra and XMM observations of three powerful radio-loud quasars

The connection between the accretion process that powers AGN and the formation of jets is still poorly understood. Here we tackle this issue using new, deep \chandra and \xmm observations of the cores of three powerful radio loud quasars: 1136--135, 1150+497 (\chandra), and 0723+679 (\xmm), in the redshift range $z$=0.3--0.8. These sources are known from our previous \chandra snapshot survey to have kpc-scale X-ray jets. In 1136--135 and 1150+497, evidence is found for the presence of diffuse thermal X-ray emission around the cores, on scales of 40--50 kpc and with luminosity L$_{0.3-2 keV} \sim 10^{43}$ \lum, suggesting thermal emission from the host galaxy or a galaxy group. The X-ray continua of the cores in the three sources are described by an upward-curved (concave) broken power law, with photon indices $Γ_{soft} \sim 1.8-2.1$ and $Γ_{hard} \sim 1.7$ below and above $\approx 2$ keV, respectively. There is evidence for an unresolved \feka line with EW $\sim$ 70 eV in the three quasars. The Spectral Energy Distributions of the sources can be well described by a mix of jet and disk emission, with the jet dominating the radio and hard X-rays (via synchrotron and external Compton) and the disk dominating the optical/UV through soft X-rays. The ratio of the jet-to-disk powers is $\sim 1$, consistent with those derived for a number of gamma-ray emitting blazars. This indicates that near equality of accretion and jet power may be common in powerful radio-loud AGN.

preprint2006arXiv

The unexpected clustering of the optical afterglow luminosities

We studied the behaviour of the optical afterglow lightcurves of a sample of 24 Gamma--Ray Bursts (GRBs) with known redshift and published estimates of the optical extinction in the source frame, detected before the SWIFT satellite launch. We found an unexpected clustering of the optical luminosities at 12 hours in the source frame. The distribution of the optical luminosities is narrower than the distribution of X-ray luminosities at the same time. Few (3) bursts stand apart from the main optical distribution, being fainter by a factor of about 15. We also analysed the optical luminosities of the SWIFT burst with known redshift finding that the luminosity distribution is similar to the pre SWIFT GRBs one, even if they have a different mean redshift. These results can suggest the existence of a family of intrinsically optically under--luminous dark GRBs.

preprint2004arXiv

Structured jets in TeV BL Lac objects and radiogalaxies

TeV BL Lacertae objects require extreme relativistic bulk motion in the gamma-ray emission region, but at the VLBI scale their radio knots hardly move. The same sources show evidence, in radio, of a structure made of a fast spine plus a slow layer. We propose that this structure exists even on the spatial scale of regions responsible for the gamma-ray emission. One component sees the (beamed) radiation produced by the other, and this enhances the inverse Compton emission of both components. In addition, this allows the magnetic field to be nearly in equipartition with the emitting particles. The inverse Compton emission of the spine is anisotropic in its frame, possibly producing the deceleration of the spine by the Compton rocket effect. In this scenario, also the slow layer is a relatively strong high energy emitter, and thus radiogalaxies become potentially detectable by GLAST.

preprint2002arXiv

A survey of extended radio jets in AGN with Chandra and HST: First Results

We present the first results from an X-ray and optical survey of a sample of AGN radio jets with Chandra and HST. We focus here on the first six sources observed at X-rays, in four of which a bright X-ray jet was detected for the first time. In three out of four cases optical emission from the jet is also detected in our HST images. We compare the X-ray morphology with the radio as derived from improved processing of archival VLA data and we construct spectral energy distributions (SED) for the most conspicuous emission knots. In most cases the SEDs, together with the similarity of the X-ray and radio morphologies, favor an inverse Compton origin of the X-rays. The most likely origin of the seed photons is the Cosmic Microwave Background, implying the jets are still relativistic on kiloparsec scales. However, in the first knot of the PKS 1136-135 jet, X-rays are likely produced via the synchrotron process. In all four cases bulk Lorentz factors of a few are required. The radio maps of the two jets not detected by either Chandra or HST suggest that they are less beamed at large scales than the other four detected sources. Our results demonstrate that, at the sensitivity and resolution of Chandra, X-ray emission from extragalactic jets is common, yielding essential information on their physical properties.

preprint2001arXiv

Extreme synchrotron BL Lac objects

We performed an observational program with the X-ray satellite BeppoSAX to study objects with extreme synchrotron peak frequencies (nu_peak > 1 keV). Of the seven sources observed, four showed peak frequencies in the range 1-5 keV, while one (1ES 1426+428) displayed a flat power law spectrum (alpha= 0.92), locating its synchrotron peak at or above 100 keV. This is the third source of this type ever found, after Mkn 501 and 1ES 2344+514. Our data confirm the large nu_peak variability of this class of sources, compared with lower peaked objects. The high synchrotron peak energies, flagging the presence of high energy electrons, make the extreme BL Lacs also good candidates for TeV emission, and therefore good probes for the IR background.

preprint2000arXiv

Gamma-ray loud quasars: a view with BeppoSAX

We present BeppoSAX observations of three gamma-ray emitting quasars, namely 0836+710, 1510-089 and 2230+114. The three objects have been detected up to 100 keV showing extremely flat power-law spectra above 2 keV (energy index alpha_{2-10}=0.3-0.5). The soft X-ray spectrum of 0836+710 implies either an absorption column density higher than the galactic one or an intrinsically very hard slope (alpha_{0.1-1}=-0.2) below 1 keV. 1510-089 shows a soft excess, with the low energy spectrum steeper (alpha_{0.1-1}=1.6) than the high energy power-law. The results are discussed in the framework of current Inverse Compton models for the high energy emission of Flat Spectrum Radio Quasars and are used to estimate the physical quantities in the jet emitting region and to shed light on the energy transport mechanisms in jets. Finally we discuss the estimates of the jet luminosity in the context of the Blandford & Znajek mechanism for jet production.

preprint1998arXiv

AGILE: a Gamma-Ray Mission for a Light Imaging Detector

AGILE is an innovative, cost-effective gamma-ray mission proposed to the ASI Program of Small Scientific Missions. It is planned to detect gamma-rays in the 30 MeV - 50 GeV energy band and operate as an {\bf Observatory} open to the international community. Primary scientific goals include the study of AGNs, gamma-ray bursts, Galactic sources, unidentified gamma-ray sources, solar flares, and diffuse gamma-ray emission. AGILE is planned to be operational during the year 2001 for a 3-year mission. It will ideally `fill the gap' between EGRET and GLAST, and support ground-based multiwavelength studies of high-energy sources.