REVIEW
Ranking Prompt--Plateau Correlations in Gamma-Ray Bursts with Symbolic Regression
High Energy Astrophysical Phenomena · Submitted 16 Oct 2026 · 295 references
Pith's verdict: Competent symbolic-regression study that systematically ranks GRB plateau correlations and recovers a prompt-plateau plane; limitations on the intrinsic claim are clearly stated and honestly scoped.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Symbolic regression reveals a tighter GRB plateau relation linking prompt and afterglow emission.
The reading
What carries the argument
Symbolic regression, implemented via the PySR framework (Cranmer 2023), searches over mathematical expressions built from $+, -, \times, /, \exp, \log, \sqrt{}$ without assuming a functional form. It evaluates each candidate on a held-out validation split, ranking feature subsets by validation mean absolute error. This allows the authors to ask not just what the best fit is, but how the classical Dainotti relation ranks among all low-dimensional alternatives. The method is used in a Dainotti-anchored search (where $\log T_a^{\rm rest}$ is always present) and an unrestricted search over all non-empty subsets of at most three of the 13 candidate variables.
Extended reading notes
Core claim
The central claim is that the X-ray plateau luminosity $L_a$ depends on both the prompt peak luminosity $L_{\rm peak}$ and the rest-frame plateau duration $T_a^{\rm rest}$, not on $T_a^{\rm rest}$ alone. Symbolic regression ranks the two-variable relation $\log L_a \approx 0.855 \log L_{\rm peak} - 0.782 \log T_a^{\rm rest}$ as the best among 377 candidate spaces, with validation $R^2 \approx 0.88$. Direct regression yields $\log L_a = 4.718 + 0.776 \log L_{\rm peak} - 0.805 \log T_a^{\rm rest}$, with $R_P = 0.91$ for 246 bursts. After removing the shared luminosity-distance dependence (which both $L_a$ and $L_{\rm peak}$ inherit from $D_L^2$), the relation retains $R^2 \approx 0.62$, consistent with an intrinsic physical component. The classical Dainotti relation ($\log L_a$ vs. $\log T_a^{\rm rest}$ alone) ranks within the top $\sim25\%$ of the search spaces, confirming it is competitive but not the tightest.
Load-bearing premise
The apparent intrinsic component of the prompt–plateau relation is estimated by removing only the linear dependence on luminosity distance $D_L$, not by a full Efron–Petrosian treatment of sample truncation and luminosity evolution; if those biases are significant, the residual correlation could be an artifact rather than a true physical relation.
Editorial extensions
If this is right
- The Dainotti relation is shown to be a two-dimensional projection of a three-dimensional prompt–plateau plane, placing earlier one-dimensional fits in a broader context.
- Plateau luminosity carries information about the prompt energy release, meaning that any model of the plateau phase must couple to the central engine activity that powers the prompt emission.
- Additional spectral, absorption, and afterglow-shape variables (photon index, column density, decay indices) do not improve the plateau luminosity prediction beyond the prompt–plateau pair, suggesting the dominant low-dimensional physics is captured by these two observables.
- The prompt–plateau plane may help standardize GRBs for cosmological distance measurements, analogous to how the Dainotti relation has been used, potentially with lower scatter.
- The residual trend after removing shared distance dependence ($R^2 \approx 0.62$) supports an intrinsic component, but a full Efron–Petrosian truncation analysis is needed to separate it from evolution and selection effects.
Reading between the lines
- Symbolic regression could be applied to other GRB correlations (e.g., the Amati relation, the Yonetoku relation) to test whether they are projections of higher-dimensional planes, or whether the apparent correlation changes when additional variables are considered.
- The fact that the ratio $L_a T_a / L_{\rm peak}$ is not constant (it increases with $T_a$) argues against simple constant-efficiency energy injection models; future theoretical work should explore prompt–plateau coupling with variable conversion efficiency.
- The class-dependent fits (XRR, LGRB, SNe-GRB, short/SEE) show the same sign pattern, hinting that the prompt–plateau link may be universal across GRB subclasses, though small-number statistics prevent a definitive test here.
Original abstract
Gamma-ray burst (GRB) correlations probe central engines, relativistic outflows, and radiation processes. The Dainotti relation is an anticorrelation between the X-ray plateau luminosity $L_a$ and rest-frame plateau duration $T_a^{\rm rest}$. We use symbolic regression to compare it with other low-dimensional correlations in 251 GRBs with Swift-XRT plateau properties and Swift-BAT prompt-emission parameters. We perform a Dainotti-guided search, in which $\log T_a^{\rm rest}$ is always present, and an unrestricted search over subsets of at most three variables. The guided search recovers the Dainotti dependence; a D'Agostini fit gives $\log L_a=(51.92\pm0.30)-(1.31\pm0.09)\log T_a^{\rm rest}$, with intrinsic scatter $\sigma_{\rm int}=0.82\pm0.05$ dex. In the unrestricted ranking of 377 feature subsets, the classical Dainotti relation remains competitive, within the top $\sim25\%$ of the tested search spaces. The highest-ranked screened expression is $\log L_a\approx0.855\log L_{\rm peak}-0.782\log T_a^{\rm rest}$, with RMSE $\approx0.46$ and $R^2\approx0.88$. The highest-ranked runs repeatedly reduce to this prompt--plateau structure, without retaining additional spectral, absorption, or afterglow-shape variables. A direct regression gives $\log L_a=4.718+0.776\log L_{\rm peak}-0.805\log T_a^{\rm rest}$, with $R_{\rm P}=0.91$ for 246 bursts. An errors-in-variables refit gives consistent coefficients. After removing shared luminosity-distance dependence, the relation retains non-negligible explanatory power ($R^2\approx0.62$), consistent with a possible intrinsic component. Restricted Monte Carlo null tests give $p\approx0.01$. Plateau luminosity therefore shows a coupled dependence on prompt luminosity and plateau duration.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
Gamma-ray burst (GRB) correlations probe central engines, relativistic outflows, and radiation processes. The Dainotti relation is an anticorrelation between the X-ray plateau luminosity $L_a$ and rest-frame plateau duration $T_a^{\rm rest}$. We use symbolic regression to compare it with other low-dimensional correlations in 251 GRBs with Swift-XRT plateau properties and Swift-BAT prompt-emission parameters. We perform a Dainotti-guided search, in which $\log T_a^{\rm rest}$ is always present, and an unrestricted search over subsets of at most three variables. The guided search recovers the Dainotti dependence; a D'Agostini fit gives $\log L_a=(51.92\pm0.30)-(1.31\pm0.09)\log T_a^{\rm rest}$, with intrinsic scatter $\sigma_{\rm int}=0.82\pm0.05$ dex. In the unrestricted ranking of 377 feature subsets, the classical Dainotti relation remains competitive, within the top $\sim25\%$ of the tested search spaces. The highest-ranked screened expression is $\log L_a\approx0.855\log L_{\rm peak}-0.782\log T_a^{\rm rest}$, with RMSE $\approx0.46$ and $R^2\approx0.88$. The highest-ranked runs repeatedly reduce to this prompt--plateau structure, without retaining additional spectral, absorption, or afterglow-shape variables. A direct regression gives $\log L_a=4.718+0.776\log L_{\rm peak}-0.805\log T_a^{\rm rest}$, with $R_{\rm P}=0.91$ for 246 bursts. An errors-in-variables refit gives consistent coefficients. After removing shared luminosity-distance dependence, the relation retains non-negligible explanatory power ($R^2\approx0.62$), consistent with a possible intrinsic component. Restricted Monte Carlo null tests give $p\approx0.01$. Plateau luminosity therefore shows a coupled dependence on prompt luminosity and plateau duration.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Gamma-ray burst (GRB) correlations probe central engines, relativistic outflows, and radiation processes. | background 1.00 | 0.05 | 0.02 | asserted 1.00 |
| s1 | The Dainotti relation is an anticorrelation between the X-ray plateau luminosity $L_a$ and rest-frame plateau duration $T_a^{\rm rest}$. | background 1.00 | 0.05 | 0.03 | asserted 1.00 |
| s2 | We use symbolic regression to compare it with other low-dimensional correlations in 251 GRBs with Swift-XRT plateau properties and Swift-BAT prompt-emission parameters. | method 0.77 | 0.94 | 0.04 | asserted 0.97 |
| s3 | We perform a Dainotti-guided search, in which $\log T_a^{\rm rest}$ is always present, and an unrestricted search over subsets of at most three variables. | method 0.99 | 0.94 | 0.03 | asserted 0.94 |
| s4 | The guided search recovers the Dainotti dependence; a D'Agostini fit gives $\log L_a=(51.92\pm0.30)-(1.31\pm0.09)\log T_a^{\rm rest}$, with intrinsic scatter $\sigma_{\rm int}=0.82\pm0.05$ dex. | result 1.00 | 0.95 | 0.11 | asserted 1.00 |
| s5 | In the unrestricted ranking of 377 feature subsets, the classical Dainotti relation remains competitive, within the top $\sim25\%$ of the tested search spaces. | result 1.00 | 0.95 | 0.09 | asserted 1.00 |
| s6 | The highest-ranked screened expression is $\log L_a\approx0.855\log L_{\rm peak}-0.782\log T_a^{\rm rest}$, with RMSE $\approx0.46$ and $R^2\approx0.88$. | result 1.00 | 0.95 | 0.17 | asserted 1.00 |
| s7 | The highest-ranked runs repeatedly reduce to this prompt--plateau structure, without retaining additional spectral, absorption, or afterglow-shape variables. | result 0.99 | 0.93 | 0.12 | asserted 1.00 |
| s8 | A direct regression gives $\log L_a=4.718+0.776\log L_{\rm peak}-0.805\log T_a^{\rm rest}$, with $R_{\rm P}=0.91$ for 246 bursts. | result 0.99 | 0.94 | 0.12 | asserted 1.00 |
| s9 | An errors-in-variables refit gives consistent coefficients. | result 0.77 | 0.89 | 0.03 | asserted 1.00 |
| s10 | After removing shared luminosity-distance dependence, the relation retains non-negligible explanatory power ($R^2\approx0.62$), consistent with a possible intrinsic component. | result 0.84 | 0.93 | 0.10 | hedged possibility 0.95 |
| s11 | Restricted Monte Carlo null tests give $p\approx0.01$. | result 0.93 | 0.92 | 0.04 | asserted 0.99 |
| s12 | Plateau luminosity therefore shows a coupled dependence on prompt luminosity and plateau duration. | interpretation 0.89 | 0.87 | 0.84 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2610.23155 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.