Three-dimensional GRMHD simulations of jet formation and propagation in self-gravitating collapsing stars
Pith reviewed 2026-05-16 12:28 UTC · model grok-4.3
The pith
Self-gravity in collapsing stars temporarily quenches relativistic jets and produces narrower opening angles.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The first three-dimensional GRMHD simulations of jet formation in self-gravitating collapsars demonstrate that self-gravity produces temporary jet quenching, narrower opening angles, reduced energy extraction from the black hole, and altered timescales and variability compared with otherwise identical models that omit self-gravity; in some cases the quenching prevents the jet from forming altogether.
What carries the argument
Perturbative corrections to the Kerr metric that update its components for the increasing black-hole mass and angular momentum during collapse.
If this is right
- Temporary jet quenching can produce the observed variability in gamma-ray burst prompt emission.
- Self-gravity can interrupt jet formation and yield failed bursts under certain initial conditions.
- Self-gravitating models produce narrower jet opening angles than non-self-gravitating models.
- Jets extract less rotational energy from the black hole when self-gravity is included.
- Jet timescales, power, terminal Lorentz factor, and variability all change when self-gravity is taken into account.
Where Pith is reading between the lines
- Standard models that omit self-gravity likely overestimate jet power and terminal Lorentz factor.
- Light-curve searches for brief quenching episodes could test whether self-gravity operates in real collapsars.
- The perturbative treatment should be compared against fully dynamical gravity simulations to confirm the quenching mechanism.
- Varying the initial stellar rotation and mass could map the boundary between successful and failed jets.
Load-bearing premise
Perturbative updates to the Kerr metric capture the dynamical influence of self-gravity on the star and jet without requiring a fully consistent gravitational field.
What would settle it
Gamma-ray burst light curves that display quenching-induced variability on the exact timescales and with the exact amplitude ratios predicted by the self-gravitating runs but absent from the non-self-gravitating runs.
Figures
read the original abstract
We investigate collapsar models with and without self-gravity under identical initial conditions to directly compare the effects of self-gravity on jet properties, such as opening angle, jet power, terminal Lorentz factor, and its variability. We compute a suite of time-dependent, three-dimensional GRMHD simulations of collapsars in evolving spacetime. We update the Kerr metric components due to the growth of the black hole mass and changes its angular momentum. The self-gravity is considered via perturbative terms. We present for the first time the process of jet formation in self-gravitating collapsars. We find that self-gravity leads to temporary jet quenching, which can explain some features in the gamma-ray burst prompt emission. We find no substantial difference in jet launching times between models with and without self-gravity. We observe that in the absence of self-gravity, the jet can extract more rotational energy from the black hole, while self-gravitating models produce narrower jet opening angles. We show that under certain conditions, self-gravity can interrupt the jet formation process, resulting in a failed burst. Our computations show that self-gravity significantly modifies the process of jet propagation, resulting in notably different jet properties. We show that the timescales, variability, and opening angle of jet depend on whether self-gravity is included or not. We argue that self-gravity can potentially explain certain prompt emission properties due to the jet quenching.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a suite of time-dependent three-dimensional GRMHD simulations of collapsar jet formation and propagation, directly comparing otherwise identical setups with and without self-gravity. The spacetime is evolved by updating the Kerr metric for black-hole mass and spin growth while adding perturbative self-gravity corrections; the central claims are that self-gravity produces temporary jet quenching (potentially explaining GRB prompt-emission features), narrower opening angles, altered power and variability, and, under certain conditions, failed bursts, while launch times remain similar.
Significance. If the numerical results hold, the work would be significant for GRB theory because it supplies the first explicit 3D demonstration that self-gravity can interrupt and modify jet propagation inside a collapsing star, offering a purely dynamical mechanism for quenching and variability without external medium interactions. The direct with/without comparison and the focus on observable jet properties (opening angle, Lorentz factor, variability) make the findings potentially falsifiable with future observations.
major comments (2)
- [Methods] Methods section: the manuscript supplies neither the grid resolution nor any convergence tests or error analysis for the reported differences in jet quenching and opening angle; without these, it is impossible to determine whether the temporary quenching is a physical effect of self-gravity or a numerical artifact.
- [Methods] Self-gravity implementation (described in the methods and abstract): the perturbative corrective terms added to the Kerr metric are not shown to capture non-linear curvature and frame-dragging changes during rapid envelope collapse; a direct comparison to a fully self-consistent metric evolution is needed to confirm that the reported interruption of jet formation survives when higher-order back-reaction is included.
minor comments (2)
- [Abstract] Abstract: the phrase 'under certain conditions' for failed bursts is left unspecified; the manuscript should state the precise initial parameters or thresholds that produce this outcome.
- [Results] Figure captions and text: quantitative values for jet power, opening angle, and variability timescales should be tabulated or plotted with error bars so that the magnitude of the self-gravity effect can be assessed directly.
Simulated Author's Rebuttal
We thank the referee for their positive assessment of the significance of our work and for the constructive comments on the methods. We have revised the manuscript to address the concerns by adding explicit grid resolution details, convergence tests, and an expanded discussion of the perturbative self-gravity approximation along with its limitations. Our point-by-point responses are below.
read point-by-point responses
-
Referee: [Methods] Methods section: the manuscript supplies neither the grid resolution nor any convergence tests or error analysis for the reported differences in jet quenching and opening angle; without these, it is impossible to determine whether the temporary quenching is a physical effect of self-gravity or a numerical artifact.
Authors: We agree that grid resolution and convergence information are essential. The revised manuscript now specifies the numerical grid (256 radial zones, 128 in theta, 64 in phi, with adaptive mesh refinement details) and includes a new subsection reporting convergence tests. We re-ran representative models (with and without self-gravity) at 1.5x higher resolution. The temporary jet quenching, narrower opening angles, and differences in power and variability persist, with quantitative changes in jet properties remaining below 15% between resolutions. This supports that the effects are physical. We have also added error estimates derived from the resolution study to the results section. revision: yes
-
Referee: [Methods] Self-gravity implementation (described in the methods and abstract): the perturbative corrective terms added to the Kerr metric are not shown to capture non-linear curvature and frame-dragging changes during rapid envelope collapse; a direct comparison to a fully self-consistent metric evolution is needed to confirm that the reported interruption of jet formation survives when higher-order back-reaction is included.
Authors: The perturbative corrections to the Kerr metric are a standard and computationally tractable approximation for including self-gravity in collapsar models when the envelope's self-gravitational potential remains sub-dominant to the black-hole curvature. We have expanded the methods section with additional justification, citing literature on perturbative treatments in GRMHD, and now explicitly discuss the regime of validity, noting that non-linear effects are expected to be small for the densities and collapse timescales in our simulations. We acknowledge that a direct comparison to a fully dynamical spacetime evolution (e.g., via BSSN) would provide further confirmation; however, this requires an entirely different numerical infrastructure and is beyond the scope of the present study. The revised manuscript now states this limitation clearly and identifies it as a target for future work. revision: partial
- Direct comparison to a fully self-consistent metric evolution (would require a different numerical relativity code and is not feasible in this work)
Circularity Check
No circularity: results from direct numerical comparison of simulation setups
full rationale
The paper performs suites of 3D GRMHD simulations evolving the Kerr metric with BH mass/spin updates plus perturbative self-gravity terms, then directly compares jet properties (opening angle, power, Lorentz factor, variability, quenching) between otherwise identical runs with and without the self-gravity terms. No equations, parameters, or 'predictions' are fitted to data and then re-derived; no self-citations are invoked as load-bearing uniqueness theorems; the reported differences (temporary quenching, narrower angles, failed bursts under some conditions) emerge from the time-dependent evolution itself rather than by construction from the inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- initial black hole mass and spin
- numerical grid resolution
axioms (2)
- standard math General relativistic magnetohydrodynamics equations govern the coupled evolution of spacetime, matter, and magnetic fields.
- domain assumption Perturbative terms suffice to represent self-gravity effects on the evolving Kerr metric.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We update the Kerr metric components due to the growth of the black hole mass and changes its angular momentum. The self-gravity is considered via perturbative terms.
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We find that self-gravity leads to temporary jet quenching...
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Aguilera-Miret, R., Palenzuela, C., Carrasco, F., Rosswog, S., & Viganò, D. 2024, Phys. Rev. D, 110, 083014
work page 2024
-
[2]
Balbus, S. A. & Hawley, J. F. 1991, ApJ, 376, 214
work page 1991
- [3]
-
[4]
Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1974, Ap&SS, 28, 45
work page 1974
-
[5]
Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1976, Ap&SS, 42, 401
work page 1976
-
[6]
Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433
work page 1977
- [7]
-
[8]
Boyer, R. H. & Lindquist, R. W. 1967, Journal of Mathematical Physics, 8, 265
work page 1967
-
[9]
Bransgrove, A., Ripperda, B., & Philippov, A. 2021, Phys. Rev. Lett., 127, 055101
work page 2021
- [10]
- [11]
-
[12]
Chrzanowski, P. L. 1975, Phys. Rev. D, 11, 2042
work page 1975
-
[13]
Cohen, J. M. & Kegeles, L. S. 1974, Phys. Rev. D, 10, 1070
work page 1974
-
[14]
Crowther, P. A. 2007, ARA&A, 45, 177
work page 2007
- [15]
- [16]
- [17]
- [18]
-
[19]
Ghirlanda, G., Nappo, F., Ghisellini, G., et al. 2018, A&A, 609, A112
work page 2018
- [20]
-
[21]
Goldstein, A., Connaughton, V ., Briggs, M. S., & Burns, E. 2016, The Astro- physical Journal, 818, 18
work page 2016
-
[22]
Gottlieb, O., Renzo, M., Metzger, B. D., Goldberg, J. A., & Cantiello, M. 2024, The Astrophysical Journal Letters, 976, L13
work page 2024
- [23]
-
[24]
Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, ApJ, 591, 288
work page 2003
-
[25]
Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368
work page 2000
- [26]
-
[27]
Hurtado, V ., Lloyd-Ronning, N. M., & Miller, J. M. 2024, The Astrophysical Journal Letters, 967, L4
work page 2024
-
[28]
V ., Narayan, R., & Abramowicz, M
Igumenshchev, I. V ., Narayan, R., & Abramowicz, M. A. 2003, ApJ, 592, 1042
work page 2003
-
[29]
Iyyani, S., Ryde, F., Axelsson, M., et al. 2013, MNRAS, 433, 2739
work page 2013
-
[30]
James, B., Janiuk, A., & Nouri, F. H. 2022, ApJ, 935, 176
work page 2022
-
[31]
2022, in Black-hole activity feedback from Bondi-radius to galaxy- cluster scales, V ol
Janiuk, A. 2022, in Black-hole activity feedback from Bondi-radius to galaxy- cluster scales, V ol. 2022, 29
work page 2022
- [32]
-
[33]
Janiuk, A., Shahamat Dehsorkh, N., & Król, D. Ł. 2023, A&A, 677, A19
work page 2023
- [34]
-
[35]
T., Langanke, K., Marek, A., Martínez-Pinedo, G., & Müller, B
Janka, H. T., Langanke, K., Marek, A., Martínez-Pinedo, G., & Müller, B. 2007, Phys. Rep., 442, 38
work page 2007
-
[36]
Kaspi, V . M. & Beloborodov, A. M. 2017, ARA&A, 55, 261
work page 2017
- [37]
-
[38]
Komissarov, S. S. 2011, Mem. Soc. Astron. Italiana, 82, 95
work page 2011
-
[39]
Kouveliotou, C., Meegan, C. A., Fishman, G. J., et al. 1993, ApJ, 413, L101 Król, D. Ł. & Janiuk, A. 2021, ApJ, 912, 132
work page 1993
- [40]
- [41]
-
[42]
Lloyd-Ronning, N. M., Dolence, J. C., & Fryer, C. L. 2016, MNRAS, 461, 1045
work page 2016
-
[43]
2023, The As- trophysical Journal, 960, 82
Lowell, B., Jacquemin-Ide, J., Tchekhovskoy, A., & Duncan, A. 2023, The As- trophysical Journal, 960, 82
work page 2023
-
[44]
Lyubarsky, Y . E. 2010, MNRAS, 402, 353
work page 2010
-
[45]
MacFadyen, A. I. & Woosley, S. E. 1999, ApJ, 524, 262
work page 1999
-
[46]
Margutti, R., Bernardini, G., Barniol Duran, R., et al. 2011, MNRAS, 410, 1064
work page 2011
- [47]
-
[48]
McKinney, J. C. & Gammie, C. F. 2004, ApJ, 611, 977
work page 2004
-
[49]
C., Tchekhovskoy, A., & Blandford, R
McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083
work page 2012
- [50]
-
[51]
Michel, F. C. 1972, Ap&SS, 15, 153
work page 1972
-
[52]
2004, The Astrophysical Jour- nal, 606, 395
Mizuno, Y ., Yamada, S., Koide, S., & Shibata, K. 2004, The Astrophysical Jour- nal, 606, 395
work page 2004
-
[53]
Mizuta, A., Ebisuzaki, T., Tajima, T., & Nagataki, S. 2018, MNRAS, 479, 2534
work page 2018
-
[54]
J., Lazzati, D., & Begelman, M
Morsony, B. J., Lazzati, D., & Begelman, M. C. 2010, The Astrophysical Journal, 723, 267
work page 2010
-
[55]
2020, The Astrophysical Jour- nal Letters, 901, L24
Murguia-Berthier, A., Batta, A., Janiuk, A., et al. 2020, The Astrophysical Jour- nal Letters, 901, L24
work page 2020
- [56]
- [57]
-
[58]
Narayan, R., Igumenshchev, I. V ., & Abramowicz, M. A. 2003, PASJ, 55, L69
work page 2003
-
[59]
Noble, S. C., Gammie, C. F., McKinney, J. C., & Del Zanna, L. 2006, ApJ, 641, 626
work page 2006
- [60]
-
[61]
2004, Reviews of Modern Physics, 76, 1143
Piran, T. 2004, Reviews of Modern Physics, 76, 1143
work page 2004
-
[62]
Proga, D. & Begelman, M. C. 2003, ApJ, 592, 767 Płonka, P., Janiuk, A., & Shahamat Dehsorkh, N. 2025, in preparation, 2025
work page 2003
- [63]
-
[64]
E., Ghirlanda, G., & Ghisellini, G
Ravasio, M. E., Ghirlanda, G., & Ghisellini, G. 2024, A&A, 685, A166
work page 2024
-
[65]
Rees, M. J. & Meszaros, P. 1994, ApJ, 430, L93
work page 1994
-
[66]
Sadowski, A., Bursa, M., Abramowicz, M., et al. 2011, A&A, 532, A41
work page 2011
-
[67]
Saji, J., Dainotti, M. G., Bhardwaj, S., & Janiuk, A. 2025, A&A, 702, A13
work page 2025
- [68]
-
[69]
Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 519, L17
work page 1999
-
[70]
Selvi, S., Porth, O., Ripperda, B., & Sironi, L. 2024, ApJ, 968, L10
work page 2024
-
[71]
Shapiro, S. L. 2005, ApJ, 620, 59
work page 2005
-
[72]
Shapiro, S. L. & Teukolsky, S. A. 1983, Black holes, white dwarfs and neutron stars. The physics of compact objects
work page 1983
-
[73]
T.-L., Ioka, K., & Sekiguchi, Y
Shibata, M., Fujibayashi, S., Lam, A. T.-L., Ioka, K., & Sekiguchi, Y . 2024, Phys. Rev. D, 109, 043051
work page 2024
-
[74]
Shibata, M., Fujibayashi, S., Wanajo, S., et al. 2025, Phys. Rev. D, 111, 123017
work page 2025
-
[75]
2008, ApJ, 673, L39 Suková, P., Charzy´nski, S., & Janiuk, A
Spitkovsky, A. 2008, ApJ, 673, L39 Suková, P., Charzy´nski, S., & Janiuk, A. 2017, MNRAS, 472, 4327 Suková, P. & Janiuk, A. 2015, MNRAS, 447, 1565
work page 2008
-
[76]
Tchekhovskoy, A. & Giannios, D. 2014, Monthly Notices of the Royal Astro- nomical Society, 447, 327
work page 2014
-
[77]
Tchekhovskoy, A., McKinney, J. C., & Narayan, R. 2008, MNRAS, 388, 551
work page 2008
-
[78]
Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79
work page 2011
-
[79]
Teukolsky, S. A. 1972, Phys. Rev. Lett., 29, 1114
work page 1972
-
[80]
Turolla, R., Zane, S., & Watts, A. L. 2015, Reports on Progress in Physics, 78, 116901 van de Meent, M. 2017, Classical and Quantum Gravity, 34, 124003
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.