• ### High-energy gamma-ray observations of the accreting black hole V404 Cygni during its June 2015 outburst(1607.06239)

July 21, 2016 astro-ph.HE
We report on Fermi/Large Area Telescope observations of the accreting black hole low-mass X-ray binary V404 Cygni during its outburst in June-July 2015. Detailed analyses reveal a possible excess of $\gamma$-ray emission on 26 June 2015, with a very soft spectrum above $100$ MeV, at a position consistent with the direction of V404 Cyg (within the $95\%$ confidence region and a chance probability of $4 \times 10^{-4}$). This emission cannot be associated with any previously-known Fermi source. Its temporal coincidence with the brightest radio and hard X-ray flare in the lightcurve of V404 Cyg, at the end of the main active phase of its outburst, strengthens the association with V404 Cyg. If the $\gamma$-ray emission is associated with V404 Cyg, the simultaneous detection of $511\,$keV annihilation emission by INTEGRAL requires that the high-energy $\gamma$ rays originate away from the corona, possibly in a Blandford-Znajek jet. The data give support to models involving a magnetically-arrested disk where a bright $\gamma$-ray jet can re-form after the occurrence of a major transient ejection seen in the radio.
• ### Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113(1509.02063)

Oct. 12, 2015 astro-ph.HE
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.
• ### Rapid TeV Gamma-Ray Flaring of BL Lacertae(1211.3073)

Nov. 13, 2012 astro-ph.HE
We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.
• ### The Spitzer Extragalactic Representative Volume Survey (SERVS): survey definition and goals(1206.4060)

Sept. 17, 2012 astro-ph.CO
We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.
• 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 {\gamma} -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 {\gamma} -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.
• ### The structure and emission model of the relativistic jet in the quasar 3C 279 inferred from radio to high-energy gamma-ray observations in 2008-2010(1206.0745)

June 4, 2012 astro-ph.CO, astro-ph.HE
We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.
• ### Simultaneous Planck, Swift, and Fermi observations of X-ray and gamma-ray selected blazars(1108.1114)

May 23, 2012 astro-ph.CO
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 <\alpha>~-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.
• ### Long-term monitoring of the high-energy gamma-ray emission from LS I +61{\deg} 303 and LS 5039(1202.1866)

Feb. 9, 2012 astro-ph.GA, astro-ph.HE
The Fermi Large Area Telescope (LAT) reported the first definitive GeV detections of the binaries LS I +61\degree 303 and LS 5039 in the first year after its launch in June, 2008. These detections were unambiguous as a consequence of the reduced positional uncertainty and the detection of modulated gamma-ray emission on the corresponding orbital periods. An analysis of new data from the LAT, comprising 30 months of observations, identifies a change in the gamma-ray behavior of LS I +61\degree 303. An increase in flux is detected in March 2009 and a steady decline in the orbital flux modulation is observed. Significant emission up to 30GeV is detected by the LAT; prior datasets led to upper limits only. Contemporaneous TeV observations no longer detected the source, or found it -in one orbit- close to periastron, far from the phases at which the source previously appeared at TeV energies. The detailed numerical simulations and models that exist within the literature do not predict or explain many of these features now observed at GeV and TeV energies. New ideas and models are needed to fully explain and understand this behavior. A detailed phase-resolved analysis of the spectral characterization of LS I +61\degree 303 in the GeV regime ascribes a power law with an exponential cutoff spectrum along each analyzed portion of the system's orbit. The on-source exposure of LS 5039 is also substantially increased with respect to our prior publication. In this case, whereas the general gamma-ray properties remain consistent, the increased statistics of the current dataset allows for a deeper investigation of its orbital and spectral evolution.
• Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.
• ### The Luminosity Function of Fermi-detected Flat-Spectrum Radio Quasars(1110.3787)

Oct. 17, 2011 astro-ph.CO, astro-ph.HE
Fermi has provided the largest sample of {\gamma}-ray selected blazars to date. In this work we use a complete sample of FSRQs detected during the first year of operation to determine the luminosity function (LF) and its evolution with cosmic time. The number density of FSRQs grows dramatically up to redshift \sim0.5-2.0 and declines thereafter. The redshift of the peak in the density is luminosity dependent, with more luminous sources peaking at earlier times; thus the LF of {\gamma}-ray FSRQs follows a luminosity-dependent density evolution similarly to that of radio-quiet AGN. Also using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution of FSRQs in the 10 keV-100 GeV band and show that there is no correlation of the peak {\gamma}-ray luminosity with {\gamma}-ray peak frequency. The coupling of the SED and LF allows us to predict that the contribution of FSRQs to the Fermi isotropic {\gamma}-ray background is 9.3(+1.6/-1.0) (\pm3% systematic uncertainty) in the 0.1-100GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of {\gamma} = 11.7(+3.3/-2.2), that most are seen within 5\circ of the jet axis, and that they represent only ~0.1 % of the parent population.
• ### Gamma-Ray and Parsec-Scale Jet Properties of a Complete Sample of Blazars From the MOJAVE Program(1107.4977)

Aug. 31, 2011 astro-ph.CO, astro-ph.HE
We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.
• ### The radio/gamma-ray connection in Active Galactic Nuclei in the era of the Fermi Large Area Telescope(1108.0501)

Aug. 2, 2011 astro-ph.CO, astro-ph.HE
We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.