Abstract

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 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 approximate to 1.2). Low-luminosity, high-synchrotron-peaked (HSP) BL Lac objects are an exception, showing strong negative evolution, with number density increasing for z less than or similar to 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 gamma = 6.1(-0.8)(+1.1), and that most are seen within 10 degrees of the jet axis.