Self-Consistent Model of Extragalactic Neutrino Flux from Evolving Blazar Population

We study constraints on the population of neutrino emitting blazars imposed by the absence of doublets in astrophysical muon neutrino signal and z similar or equal to 0.3 redshift of nearest identified neutrino-emitting blazar (an order of magnitude further away than the nearest gamma-ray emitting blazar). We show that in spite of the absence of correlation of neutrino arrival directions with positions of gamma-ray emitting blazars, cumulative blazar flux could explain most of astrophysical neutrino flux measured in muon neutrino channel. This is possible if the population of neutrino emitting blazars has experienced rapid positive evolution at least as (1 + z)(5) at z less than or similar to 1. Such a model avoids previously derived constraint on the low level of blazar contribution to extragalactic neutrino flux because gamma-ray and neutrino fluxes are dominated by different sets of blazars. Rapid evolution of neutrino emitting blazars could be explained by the fact that only rapidly evolving sub-population of blazars, which can include parts of Flat Spectrum Radio Quasar and bright BL Lac populations, are efficient neutrino sources, although their neutrino luminosity has to be systematically lower than the gamma-ray luminosity.