Dainotti Relation in GRB Afterglows

• Dainotti relation is an empirical correlation linking GRB plateau luminosity and rest-frame end time, with slopes near -1 indicating key energy dissipation mechanisms.
• The extension to a 3D fundamental plane by including prompt peak luminosity reduces intrinsic scatter and better constrains the prompt-afterglow energetics.
• Robust calibration using methods like the Efron–Petrosian technique and MCMC frameworks makes it a promising tool for standardizing GRBs in cosmological distance measurements.

The Dainotti relation is an empirical gamma-ray-burst (GRB) afterglow correlation linking the luminosity at the end of the plateau phase to the rest-frame time at which that plateau ends. In its bidimensional form it is a luminosity–time anti-correlation, commonly written as $\log L_a = c + a \log T_a^*$, while in much of the recent literature the term also encompasses the three-parameter “fundamental plane” obtained by adding the prompt peak luminosity, $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$. It occupies a central position in attempts to standardize plateau GRBs, connect prompt and afterglow energetics, and extend cosmological distance measurements beyond the supernova regime, but its practical use is inseparable from questions of intrinsic scatter, redshift evolution, sample selection, and calibration strategy (Dainotti et al., 2023, Favale et al., 2024).

1. Definition, notation, and observational meaning

The relation is used in two closely related senses. The original form is the anti-correlation between the X-ray luminosity at the end of the plateau and the plateau end time in the burst rest frame. The literature alternates between $L_a$ and $L_X$ for the plateau-end luminosity, and between $T_a^*$, $T_X^*$, and $T_{\rm brk}$ for the rest-frame plateau-end time. In this notation, the standard form is

$\log L_a = c + a \log T_a^* ,$

with $a<0$, typically close to $-1$. The three-dimensional extension adds the prompt peak luminosity $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$0 or $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$1,

$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$2

and is widely termed the GRB fundamental plane; some reviews explicitly identify this fundamental plane with the Dainotti relation in the broader sense (Dainotti et al., 2023).

The quantities entering the relation are rest-frame observables. The time variable is the observed plateau end time corrected by cosmological time dilation, $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$3. In the X-ray construction, the plateau luminosity is computed from the observed flux through

$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$4

or equivalently $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$5, where $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$6 or $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$7 is the bandpass correction and $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$8 is the luminosity distance. When the 3D relation is used, the prompt term is likewise a luminosity derived from the observed 1 s peak prompt flux. The relation is therefore not a direct observer-frame flux–time law; it is a rest-frame luminosity correlation whose construction already depends on redshift handling and on a luminosity-distance prescription (Dainotti et al., 2016, Cao et al., 2022).

A persistent source of notational ambiguity is that the same physical relation appears under different variable labels in different subfields. In afterglow plateau work the pair $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$9 is standard; in some cosmological and prompt–afterglow comparison studies the same end-time quantity is denoted $L_a$0 or $L_a$1. These are not distinct observables in substance, but distinct conventions around the end of the shallow or plateau phase (Sultana et al., 2012).

2. Measurement pipeline and statistical form

The plateau observables are model-derived rather than purely visual. In the main Swift-XRT literature, the light curves are commonly fit with the Willingale et al. phenomenological prompt-plus-afterglow representation, and the plateau end time is identified with the afterglow transition time $L_a$2, where the afterglow component changes from an exponential-like phase to a power-law decline. In optical work, a broken power law is often used instead, with an optical plateau operationally defined by $L_a$3, so that $L_a$4 is the break time and $L_a$5 is the optical luminosity at that epoch. These constructions make the Dainotti relation sensitive to the adopted light-curve model, data coverage, flare contamination, and plateau morphology (Dainotti et al., 2016, Dainotti et al., 2022).

Because both variables carry measurement uncertainty and the astrophysical dispersion is substantial, modern fits usually include an intrinsic-scatter term $L_a$6. The standard regression framework in this literature is the D’Agostini method, implemented either directly or within MCMC samplers. Recent analyses use this framework with cobaya, MontePython interfaced with [CLASS](https://www.emergentmind.com/topics/colorado-learning-attitudes-about-science-survey-class), and emcee, depending on whether the aim is correlation fitting, joint cosmology–correlation inference, or model-independent calibration. In this setting the Dainotti relation is treated not simply as a line or plane in log space, but as a stochastic relation with measurement errors in all coordinates and an extra variance term absorbing unresolved astrophysical and instrumental systematics (Dainotti et al., 2022, Cao et al., 2022).

The dependence on $L_a$7 is the key technical obstacle. Plateau luminosity and prompt peak luminosity are not directly observed; they are inferred from fluxes and a cosmological distance law. This is why any cosmological use of the relation faces the circularity problem: a cosmology is needed to construct the relation, yet the calibrated relation is then intended to constrain cosmology. Much of the methodological development of the subject consists of different attempts to weaken or bypass this dependence (Cao et al., 2021).

3. From the 2D anti-correlation to the 3D fundamental plane

The 3D extension emerged from the recognition that the 2D luminosity–time anti-correlation leaves substantial residual scatter. By adding the prompt peak luminosity, the relation becomes a plane in $L_a$8 space. This construction is motivated empirically by the prompt–afterglow coupling and statistically by a measurable reduction of intrinsic scatter (Dainotti et al., 2017).

A detailed comparative analysis of three long-GRB samples—Platinum (50 GRBs), LGRB95, and the combined LGRB145—found that the 3D form is very strongly favored over the 2D form by AIC, BIC, and DIC in every cosmological model tested. In flat $L_a$9CDM, the quoted information-criterion differences are $L_X$0, $L_X$1, and $L_X$2 for Platinum, LGRB95, and LGRB145, respectively, all far above the threshold used there for “very strong” evidence. The same study reported an intrinsic-scatter reduction of about $L_X$3–$L_X$4 when moving from the 2D relation to the 3D plane, with representative flat-$L_X$5CDM values changing from $L_X$6 for Platinum, $L_X$7 for LGRB95, and $L_X$8 for LGRB145 (Cao et al., 2022).

The high-quality “gold” and “platinum” subsamples are central to this refinement. In an updated Swift-XRT study of 183 plateau GRBs, the gold sample of 45 GRBs—selected for good coverage and relatively flat plateaus—defined the tightest plane,

$L_X$9

with $T_a^*$0. In that analysis, most GRB categories were statistically compatible with the gold plane, while short bursts with extended emission were the notable exception and were interpreted as a physically distinct class (Dainotti et al., 2017).

The same transition from 2D to 3D has been extended beyond X-rays. Optical and X-ray fundamental planes have been treated in parallel as cosmological distance indicators, with full optical samples and carefully trimmed subsamples showing that the optical 3D Dainotti correlation can be as efficacious as the X-ray one in constraining $T_a^*$1 when combined with Pantheon supernovae. After correcting for redshift evolution, the X-ray full sample reaches $T_a^*$2, while the optical full sample reaches $T_a^*$3, reinforcing the view that the prompt term carries genuine explanatory power rather than acting as a purely phenomenological nuisance variable (Dainotti et al., 2022).

4. Selection effects, redshift evolution, and calibration strategies

The relation cannot be interpreted or used without systematic control. A recurring theme across the literature is that apparent luminosity–time correlations can be distorted by detector thresholds, truncation, and redshift evolution. The most widely adopted corrective framework is the Efron–Petrosian method, used to de-evolve observables and recover intrinsic correlations. In a targeted comparison of long GRBs with and without associated supernovae, the Efron–Petrosian analysis of the LONG-NO-SNe sample yielded very weak evolution,

$T_a^*$4

and an intrinsic slope

$T_a^*$5

essentially identical to the observed LONG-NO-SNe slope $T_a^*$6. In that case, the observed Dainotti slope was therefore not appreciably steepened by redshift evolution or threshold bias (Dainotti et al., 2016).

Selection-bias correction has also been treated as a prerequisite for using the 3D plane in cosmology. Review-style discussions of GRBs as high-$T_a^*$7 probes emphasize that $T_a^*$8, $T_a^*$9, and $T_X^*$0 should be corrected for selection biases and redshift evolution through the Efron–Petrosian method, and that reliable GRB standardization requires correlations independent of the cosmological model. This perspective is especially explicit in the plateau-based “platinum sample” program, where bias correction and morphological sample cleaning are treated as coupled requirements rather than as separate refinements (Dainotti et al., 2023).

A second strategy is joint inference of cosmology and correlation parameters. In this approach, the Dainotti coefficients and intrinsic scatter are fitted simultaneously with cosmological parameters so that the relation is not pre-calibrated under a fixed background model. This is the logic behind the claim that some Dainotti-correlated samples are “standardizable”: across flat and non-flat $T_X^*$1CDM, XCDM, and $T_X^*$2CDM, the inferred correlation parameters remain mutually consistent within errors. Both the 50-burst Platinum sample and the three Dainotti-correlated data sets analyzed in later cosmological studies were presented in this operational sense as cosmological-model-independent standardizable samples, even though their GRB-only cosmological constraints remain weaker than those from $T_X^*$3+BAO (Cao et al., 2022, Cao et al., 2021).

A third strategy is external but model-independent calibration. A recent calibration program reconstructs $T_X^*$4 from 33 Cosmic Chronometer measurements over $T_X^*$5 using Gaussian Processes, converts the reconstruction into luminosity distances, and calibrates the Dainotti relation directly in distance space. To stay within the Cosmic Chronometer redshift range, 20 Platinum GRBs in $T_X^*$6 are used as anchors. The resulting model-independent fit gives $T_X^*$7, $T_X^*$8 at $T_X^*$9 C.L., $T_{\rm brk}$0, and $T_{\rm brk}$1 for the 3D relation, while the corresponding 2D fit gives $T_{\rm brk}$2, $T_{\rm brk}$3, and the same central scatter. When evolutionary effects are included, the intrinsic scatter tightens slightly further to $T_{\rm brk}$4, and the calibrated distance ladder can be extended to $T_{\rm brk}$5 (Favale et al., 2024).

An important interpretive caution is that not every paper that discusses the Dainotti relation actually uses it in its inference pipeline. A prominent example is a Hubble-parameter reconstruction study in which the Dainotti relation appears only in a methodological comparison section. The GRB cosmology in that work is based on an Amati-type relation for 162 long GRBs over $T_{\rm brk}$6; no Dainotti coefficients, intrinsic scatter, distance moduli, or Dainotti-specific likelihood are fitted. The Dainotti relation is there a contextual bridge to recent likelihood and calibration work, not the engine of the reported cosmological constraints (Sudharani et al., 2023).

5. Physical interpretations and population dependence

The empirical slope near $T_{\rm brk}$7 has motivated several physical interpretations. The most developed recent account is the magnetar spin-down picture generalized to multipolar magnetic fields. In this framework, the observed plateau luminosity $T_{\rm brk}$8 is identified with the spin-down plateau luminosity and the observed rest-frame plateau end time $T_{\rm brk}$9 with the spin-down timescale. For a single dominant multipole of order $\log L_a = c + a \log T_a^* ,$0, the spin-down luminosity obeys

$\log L_a = c + a \log T_a^* ,$1

so that $\log L_a = c + a \log T_a^* ,$2 and therefore

$\log L_a = c + a \log T_a^* ,$3

This yields a theoretical Dainotti slope $\log L_a = c + a \log T_a^* ,$4 independent of multipole order, while the post-plateau decay index becomes $\log L_a = c + a \log T_a^* ,$5, spanning $\log L_a = c + a \log T_a^* ,$6 for a dipole and approaching $\log L_a = c + a \log T_a^* ,$7 for higher orders. In a 238-burst Swift-XRT sample updated to the end of December 2024, $\log L_a = c + a \log T_a^* ,$8 of decay indices lie between $\log L_a = c + a \log T_a^* ,$9 and $a<0$0, with median $a<0$1 and mean $a<0$2, a distribution used to argue that multipolar magnetars can explain both the luminosity–time slope and the diversity of post-plateau decays (Yorgancioglu et al., 12 Jul 2025).

That interpretation is not unique. The same multipolar study explicitly notes that black-hole spin-down through the Blandford–Znajek process in a MAD state can also yield an $a<0$3 scaling, so the Dainotti relation by itself does not prove a magnetar origin. Earlier comparative work also proposed a kinematic interpretation by placing the Dainotti relation beside the prompt lag–luminosity relation. In a common Swift subset, the afterglow relation

$a<0$4

with $a<0$5 was found to align strikingly with the prompt lag–luminosity trend, leading to the suggestion that both may be manifestations of Doppler and viewing-angle effects rather than entirely separate prompt and afterglow microphysics (Sultana et al., 2012).

Population dependence is another major theme. Long GRBs associated with supernovae define a steeper and tighter relation than long GRBs without observed supernova association. In a Swift sample of 176 plateau GRBs, the LONG-NO-SNe subset gives $a<0$6 with $a<0$7, the LONG-SNe subset gives $a<0$8 with $a<0$9, and the most secure spectroscopic A+B subset of seven GRB-SNe gives $-1$0 with $-1$1. The slope difference between LONG-NO-SNe and the A+B subset yields $-1$2 in the reported Student’s $-1$3-test, suggesting that the plateau energetics of securely SN-associated bursts may differ from those of the broader long-GRB population. Yet this conclusion is not immune to modeling choices: after rough beaming corrections the A+B slope becomes $-1$4, and the slope-difference significance weakens to $-1$5 (Dainotti et al., 2016).

Attempts to define more physically homogeneous optical subsamples have given a more mixed outcome. Closure-relation-selected optical plateau GRBs in the favored $-1$6 regime still give Dainotti slopes consistent with $-1$7, for example $-1$8 and $-1$9 before correction for the ISM-like and wind-like subsets, respectively. However, the intrinsic scatter is not significantly lower than in the larger parent optical sample. In this sense, closure-relation selection supports the robustness of the slope but has not yet isolated a dramatically tighter optical standardization subclass (Dainotti et al., 2022).

6. Cosmological role, present performance, and research frontier

The cosmological appeal of the relation is straightforward: GRBs extend the Hubble diagram far beyond Type Ia supernovae. Review papers highlight that Pantheon supernovae reach $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$00, quasars can reach $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$01, and GRBs reach $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$02, making plateau GRBs potentially important probes of the $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$03–$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$04 regime that is otherwise sparsely constrained (Dainotti et al., 2023).

In current data, the Dainotti relation functions as a standardizable-candle relation rather than a precision standard candle. Joint analyses of Pantheon SNe Ia and GRB fundamental-plane samples give cosmological parameters consistent with supernova-only results but do not materially improve them. A representative result is $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$05 from SNe Ia plus the X-ray Platinum fundamental-plane sample, with the same central value recovered for optical and trimmed variants; after evolution corrections the uncertainties remain in the $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$06–$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$07 range and still track the Pantheon constraint rather than surpass it (Dainotti et al., 2022).

GRB-only cosmology remains substantially weaker. Studies of the Platinum sample alone, or of jointly analyzed Dainotti- and Amati-correlated samples, report constraints broadly consistent with $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$08+BAO but markedly less precise. The significance of these analyses lies less in present parameter precision than in the demonstration that some plateau-GRB samples can be treated as standardizable across multiple dark-energy models, and in the fact that they probe redshift intervals not covered by the standard late-time probes with comparable leverage (Cao et al., 2022, Cao et al., 2021).

Forecasting studies are correspondingly future-oriented. Simulations based on optical and X-ray fundamental planes suggest that the optical plane may become competitive faster than the X-ray plane if errors shrink and the sample increases. With halved errors, 142 and 284 simulated optical plateau GRBs are reported as sufficient to match the SN Ia precision levels of 2011 and 2014, respectively, while 390 optical GRBs would be needed to match current supernova precision, with the corresponding date estimate extending to 2054 under the stated assumptions (Dainotti et al., 2022).

The relation has also begun to support applications beyond standard parameter fitting. A recent anisotropy analysis standardized 176 long Swift GRBs with the redshift-corrected bidimensional X-ray Dainotti relation and searched for a dipolar modulation in the resulting GRB Hubble diagram. The main combined-sample fit gives

$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$09

together with a dipole amplitude $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$10 pointing toward $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$11. In that analysis, residual directional correlations present under isotropic $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$12CDM disappear once the dipole term is included (Santiago et al., 23 Oct 2025).

At the same time, the literature is explicit that the Dainotti relation has not yet provided a direct resolution of the Hubble-constant tension. In conference-style discussions of the tension, GRBs and the Dainotti relation are presented as promising high-$\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$13 probes and as part of a future toolkit, but the main $\log L_a = a \log T_a^* + b \log L_{\rm peak} + c$14-trend analyses still rely on binned SNe Ia and SNe Ia+BAO rather than on a Dainotti-calibrated GRB Hubble diagram (Dainotti et al., 2023). The current state of the subject is therefore dual: the Dainotti relation is one of the most developed empirical routes toward GRB standardization, yet its full cosmological utility still depends on reducing intrinsic scatter, controlling evolution and selection effects, and clarifying which GRB subclasses obey the relation most fundamentally.