[en] Fermi has provided the largest sample of γ-ray-selected blazars to date. In this work we use a complete sample of flat spectrum radio quasars (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 ∼0.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 γ-ray FSRQs follows a luminosity-dependent density evolution similar to that of radio-quiet active galactic nuclei. Also, using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution (SED) of FSRQs in the 10 keV-300 GeV band and show that there is no correlation between the luminosity at the peak of the γ-ray emission component and its peak frequency. Using this luminosity-independent SED with the derived LF allows us to predict that the contribution of FSRQs to the Fermi isotropic γ-ray background is 9.3^+1.6_–1.0% (±3% systematic uncertainty) in the 0.1-100 GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of γ = 11.7^+3.3_–2.2, that most are seen within 5° of the jet axis, and that they represent only ∼0.1% of the parent population.

[en] Whether or not 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: although it may fail to confirm an existing intrinsic correlation in a small sample that cannot be split into many redshift bins, it will 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.

[en] Since mid-2007, we have been monitoring ∼1200 sources at 15 GHz with the 40 M Telescope at the Owens Valley Radio Observatory. Our sample, mostly blazars, is monitored at least twice per week, yielding densely-sampled light curves. A large fraction of the sources in our sample exhibit significant variation in 15 GHz flux density, enabling variability studies and cross-correlations with other bands. Additionally, many have been detected by the Fermi Gamma-ray Space Telescope. We compare our data with gamma-ray data from Fermi and find a statistically significant flux density correlation after accounting for red shift and selection biases using a new Monte Carlo method. The OVRO program is a part of the F-GAMMA project, which also obtains monthly 2.6-270 GHz radio spectra for a smaller, overlapping blazar sample.

[en] Fermi has provided the largest sample of γ-ray-selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by Fermi during its first year of operation. We obtained redshift constraints for 206 out of the 211 BL Lac objects in our sample, making it the largest and most complete sample of BL Lac objects available in the literature. We use this sample to determine the luminosity function of BL Lac objects 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). Low-luminosity, high-synchrotron-peaked (HSP) BL Lac objects are an exception, showing strong negative evolution, with number density increasing for z ≲ 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, by estimating the beaming corrections to the luminosity function, we find that BL Lac objects have an average Lorentz factor of γ=6.1_−0.8^+1.1, and that most are seen within 10° of the jet axis.

[en] The Fermi Large Area Telescope (LAT) discovered a new gamma-ray source near the Galactic plane, Fermi J0109+6134, when it flared brightly in 2010 February. The low Galactic latitude (b = -1.⁰2) 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 Swift and the Chandra X-ray Observatory. We determined the source redshift, z = 0.783, using a Keck Low-Resolution Imaging Spectrometer 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 A, this blazar belongs in the flat-spectrum radio quasar category.