Source author record

F. D'Ammando

F. D'Ammando 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

138works
6topics
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

138 published item(s)

preprint2025arXiv

Optical spectroscopy of blazars for the Cherenkov Telescope Array Observatory- IV

Blazars, including BL Lacs and FSRQs, are the most luminous extragalactic γ-ray sources. They account for about 56% of the sources listed in the recent Fermi-LAT catalog (4FGL-DR4). The optical and UV spectra of BL Lacs are nearly featureless, making it difficult to precisely determine their redshifts. Consequently, nearly half of the γ-ray BL Lacs lack reliable redshift measurements. This poses a major challenge, since redshift is crucial for studying the cosmic evolution of the blazar population and γ-ray propagation studies such as indirect evidence of EBL, placing constraints on IGMF and searches for LIV and ALPs. This paper is the fourth in a series dedicated to determining the redshift of a sample of blazars identified as key targets for future observations with the Cherenkov Telescope Array Observatory (CTAO). We performed Monte Carlo simulations to select γ-ray blazars detected by Fermi-LAT with hard spectra, that lack redshift measurements. These blazars are expected to be detectable by CTAO within 30 hours or less of exposure assuming an average flux state. In this fourth paper, we report the results of detailed spectroscopic observations of 29 blazars using the ESO/VLT, Keck II, and SALT telescopes. Our analysis involved a thorough search for spectral lines in the spectra of each blazar, and when features of the host galaxy were identified, we modeled its properties. We also compared the magnitudes of the targets during the observations to their long-term light curves. In the sample studied, 9 of 29 sources were observed with a high signal-to-noise ratio (S/N > 100), while the remaining 20 were observed with a moderate or low S/N. We successfully determined firm redshifts for 12 blazars, ranging from 0.1636 to 1.1427, and identified two lower limit redshifts at z > 1.0196 and z > 1.4454. The remaining 15 BL Lac objects exhibited featureless spectra.

preprint2022arXiv

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

preprint2022arXiv

A truncated inner disc in the Seyfert 1 galaxy WKK 4438

Understanding if and when the accretion disc extends down to the innermost stable circular orbit (ISCO) is important since it is the fundamental assumption behind measuring black hole spin. Here, we examine the 2013 and 2018 NuSTAR and Swift data (0.5 - 50 keV) of the narrow-line Seyfert 1 galaxy, WKK 4438. The X-ray emission can be fitted well with models depicting a corona and blurred reflection originating from a disc around a low spin (a* ~ 0) black hole. However, such models result in unconventional values for some of the parameters (e.g. inverted emissivity profile and high coronal height). Alternatively, equally good fits can be achieved if the disc is truncated at ~10 rg and the black hole is spinning at the Thorne limit (a* = 0.998). In these cases, the model parameters are consistent with the interpretation that the corona is centrally located close to the black hole and illuminating the disc at a larger distance.

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

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

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

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 ]"

preprint2021arXiv

The First Fermi-LAT Solar Flare Catalog

We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics.

preprint2020arXiv

A parsec-scale wobbling jet in the high-synchrotron peaked blazar PG 1553+113

PG 1553+113 is the first blazar showing an approximately two-year quasi-periodic pattern in its gamma-ray light curve. Such quasi-periodicity might have a geometrical origin, possibly related to the precessing nature of the jet, or could be intrinsic to the source and related to pulsational accretion flow instabilities. By means of a ~2yr very long baseline array (VLBA) monitoring at 15, 24, and 43 GHz we investigate the source pc-scale properties during an entire cycle of gamma-ray activity in the period 2015-2017. In contrast to the well-defined periodicity in the gamma-ray emission, at radio frequencies no clear periodic pattern can be recognized. The jet position angle, constrained by means of the total intensity ridge line, varies across the different observing epochs in the range 40-60 deg. We also investigate the time evolution of the source polarization properties, including the rotation measure. The brightness temperature is found to decrease as the frequency increases with an intrinsic value of ~1.5 x 10^10 K and the estimated Doppler factor is ~1.4.

preprint2020arXiv

Fermi and Swift Observations of GRB 190114C: Tracing the Evolution of High-Energy Emission from Prompt to Afterglow

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

preprint2020arXiv

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

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

preprint2020arXiv

Jansky VLA observations of synchrotron emitting optical hotspots of 3C 227 and 3C 445 radio galaxies

We report results on deep Jansky Very Large Array A-configuration observations at 22 GHz of the hotspots of the radio galaxies 3C227 and 3C445. Synchrotron emission in the optical on scales up to a few kpc was reported for the four hotspots. Our VLA observations point out the presence of unresolved regions with upper limit to their linear size of about 100 pc. This is the first time that such compact components in hotspots have been detected in a mini-sample, indicating that they are not a peculiar characteristic of a few individual hotspots. The polarization may reach values up to 70 per cent in compact (about 0.1 kpc scale) regions within the hotspot, indicating a highly ordered magnetic field with size up to a hundred parsecs. On larger scales, the average polarization of the hotspot component is about 30-45 per cent, suggesting the presence of a significant random field component, rather than an ordered magnetic field. This is further supported by the displacement between the peaks in polarized intensity and in total intensity images that is observed in all the four hotspots. The electric vector position angle is not constant, but changes arbitrarily in the central part of the hotspot regions, whereas it is usually perpendicular to the total intensity contours of the outermost edge of the hotspot structure, likely marking the large-scale shock front. The misalignment between X-ray and radio-to-optical emission suggests that the former is tracing the current particle acceleration, whereas the latter marks older shock fronts.

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

NGC 3894: a young radio galaxy seen by Fermi-LAT

Context. According to radiative models, radio galaxies may produce gamma-ray emission from the first stages of their evolution. However, very few such galaxies have been detected by the Fermi Large Area Telescope (LAT) so far. Aims. NGC 3894 is a nearby (z = 0.0108) object that belongs to the class of compact symmetric objects (CSOs, i.e., the most compact and youngest radio galaxies), which is associated with a gamma-ray counterpart in the Fourth Fermi-LAT source catalog. Here we present a study of the source in the gamma-ray and radio bands aimed at investigating its high-energy emission and assess its young nature. Methods. We analyzed 10.8 years of Fermi-LAT data between 100 MeV and 300 GeV and determined the spectral and variability characteristics of the source. Multi-epoch very long baseline array (VLBA) observations between 5 and 15 GHz over a period of 35 years were used to study the radio morphology of NGC 3894 and its evolution. Results. NGC 3894 is detected in gamma-rays with a significance >9 sigma over the full period, and no significant variability has been observed in the gamma-ray flux on a yearly time-scale. The spectrum is modeled with a flat power law ($Γ$ = 2.0$\pm$0.1) and a flux on the order of 2.2 $\times$ 10$^{-9}$ ph cm$^{-2}$ s$^{-1}$. For the first time, the VLBA data allow us to constrain with high precision the apparent velocity of the jet and counter-jet side to be $β_{\mathrm{app,NW}}$ = 0.132$\pm$0.004 and $β_{\mathrm{app,SE}}$ = 0.065$\pm$0.003, respectively. Conclusions. Fermi-LAT and VLBA results favor the youth scenario for the inner structure of this object, with an estimated dynamical age of 59$\pm$5 years. The estimated range of viewing angle (10° < $θ$ < 21°) does not exclude a possible jet-like origin of the gamma-ray emission.

preprint2020arXiv

Particle acceleration in low-power hotspots: modelling the broad-band spectral energy distribution

The acceleration and radiative processes active in low-power radio hotspots are investigated by means of new deep near-infrared (NIR) and optical VLT observations, complemented with archival, high-sensitivity VLT, radio VLA and X-ray Chandra data. For the three studied radio galaxies (3C 105, 3C 195 and 3C 227), we confirm the detection of NIR/optical counterparts of the observed radio hotspots. We resolve multiple components in 3C 227 West and in 3C 105 South and characterize the diffuse NIR/optical emission of the latter. We show that the linear size of this component ($\gtrsim$4 kpc) makes 3C 105 South a compelling case for particles' re-acceleration in the post-shock region. Modeling of the radio-to-X-ray spectral energy distribution (SED) of 3C 195 South and 3C 227 W1 gives clues on the origin of the detected X-ray emission. In the context of inverse Compton models, the peculiarly steep synchrotron curve of 3C 195 South sets constraints on the shape of the radiating particles' spectrum that are testable with better knowledge of the SED shape at low ($\lesssim$GHz) radio frequencies and in X-rays. The X-ray emission of 3C 227 W1 can be explained with an additional synchrotron component originating in compact ($<$100 pc) regions, such those revealed by radio observations at 22 GHz, provided that efficient particle acceleration ($γ\gtrsim$10$^7$) is ongoing. The emerging picture is that of systems in which different acceleration and radiative processes coexist.

preprint2020arXiv

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

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

preprint2020arXiv

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.

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

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

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

preprint2016arXiv

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

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

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

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

preprint2016arXiv

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

preprint2016arXiv

NuSTAR, Swift, and GROND observations of the flaring MeV blazar: PMN J0641$-$0320

MeV blazars are a sub--population of the blazar family, exhibiting larger--than--average jet powers, accretion luminosities and black hole masses. Because of their extremely hard X--ray continua, these objects are best studied in the X-ray domain. Here, we report on the discovery by the $Fermi$ Large Area Telescope and subsequent follow-up observations with $NuSTAR$, $Swift$ and GROND of a new member of the MeV blazar family: PMN J0641$-$0320. Our optical spectroscopy provides confirmation that this is a flat--spectrum radio quasar located at a redshift of $z=1.196$. Its very hard $NuSTAR$ spectrum (power--law photon index of $\sim$1 up to $\sim$80 keV) indicates that the emission is produced via inverse Compton scattering off photons coming from outside the jet.The overall spectral energy distribution of PMN J0641$-$0320 is typical of powerful blazars and by reproducing it with a simple one-zone leptonic emission model we find the emission region to be located either inside the broad line region or within the dusty torus.

preprint2016arXiv

Radio follow-up of the gamma-ray flaring gravitational lens JVAS B0218+357

We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the gamma-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density variability, suggesting that no clear correlation between the high energy and low-energy emission is present. This behaviour was confirmed also by the long-term Owens Valley Radio Observatory monitoring data at 15 GHz. The milliarcsecond-scale resolution provided by the VLBA observations allowed us to resolve the two images of the lensed blazar, which have a core-jet structure. No significant morphological variation is found by the analysis of the multi-epoch data, suggesting that the region responsible for the gamma-ray variability is located in the core of the AGN, which is opaque up to the highest observing frequency of 22 GHz.

preprint2016arXiv

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

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

preprint2016arXiv

The awakening of the gamma-ray narrow-Line Seyfert 1 galaxy PKS 1502+036

After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-ray observations of PKS 1502+036 during 2008 August-2016 March by Fermi-LAT, Swift, XMM-Newton, Catalina Real-Time Transient Survey, and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV, and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the gamma-ray and 15 GHz emission peaks. This suggests that the gamma-ray emitting region is located beyond the broad line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10$^8$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar.

preprint2016arXiv

VLBA observations of radio faint Fermi-LAT sources above 10 GeV

The first Fermi-LAT High-energy source catalog (1FHL), containing gamma-ray sources detected above 10 GeV, is an ideal sample to characterize the physical properties of the most extreme gamma-ray sources. We investigate the pc scale properties of a sub-sample of radio faint 1FHL sources with the aim to confirm the proposed blazar associations, by revealing a compact high brightness temperature radio core, and we propose new low-frequency counterparts for the unassociated gamma-ray sources (UGS). Moreover, we increase the number of 1FHL sources with high resolution observations to explore the possible connection between radio and gamma rays at E >10 GeV. We observed 84 1FHL sources, mostly blazars of High Synchrotron Peaked (HSP) type, in the northern sky with the Very Long Baseline Array (VLBA) at 5 GHz. These sources lack high resolution radio observations and have at least one NVSS counterpart within the 95% confidence radius. For those sources without a well identified radio counterpart we exploit the VLBA multiple phase-center correlation capability to discern among the possible low-frequency candidates. For about 93% of the sources of our sample we reveal a compact high brightness temperature radio core, supporting their proposed blazar association. The vast majority of the detected sources are radio weak, with a median VLBI flux density value of 16.3 mJy. For the detected sources we obtain an average brightness temperature of the order of $2\times10^{10} \, \rm{K}$. We find a compact component for 16 UGS, for which we propose a new low-frequency association. We find brightness temperature values which do not require high Doppler factors, and are in agreement with the expected values for the equipartition of energy between particles and magnetic field. We find strong indications about the blazar nature of all of the detected UGS, for which we propose new low-frequency associations.

preprint2015arXiv

Flaring gamma-ray emission from high redshift blazars

High redshift blazars are among the most powerful objects in the Universe. Although they represent a significant fraction of the extragalactic hard X-ray sky, they are not commonly detected in gamma-rays. High redshift (z>2) objects represent <10 per cent of the AGN population observed by Fermi so far, and gamma-ray flaring activity from these sources is even more uncommon. The characterization of the radio-to-gamma-ray properties of high redshift blazars represent a powerful tool for the study of both the energetics of such extreme objects and the Extragalactic Background Light. We present results of a multi-band campaign on TXS 0536+145, which is the highest redshift flaring gamma-ray blazar detected so far. At the peak of the flare the source reached an apparent isotropic gamma-ray luminosity of 6.6x10^49 erg/s, which is comparable with the luminosity observed from the most powerful blazars. The physical properties derived from the multi-wavelength observations are then compared with those shown by the high redshift population. In addition preliminary results from the high redshift flaring blazar PKS 2149-306 will be discussed.

preprint2015arXiv

Gamma-ray flaring activity from the gravitationally lensed blazar PKS 1830-211 observed by Fermi LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

preprint2015arXiv

High-energy properties of the high-redshift flat spectrum radio quasar PKS 2149-306

We investigate the gamma-ray and X-ray properties of the flat spectrum radio quasar PKS 2149-306 at redshift z = 2.345. A strong gamma-ray flare from this source was detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope satellite in 2013 January, reaching on January 20 a daily peak flux of (301$\pm$36)$\times$10$^{-8}$ ph/cm$^2$/s in the 0.1-100 GeV energy range. This flux corresponds to an apparent isotropic luminosity of (1.5$\pm$0.2)$\times$10$^{50}$ erg/s, comparable to the highest values observed by a blazar so far. During the flare the increase of flux was accompanied by a significant change of the spectral properties. Moreover significant flux variations on a 6-h time-scale were observed, compatible with the light crossing time of the event horizon of the central black hole. The broad band X-ray spectra of PKS 2149-306 observed by Swift-XRT and NuSTAR are well described by a broken power-law model, with a very hard spectrum ($Γ$$_1$ $\sim$ 1) below the break energy, at E$_{\rm\,break}$ = 2.5-3.0 keV, and $Γ$$_2$ $\sim$ 1.4-1.5 above the break energy. The steepening of the spectrum below $\sim$ 3 keV may indicate that the soft X-ray emission is produced by the low-energy relativistic electrons. This is in agreement with the small variability amplitude and the lack of spectral changes in that part of the X-ray spectrum observed between the two NuSTAR and Swift joint observations. As for the other high-redshift FSRQ detected by both Fermi-LAT and Swift-BAT, the photon index of PKS 2149-306 in hard X-ray is 1.6 or lower and the average gamma-ray luminosity higher than 2$\times$10$^{48}$ erg/s.

preprint2015arXiv

Investigating powerful jets in radio-loud Narrow Line Seyfert 1s

We report results on multiband observations from radio to gamma-rays of the two radio-loud narrow-line Seyfert 1 (NLSy1) galaxies PKS 2004-447 and J1548+3511. Both sources show a core-jet structure on parsec scale, while they are unresolved at the arcsecond scale. The high core dominance and the high variability brightness temperature make these NLSy1 galaxies good gamma-ray source candidates. Fermi-LAT detected gamma-ray emission only from PKS 2004-447, with a gamma-ray luminosity comparable to that observed in blazars. No gamma-ray emission is observed for J1548+3511. Both sources are variable in X-rays. J1548+3511 shows a hardening of the spectrum during high activity states, while PKS 2004-447 has no spectral variability. A spectral steepening likely related to the soft excess is hinted below 2 keV for J1548+3511, while the X-ray spectra of PKS 2004-447 collected by XMM-Newton in 2012 are described by a single power-law without significant soft excess. No additional absorption above the Galactic column density or the presence of an Fe line is detected in the X-ray spectra of both sources.

preprint2015arXiv

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

preprint2015arXiv

Parsec scale polarization properties of the TeV blazar Markarian 421

In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.

preprint2015arXiv

PSR J1906+0722: An Elusive Gamma-ray Pulsar

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

preprint2015arXiv

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).

preprint2015arXiv

Search for extended gamma-ray emission from the Virgo galaxy cluster with Fermi-LAT

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3°that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $γ$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle σv \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $γ$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

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 1st Fermi Lat Supernova Remnant Catalog

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

preprint2015arXiv

The first gamma-ray detection of the narrow-line Seyfert 1 FBQS J1644+2619

We report the discovery of gamma-ray emission from the narrow-line Seyfert 1 (NLSy1) galaxy FBQS J1644+2619 by the Large Area Telescope on board the Fermi satellite. The Third Fermi LAT Source catalogue reports an unidentified gamma-ray source, detected over the first four years of Fermi operation, 0.23 deg from the radio position of the NLSy1. Analysing 76 months of gamma-ray data (2008 August 4 - 2014 December 31) we are able to better constrain the localization of the gamma-ray source. The new position of the gamma-ray source is 0.05 deg from FBQS J1644+2619, suggesting a spatial association with the NLSy1. This is the sixth NLSy1 detected at high significance by Fermi-LAT so far. Notably, a significant increase of activity was observed in gamma-rays from FBQS J1644+2619 during 2012 July-October, and an increase of activity in V-band was detected by the Catalina Real-Time Sky Survey in the same period.

preprint2015arXiv

The Third Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10°), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

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.

preprint2014arXiv

A strong radio brightening at the jet base of M87 during the elevated very-high-energy gamma-ray state in 2012

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure at radio. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central supermassive black hole.

preprint2014arXiv

A strong radio brightening at the jet base of M87 during the elevated very-high-energy gamma-ray state in 2012

We report our intensive radio monitoring observations of the jet in M87 with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi-LAT. During this period, an elevated level of the M87 flux is reported at VHE gamma rays. We detected a remarkable increase of the radio flux density from the unresolved jet base (radio core) with VERA at 22 and 43GHz coincident with the VHE activity. Meanwhile, we confirmed with EVN at 5GHz that HST-1 (an alternative gamma-ray production candidate site) remained quiescent in terms of its flux density and structure. These results in the radio bands strongly suggest that the VHE gamma-ray activity in 2012 originates in the jet base within 0.03pc or 56 Schwarzschild radii from the central supermassive black hole. We further conducted VERA astrometry for the M87 core during the flaring period, and detected core shifts between 22 and 43GHz. We also discovered a clear frequency-dependent evolution of the radio core flare at 43, 22 and 5GHz; the radio flux density increased more rapidly at higher frequencies with a larger amplitude, and the light curves clearly showed a time-lag between the peaks at 22 and 43GHz. This indicates that a new radio-emitting component was created near the black hole in the period of the VHE event, and then propagated outward with progressively decreasing synchrotron opacity. By combining these results, we estimated an apparent speed of the newborn component, and derived a sub-luminal speed of less than ~0.2c. This value is significantly slower than the super-luminal (~1.1c) features that appeared from the core during the prominent VHE flaring event in 2008, suggesting that the stronger VHE activity can be associated with the production of the higher Lorentz factor jet.

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

Dark Matter Constraints from Observations of 25 Milky Way Satellite Galaxies with the Fermi Large Area Telescope

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

preprint2014arXiv

Deep Broadband Observations of the Distant Gamma-ray Blazar PKS 1424+240

We present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hours of VERITAS observations from three years, a multiwavelength light curve and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured differential very high energy (VHE; $E\ge100$ GeV) spectral indices are $Γ=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than $τ=2$, where it is postulated that any variability would be small and occur on longer than year timescales if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.

preprint2014arXiv

Exploring the multiband emission of TXS 0536+145: the most distant gamma-ray flaring blazar

We report results of a multiband monitoring campaign of the flat spectrum radio quasar TXS 0536+145 at redshift 2.69. This source was detected during a very high gamma-ray activity state in 2012 March by the Large Area Telescope on board Fermi, becoming the gamma-ray flaring blazar at the highest redshift detected so far. At the peak of the flare the source reached an apparent isotropic gamma-ray luminosity of 6.6 x 10^49 erg/s which is comparable to the values achieved by the most luminous blazars. This activity triggered radio-to-X-rays monitoring observations by Swift, Very Long Baseline Array, European VLBI Network, and Medicina single-dish telescope. Significant variability was observed from radio to X-rays supporting the identification of the gamma-ray source with TXS 0536+145. Both the radio and gamma-ray light curves show a similar behaviour, with the gamma-rays leading the radio variability with a time lag of about 4-6 months. The luminosity increase is associated with a flattening of the radio spectrum. No new superluminal component associated with the flare was detected in high resolution parsec-scale radio images. During the flare the gamma-ray spectrum seems to deviate from a power law, showing a curvature that was not present during the average activity state. The gamma-ray properties of TXS 0536+145 are consistent with those shown by the high-redshift gamma-ray blazar population.

preprint2014arXiv

Fermi-LAT Detection of Gravitational Lens Delayed Gamma-ray Flares from Blazar B0218+357

Using data from the Fermi Large Area Telescope (LAT), we report the first clear gamma-ray measurement of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected in B0218+357, a known double-image lensed system, during a period of enhanced gamma-ray activity with peak fluxes consistently observed to reach >20-50 times its previous average flux. An auto-correlation function analysis identified a delay in the gamma-ray data of 11.46 +/- 0.16 days (1 sigma) that is ~1 day greater than previous radio measurements. Considering that it is beyond the capabilities of the LAT to spatially resolve the two images, we nevertheless decomposed individual sequences of superposing gamma-ray flares/delayed emissions. In three such ~8-10 day-long sequences within a ~4-month span, considering confusion due to overlapping flaring emission and flux measurement uncertainties, we found flux ratios consistent with ~1, thus systematically smaller than those from radio observations. During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we observed flux doubling timescales of ~3-6 hrs implying as well extremely compact gamma-ray emitting regions.

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

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $γ$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

preprint2014arXiv

Infrared properties of blazars: putting the GASP-WEBT sources into context

The infrared properties of blazars can be studied from the statistical point of view with the help of sky surveys, like that provided by the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All Sky Survey (2MASS). However, these sources are known for their strong and unpredictable variability, which can be monitored for a handful of objects only. In this paper we consider the 28 blazars (14 BL Lac objects and 14 flat-spectrum radio quasars, FSRQs) that are regularly monitored by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) since 2007. They show a variety of infrared colours, redshifts, and infrared-optical spectral energy distributions (SEDs), and thus represent an interesting mini-sample of bright blazars that can be investigated in more detail. We present near-IR light curves and colours obtained by the GASP from 2007 to 2013, and discuss the infrared-optical SEDs. These are analysed with the aim of understanding the interplay among different emission components. BL Lac SEDs are accounted for by synchrotron emission plus an important contribution from the host galaxy in the closest objects, and dust signatures in 3C 66A and Mkn 421. FSRQ SEDs require synchrotron emission with the addition of a quasar-like contribution, which includes radiation from a generally bright accretion disc, broad line region, and a relatively weak dust torus.

preprint2014arXiv

Limb-Brightened Jet of 3C 84 Revealed by the 43-GHz Very-Long-Baseline-Array Observation

We present a study of sub-pc scale radio structure of the radio galaxy 3C 84/NGC 1275 based on the Very Long Baseline Array (VLBA) data at 43 GHz. We discover a limb-brightening in the "restarted" jet associated with the 2005 radio outburst. In the 1990s, the jet structure was ridge-brightening rather than limb-brightening, despite the observations being done with similar angular resolution. This indicates that the transverse jet structure has changed recently. This change in the morphology shows an interesting agreement with the $γ$-ray flux increase, i.e., the $γ$-ray flux in 1990s was at least seven times lower than the current one. One plausible explanation for the limb-brightening is the velocity structure of the jet in the context of the stratified jet, which is a successful scenario to explain the $γ$-ray emission in some active galactic nuclei (AGNs). If this is the case, the change in apparent transverse structure might be caused by the change in the transverse velocity structure. We argue the possibility that the transition from ridge-brightening to limb-brightening is related to the $γ$-ray time variability on the timescale of decades. We also discuss the collimation profile of the jet.

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

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

Radio-Gamma-ray connection and spectral evolution in 4C +49.22 (S4 1150+49): the Fermi, Swift and Planck view

The Large Area Telescope on board the Fermi Gamma-ray Space Telescope detected a strong gamma-ray flare on 2011 May 15 from a source identified as 4C 49.22, a flat spectrum radio quasar also known as S4 1150+49. This blazar, characterised by a prominent radio-optical-X-ray jet, was in a low gamma-ray activity state during the first years of Fermi observations. Simultaneous observations during the quiescent, outburst and post-flare gamma-ray states were obtained by Swift, Planck and optical-IR-radio telescopes (INAOE, Catalina CSS, VLBA, Metsahovi). The flare is observed from microwave to X-ray bands with correlated variability and the Fermi, Swift and Planck data for this FSRQ show some features more typical of BL Lac objects, like the synchrotron peak in the optical band that outshines the thermal blue-bump emission, and the X-ray spectral softening. Multi-epoch VLBA observations show the ejection of a new component close in time with the GeV gamma-ray flare. The radio-to-gamma-ray spectral energy distribution is modeled and fitted successfully for the outburst and the post-flare epochs using either a single flaring blob with two emission processes (synchrotron self Compton, and external-radiation Compton), and a two-zone model with SSC-only mechanism.

preprint2014arXiv

Search for cosmic-ray induced gamma-ray emission in Galaxy Clusters

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $γ$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $γ$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7σ$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $γ$-ray flux from individual clusters in our sample.

preprint2014arXiv

The most powerful flaring activity from the NLSy1 PMN J0948+0022

We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

preprint2014arXiv

The physics of the radio emission in the quiet side of the AGN population with the SKA

Despite targets of many multiwavelength campaigns, the main physical processes at work in AGN are still under debate. In particular the origin of the radio emission and the mechanisms involved are among the open questions in astrophysics. In the radio-loud AGN population the radio emission is linked to the presence of bipolar outflows of relativistic jets. However, the large majority of the AGN population do not form powerful highly-relativistic jets on kpc scales and are characterized by radio luminosity up to 10^23 W/Hz at 1.4 GHz, challenging our knowledge on the physical processes at the basis of the radio emission in radio-quiet objects. The main mechanisms proposed so far are synchrotron radiation from mildly relativistic mini-jets, thermal cyclo-synchrotron emission by low-efficiency accretion flow (like ADAF or ADIOS), or thermal free-free emission from the X-ray heated corona or wind. The difficulty in understanding the main mechanism involved is related to the weakness of these objects, which precludes the study of non-local radio-quiet AGN. Multifrequency, high-sensitivity radio observations are crucial to constrain the nature of the power engine, and they may help in distinguishing between the contribution from star formation and AGN activity. The advent of the SKA, with its sub-arcsecond resolution and unprecedented sensitivity will allow us to investigate these processes in radio-quiet AGN, even at high redshift for the first time. Both the broad-band radio spectrum and the polarization information will help us in disentangling between non-thermal and thermal origin of the radio emission. The jump in sensitivity of a few order of magnitudes at the (sub-)uJy level will enable us to detect radio emission from a large number of radio-quiet AGN at high redshift, providing a fundamental step in our understanding of their cosmological evolution. (Abridged)

preprint2014arXiv

The spectrum of isotropic diffuse gamma-ray emission between 100 MeV and 820 GeV

The γ-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse γ-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

preprint2014arXiv

Two active states of the narrow-line gamma-ray-loud AGN GB 1310+487

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

preprint2014arXiv

Unusual Flaring Activity in the Blazar PKS 1424-418 during 2008-2011

Context. Blazars are a subset of active galactic nuclei (AGN) with jets that are oriented along our line of sight. Variability and spectral energy distribution (SED) studies are crucial tools for understanding the physical processes responsible for observed AGN emission. Aims. We report peculiar behaviour in the bright gamma-ray blazar PKS 1424-418 and use its strong variability to reveal information about the particle acceleration and interactions in the jet. Methods. Correlation analysis of the extensive optical coverage by the ATOM telescope and nearly continuous gamma-ray coverage by the Fermi Large Area Telescope is combined with broadband, time-dependent modeling of the SED incorporating supplemental information from radio and X-ray observations of this blazar. Results. We analyse in detail four bright phases at optical-GeV energies. These flares of PKS 1424-418 show high correlation between these energy ranges, with the exception of one large optical flare that coincides with relatively low gamma-ray activity. Although the optical/gamma-ray behaviour of PKS 1424-418 shows variety, the multiwavelength modeling indicates that these differences can largely be explained by changes in the flux and energy spectrum of the electrons in the jet that are radiating. We find that for all flares the SED is adequately represented by a leptonic model that includes inverse Compton emission from external radiation fields with similar parameters. Conclusions. Detailed studies of individual blazars like PKS 1424-418 during periods of enhanced activity in different wavebands are helping us identify underlying patterns in the physical parameters in this class of AGN.

preprint2014arXiv

Very Long Baseline polarimetry and the gamma-ray connection in Markarian 421 during the broadband campaign in 2011

We analyzed data in polarized intensity obtained with the Very Long Baseline Array (VLBA) at twelve epochs (one observation per month from January to December 2011) at 15, 24, and 43 GHz. For the absolute orientation of the electric vector position angles (EVPA) we used the D-terms method. We also used gamma-ray data from the Fermi Large Area Telescope on weekly time bins throughout 2011. The source shows polarized emission, and its properties vary with time, frequency, and location along the jet. The core mean polarization fraction is generally between 1% and 2%, with a 4% peak at 43 GHz in March; the polarization angle is variable, mainly at 15 GHz, where it changes frequently, and less so at 43 GHz, where it oscillates in the range 114 - 173 deg. The jet polarization properties are more stable, with a fractional polarization of around 16% and a polarization angle nearly perpendicular to the jet axis. The average flux and photon index at gamma-ray energies are (17.7+-0.5)x10^8 ph cm^-2 s^-1 and the photon index is 1.77+-0.02. The gamma-ray light curve shows variability, with a main peak that appears to be associated with the peak in the core polarized emission at 43 GHz, as well as with the total intensity light curve. A discrete correlation function analysis yields a correlation coefficient of 0.54 at zero delay, with a significance level above 99.7%. We accurately determine the polarization properties of Mrk 421, both in the core and in the jet region. The radio and gamma-ray light curves are correlated. The observed EVPA variability at 15 GHz is partly due to opacity and partly to a variable Faraday rotation effect. To explain the residual variability of the intrinsic polarization angle and the low degree of polarization in the core region, we invoke a blend of variable cross-polarized subcomponents with different polarization properties within the beam.

preprint2013arXiv

Constraints on the Galactic Halo Dark Matter from Fermi-LAT Diffuse Measurements

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

preprint2013arXiv

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

preprint2013arXiv

Detection of the Characteristic Pion-Decay Signature in Supernova Remnants

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

preprint2013arXiv

Flaring gamma-ray AGNs seen with the FERMI LAT

In the first 3.5 years of operations the Fermi LAT detected several sources with daily fluxes brighter than F(E>100MeV) > 10^{-6} ph/cm^2/s, the threshold set by the Fermi Collaboration for issuing an Astronomer Telegram. We focus the attention on these flaring sources, most of which are blazars, known to be extremely variable over the whole electromagnetic spectrum, from radio to gamma-ray energies and we investigate the properties of the selected sample comparing them to the general characteristics of the Fermi catalog sources.

preprint2013arXiv

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

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

preprint2013arXiv

Multifrequency studies of the narrow-line Seyfert 1 galaxy SBS 0846+513

The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.

preprint2013arXiv

Multiwavelength observations of the gamma-ray emitting narrow-line Seyfert 1 PMN J0948+0022 in 2011

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

preprint2013arXiv

On the connection between radio and gamma rays. Variability and polarization properties in relativistic jets

Relativistic jets are one of the most powerful manifestations of the release of energy produced around supermassive black holes at the centre of active galactic nuclei (AGN). Their emission is observed across the entire electromagnetic spectrum, from the radio band to gamma rays. Despite decades of efforts, many aspects of the physics of relativistic jets remain elusive. In particular, the location and the mechanisms responsible for the high-energy emission and the connection of the variability at different wavelengths are among the greatest challenges in the study of AGN. From the comparison of the radio and gamma-ray light curves of gamma-ray flaring objects, there is evidence that some flares, either in radio or in gamma rays, have not an obvious connection at the other extreme of the electromagnetic spectrum, like in the case of the Narrow-Line Seyfert 1 SBS 0846+513. An intriguing aspect pointed out by high resolution radio observations is the change of the polarization properties close in time with some high energy flares. In particular, in PKS 1510-089 and 3C 454.3 a rotation of almost 90 degrees has been observed after strong gamma-ray flares. The swing of the polarization angle may be related either to the propagation of a shock along the jet that orders the magnetic field, or a change of the opacity regime.

preprint2013arXiv

The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

preprint2013arXiv

The ordinary life of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036

We report on multifrequency observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036 performed from radio to gamma-rays during 2008 August-2012 November by Fermi-LAT, Swift (XRT and UVOT), OVRO, VLBA, and VLA. No significant variability has been observed in gamma-rays, with 0.1-100 GeV flux that ranged between (3-7)x10^-8 ph/cm^2/s using 3-month time bins. The photon index of the LAT spectrum (Gamma=2.60+/-0.06) and the apparent isotropic gamma-ray luminosity, L(0.1-100 GeV)= 7.8x10^45 erg/s, over 51 months are typical of a flat spectrum radio quasar. The radio spectral variability and the one-sided structure, in addition to the observed gamma-ray luminosity, suggest a relativistic jet with a high Doppler factor. In contrast to SBS 0846+513, the VLBA at 15 GHz did not observe superluminal motion for PKS 1502+036. Despite having the optical characteristics typical of a narrow-line Seyfert 1 galaxy, radio and gamma-ray properties of PKS 1502+036 are found to be similar to those of a blazar at the low end of the black hole mass distribution for blazars. This is in agreement with what has been found in the case of the other gamma-ray emitting narrow-line Seyfert 1 SBS 0846+513.

preprint2013arXiv

Very Rapid High-Amplitude Gamma-ray Variability in Luminous Blazar PKS 1510-089 Studied with Fermi-LAT

Here we report on the detailed analysis of the gamma-ray light curve of a luminous blazar PKS1510-089 observed in the GeV range with the Large Area Telescope (LAT) onboard the Fermi satellite during the period 2011 September -- December. By investigating the properties of the detected three major flares with the shortest possible time binning allowed by the photon statistics, we find a variety of temporal characteristics and variability patterns. This includes a clearly asymmetric profile (with a faster flux rise and a slower decay) of the flare resolved on sub-daily timescales, a superposition of many short uncorrelated flaring events forming the apparently coherent longer-duration outburst, and a huge single isolated outburst unresolved down to the timescale of three-hours. In the latter case we estimate the corresponding gamma-ray flux doubling timescale to be below one hour, which is extreme and never previously reported for any active galaxy in the GeV range. The other unique finding is that the total power released during the studied rapid and high-amplitude flares constitute the bulk of the power radiatively dissipated in the source, and a significant fraction of the total kinetic luminosity of the underlying relativistic outflow. Our analysis allows us to access directly the characteristic timescales involved in shaping the energy dissipation processes in the source, and to provide constraints on the location and the structure of the blazar emission zone in PKS1510-089.

preprint2012arXiv

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

preprint2012arXiv

Characterization of a tagged $γ$-ray beam line at the DA$Φ$NE Beam Test Facility

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

preprint2012arXiv

Fermi LAT Search for Dark Matter in Gamma-ray Lines and the Inclusive Photon Spectrum

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

preprint2012arXiv

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

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

preprint2012arXiv

Gamma-ray Activity in the Crab Nebula: The Exceptional Flare of April 2011

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

preprint2012arXiv

Gamma-ray observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1° \times1° pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212°), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

preprint2012arXiv

Multi-wavelength observations of blazar AO 0235+164 in the 2008-2009 flaring state

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the γ -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

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

PKS 2123-463: a confirmed gamma-ray blazar at high redshift

The flat spectrum radio quasar (FSRQ) PKS 2123-463 was associated in the First Fermi-LAT source catalog with the gamma-ray source 1FGL J2126.1-4603, but when considering the full first two years of Fermi observations, no gamma-ray source at a position consistent with this FSRQ was detected, and thus PKS 2123-463 was not reported in the Second Fermi-LAT source catalog. On 2011 December 14 a gamma-ray source positionally consistent with PKS 2123-463 was detected in flaring activity by Fermi-LAT. This activity triggered radio-to-X-ray observations by the Swift, GROND, ATCA, Ceduna, and KAT-7 observatories. Results of the localization of the gamma-ray source over 41 months of Fermi-LAT operation are reported here in conjunction with the results of the analysis of radio, optical, UV and X-ray data collected soon after the gamma-ray flare. The strict spatial association with the lower energy counterpart together with a simultaneous increase of the activity in optical, UV, X-ray and gamma-ray bands led to a firm identification of the gamma-ray source with PKS 2123-463. A new photometric redshift has been estimated as z = 1.46+/-0.05 using GROND and Swift/UVOT observations, in rough agreement with the disputed spectroscopic redshift of z = 1.67. We fit the broadband spectral energy distribution with a synchrotron/external Compton model. We find that a thermal disk component is necessary to explain the optical/UV emission detected by Swift/UVOT. This disk has a luminosity of about 1.8x10^46 erg/s, and a fit to the disk emission assuming a Schwarzschild (i.e., nonrotating) black hole gives a mass of about 2x10^9 solar masses. This is the first black hole mass estimate for this source.

preprint2012arXiv

Radio and gamma-ray follow-up of the exceptionally high activity state of PKS 1510-089 in 2011

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

preprint2012arXiv

SBS 0846+513: a new gamma-ray emitting Narrow-Line Seyfert 1 galaxy

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

preprint2012arXiv

Simultaneous Planck, Swift, and Fermi observations of X-ray and gamma-ray selected blazars

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <α>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

preprint2012arXiv

Variability of the blazar 4C 38.41 (B3 1633+382) from GHz frequencies to GeV energies

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and γ-ray flux apparently fades in time, likely because of an increasing optical to γ-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor δ. Under the hypothesis that in the period 2008-2011 all the γ-ray and optical variability on a one-week timescale were due to changes in δ, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle θ only, then θ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

preprint2011arXiv

AGILE detection of Cygnus X-3 γ-ray active states during the period mid-2009/mid-2010

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

preprint2011arXiv

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

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

preprint2011arXiv

Broad Line Radio Galaxies Observed with Fermi-LAT: The Origin of the GeV Gamma-Ray Emission

We report on a detailed investigation of the gamma-ray emission from 18 broad line radio galaxies (BLRGs) based on two years of Fermi Large Area Telescope (LAT) data. We confirm the previously reported detections of 3C 120 and 3C 111 in the GeV photon energy range; a detailed look at the temporal characteristics of the observed gamma-ray emission reveals in addition possible flux variability in both sources. No statistically significant gamma-ray detection of the other BLRGs was however found in the considered dataset. Though the sample size studied is small, what appears to differentiate 3C 111 and 3C 120 from the BLRGs not yet detected in gamma-rays is the particularly strong nuclear radio flux. This finding, together with the indications of the gamma-ray flux variability and a number of other arguments presented, indicate that the GeV emission of BLRGs is most likely dominated by the beamed radiation of relativistic jets observed at intermediate viewing angles. In this paper we also analyzed a comparison sample of high accretion-rate Seyfert 1 galaxies, which can be considered radio-quiet counterparts of BLRGs, and found none were detected in gamma-rays. A simple phenomenological hybrid model applied for the broad-band emission of the discussed radio-loud and radio-quiet type 1 active galaxies suggests that the relative contribution of the nuclear jets to the accreting matter is > 1 percent on average for BLRGs, whilst <0.1 percent for Seyfert 1 galaxies.

preprint2011arXiv

Discovery of Powerful Gamma-Ray Flares from the Crab Nebula

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

preprint2011arXiv

Fermi Large Area Telescope Detection of Bright Gamma-ray Outbursts from a Peculiar Quasar 4C +21.35

In this paper we report on the two-year-long Fermi-LAT observation of the peculiar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from the start of science operations of the Fermi Gamma-ray Space Telescope in 2008 August until 2009 September, and then became more active, with gradually increasing flux and some moderately-bright flares. In 2010 April and June, 4C +21.35 underwent a very strong GeV outburst composed of several major flares characterized by rise and decay timescales of the order of a day. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form with spectral breaks observed near 1-3 GeV photon energies. We demonstrate that, at least during the major flares, the jet in 4C +21.35 carried a total kinetic luminosity comparable to the total accretion power available to feed the outflow. We also discuss the origin of the break observed in the flaring spectra of 4C +21.35. We show that, in principle, a model involving annihilation of the GeV photons on the He II Lyman recombination continuum and line emission of "broad line region" clouds may account for such. However, we also discuss the additional constraint provided by the detection of 4C +21.35 at 0.07-0.4 TeV energies by the MAGIC telescope, which coincided with one of the GeV flares of the source. We argue that there are reasons to believe that the $lesssim$,TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely to be produced inside the broad line region zone of highest ionization ($sim 10^{17}$,cm from the nucleus), but instead originates further away from the active center, namely around the characteristic scale of the hot dusty torus surrounding the 4C +21.35 nucleus ($sim 10^{19}$,cm).

preprint2011arXiv

First results about on-ground calibration of the Silicon Tracker for the AGILE satellite

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

preprint2011arXiv

Identification of gamma-ray emission from 3C345 and NRAO512

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

preprint2011arXiv

Modeling the flaring activity of the high z, hard X-ray selected blazar IGR J22517+2217

We present new Suzaku and Fermi data, and re-analyzed archival hard X-ray data from INTEGRAL and Swift-BAT survey, to investigate the physical properties of the luminous, high-redshift, hard X-ray selected blazar IGR J22517+2217, through the modelization of its broad band spectral energy distribution (SED) in two different activity states. Through the analysis of the new Suzaku data and the flux selected data from archival hard X-ray observations, we build the source SED in two different states, one for the newly discovered flare occurred in 2005 and one for the following quiescent period. Both SEDs are strongly dominated by the high energy hump peaked at 10^20 -10^22 Hz, that is at least two orders of magnitude higher than the low energy (synchrotron) one at 10^11 -10^14 Hz, and varies by a factor of 10 between the two states. In both states the high energy hump is modeled as inverse Compton emission between relativistic electrons and seed photons produced externally to the jet, while the synchrotron self-Compton component is found to be negligible. In our model the observed variability can be accounted for by a variation of the total number of emitting electrons, and by a dissipation region radius changing from within to outside the broad line region as the luminosity increases. In its flaring activity, IGR J22517+2217 shows one of the most powerful jet among the population of extreme, hard X-ray selected, high redshift blazar observed so far.

preprint2011arXiv

Multi-epoch parsec-scale observations of the blazar PKS 1510-089

(Abridged) We investigate the flux density variability and changes in the parsec-scale radio structure of the flat spectrum radio quasar PKS 1510-089. This source was target of multi-epoch VLBI and Space-VLBI observations at 4.8, 8.4 and 22 GHz carried out between 1999 and 2001. The comparison of the parsec-scale structure observed at different epochs shows the presence of a non-stationary jet feature moving with a superluminal apparent velocity of 16.2c+-0.7c. Over three epochs at 8.4 GHz during this period the core flux density varies of about 50%, while the scatter in the jet flux density is within 10%. The polarization percentage of both core and jet components significantly change from 2 to 9 per cent, while the polarization angle of the core shows an abrupt change of about 90 degrees becoming roughly perpendicular to the jet direction, consistent with a change in the opacity. To complete the picture of the physical processes at work, we complemented our observations with multi-epoch VLBA data at 15 GHz from the MOJAVE programme spanning a time baseline from 1995 to 2010. Since 1995 jet components are ejected roughly once per year with the same position angle and an apparent speed between 15c and 20c, indicating that no jet precession is taking place on a timescale longer than a decade in our frame. The variability of the total intensity flux density together with variations in the polarization properties may be explained assuming either a change between the optically-thick and -thin regimes produced by a shock that varies the opacity, or a highly ordered magnetic field produced by the compression of the relativistic plasma by a shock propagating along the jet. Taking into account the high gamma-ray emission observed from this source by the AGILE and Fermi satellites we investigated the connection between the radio and gamma-ray activity during 2007-2010.

preprint2011arXiv

Spectral Energy Distribution variation in BL Lacs and FSRQs

We present the results of our study of spectral energy distributions (SEDs) of a sample of ten low- to intermediate-synchrotron-peaked blazars. We investigate some of the physical parameters most likely responsible for the observed short-term variations in blazars. To do so, we focus on the study of changes in the SEDs of blazars corresponding to changes in their respective optical fluxes. We model the observed spectra of blazars from radio to optical frequencies using a synchrotron model that entails a log-parabolic distribution of electron energies. A significant correlation among the two fitted spectral parameters ($a$, $b$) of log-parabolic curves and a negative trend among the peak frequency and spectral curvature parameter, $b$, emphasize that the SEDs of blazars are fitted well by log-parabolic curves. On considering each model parameter that could be responsible for changes in the observed SEDs of these blazars, we find that changes in the jet Doppler factors are most important.

preprint2011arXiv

The AGILE observations of the hard and bright GRB 100724B

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

preprint2011arXiv

The 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 flaring blazars of the first 1.5 years of the AGILE mission

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

preprint2011arXiv

The 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

Young radio sources: the duty-cycle of the radio emission and prospects for gamma-ray emission

The evolutionary stage of a powerful radio source originated by an AGN is related to its linear size. In this context, compact symmetric objects (CSOs), which are powerful and intrinsically small objects, should represent the young stage in the individual radio source life. However, the fraction of young radio sources in flux density-limited samples is much higher than what expected from the number counts of large radio sources.This indicates that a significant fraction of young radio sources does not develop to the classical Fanaroff-Riley radio galaxies,suggesting an intermittent jet activity. As the radio jets are expanding within the dense and inhomogeneous interstellar medium,the ambient may play a role in the jet growth, for example slowing down or even disrupting its expansion when a jet-cloud interaction takes place. Moreover, this environment may provide the thermal seed photons that scattered by the lobes' electrons may be responsible for high energy emission, detectable by Fermi-LAT.

preprint2010arXiv

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

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

preprint2010arXiv

AGILE detection of GeV gamma-ray emission from the SNR W28

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

preprint2010arXiv

Agile Observations of the "Soft" Gamma-Ray Pulsar PSR B1509-58

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

preprint2010arXiv

Detection of Gamma-Ray Emission from the Vela Pulsar Wind Nebula with AGILE

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

Flux variability in the FSRQ 1510-089. A radio-gamma perspective

The knowledge of the physical conditions occurring in the relativistic jet of radio-loud active galactic nuclei (AGNs) is important to understand the mechanisms at the basis of their multiband emission. From parsec-scale radio observations of blazar objects it has been suggested a connection between the ejection of new jet features and strong gamma-ray flares. We present results from multi-epoch Very Long Baseline Interferometer (VLBI) and Space-VLBI observations of the Flat Spectrum Radio Quasar (FSRQ) PKS 1510-089. The comparison of the radio structure observed at different epochs shows the presence of jet features moving with highly superluminal apparent velocity. Radio flux density variability and changes in the source structure and in the polarization properties are then compared with the information on the gamma-ray emission in order to find a possible connection between radio and gamma-ray emission.

preprint2010arXiv

Gamma-Ray Localization of Terrestrial Gamma-Ray Flashes

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

preprint2010arXiv

Multiwavelength observations of 3C 454.3. III. Eighteen months of AGILE monitoring of the "Crazy Diamond"

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

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

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

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

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

preprint2010arXiv

The December 2009 gamma-ray flare of 3C 454.3: the multifrequency campaign

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

preprint2010arXiv

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

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

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2009arXiv

Discovery of extreme particle acceleration in the microquasar Cygnus X-3

The study of relativistic particle acceleration is a major topic of high-energy astrophysics. It is well known that massive black holes in active galaxies can release a substantial fraction of their accretion power into energetic particles, producing gamma-rays and relativistic jets. Galactic microquasars (hosting a compact star of 1-10 solar masses which accretes matter from a binary companion) also produce relativistic jets. However, no direct evidence of particle acceleration above GeV energies has ever been obtained in microquasar ejections, leaving open the issue of the occurrence and timing of extreme matter energization during jet formation. Here we report the detection of transient gamma-ray emission above 100 MeV from the microquasar Cygnus X-3, an exceptional X-ray binary which sporadically produces powerful radio jets. Four gamma-ray flares (each lasting 1-2 days) were detected by the AGILE satellite simultaneously with special spectral states of Cygnus X-3 during the period mid-2007/mid-2009. Our observations show that very efficient particle acceleration and gamma-ray propagation out of the inner disk of a microquasar usually occur a few days before major relativistic jet ejections. Flaring particle energies can be thousands of times larger than previously detected maximum values (with Lorentz factors of 105 and 102 for electrons and protons, respectively). We show that the transitional nature of gamma-ray flares and particle acceleration above GeV energies in Cygnus X-3 is clearly linked to special radio/X-ray states preceding strong radio flares. Thus gamma-rays provide unique insight into the nature of physical processes in microquasars.

preprint2009arXiv

First AGILE Catalog of High Confidence Gamma-Ray Sources

We present the first catalog of high-confidence gamma-ray sources detected by the AGILE satellite during observations performed from July 9, 2007 to June 30, 2008. Catalogued sources are detected by merging all the available data over the entire time period. AGILE, launched in April 2007, is an ASI mission devoted to gamma-ray observations in the 30 MeV - 50 GeV energy range, with simultaneous X-ray imaging capability in the 18-60 keV band. This catalog is based on Gamma-Ray Imaging Detector (GRID) data for energies greater than 100 MeV. For the first AGILE catalog we adopted a conservative analysis, with a high-quality event filter optimized to select gamma-ray events within the central zone of the instrument Field of View (radius of 40 degrees). This is a significance-limited (4 sigma) catalog, and it is not a complete flux-limited sample due to the non-uniform first year AGILE sky coverage. The catalog includes 47 sources, 21 of which are associated with confirmed or candidate pulsars, 13 with Blazars (7 FSRQ, 4 BL Lacs, 2 unknown type), 2 with HMXRBs, 2 with SNRs, 1 with a colliding-wind binary system, 8 with unidentified sources.

preprint2009arXiv

Monitoring the hard X-ray sky with SuperAGILE

SuperAGILE is the hard X-ray monitor of the AGILE gamma ray mission, in orbit since 23$^{rd}$ April 2007. It is an imaging experiment based on a set of four independent silicon strip detectors, equipped with one-dimensional coded masks, operating in the nominal energy range 18-60 keV. The main goal of SuperAGILE is the observation of cosmic sources simultaneously with the main gamma-ray AGILE experiment, the Gamma Ray Imaging Detector (GRID). Given its $\sim$steradian-wide field of view and its $\sim$15 mCrab day-sensitivity, SuperAGILE is also well suited for the long-term monitoring of Galactic compact objects and the detection of bright transients. The SuperAGILE detector properties and design allow for a 6 arcmin angular resolution in each of the two independent orthogonal projections of the celestial coordinates. Photon by photon data are continuously available by the experiment telemetry, and are used to derive images and fluxes of individual sources, with integration times depending on the source intensity and position in the field of view. In this paper we report on the main scientific results achieved by SuperAGILE over its first two years in orbit, until April 2009.

preprint2009arXiv

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

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

preprint2009arXiv

Multiwavelength observations of a TeV-Flare from W Comae

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

preprint2009arXiv

Powerful high energy emission of the remarkable BL Lac object S5 0716+714

BL Lac objects of the intermediate subclass (IBLs) are known to emit a substantial fraction of their power in the energy range 0.1--10 GeV. Detecting gamma-ray emission from such sources provides therefore a direct probe of the emission mechanisms and of the underlying powerhouse. The AGILE gamma-ray satellite detected the remarkable IBL S5 0716+714 (z \simeq 0.3) during a high state in the period from 2007 September - October, marked by two very intense flares reaching peak fluxes of 200\times10^{-8} ph / cm^2 s above 100 MeV, with simultaneous optical and X-ray observations. We present here a theoretical model for the two major flares and discuss the overall energetics of the source. We conclude that 0716+714 is among the brightest BL Lac's ever detected at gamma-ray energies. Because of its high power and lack of signs for ongoing accretion or surrounding gas, the source is an ideal candidate to test the maximal power extractable from a rotating supermassive black hole via the pure Blandford-Znajek (BZ) mechanism. We find that during the 2007 gamma-ray flares our source approached or just exceeded the upper limit set by BZ for a black hole of mass 10^9 M_sun

preprint2009arXiv

SSC radiation in BL Lac sources, the end of the tether

The synchrotron-self Compton (SSC) radiation process is widely held to provide a close representation of the double peaked spectral energy distributions from BL Lac Objects (BL Lacs), which are marked by non-thermal beamed radiations, highly variable on timescales of days or less. Their outbursts in the gamma ray relative to the optical/X rays might be surmised to be enhanced in BL Lacs as these photons are upscattered via the inverse Compton (IC) process. From the observed correlations among the spectral parameters during optical/X-ray variations we aim at predicting corresponding correlations in the gamma-ray band, and the actual relations between the gamma-ray and the X-ray variability consistent with the SSC emission process. We start from the homogeneous single-zone SSC source model, with log-parabolic energies distributions of emitting electron as required by the X-ray data of many sources. We find relations among spectral parameters of the IC radiation in both the Thomson (for Low energy BL Lacs) and the Klein-Nishina regimes (mainly for High energy BL Lacs) and we compute how variability is driven by a smooth increase of key source parameters, primarily the root mean square electron energy. The single component SSC source model in the Thomson regime turns out to be adequate for many LBL sources. However, the simple model meets its limits with the fast/strong flares recently reported for a few sources in the TeV range; these require sudden accelerations of emitting electrons in a second source component.

preprint2009arXiv

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

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

preprint2008arXiv

AGILE detection of delayed gamma-ray emission from GRB 080514B

GRB 080514B is the first gamma ray burst (GRB), since the time of EGRET, for which individual photons of energy above several tens of MeV have been detected with a pair-conversion tracker telescope. This burst was discovered with the Italian AGILE gamma-ray satellite. The GRB was localized with a cooperation by AGILE and the interplanetary network (IPN). The gamma-ray imager (GRID) estimate of the position, obtained before the SuperAGILE-IPN localization, is found to be consistent with the burst position. The hard X-ray emission observed by SuperAGILE lasted about 7 s, while there is evidence that the emission above 30 MeV extends for a longer duration (at least ~13 s). Similar behavior was seen in the past from a few other GRBs observed with EGRET. However, the latter measurements were affected, during the brightest phases, by instrumental dead time effects, resulting in only lower limits to the burst intensity. Thanks to the small dead time of the AGILE/GRID we could assess that in the case of GRB 080514B the gamma-ray to X-ray flux ratio changes significantly between the prompt and extended emission phase.

preprint2008arXiv

AGILE observation of a gamma-ray flare from the blazar 3C 279

Context. We report the detection by the AGILE satellite of an intense gamma-ray flare from the gamma-ray source 3EG J1255-0549, associated to the Flat Spectrum Radio Quasar 3C 279, during the AGILE pointings towards the Virgo Region on 2007 July 9-13. Aims. The simultaneous optical, X-ray and gamma-ray covering allows us to study the spectral energy distribution (SED) and the theoretical models relative to the flaring episode of mid-July. Methods. AGILE observed the source during its Science Performance Verification Phase with its two co-aligned imagers: the Gamma- Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE) sensitive in the 30 MeV - 50 GeV and 18 - 60 keV respectively. During the AGILE observation the source was monitored simultaneously in optical band by the REM telescope and in the X-ray band by the Swift satellite through 4 ToO observations. Results. During 2007 July 9-13 July 2007, AGILE-GRID detected gamma-ray emission from 3C 279, with the source at ~2 deg from the center of the Field of View, with an average flux of (210+-38) 10^-8 ph cm^-2 s^-1 for energy above 100 MeV. No emission was detected by Super-AGILE, with a 3-sigma upper limit of 10 mCrab. During the observation lasted about 4 days no significative gamma-ray flux variation was observed. Conclusions. The Spectral Energy Distribution is modelled with a homogeneous one-zone Synchrotron Self Compton emission plus the contributions by external Compton scattering of direct disk radiation and, to a lesser extent, by external Compton scattering of photons from the Broad Line Region.

preprint2008arXiv

Gamma-ray burst detection with the AGILE mini-calorimeter

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

preprint2008arXiv

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

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

preprint2008arXiv

High-Resolution Timing Observations of Spin-Powered Pulsars with the AGILE Gamma-Ray Telescope

AGILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to high-energy astrophysics launched in 2007 April. Its 1 microsecond absolute time tagging capability coupled with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the CGRO/EGRET heritage. In this paper we present the first AGILE timing results on the known gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a large number of gamma-ray photons from these pulsars (about 10,000 pulsed counts for Vela) in only few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves.