The origin of the cosmic gamma-ray background (CGB) is a longstanding mystery in high-energy astrophysics. Possible candidates include ordinary astrophysical objects such as unresolved blazars, as well as more exotic processes such as dark matter annihilation. While it would be difficult to distinguish them from the mean intensity data alone, one can use anisotropy data instead. We investigate the CGB anisotropy both from unresolved blazars and dark matter annihilation (including contributions from dark matter substructures), and find that the angular power spectra from these sources are very different. We then focus on detectability of dark matter annihilation signals using the anisotropy data, by treating the unresolved blazar component as a known background. We find that the dark matter signature should be detectable in the angular power spectrum of the CGB from two-year all-sky observations with the Gamma Ray Large Area Space Telescope (GLAST), as long as the dark matter annihilation contributes to a reasonable fraction, e.g., >~ 0.3, of the CGB at around 10 GeV. We conclude that the anisotropy measurement of the CGB with GLAST should be a powerful tool for revealing the CGB origin, and potentially for the first detection of dark matter annihilation.

Blazars, a beamed population of active galactic nuclei, radiate high-energy gamma-rays, and thus are a good target for the Gamma Ray Large Area Space Telescope (GLAST). As the blazars trace the large-scale structure of the universe, one may observe spatial clustering of blazars. We calculate the angular power spectrum of blazars that would be detected by GLAST. We show that we have the best chance of detecting their clustering at large angular scales, \theta >~ 10 deg, where shot noise is less important, and the dominant contribution to the correlation comes from relatively low redshift, z <~ 0.1. The GLAST can detect the correlation signal, if the blazars detected by GLAST trace the distribution of low-z quasars observed by optical galaxy surveys, which have the bias of unity. If the bias of blazars is greater than 1.5, GLAST will detect the correlation signal unambiguously. We also find that GLAST may detect spatial clustering of clusters of galaxies in gamma-rays. The shape of the angular power spectrum is different for blazars and clusters of galaxies; thus, we can separate these two contributions on the basis of the shape of the power spectrum.

We present a comprehensive study for the gamma-ray luminosity function (GLF) of blazars and their contribution to the extragalactic diffuse gamma-ray background (EGRB). Radio and gamma-ray luminosity correlation is introduced to take into account the radio detectability which is important for the blazar identification. Previous studies considered only pure luminosity evolution (PLE) or pure density evolution, but we introduce the luminosity dependent density evolution (LDDE) model, which is favored from the evolution of X-ray luminosity function (XLF) of AGNs. The model parameters are constrained by likelihood analyses about the observed redshift and gamma-ray flux distributions of the EGRET blazars. We find that the LDDE model gives a better fit to the observed distributions than the PLE model, indicating that the LDDE model is also appropriate for gamma-ray blazars, and that the jet activity is universally correlated with the accretion history of AGNs. The normalization between the GLF and XLF is consistent with the unified picture of AGNs, when the beaming and a reasonable duty cycle of jet activity are taken into account. We then find that only 25--50% of the EGRB can be explained by unresolved blazars with the best-fit LDDE parameters. Unresolved blazars can account for all the EGRB only with a steeper index of the faint-end slope of the GLF, which is marginally consistent with the EGRET data but inconsistent with that of the XLF. Therefore unresolved AGNs cannot be the dominant source of the EGRB, unless there is a new population of gamma-ray emitting AGNs that evolves differently from the XLF of AGNs. Predictions for the GLAST mission are made, and we find that the best-fit LDDE model predicts about 3000 blazars in the entire sky, which is considerably fewer than a previous estimate.