Possible Origins of Dispersion of the Peak Energy--Brightness Correlations of Gamma-Ray Bursts

We collect and reanalyze about 200 GRB data of prompt-emission with known redshift observed until the end of 2009, and select 101 GRBs which were well observed to have good spectral parameters to determine the spectral peak energy ($E_p$), 1-second peak luminosity ($L_p$) and isotropic energy ($E_{\rm iso}$). Using our newly-constructed database with 101 GRBs, we first revise the $E_p$–$L_p$ and $E_p$–$E_{\rm iso}$ correlations. The correlation coefficients of the revised correlations are 0.889 for 99 degree of freedom for the $E_p$–$L_p$ correlation and 0.867 for 96 degree of freedom for the $E_p$–$E_{\rm iso}$ correlation. These values correspond to the chance probability of $2.18 \times 10^{-35}$ and $4.27 \times 10^{-31}$, respectively. It is a very important issue whether these tight correlations are intrinsic property of GRBs or caused by some selection effect of observations. In this paper, we examine how the truncation of the detector sensitivity affects the correlations, and we conclude they are surely intrinsic properties of GRBs. Next we investigate origins of the dispersion of the correlations by studying their brightness and redshift dependence. Here the brightness (flux or fluence) dependence would be regarded as an estimator of the bias due to the detector threshold. We find a weak fluence-dependence in the $E_p$–$E_{\rm iso}$ correlations and a redshift dependence in the $E_p$–$L_p$ correlation both with 2 $\sigma$ statistical level. These two effects may contribute to the dispersion of the correlations which is larger than the statistical uncertainty. We discuss a possible reason of these dependence and give a future prospect to improve the correlations.

💡 Research Summary

The authors set out to determine whether the well‑known empirical relations linking the rest‑frame spectral peak energy (Eₚ) of gamma‑ray bursts (GRBs) to their brightness indicators—1‑second peak luminosity (Lₚ) and isotropic equivalent energy (E_iso)—are intrinsic physical properties or artifacts of observational selection. They assembled a homogeneous database of 101 GRBs with measured redshifts and reliable spectral parameters, drawn from roughly 200 events observed up to the end of 2009 by a variety of missions (BATSE, Swift/BAT, Konus‑Wind, Fermi/GBM, etc.). For each burst they fitted the prompt spectrum with the Band function (or a fixed high‑energy index β = –2.25 when β could not be constrained), derived the rest‑frame peak energy Eₚ = Eₚ,obs (1+z), and performed k‑corrections to compute bolometric fluence (S_bol) and bolometric 1‑s peak flux (F_p,bol) in the 1–10 000 keV rest‑frame band. Using a standard ΛCDM cosmology (Ω_m = 0.3, Ω_Λ = 0.7, H₀ = 70 km s⁻¹ Mpc⁻¹) they obtained E_iso = 4πd_L² S_bol/(1+z) and Lₚ = 4πd_L² F_p,bol.

Linear regression in log–log space yields very tight correlations: r = 0.889 (99 d.o.f.) for Eₚ–Lₚ and r = 0.867 (96 d.o.f.) for Eₚ–E_iso, corresponding to chance probabilities of 2 × 10⁻³⁵ and 4 × 10⁻³¹, respectively. These results confirm that the correlations persist even with a larger, uniformly processed sample.

The central question, however, is why the scatter around the best‑fit lines exceeds the statistical uncertainties. The authors explore two possible contributors:

1. Detector‑threshold (truncation) effects – Low‑flux or low‑fluence bursts lie close to the instrumental sensitivity limit, making spectral parameter estimation less reliable. By dividing the sample into bright and dim subsets based on observed flux/fluence, they find a modest fluence dependence in the Eₚ–E_iso relation: the dimmer half shows a slightly shallower slope, significant at the ~2σ level. This suggests that the spectral analysis threshold (rather than the trigger threshold) can introduce systematic offsets for faint events.

2. Redshift evolution – The sample is split into three redshift bins (z < 1, 1 ≤ z < 2, z ≥ 2). The Eₚ–Lₚ correlation exhibits a noticeable change in slope for the highest‑z bin, again at ~2σ significance. This may reflect genuine cosmological evolution of the GRB emission mechanism or selection biases that become more pronounced at large distances.

To assess whether these effects could fully account for the observed correlations, the authors perform Monte‑Carlo simulations that impose realistic flux/fluence thresholds on an underlying GRB population. The resulting “truncated” datasets reproduce the observed scatter but do not substantially alter the correlation slopes or intercepts, indicating that the intrinsic Eₚ–Lₚ and Eₚ–E_iso relations are robust against the examined selection biases.

Methodologically, the paper pays careful attention to cross‑instrument consistency: peak fluxes measured on different time scales (64 ms, 256 ms, 1024 ms) are rescaled to a common 1‑s interval using empirical correction factors; when necessary, the authors re‑analyze raw light curves (e.g., Konus/Wind) to obtain uniform peak fluxes. They also discuss the impact of fixing β, the treatment of cutoff‑power‑law spectra, and the choice of the 1–10 000 keV bolometric band.

In conclusion, the study provides strong evidence that the Eₚ–Lₚ and Eₚ–E_iso correlations are genuine, intrinsic properties of GRBs rather than mere artifacts of detector thresholds. Nonetheless, the residual dispersion is partially attributable to (i) a weak fluence dependence (likely reflecting spectral analysis limits for faint bursts) and (ii) a modest redshift dependence (potentially hinting at evolutionary effects). The authors recommend future work to (a) expand the sample with more broadband (keV–GeV) observations to better constrain β, (b) develop rigorous selection‑function models to construct truly “threshold‑free” samples, and (c) explore physical models (e.g., internal shock efficiency, jet composition) that could naturally produce the observed scaling laws. By reducing systematic scatter, GRBs could become more reliable standardizable candles for cosmology, extending distance measurements well beyond the reach of Type Ia supernovae.