A Multivariate Fit Luminosity Function and World Model for Long GRBs

It is proposed that the luminosity function, the rest-frame spectral correlations and distributions of cosmological Long-duration (Type-II) Gamma-Ray Bursts (LGRBs) may be very well described as multivariate log-normal distribution. This result is based on careful selection, analysis and modeling of LGRBs’ temporal and spectral variables in the largest catalog of Gamma-Ray Bursts available to date: 2130 BATSE GRBs, while taking into account the detection threshold and possible selection effects. Constraints on the joint rest-frame distribution of the isotropic peak luminosity (Liso), total isotropic emission (Eiso), the time-integrated spectral peak energy (Epkz) and duration (T90z) of LGRBs are derived. The presented analysis provides evidence for a relatively large fraction of LGRBs that have been missed by BATSE detector with Eiso extending down to ~ 10^49 [erg] and observed spectral peak energies (Epk) as low as ~ 5 [keV]. LGRBs with rest-frame duration T90z < 1 [s] or observer-frame duration T90 < 2 [s] appear to be rare events (<0.1% chance of occurrence). The model predicts a fairly strong but highly significant correlation (rho=0.58 \pm 0.04) between Eiso & Epkz of LGRBs. Also predicted are strong correlations of Liso & Eiso with T90z and moderate correlation between Liso & Epkz. The strength and significance of the correlations found, encourage the search for underlying mechanisms, though undermine their capabilities as probes of Dark Energy’s equation of state at high redshifts. The presented analysis favors – but does not necessitate – a cosmic rate for BATSE LGRBs tracing metallicity evolution consistent with a cutoff ~ 0.2-0.5, assuming no luminosity-redshift evolution.

💡 Research Summary

This paper presents a comprehensive statistical analysis of long‑duration gamma‑ray bursts (LGRBs) using the full BATSE catalog of 2,130 events, of which 1,366 are classified as LGRBs after a careful multivariate clustering procedure. The authors propose that the rest‑frame isotropic peak luminosity (L_iso), total isotropic energy (E_iso), spectral peak energy corrected for redshift (E_pkz), and rest‑frame duration (T_90z) are jointly described by a four‑dimensional log‑normal distribution. By explicitly modeling the BATSE detection threshold as a function of both peak flux and spectral peak energy, they incorporate selection effects directly into a maximum‑likelihood (or Bayesian) framework, allowing unbiased inference of the underlying distribution parameters.

The clustering step employs fuzzy C‑means on the observed peak energy (E_p) and duration (T_90) to assign each burst a probability of belonging to the long‑duration class, mitigating the well‑known misclassification that arises when a simple T_90 cut is used. Visual inspection of borderline cases further refines the sample, yielding a highly pure LGRB set.

Fitting the log‑normal model yields a mean vector and covariance matrix that encode several strong, statistically significant correlations: the most notable is a Pearson coefficient ρ ≈ 0.58 ± 0.04 between E_iso and E_pkz, reminiscent of the Amati relation but derived from a much larger, bias‑corrected sample. L_iso and E_iso are even more tightly linked (ρ ≈ 0.70), and both correlate strongly with T_90z (ρ ≈ 0.55–0.60). A moderate correlation (ρ ≈ 0.30) exists between L_iso and E_pkz. These relationships suggest that the total radiated energy, the spectral hardness, and the burst duration are governed by common physical processes, likely related to the central engine activity time and the efficiency of energy conversion in the relativistic jet.

Importantly, the inferred distribution extends down to E_iso ≈ 10^49 erg and observed peak energies as low as ~5 keV, indicating that BATSE missed a substantial population of low‑luminosity, soft bursts. The authors argue that this hidden population is consistent with a cosmic LGRB rate that tracks metallicity evolution, featuring a cutoff at Z/Z_⊙ ≈ 0.2–0.5. No additional luminosity‑redshift evolution is required to reproduce the observed data, challenging models that invoke strong evolutionary effects.

The paper also discusses the implications for using LGRBs as cosmological probes. While the multivariate correlations are robust, the intrinsic scatter and residual selection biases limit the precision of LGRBs as standard candles for constraining dark energy equation‑of‑state parameters at high redshift. Nonetheless, the strong E_iso–E_pkz and L_iso–T_90z correlations provide valuable constraints on theoretical models of prompt emission, such as internal shock versus photospheric scenarios.

In summary, the study demonstrates that a simple multivariate log‑normal model, when combined with rigorous treatment of detector thresholds and a sophisticated classification scheme, can successfully reproduce the full set of BATSE LGRB observables. The results highlight the existence of a previously unrecognized low‑energy tail of the LGRB population, support a metallicity‑dependent cosmic rate, and provide a statistically solid foundation for future investigations with more recent instruments (Swift, Fermi, and upcoming missions) that will include redshift measurements and extend sensitivity to softer, fainter bursts.