2004ApJ...611.1033F

We constrain the isotropic luminosity function (LF) and formation rate of long γ-ray bursts (GRBs) by fitting models jointly to both the observed differential peak-flux and redshift distributions. We find evidence supporting an evolving LF, where the luminosity scales as (1+z)^δ, with an optimal δ of 1.0±0.2. For a single-power law LF, the best slope is approximately -1.57 with an upper luminosity of 10^53.3 ergs/s, while the best slopes for a double-power law LF are approximately -1.6 and -2.6, with a break luminosity of 10^52.7 ergs/s. Our finding implies a jet model intermediate between the universal structured ε(θ)∝θ^–2 model and the quasi-universal Gaussian structured model. For the uniform-jet model our result is compatible with an angle distribution between 2° and 15°. Our best-constrained GRB formation rate histories increase from z=0 to 2 by a factor of ∼30 and then continue increasing slightly. We connect these histories to the cosmic star formation history and compare with observational inferences up to z∼6. GRBs could be tracing the cosmic rates of both the normal and obscured star formation regimes. We estimate a current GRB event rate in the Milky Way of ∼5x10^–5/yr and compare it with the birthrate of massive close Wolf-Rayet + black hole binaries with orbital periods of hours. The agreement is rather good, suggesting that these systems could be the progenitors of the long GRBs.