---
title: Dainotti Relation in GRB Afterglows
url: https://www.emergentmind.com/topics/dainotti-relation
type: topic
---

# Dainotti Relation in GRB Afterglows

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 [2301.10572], [2402.13115].

## 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 \(L_{\rm peak}\) or \(L_p\),
\[
\log L_a = a \log T_a^* + b \log L_{\rm peak} + c ,
\]
and is widely termed the GRB fundamental plane; some reviews explicitly identify this fundamental plane with the Dainotti relation in the broader sense [2301.10572].

The quantities entering the relation are rest-frame observables. The time variable is the observed plateau end time corrected by cosmological time dilation, \(T_a^*=T_a/(1+z)\). In the X-ray construction, the plateau luminosity is computed from the observed flux through
\[
L_a = 4\pi D_L^2(z)\,F_X\,K ,
\]
or equivalently \(L_X=4\pi D_L^2F_XK_{\rm plateau}\), where \(K\) or \(K_{\rm plateau}\) is the bandpass correction and \(D_L\) 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 [1612.02917], [2204.08710].

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 \((L_a,T_a^*)\) is standard; in some cosmological and prompt–afterglow comparison studies the same end-time quantity is denoted \(T_X^*\) or \(T_{\rm brk}\). These are not distinct observables in substance, but distinct conventions around the end of the shallow or plateau phase [1208.1680].

## 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 \(T_a\), 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 \(|\alpha_1|<0.5\), so that \(T_a\) is the break time and \(L_a\) 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 [1612.02917], [2210.03870].

Because both variables carry measurement uncertainty and the astrophysical dispersion is substantial, modern fits usually include an intrinsic-scatter term \(\sigma_{\rm int}\). 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`, 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 [2203.15538], [2201.05245].

The dependence on \(D_L(z)\) 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 [2110.14840].

## 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 \((\log T_a^*,\log L_{\rm peak},\log L_a)\) space. This construction is motivated empirically by the prompt–afterglow coupling and statistically by a measurable reduction of intrinsic scatter [1704.04908].

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 \(\Lambda\)CDM, the quoted information-criterion differences are \(\Delta{\rm AIC}'=25.86\), \(51.97\), and \(83.58\) 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 \(27\%\)–\(29\%\) when moving from the 2D relation to the 3D plane, with representative flat-\(\Lambda\)CDM values changing from \(0.515\to0.369\) for Platinum, \(0.747\to0.543\) for LGRB95, and \(0.673\to0.479\) for LGRB145 [2204.08710].

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,
\[
\log L_a=(17.65\pm5.7)-(0.83\pm0.10)\log T_a+(0.64\pm0.11)\log L_{\rm peak},
\]
with \(\sigma_{\rm int}=0.32\pm0.04\). 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 [1704.04908].

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 \(\Omega_{\rm M}\) when combined with Pantheon supernovae. After correcting for redshift evolution, the X-ray full sample reaches \(\sigma_{\rm int}=0.20\pm0.06\), while the optical full sample reaches \(\sigma_{\rm int}=0.41\pm0.06\), reinforcing the view that the prompt term carries genuine explanatory power rather than acting as a purely phenomenological nuisance variable [2203.15538].

## 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,
\[
k_{L_a}=-0.40^{+0.89}_{-0.83},\qquad k_{T_a^*}=-0.17^{+0.41}_{-0.37},
\]
and an intrinsic slope
\[
b_{\rm int}=-1.02\pm0.12,
\]
essentially identical to the observed LONG-NO-SNe slope \(b=-1.0\pm0.1\). In that case, the observed Dainotti slope was therefore not appreciably steepened by redshift evolution or threshold bias [1612.02917].

Selection-bias correction has also been treated as a prerequisite for using the 3D plane in cosmology. Review-style discussions of GRBs as high-\(z\) probes emphasize that \(T_a^*\), \(L_a\), and \(L_p\) 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 [2301.10572].

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 \(\Lambda\)CDM, XCDM, and \(\phi\)CDM, 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 \(H(z)\)+BAO [2201.05245], [2110.14840].

A third strategy is external but model-independent calibration. A recent calibration program reconstructs \(H(z)\) from 33 Cosmic Chronometer measurements over \(0.07<z<1.965\) 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 \(0.553\le z\le1.96\) are used as anchors. The resulting model-independent fit gives \(a=-1.00\pm0.16\), \(b<0.21\) at \(95\%\) C.L., \(C_o=47.05^{+4.21}_{-1.35}\), and \(\sigma_{\rm int}=0.21^{+0.03}_{-0.05}\) for the 3D relation, while the corresponding 2D fit gives \(a=-1.04\pm0.16\), \(C_o=51.16\pm0.53\), and the same central scatter. When evolutionary effects are included, the intrinsic scatter tightens slightly further to \(\sigma_{\rm int}=0.20^{+0.03}_{-0.05}\), and the calibrated distance ladder can be extended to \(z=5\) [2402.13115].

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 \(0.03\le z\le9.3\); 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 [2309.00077].

## 5. Physical interpretations and population dependence

The empirical slope near \(-1\) 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 \(L_X\) is identified with the spin-down plateau luminosity and the observed rest-frame plateau end time \(T_a^*\) with the spin-down timescale. For a single dominant multipole of order \(l\), the spin-down luminosity obeys
\[
L_l(t)=L_{l,0}\left(1+\frac{t}{\tau_l}\right)^{-(1+1/l)},
\]
so that \(L_{l,0}\propto \tau_l^{-1}\) and therefore
\[
L_X\propto (T_a^*)^{-1}.
\]
This yields a theoretical Dainotti slope \(b=-1\) independent of multipole order, while the post-plateau decay index becomes \(\alpha=-(1+1/l)\), spanning \(-2\) for a dipole and approaching \(-1\) for higher orders. In a 238-burst Swift-XRT sample updated to the end of December 2024, \(79\%\) of decay indices lie between \(-2\) and \(-1\), with median \(-1.39\) and mean \(-1.54\), a distribution used to argue that multipolar magnetars can explain both the luminosity–time slope and the diversity of post-plateau decays [2507.09292].

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 \(L_X\propto T_a^{*-1}\) 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
\[
\log L_X=(51.57\pm0.10)-(1.10\pm0.03)\log T_{\rm brk}
\]
with \(\rho=-0.88\) 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 [1208.1680].

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 \(b=-1.0\pm0.1\) with \(\rho_{LT}=-0.74\), the LONG-SNe subset gives \(b=-1.5\pm0.3\) with \(\rho_{LT}=-0.83\), and the most secure spectroscopic A+B subset of seven GRB-SNe gives \(b=-1.9\pm0.3\) with \(\rho_{LT}=-0.96\). The slope difference between LONG-NO-SNe and the A+B subset yields \(P=0.005\) in the reported Student’s \(t\)-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 \(b=-1.35\pm0.24\), and the slope-difference significance weakens to \(P=0.10\) [1612.02917].

Attempts to define more physically homogeneous optical subsamples have given a more mixed outcome. Closure-relation-selected optical plateau GRBs in the favored \(\nu>\max(\nu_c,\nu_m)\) regime still give Dainotti slopes consistent with \(-1\), for example \(-1.05\pm0.30\) and \(-1.10\pm0.28\) 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 [2210.03870].

## 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 \(z\le2.26\), quasars can reach \(z=7.642\), and GRBs reach \(z=9.4\), making plateau GRBs potentially important probes of the \(z\sim2.3\)–\(9\) regime that is otherwise sparsely constrained [2301.10572].

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 \(\Omega_{\rm M}=0.299\pm0.009\) 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 \(\pm0.008\)–\(\pm0.009\) range and still track the Pantheon constraint rather than surpass it [2203.15538].

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 \(H(z)\)+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 [2201.05245], [2110.14840].

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 [2203.15538].

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
\[
a=-1.046\pm0.073,\qquad c=50.56\pm0.28,\qquad \sigma_{\rm int}=0.251^{+0.017}_{-0.020},
\]
together with a dipole amplitude \(A_d\simeq0.6\pm0.2\) pointing toward \(({\rm RA},{\rm DEC})\approx(134^\circ\pm30^\circ,-36^\circ\pm21^\circ)\). In that analysis, residual directional correlations present under isotropic \(\Lambda\)CDM disappear once the dipole term is included [2510.20705].

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-\(z\) probes and as part of a future toolkit, but the main \(H_0\)-trend analyses still rely on binned SNe Ia and SNe Ia+BAO rather than on a Dainotti-calibrated GRB Hubble diagram [2301.10572]. 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.

Source: https://www.emergentmind.com/topics/dainotti-relation