Statistics of cosmological black hole jet sources:: Blazar predictions for the Gamma-Ray Large Area Space Telescope

A study of the statistics of cosmological black hole jet sources is applied to EGRET blazar data and predictions are made for GLAST. Black hole jet sources are modeled as collimated relativistic plasma outflows with radiation beamed along the jet axis due to strong Doppler boosting. The comoving rate density of blazar flares is assumed to follow a blazar formation rate (BFR), modeled by analytic functions based on astronomical observations and fits to EGRET data. The redshift and size distributions of gamma-ray blazars observed with EGRET, separated into BL Lac objects (BLs) and flat spectrum radio quasar (FSRQ) distributions, are fit with monoparametric functions for the distributions of the jet Lorentz factor Gamma, comoving directional power l 0 e, and spectral slope. A BFR factor approximate to 10; greater at z greater than or similar to 1 than at present is found to fit the FSRQ data. A smaller comoving rate density and greater luminosity of BL flares at early times compared to the present epoch fits the BL data. Based on the EGRET observations, approximate to 1000 blazars consisting of approximate to 800 FSRQs and FR 2 radio galaxies and approximate to 200 BL Lac objects and FR 1 radio galaxies will be detected with GLAST during the first year of the mission. Additional AGN classes, such as hard-spectrum BL Lac objects that were mostly missed with EGRET, could add more GLAST sources. The FSRQ and BL contributions to the EGRET gamma-ray background at 1 GeVare estimated at the level of approximate to 10%-15% and approximate to 2%-4%, respectively. EGRET and GLAST sensitivities to blazar flares are considered in the optimal case, and a GLAST analysis method for blazar detection is outlined.