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.

Fermi has provided the largest sample of gamma-ray selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by it Fermi during its first year of operation. We have obtained redshift constraints for 206 out of the 211 BL Lacs in our sample making it the largest and most complete sample of BL Lacs available in the literature. We use this sample to determine the luminosity function of BL Lacs and its evolution with cosmic time. We find that for most BL Lac classes, the evolution is positive with a space density peaking at modest redshift (z~1.2). The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception, showing strong negative evolution, with number density increasing for z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of flat-spectrum radio quasars (FSRQs), a possible interpretation is that these HSPs represent an accretion-starved end-state of an earlier merger-driven gas-rich phase. We additionally find that the known BL Lac correlation between luminosity and photon spectral index persists after correction for the substantial observational selection effects with implications for the so called `blazar sequence'. Finally, estimating the beaming corrections to the luminosity function, we find that BL Lacs have an average Lorentz factor of $\gamma=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the jet axis.

We report on spectroscopic observations covering most of the 475 BL Lacs in the 2nd Fermi LAT catalog of AGN. Including archival measurements (correcting several erroneous literature values) we now have spectroscopic redshifts for 44% of the BL Lacs. We establish firm lower redshift limits via intervening absorption systems and statistical lower limits via searches for host galaxies for an additional 51% of the sample leaving only 5% of the BL Lacs unconstrained. The new redshifts raise the median spectroscopic z from 0.23 to 0.33 and include redshifts as large as z=2.471. Spectroscopic redshift minima from intervening absorbers have ~ z= 0.70, showing a substantial fraction at large z and arguing against strong negative evolution. We find that detected BL Lac hosts are bright ellipticals with black hole masses M_\bullet ~ 10^{8.5-9}, substantially larger than the mean of optical AGN and LAT Flat Spectrum Radio Quasar samples. A slow increase in M_\bullet with z may be due to selection bias. We find that the power-law dominance of the optical spectrum extends to extreme values, but this does not strongly correlate with the gamma-ray properties, suggesting that strong beaming is the primary cause of the range in continuum dominance.

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

We report on optical spectroscopy of 165 Flat Spectrum Radio Quasars (FSRQs) in the Fermi 1LAC sample, which have helped allow a nearly complete study of this population. Fermi FSRQ show significant evidence for non-thermal emission even in the optical; the degree depends on the gamma-ray hardness. They also have smaller virial estimates of hole mass than the optical quasar sample. This appears to be largely due to a preferred (axial) view of the gamma-ray FSRQ and non-isotropic (H/R ~ 0.4) distribution of broad-line velocities. Even after correction for this bias, the Fermi FSRQ show higher mean Eddington ratios than the optical population. A comparison of optical spectral properties with Owens Valley Radio Observatory radio flare activity shows no strong correlation.

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.

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 (declination>-20 deg) sources from the Candidate Gamma-ray Blazar Survey (CGRaBS) and now encompasses over 1500 sources, each observed twice per week with a ~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-sigma), 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-sigma) 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) FSRQs are found to vary more strongly than high-redshift FSRQs, with 3-sigma significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

The radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). This large, flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux. We find that the LAT-detected BL Lacs tend to be similar to the non-LAT BL Lacs, but that the LAT-detected FSRQs are often significantly different from the non-LAT FSRQs. The differences between the gamma-ray loud and quiet FSRQs can be explained by Doppler boosting; these objects appear to require larger Doppler factors than those of the BL Lacs. It is possible that the gamma-ray loud FSRQs are fundamentally different from the gamma-ray quiet FSRQs. Strong polarization at the base of the jet appears to be a signature for gamma-ray loud AGN.

X-shaped radio galaxies are defined by their peculiar large-scale radio morphology. In addition to the classical double-lobed structure they have a pair of low-luminosity wings that straddles the nucleus at almost right angles to the active lobes, thus giving the impression of an 'X'. In this paper we study for the first time the optical spectral properties of this object class using a large sample (~50 sources). We find that the X-shaped radio population is composed roughly equally of sources with weak and strong emission line spectra, which makes them, in combination with the well-known fact that they preferentially have radio powers intermediate between those of Fanaroff-Riley type I (FR I) and type II (FR II) radio galaxies, the archetypal transition population. We do not find evidence in support of the proposition that the X-shape is the result of a recent merger: X-shaped radio sources do not have unusually broad emission lines, their nuclear environments are in general not dusty, and their host galaxies do not show signs of enhanced star formation. Instead, we observe that the nuclear regions of X-shaped radio sources have relatively high temperatures. This finding favours models, which propose that the X-shape is the result of an overpressured environment.

The \emph{Fermi} Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic plane, \object{Fermi J0109+6134}, when it flared brightly in 2010 February. The low Galactic latitude (b =-1.2\degr) indicated that the source could be located within the Galaxy, which motivated rapid multi-wavelength follow-up including radio, optical, and X-ray observations. We report the results of analyzing all 19 months of LAT data for the source, and of X-ray observations with both \emph{Swift} and the \emph{Chandra X-ray Observatory}. We determined the source redshift, z =0.783, using a Keck LRIS observation. Finally, we compiled a broadband spectral energy distribution (SED) from both historical and new observations contemporaneous with the 2010 February flare. The redshift, SED, optical line width, X-ray absorption, and multi-band variability indicate that this new GeV source is a blazar seen through the Galactic plane. Because several of the optical emission lines have equivalent width >5\AA, this blazar belongs in the flat-spectrum radio quasar category.

We present multifrequency Very Long Baseline Array (VLBA) observations of J16021+3326. These observations, along with variability data obtained from the Owens Valley Radio Observatory (OVRO) candidate gamma-ray blazar monitoring program, clearly indicate this source is a blazar. The peculiar characteristic of this blazar, which daunted previous classification attempts, is that we appear to be observing down a precessing jet, the mean orientation of which is aligned with us almost exactly.

We present results from the first two years of our fast-cadence 15 GHz gamma-ray blazar monitoring program, part of the F-GAMMA radio monitoring project. Our sample includes the 1158 blazars north of -20 degrees declination from the Candidate Gamma-Ray Blazar Survey (CGRaBS), which encompasses a significant fraction of the extragalactic sources detected by the Fermi Gamma-ray Space Telescope. We introduce a novel likelihood analysis for computing a time series variability amplitude statistic that separates intrinsic variability from measurement noise and produces a quantitative error estimate. We use this method to characterize our radio light curves. We also present results indicating a statistically significant correlation between simultaneous average 15 GHz radio flux density and gamma-ray photon flux.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three 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 gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band 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 SEDs 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 broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-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 SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

We report on HET and Palomar 5 m spectroscopy of recently identified $\gamma$-ray blazars in the {\it Fermi} LAT Bright Source List. These data provide identifications for 10 newly discovered $\gamma$-ray flat spectrum radio quasars (FSRQ) and six new BL Lacs plus improved spectroscopy for six additional BL Lacs. We substantially improve the identification completeness of the bright LAT blazars and give new redshifts and $z$ constraints, new estimates of the black hole masses and new measurements of the optical SED.

We describe a 4.85 GHz survey of bright, flat-spectrum radio sources conducted with the Effelsberg 100 m telescope in an attempt to improve the completeness of existing surveys, such as CRATES. We report the results of these observations and of follow-up 8.4 GHz observations with the VLA of a subset of the sample. We comment on the connection to the WMAP point source catalog and on the survey's effectiveness at supplementing the CRATES sky coverage.

We report optical spectroscopic observations of X-shaped radio sources with the Hobby-Eberly Telescope and Multiple-Mirror Telescope, focused on the sample of candidates from the FIRST survey presented in Paper I (Cheung 2007). A total of 27 redshifts were successfully obtained, 21 of which are new, including that of a newly identified candidate source of this type which is presented here. With these observations, the sample of candidates from Paper I is over 50% spectroscopically identified. Two new broad emission-lined X-shaped radio sources are revealed, while no emission lines were detected in about one third of the observed sources; a detailed study of the line properties is deferred to a future paper. Finally, to explore their relation to the Fanaroff-Riley division, the radio luminosities and host galaxy absolute magnitudes of a spectroscopically identified sample of 50 X-shaped radio galaxies are calculated to determine their placement in the Owen-Ledlow plane.

We describe a uniform all-sky survey of bright blazars, selected primarily by their flat radio spectra, that is designed to provide a large catalog of likely gamma-ray AGN. The defined sample has 1625 targets with radio and X-ray properties similar to those of the EGRET blazars, spread uniformly across the |b| > 10 deg sky. We also report progress toward optical characterization of the sample; of objects with known R < 23, 85% have been classified and 81% have measured redshifts. One goal of this program is to focus attention on the most interesting (e.g., high redshift, high luminosity, ...) sources for intensive multiwavelength study during the observations by the Large Area Telescope (LAT) on GLAST.

We examine the radio properties of EGRET-detected blazars observed as part of the VLBA Imaging and Polarimetry Survey (VIPS). VIPS has a flux limit roughly an order of magnitude below the MOJAVE survey and most other samples that have been used to study the properties of EGRET blazars. At lower flux levels, radio flux density does not directly correlate with gamma-ray flux density. We do find that the EGRET-detected blazars tend to have higher brightness temperatures, greater core fractions, and possibly larger than average jet opening angles. A weak correlation is also found with jet length and with polarization. All of the well-established trends can be explained by systematically larger Doppler factors in the gamma-ray loud blazars, consistent with the measurements of higher apparent velocities found in monitoring programs carried out at radio frequencies above 10 GHz.

We have assembled an 8.4 GHz survey of bright, flat-spectrum (alpha > -0.5) radio sources with nearly uniform extragalactic (|b|>10 deg) coverage for sources brighter than S_{4.8 GHz} = 65 mJy. The catalog is assembled from existing observations (especially CLASS and the Wright et al. PMN-CA survey), augmented by reprocessing of archival VLA and ATCA data and by new observations to fill in coverage gaps. We refer to this program as CRATES, the Combined Radio All-sky Targeted Eight GHz Survey. The resulting catalog provides precise positions, sub-arcsecond structures, and spectral indices for some 11,000 sources. We describe the morphology and spectral index distribution of the sample and comment on the survey's power to select several classes of interesting sources, especially high energy blazars. Comparison of CRATES with other high-frequency surveys also provides unique opportunities for identification of high-power radio sources.

We present the first results of the VLBA Imaging and Polarimetry Survey (VIPS), a 5 GHz VLBI survey of 1,127 sources with flat radio spectra. Through automated data reduction and imaging routines, we have produced publicly available I, Q, and U images and have detected polarized flux density from 37% of the sources. We have also developed an algorithm to use each source's I image to automatically classify it as a point-like source, a core-jet, a compact symmetric object (CSO) candidate, or a complex source. The mean ratio of the polarized to total 5 GHz flux density for VIPS sources with detected polarized flux density ranges from 1% to 20% with a median value of about 5%. We have also found significant evidence that the directions of the jets in core-jet systems tend to be perpendicular to the electric vector position angles (EVPAs). The data is consistent with a scenario in which ~24% of the polarized core-jets have EVPAs that are anti-aligned with the directions of their jet components and which have a substantial amount of Faraday rotation. In addition to these initial results, plans for future follow-up observations are discussed.

In preparation for GLAST, we have compiled a sample of blazar candidates to increase the pool of well studied AGN from which GLAST counterparts will be drawn. Sources were selected with our Figure of Merit (FoM) ranking; thus, they have radio and X-ray properties very similar to the EGRET blazars. Spectroscopic confirmation of these candidates is in progress, and more than 70% of these objects have been identified as flat spectrum radio quasars and BL Lac objects. We present ~250 new optical blazar identifications based on McDonald Observatory spectroscopy, 224 with redshifts. Of these, 167 are in our FoM-selected set. To motivate the Gamma-ray nature of these objects, we analyzed the current release of the EGRET data for possible point sources at their radio positions. We develop two distinct methods to combine multiple EGRET observations of a sky position into a single detection significance. We report a detection of the signal of the set of blazar candidates in the EGRET data at the > 3 sigma level by both techniques. We predict that the majority of these blazar candidates will be found by GLAST due to its increased sensitivity, duty cycle and resolving power.