Empirical Iso-Energy Correlation in GRBs

• The paper presents a detailed empirical calibration of the Amati correlation, linking the GRB rest-frame peak photon energy to the isotropic-equivalent radiated energy.
• A combination of regression analysis, maximum-likelihood estimation, and Monte Carlo simulations yields robust parameters (e.g., slope m ~0.5 down to 0.25) with quantified intrinsic scatter.
• Findings indicate that selection effects and viewing-angle dependencies significantly shape the correlation, enhancing its role in probing jet physics and cosmological applications.

Empirical Iso-Energy Correlation in Gamma-Ray Bursts

The empirical iso-energy correlation in gamma-ray bursts (GRBs), most famously expressed as the Amati relation—linking the rest-frame peak photon energy ($E_{p,i}$) of the prompt emission spectrum to the isotropic-equivalent radiated energy ($E_{\rm iso}$)—has emerged as a central tool for both theoretical modeling and cosmological application of GRBs. This correlation, discovered empirically in long-duration GRBs, provides a quantitative bridge between microphysical parameters of relativistic jet emission processes and the global energetics of explosive stellar phenomena observable up to redshifts $z > 9$.

1. Mathematical Formulation and Calibration

The iso-energy correlation is generally formulated as a power law in rest-frame quantities,

$$E_{p,i} = K\ \left(\frac{E_{\rm iso}}{10^{52}~{\rm erg}}\right)^{m}$$

where $E_{p,i}$ is in keV, $K$ is the normalization, and $m$ the slope. Both "forward" ($E_{p,i}$ as function of $E_{\rm iso}$) and "inverse" (the logarithmically inverted) forms are used depending on regression methodology.

Calibrated best-fit parameters vary by dataset and methodology but cluster near:

• $K \simeq 100~{\rm keV}$;
• $m = 0.5 \pm 0.05$ (canonical sample, e.g. Amati et al. 2008; Ghirlanda et al. 2010 (Amati, 2010, Ghirlanda, 2010, Demianski et al., 2016)).

The correlation is log-linear and typically incorporates an intrinsic scatter term (e.g., $\sigma_{\log E_{p,i}} \sim 0.2$–$0.4$ dex), estimated via maximum-likelihood or D’Agostini-style methods accounting for heteroscedastic measurement errors and astrophysical/extrinsic dispersion (Demianski et al., 2011, Demianski et al., 2016).

Recent large-sample recalibrations using Fermi-GBM with joint GBM+LAT-LLE spectral analysis yield slopes in the flattened regime $m = 0.25$–$0.30$, particularly when selection effects and instrumental thresholds are rigorously addressed (Aldowma et al., 18 Oct 2025, Shahmoradi, 2013).

2. Physical Interpretation and Theoretical Origin

The empirical power-law form of the Amati relation finds natural theoretical grounding within standard prompt-emission models involving internal shocks or magnetic-dissipation in ultra-relativistic outflows. In these scenarios, the observed $E_{p,i}$ is linked to the characteristic electron Lorentz factor and comoving magnetic field, while $E_{\rm iso}$ traces the total radiated energy via bulk kinetic dissipation.

Under standard assumptions: $E_{p,i} \propto \gamma B' \gamma_e^2$

$$E_{\rm iso} \propto B'^2 \gamma^4 m_{\rm sh} \gamma_e^2$$

Eliminating unknowns recovers $E_{p,i}\propto E_{\rm iso}^{0.5}$ for on-axis observers, an analytic result confirmed by Monte Carlo population synthesis (Xu et al., 2022).

For off-axis bursts or those with significant viewing-angle effects, the slope flattens to $k\sim0.25$–$0.3$, providing a theoretical rationale for the observed locus of low-luminosity GRBs in the $E_{p,i}$–$E_{\rm iso}$ plane.

Intrinsic dispersion in the Amati plane is largely attributed to variation in the bulk Lorentz factor $\Gamma_0$ (i.e., a "sequence in $\Gamma_0$" (Ghirlanda et al., 2013)), with secondary contributions from microphysical parameter scatter.

3. Observational Datasets, Analysis Methodology, and Evolution Tests

Empirical calibration of the iso-energy correlation relies on long-GRB samples with secure redshift, time-integrated Band-function spectral fits, and broadband fluence measurements; typical sample sizes now exceed $N=150$ (Demianski et al., 2016). Calibration procedures rigorously avoid cosmological “circularity” by using SN Ia–anchored distance moduli at low $z$ (e.g., local regression LOESS (Demianski et al., 2016, Demianski et al., 2011)), or via approximate cosmology-independent luminosity distance estimators.

Robustness checks for redshift evolution involve:

• Sample splitting at critical $z$ (e.g., $z=2$) and joint cosmological fits,
• Rank-correlation analysis (Spearman $C(z, E_{p,i})$, $C(z, E_{\rm iso})$),
• Likelihood-based 3D evolution fitting of the form $E_{p,i}' = E_{p,i}/(1+z)^{k_p}$, $E_{\rm iso}' = E_{\rm iso}/(1+z)^{k_{\rm iso}}$, and explicit testing for nonzero evolution coefficients.

Across all techniques, results consistently show negligible redshift-dependent evolution of the Amati slope or zero-point within current measurement precision (e.g., $k_{\rm iso} = -0.04\pm0.10$, $a k_p = -0.02\pm0.20$ (Demianski et al., 2016, Demianski et al., 2011)).

4. Selection Effects, Systematics, and the Debate on Physicality

Critical investigation of the iso-energy relation demonstrates that the observed correlation arises from an interplay of intrinsic physical boundaries and instrumental/observational effects.

Key findings:

• The lower-right boundary (high $E_{\rm iso}$, low $E_{p,i}$) reflects a real absence of bright, soft GRBs—interpreted as a physical limit (Heussaff et al., 2013, Kocevski, 2011).
• The upper-left boundary (low $E_{\rm iso}$, high $E_{p,i}$) primarily arises from trigger threshold and redshift-measurement bias, as the effective area declines steeply outside the central detector band (Shahmoradi, 2013, Heussaff et al., 2013).
• Population synthesis and forward-modeling confirm that much of the apparent tightness and steepness of the Amati relation is sculpted by selection: after correction, the true population-level slope reduces to $m\sim0.25$–$0.3$ and the intrinsic scatter increases to $\sigma_{\log E_{p,i}}\sim0.4$–$0.45$ dex (Shahmoradi, 2013, Aldowma et al., 18 Oct 2025, Kocevski, 2011).
• High-quality, spectroscopically complete ("gold sample") datasets yield systematically tighter and less biased relations (Tsutsui et al., 2010, Demianski et al., 2011).

A consensus emerges that the iso-energy correlation is partly intrinsic but substantially narrowed and steepened by selection effects—resolving previous debates regarding the presence of physical or artifact origins (Heussaff et al., 2013, Kocevski, 2011).

5. Multidimensional Extensions and Pulse-wise Correlations

Beyond the time-integrated Amati relation, multidimensional correlations incorporate peak isotropic luminosity ($L_{\rm iso}$, Yonetoku relation), afterglow energetics (e.g., $E_{X,{\rm iso}}$), and duration ($T_{90}$). Statistical treatments extending to 3D or higher (e.g., Tsutsui relation: $L_p\propto E_p^{1.7} T_L^{-0.4}$ (Tsutsui et al., 2010)) yield lower intrinsic scatter and better standardization.

Time-resolved and pulse-wise analyses reveal even tighter physical correlations. The "zero-fluence" $E_{\rm peak,0}$–$E_{\rm iso}$ relation yields Spearman $r=0.96$ and reduced scatter ($\sigma_{\log E_{p,0}}\sim0.29$ dex) at the single-pulse level (Basak et al., 2012, Basak et al., 2012). These findings suggest that the physical processes governing each pulse may be even more universal than those for the integrated burst.

Three-parameter correlations (e.g., $E_{X,{\rm iso}}$–$E_{\gamma,{\rm iso}}$–$E_{pk}$ (Zaninoni et al., 2015)) extend the framework, exhibit $\sim$0.3 dex scatter, and are robust to redshift, morphology, and GRB class.

6. Cosmological Applications and Standard Candle Prospects

Once empirically calibrated, the iso-energy correlation enables construction of a GRB Hubble diagram extending to $z>9$, beyond the reach of SN Ia and baryon acoustic oscillation probes (Demianski et al., 2016, Demianski et al., 2011). Bayesian simultaneous fits of both correlation parameters and cosmological densities (e.g., $\Omega_m$, $w_0$, $w_1$) have been performed, yielding cosmological constraints consistent with Planck and SN-based values, though with larger uncertainties due to intrinsic correlation scatter (Demianski et al., 2016, Aldowma et al., 18 Oct 2025).

The current limitation is the larger intrinsic scatter of the Amati relation (0.2–0.45 dex) compared to SNe Ia ($\approx$0.15 mag). However, the high-redshift reach and prospects for future improvements (in sample size and measurement precision, particularly from Fermi, SVOM, and THESEUS) give the empirical iso-energy relation strategic importance for cosmology and for probing the dark energy equation of state at early epochs.

---

Table: Representative Best-fit Parameters for the Empirical Iso-Energy (Amati) Correlation

Dataset/Analysis	Slope $m$	Normalization $K$ (keV)	Intrinsic Scatter (dex)
Amati+2008 (long GRBs)	$0.51 \pm 0.03$	$100 \pm 9$	$\sim0.20$
Fermi–GBM (joint fits)	$0.25 \pm 0.07$	$360$	$0.42 \pm 0.06$
BATSE (debiased pop.)	$0.25 \pm 0.03$	$100$	$\sim0.4$
Pulse-wise $E_{p,0}$	$0.555 \pm 0.050$	$437 \pm 50$	$0.29$

(Corresponding sources: (Amati, 2010, Aldowma et al., 18 Oct 2025, Shahmoradi, 2013, Basak et al., 2012))

---

In sum, the empirical iso-energy correlation is a robust, physically grounded, yet selection-effect-prone statistical relationship in GRB astrophysics. It serves as a fundamental probe of relativistic outflow physics, jet structure, and cosmology. Ongoing and future work will further clarify its physical drivers, reduce systematic uncertainties, and refine its application as a cosmological standard candle.