We have conducted a detailed investigation of the broadband spectral properties of the γ-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi γ-ray spectra with Swift, radio, infra-red, optical, and other hard X-ray/γ-ray data, collected within 3 months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous spectral energy distributions (SED) for 48 LBAS blazars. The SED of these γ-ray sources is similar to that of blazars discovered at other ABDO ET AL. Vol. 716 wavelengths, clearly showing, in the usual log ν-log ν F ν representation, the typical broadband spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SED to characterize the peak intensity of both the low-and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broadband colors (i.e., the radio to optical, α ro , and optical to X-ray, α ox , spectral slopes) and from the γ-ray spectral index. Our data show that the synchrotron peak frequency (ν S peak) is positioned between 10 12.5 and 10 14.5 Hz in broad-lined flat spectrum radio quasars (FSRQs) and between 10 13 and 10 17 Hz in featureless BL Lacertae objects. We find that the γ-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron-inverse Compton scenarios. However, simple homogeneous, one-zone, synchrotron self-Compton (SSC) models cannot explain most of our SED, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. More complex models involving external Compton radiation or multiple SSC components are required to reproduce the overall SED and the observed spectral variability. While more than 50% of known radio bright high energy peaked (HBL) BL Lacs are detected in the LBAS sample, only less than 13% of known bright FSRQs and LBL BL Lacs are included. This suggests that the latter sources, as a class, may be much fainter γ-ray emitters than LBAS blazars, and could in fact radiate close to the expectations of simple SSC models. We categorized all our sources according to a new physical classification scheme based on the generally accepted paradigm for Active Galactic Nuclei and on the results of this SED study. Since the LAT detector is more sensitive to flat spectrum γ-ray sources, the correlation between ν S peak and γ-ray spectral index strongly favors the detection of high energy peaked blazars, thus explaining the Fermi overabundance of this type of sources compared to radio and EGRET samples. This selection effect is similar to that experienced in the soft X-ray band where HBL BL Lacs are the dominant type of blazars.

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope provides an unprecedented opportunity to study gamma-ray blazars. To capitalize on this opportunity, beginning in late 2007, about a year before the start of LAT science operations, we began a large-scale, fast-cadence 15 GHz radio monitoring program with the 40-m telescope at the Owens Valley Radio Observatory (OVRO). This program began with the 1158 northern (δ > −20 • ) sources from the Candidate Gammaray Blazar Survey (CGRaBS) and now encompasses over 1500 sources, each observed twice per week with about 4 mJy (minimum) and 3% (typical) uncertainty. Here, we describe this monitoring program and our methods, and present radio light curves from the first two years (2008 and 2009). As a first application, we combine these data with a novel measure of light curve variability amplitude, the intrinsic modulation index, through a likelihood analysis to examine the variability properties of subpopulations of our sample. We demonstrate that, with high significance (7-σ), gamma-ray-loud blazars detected by the LAT during its first 11 months of operation vary with about a factor of two greater amplitude than do the gamma-ray quiet blazars in our sample. We also find a significant (3-σ) difference between variability amplitude in BL Lacertae objects and flat-spectrum radio quasars (FSRQs), with the former exhibiting larger variability amplitudes. Finally, low-redshift (z < 1) FS-RQs are found to vary more strongly than high-redshift FSRQs, with 3-σ significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

We have conducted a multiwavelength survey of 42 radio loud narrow-1ine Seyfert 1 galaxies (RLNLS1s), selected by searching among all the known sources of this type and omitting those with steep radio spectra. We analyse data from radio frequencies to X-rays, and supplement these with information available from online catalogues and the literature in order to cover the full electromagnetic spectrum. This is the largest known multiwavelength survey for this type of source. We detected 90% of the sources in X-rays and found 17% at γ rays. Extreme variability at high energies was also found, down to timescales as short as hours. In some sources, dramatic spectral and flux changes suggest interplay between a relativistic jet and the accretion disk. The estimated masses of the central black holes are in the range ∼ 10 6−8 M ⊙ , lower than those of blazars, while the accretion luminosities span a range from ∼ 0.01 to ∼ 0.49 times the Eddington limit, with an outlier at 0.003, similar to those of quasars. The distribution of the calculated jet power spans a range from ∼ 10 42.6 to ∼ 10 45.6 erg s −1 , generally lower than quasars and BL Lac objects, but partially overlapping with the latter. Once normalised by the mass of the central black holes, the jet power of the three types of active galactic nuclei are consistent with each other, indicating that the jets are similar and the observational differences are due to scaling factors. Despite the observational differences, the central engine of RLNLS1s is apparently quite similar to that of blazars. The historical difficulties in finding radio-loud narrow-line Seyfert 1 galaxies might be due to their low power and to intermittent jet activity.

We report on the γ -ray activity of the high-synchrotron-peaked BL Lacertae object Markarian 421 (Mrk 421) during the first 1.5 years of Fermi operation, from 2008 August 5 to 2010 March 12. We find that the Large Area Telescope (LAT) γ -ray spectrum above 0.3 GeV can be well described by a power-law function with photon index Γ = 1.78 ± 0.02 and average photon flux F (>0.3 GeV) = (7.23 ± 0.16) × 10 −8 ph cm −2 s −1 . Over this time period, the Fermi-LAT spectrum above 0.3 GeV was evaluated on seven-day-long time intervals, showing significant variations in the photon flux (up to a factor ∼3 from the minimum to the maximum flux) but mild spectral variations. The variability amplitude at X-ray frequencies measured by RXTE/ASM and Swift/BAT is substantially larger than that in γ -rays measured by Fermi-LAT, and these two energy ranges are not significantly correlated. We also present the first results from the 4.5 month long multifrequency campaign on Mrk 421, which included the VLBA, Swift, RXTE, MAGIC, the F-GAMMA, GASP-WEBT, and other collaborations and instruments that provided excellent temporal and energy coverage of the source throughout the entire campaign During this campaign, Mrk 421 showed a low activity at all wavebands. The extensive multi-instrument (radio to TeV) data set provides an unprecedented, complete look at the quiescent spectral energy distribution (SED) for this source. The broadband SED was reproduced with a leptonic (one-zone synchrotron self-Compton) and a hadronic model (synchrotron proton blazar). Both frameworks are able to describe the average SED reasonably well, implying comparable jet powers but very different characteristics for the blazar emission site.

We report on the γ-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) γ-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 ± 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size 0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (10 44 erg s −1) constitutes only a small fraction (∼10 −3) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude.

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

The exact location of the γ-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and γ-ray (0.1-300 GeV) ∼ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and γ-ray correlations (radio lagging γ rays) when averaging over the whole sample, (ii) average time delays (source frame: 76 ± 23 to 7 ± 9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence τ r,γ (ν) ∝ ν −1 , in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk γ-ray production region typically located within/upstream of the 3 mm core region (τ 3mm,γ = 12 ± 8 days), (iv) mean distances between the region of γ-ray peak emission and the radio "τ = 1 photosphere" decreasing from 9.8 ± 3.0 pc (11 cm) to 0.9 ± 1.1 pc (2 mm) and 1.4 ± 0.8 pc (0.8 mm), (v) 3 mm/γ-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 × 10 −6 ), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk γ-ray production region from the supermassive black hole (SMBH) of ∼ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm τ = 1 surface at ∼ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the data set obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations during non-rotating periods have significantly larger amplitude and faster variations of polarization angle than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability ≤ 1.5 × 10 −2 ) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely (∼ 5 × 10 −5 ) that none of our rotations is physically connected with an increase in gamma-ray activity.