{"id":"23b637ed-4403-4122-b14a-90ddbbf7f26c","arxiv_id":"2502.02077","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A collapsar disk can grow its own poloidal magnetic field through a dynamo, launching a gamma-ray-burst jet and exploding the star without a pre-existing strong poloidal field.","lead":"This paper simulates a massive collapsing star around a spinning black hole and finds that the surrounding disk can build up its own magnetic field, launching a jet and exploding the star. It offers a self-consistent path from ordinary initial conditions to long gamma-ray bursts and the bright supernovae that accompany them.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central scenario depends on an uncalibrated axisymmetric alpha-Omega dynamo closure; the jet/no-jet result itself flips with resolution (B12.3.8l vs B12.3.8l-H), so the claim remains conditional until a 3D MRI-based closure check is done.","rationale":"The reader's weakest assumption is the right one. For the advertised scenario to be true, the mean-field alpha-Omega closure must be a reliable stand-in for the real MRI dynamo in a collapsar torus. This is not established by the paper: alpha_d, sigma_c, and rho_cut are hand-set (Table I), and the paper's own runs show a high-resolution jet/no-jet flip (B12.3.8l vs B12.3.8l-H, Sec. III.C) and a magnetosphere that is explicitly described as possibly mis-modeled (Sec. IV). This makes the concern concrete: the central phenomenon, jet launch, is sensitive to numerical and closure choices at the factor-of-few level, while the claimed agreement with observed Ic-BL SNe uses those same runs. I would not call this a reason to reject the paper. The authors are transparent about the limitation, the long-timescale general-relativistic neutrino-MHD framework is substantial, the parameter dependence is explored, and the viscous-hydro comparison gives independent support for the explosion-energy part of the story. The jet-launch part, however, is not independently supported. A 3D MRI-resolving simulation, or at least measured closure coefficients from such a simulation, is the check that would settle whether the dynamo-generated poloidal flux is physically real. Because the concern is an explicit, acknowledged limitation rather than an internal contradiction, the appropriate verdict remains CONDITIONAL, and my stress-test does not move the reader's verdict.","tokens_in":30924,"tokens_out":7741,"duration_ms":77335,"concrete_test":"Run a 3D GRMHD simulation of the same AD35 setup with a weak toroidal seed through MRI saturation in the inner torus; measure the mean-field alpha and turbulent resistivity profiles from the EMF and insert them into the 2D axisymmetric code in place of the Table I constants and Eq. (2) cutoff. Recompute B12.1.8l. The central claim is supported only if the horizon poloidal flux, LBZ, and Eexp remain within the stated factor-of-few uncertainty. Auxiliary: rerun B12.1.8l with alpha_d = 0 and with alpha_d of opposite sign; a no-jet outcome in both cases would show that the result is controlled by the imposed closure rather than by the dynamo physics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that a weak toroidal seed field grows through a dynamo, builds a poloidal field threading the black hole, and thereby launches a Blandford-Znajek jet and an Ic-BL-like explosion—requires that the mean-field alpha-Omega closure in Sec. II (Eq. 2 and Table I) faithfully represents the turbulent dynamo that a real 3D MRI-unstable disk would produce. This is the least secure condition in the argument. The coefficients alpha_d, sigma_c, and rho_cut are prescribed rather than derived from or calibrated to 3D MRI simulations of this system; they set growth rate, dissipation, and spatial cutoff of the dynamo, and the key results vary strongly across Table I. Model B12.3.8l launches a jet at Delta x = 360 m, while the higher-resolution run B12.3.8l-H does not (Table I and Sec. III.C), so jet launch is not numerically converged even within the closure. The authors state that the post-launch magnetosphere may be poorly modeled and that no quasi-steady magnetosphere forms, possibly because of the simple dynamo treatment (Sec. III.D and Sec. IV). The imposed alpha term is spatially uniform apart from the rho_cut factor and has a fixed positive sign, whereas a real MRI dynamo alpha is fluctuating and spatially structured; the polarity reversals seen in Phi_AH (Fig. 8) indicate how sensitive the horizon flux is to stochastic polarity. If the true alpha is weaker, differently structured, or sign-reversing, the ~1e14 G horizon poloidal field and BZ jet may never be established. This is not an internal inconsistency, but it is the load-bearing unvalidated assumption behind the paper's central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports fully general-relativistic, axisymmetric, neutrino-radiation resistive MHD simulations of a collapsar, with a phenomenological mean-field alpha-Omega dynamo term added to mimic MRI-driven turbulence. Starting from a spinning 16 solar-mass black hole and infalling stellar matter with a weak toroidal seed magnetic field, the simulations follow 30-50 s of evolution. In several models a turbulent torus develops, poloidal field threading the black hole builds up, and a Blandford-Znajek jet plus a stellar explosion result, with explosion energies 1e51-1e52 erg, ejecta masses 1.6-4.9 solar masses, and 56Ni masses 0.1-1.1 solar masses, in rough agreement with observed broad-lined type Ic supernovae. The authors also report large Zn production and weak r-process nucleosynthesis in jet-launching models, and they explicitly state that the modeling is qualitative/semi-quantitative and that full 3D simulations are needed for quantitative conclusions.","tokens_in":31314,"tokens_out":5918,"duration_ms":57722,"significance":"If correct, the paper would demonstrate a self-consistent path from a weak toroidal seed field to a Blandford-Znajek jet and an Ic-BL-like explosion, removing the need to assume a pre-existing ~1e14 G poloidal field in collapsar models. This is an important step for the collapsar scenario. The paper's strengths include long-duration full-GR simulations with neutrino radiation transport, a systematic multi-model parameter study, transparent reporting of resolution dependence, and post-processing nucleosynthesis with comparisons to observed supernovae and metal-poor stellar abundances. The central mechanism, however, is conditioned on an uncalibrated mean-field dynamo closure and on resolution-dependent jet/no-jet outcomes; the paper itself acknowledges these limitations. The result should therefore be read as a plausible scenario rather than a demonstrated mechanism, and the observational agreement, while encouraging, does not independently validate the dynamo prescription.","major_comments":[{"comment":"The dynamo closure is load-bearing for the central scenario, yet it is introduced phenomenologically. The coefficients alpha_d, sigma_c, and rho_cut set the growth rate, dissipation, and spatial cutoff of the dynamo, and the results in Table II vary substantially across the small parameter grid (Eexp from 1.33 to 11.6 x 10^51 erg, Mej from 1.56 to 4.90 solar masses, EBZ from 0.22 to 1.59 x 10^51 erg). Because the simulations are axisymmetric, MRI turbulence is not resolved and the closure cannot be validated internally; the references to neutron-star merger simulations (Refs. [30-33]) establish that a dynamo can operate in disks, but not that this specific alpha_d(rho) prescription with fixed positive sign and rho_cut is quantitatively faithful for a collapsar torus. Please provide a sensitivity study of the qualitative conclusions to the sign and spatial structure of alpha_d, or explicitly reframe the central claim as contingent on the closure rather than 'self-consistent'.","section":"Section II, Eq. (2) and Table I"},{"comment":"The jet/no-jet outcome is not converged for a key model. Model B12.3.8l launches a jet at Delta x = 360 m, while the higher-resolution run B12.3.8l-H at Delta x = 300 m does not launch a jet even though the explosion energy is similar (Table II). The authors state that convergence is 'fair' after turbulence develops and attribute the difference to stochastic polarity of the horizon flux (Sec. III.A), but a binary flip from jet to no-jet with resolution means that the central claim 'a jet is launched' is not established for this model. The paper should show that at least one jet-launching model is robust at higher resolution, or identify a resolution-converged diagnostic (e.g., a threshold in horizon poloidal flux) that determines jet formation.","section":"Section III.C and Table I"},{"comment":"The horizon poloidal flux oscillates with sign reversals, and the Poynting luminosity decays on a timescale of about 10 s after the peak; the authors note that no quasi-steady magnetosphere forms, possibly because of the simple dynamo modeling. Since long GRBs typically last 10-100 s, the paper should discuss whether an intermittent, polarity-flipping BZ jet with LBZ of order 10^50 erg/s for only about 10 s is compatible with the duration and energetics of typical long GRBs, or whether the late-time suppression is an artifact of the closure. This is directly relevant to the abstract's claim that the jet luminosity is 'suitable for explaining typical long gamma-ray bursts.'","section":"Section III.D and Fig. 8"},{"comment":"The nucleosynthesis yields, including the claimed large Zn mass and the weak r-process, depend on low-Ye and high-entropy ejecta components that are not converged. For example, model B12.1.8l ejects no material with Ye less than about 0.35, while the higher-resolution run B12.1.8l-H ejects an appreciable low-Ye component, and the 56Ni mass drops from 1.08 to 0.50 solar masses between these two runs. The authors acknowledge this in the text ('we have to keep in mind that the convergence of the numerical results might be poor for other models as well'), but the abstract and conclusions present the 56Ni, Zn, and r-process results as part of the successful scenario. Please either soften these claims to reflect the resolution sensitivity or demonstrate convergence of the nucleosynthesis-relevant ejecta properties.","section":"Section III.E, Fig. 9, and Table III"}],"minor_comments":[{"comment":"In Eq. (8), the sound speed c_s is a velocity, but the denominator in the second factor is written with units of g/cm^3; it should be cm/s (the text preceding the equation gives c_s = 0.1c = 3 x 10^9 cm/s).","section":"Eq. (8)"},{"comment":"Typo: 'Poynging flux' should be 'Poynting flux' in the paragraph describing the funnel structure and jet propagation.","section":"Section III.A"},{"comment":"Typo: 'B12.1.8l abd B12.1.8l-H' should be 'B12.1.8l and B12.1.8l-H' in the first paragraph discussing the Poynting flux evolution.","section":"Section III.D"},{"comment":"Reference [9] and reference [47] are the same paper (Fujibayashi et al., Phys. Rev. D 109, 023031) and should be consolidated.","section":"References"},{"comment":"The model-name key '104alpha_d' is ambiguous; it should be written as '10^4 alpha_d' for clarity.","section":"Table I caption"},{"comment":"The caption says 'Poynting luminosity integrated over simulation time,' but the column EBZ is an energy, not a time-integrated luminosity; consider renaming it 'Poynting energy.'","section":"Table II caption"}],"recommendation":"major_revision","confidential_remarks":"This is a ambitious scenario paper from an experienced group, and the simulations are expensive and carefully described. The main risk is that the abstract's 'self-consistent' and 'in good agreement' framing outruns the evidence, especially given the resolution flip for jet launch and the uncalibrated nature of the dynamo closure. I would support publication after the claims are recalibrated to the closure's limitations and the resolution sensitivity is reflected in the abstract and conclusions, or after additional runs demonstrate robustness of at least one jet-launching model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nBottom line: this is the first collapsar simulation I know of that starts from a weak toroidal seed and ends with a BZ jet plus an Ic-BL-like explosion without putting a strong poloidal field in by hand. That is genuinely new, and the authors are straight about what is and isn't modeled. I think the paper deserves a serious referee; it just shouldn't be read as the last word.\n\nWhat it does well: the setup is careful—stellar evolution model, full GR, neutrino transport, 30–50s evolutions, a small parameter study across alpha_d, sigma_c, rho_cut, and resolution, plus nucleosynthesis post-processing. The Eexp–Mej and MNi correlations land on the observed Ic-BL points, and the authors stress this is scatter from dynamo parameters, not fits. The discussion of why the flux gets swallowed early and only later can a magnetosphere form is plausible and useful. Credit where due: the long runtimes matter, and the qualitative scenario—torus dynamo, flux accretion, late-time magnetosphere, joint jet/explosion—is coherent.\n\nThe soft spots are real and the stress-test note has the right emphasis. The load-bearing piece is the alpha-Omega closure: alpha_d, sigma_c, and rho_cut are hand-picked, the alpha term is spatially uniform and fixed-sign, and the whole scheme is axisymmetric, so there is no resolved MRI to validate it. The resolution flip is the most concrete problem: B12.3.8l launches a jet at 360 m and does not at 300 m. The authors report it honestly, but it means jet/no-jet is not converged even within the closure. They also note the magnetosphere after launch is likely poorly modeled and that no quasi-steady magnetosphere forms; the polarity reversals in Phi_AH show how fragile the horizon flux is. I don't see an internal contradiction, and the authors flag most of these limitations themselves. But the central claim should be labeled conditional: the dynamo closure could be wrong in sign or structure, and then the 10^14 G horizon field and the jet vanish. The quantitative SN comparisons also lack error bars; factor-of-a-few uncertainties are acknowledged.\n\nWho benefits: people working on collapsar engines, GRB-SN connections, and long-term BH-torus evolution. A referee should push for a clearer statement that the jet result is not converged, and ideally for a 3D check with resolved turbulence, even if shorter. Code and data are not released, which is a minus but not disqualifying.\n\nRecommendation: send it to review. It is a serious, honest simulation study with a genuinely new scenario; it just needs to be read as a proof-of-feasibility with a phenomenological closure.","headline":"A credible but not yet converged collapsar engine: the dynamo closure does real work, and the jet/no-jet split across resolution keeps the headline result conditional.","tokens_in":31904,"tokens_out":2173,"would_cite":true,"duration_ms":22370,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.30.Qd","97.60.Bw","98.70.Rz"],"model":"deepseek-v4-flash","headline":"A collapsar can build the magnetic field for its own gamma-ray-burst jet, starting from a weak seed.","keywords":["collapsar","long gamma-ray bursts","broad-lined type Ic supernovae","general relativistic magnetohydrodynamics","mean-field dynamo","Blandford-Znajek mechanism","nucleosynthesis","black hole accretion"],"falsifier":"A three-dimensional general-relativistic magnetohydrodynamics simulation of the same progenitor starting from the same weak toroidal seed, with resolution high enough to resolve the magnetorotational instability, would settle the claim: if no $\\sim10^{14}$ G poloidal field threads the black hole and no jet appears within about 30 seconds, the dynamo closure is carrying the result.","tokens_in":30705,"feed_emoji":"💥","tokens_out":13916,"duration_ms":137459,"temperature":0.7,"pith_summary":"Long gamma-ray bursts and the broad-lined type Ic supernovae that accompany them are usually modeled by assuming a spinning black hole already threaded by a very strong poloidal magnetic field. This paper argues that no such assumption is needed. In a fully general-relativistic simulation of a collapsar—a rotating massive star whose core collapses directly to a black hole surrounded by a torus—a mean-field dynamo (a prescription for how turbulence regenerates magnetic fields) amplifies a weak toroidal seed field, builds a poloidal field through the black hole, and launches a Blandford-Znajek jet, a jet powered by extraction of the black hole's spin energy through magnetic fields, with the luminosity of typical long GRBs, while the same torus drives a stellar explosion. The simulated explosion energies, ejecta masses, and $^{56}$Ni masses fall on the correlations measured for broad-lined type Ic supernovae, so the strong field that powers the jet is generated by the disk itself rather than inserted by hand.","feed_headline":"Weak seed field grows into a GRB jet in collapsar simulation","feed_subtitle":"A dynamo in the black hole's disk builds the 10^14-gauss field earlier models assumed by hand.","key_machinery":"The load-bearing object is the phenomenological mean-field $\\alpha\\Omega$ dynamo, inserted into general-relativistic resistive magnetohydrodynamics as a source term with amplitude $\\alpha_d$ and conductivity $\\sigma_c$, suppressed by a factor $1-\\exp(-\\rho/\\rho_{\\rm cut})$ in low-density regions. It stands in for the turbulent dynamo that a fully resolved three-dimensional magnetorotational-instability simulation would produce, sustaining both field amplification and the turbulent angular-momentum transport that makes the torus 'viscous.' The second essential element is the jet-launch criterion $B^2/8\\pi > \\rho_{\\rm inf}v_{\\rm inf}^2$: once magnetic pressure in the polar region beats the ram pressure of infalling matter, a magnetosphere inflates and the Blandford-Znajek mechanism can operate. Because the dynamo reverses field polarity quasi-periodically, reconnection in the magnetosphere limits the extracted energy to roughly $10^{51}$ erg, giving the jet a natural lifetime of order 100 s.","core_discovery":"The central claim is that a collapsar with only a weak toroidal magnetic seed self-consistently produces both a long-GRB jet and a supernova explosion. Once a massive torus forms, the dynamo amplifies the field toward equipartition, giving $\\sim10^{14}$ G near the black hole. At early times the ram pressure of infalling matter swallows magnetic flux into the horizon before a magnetosphere can grow; only when the ram pressure drops does the condition $B^2/8\\pi > \\rho_{\\rm inf}v_{\\rm inf}^2$ hold. Then a horizon-threading poloidal field is established, the Blandford-Znajek mechanism extracts black-hole spin energy as a Poynting-flux jet with $L_{BZ}\\sim10^{50}$–$10^{51}$ erg/s, and the turbulent torus plus magnetocentrifugal effects drive an explosion with $E_{\\rm exp}\\sim10^{51}$–$10^{52}$ erg, $M_{\\rm ej}\\sim1.6$–$4.9\\,M_\\odot$, and $M_{\\rm Ni}\\sim0.1$–$1.1\\,M_\\odot$, matching observed broad-lined type Ic supernovae. The paper also finds that jet activity is required for even a weak $r$-process and that large amounts of zinc are synthesized.","pith_inferences":["Inference: if the dynamo closure is approximately faithful, the observed scatter in $E_{\\rm exp}$ and $M_{\\rm ej}$ among type Ic-BL supernovae could be driven as much by dynamo activity as by progenitor angular momentum; varying $\\alpha_d$ and $\\sigma_c$ in the simulation already produces a spread similar to the observations.","Inference: the paper's negligible lanthanide and actinide yields imply that red kilonova-like transients associated with long GRBs, such as the reported GRB 230307A event, more likely come from neutron-star mergers; future events with tellurium features can test this.","Inference: if the local magnetosphere with $|\\phi_{\\rm AH}|\\sim5$ suffices for jet launch, the higher global magnetically arrested disk (MAD) threshold often used in simulations may overstate the field required for GRB jets; this can be checked with higher-resolution runs.","Inference: because the dynamo's sign, magnitude, and spatial structure are chosen rather than derived, a three-dimensional simulation resolving the magnetorotational instability is the decisive check, and until then the quantitative Poynting luminosity carries a factor-of-a-few uncertainty."],"forward_implications":["Long gamma-ray bursts do not require a pre-existing $\\sim10^{14}$ G poloidal field; the torus dynamo builds the jet-launching field from a weak seed in roughly 10–20 s.","The same torus that feeds the jet drives the stellar explosion, so the $E_{\\rm exp}$–$M_{\\rm ej}$–$M_{\\rm Ni}$ correlations of GRB-associated broad-lined type Ic supernovae emerge from a single calculation with no fine-tuning.","Jet activity is the switch for trans-iron nucleosynthesis: models with jets make modest $r$-process nuclei and sizable zinc masses, while non-jet models make essentially none.","The jet energy is capped at the order of $10^{51}$ erg because dynamo polarity flips and reconnection erode the magnetosphere, providing a natural timescale of order 100 s for the jet lifetime."],"supporting_citations":[{"why":"Previous ideal-MHD collapsar simulations that this work directly extends; they motivated the need for a dynamo-built poloidal field.","marker":"[22]"},{"why":"Introduces the resistive-MHD mean-field dynamo formulation whose parameters and source term are adopted here.","marker":"[29]"},{"why":"Viscous-hydrodynamics collapsar models with the same progenitor setup that serve as the comparison explosion engine.","marker":"[9]"},{"why":"Stellar-evolution progenitor model used to construct the initial spinning black hole plus infalling matter.","marker":"[23]"},{"why":"The Blandford-Znajek mechanism that converts black-hole spin into the Poynting-flux jet.","marker":"[6]"},{"why":"Observational sample of broad-lined type Ic supernovae used to compare ejecta mass, explosion energy, and nickel mass.","marker":"[46]"},{"why":"Definition of the magnetically arrested state and the reference for the MADness parameter threshold.","marker":"[16]"},{"why":"Long-term black hole-neutron star merger simulations with MRI dynamo-built magnetospheres and jets that motivate the dynamo treatment.","marker":"[30]"}],"fun_headline_variants":["Self-consistent collapsar: weak seed field yields jet and supernova","Dynamo turns weak seed into GRB jet and supernova explosion","Weak seed magnetism grows a jet and explosion in collapsars","Collapsar dynamo builds jet and supernova from weak seed field","Self-consistent collapsar scenario: jet and supernova from dynamo"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scenario stands or falls on the assumption that the hand-tuned two-dimensional dynamo term reproduces how real three-dimensional turbulence amplifies magnetic fields in the torus.","fun_headline_variants_meta":{"raw":{"variants":["Self-consistent collapsar: weak seed field yields jet and supernova","Dynamo turns weak seed into GRB jet and supernova explosion","Weak seed magnetism grows a jet and explosion in collapsars","Collapsar dynamo builds jet and supernova from weak seed field","Self-consistent collapsar scenario: jet and supernova from dynamo"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000621,"raw_usage":{"total_tokens":2998,"prompt_tokens":1182,"completion_tokens":1816,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":798,"completion_tokens_details":{"reasoning_tokens":1720}},"tokens_in":798,"tokens_out":1816,"duration_ms":12633,"temperature":1.0,"reasoning_tokens":1720,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T13:25:25.349964+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A three-dimensional general-relativistic magnetohydrodynamics simulation of the same progenitor starting from the same weak toroidal seed, with resolution high enough to resolve the magnetorotational instability, would settle the claim: if no $\\sim10^{14}$ G poloidal field threads the black hole and no jet appears within about 30 seconds, the dynamo closure is carrying the result.","supporting_citations":[{"cited_title":"The Connection between Gamma-Ray Bursts and Extremely Metal-Poor Stars: Black Hole-forming Supernovae with Relativistic Jets","cited_arxiv_id":"astro-ph/0702471","evidence_quote":"Previous ideal-MHD collapsar simulations that this work directly extends; they motivated the need for a dynamo-built poloidal field."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Blandford-Znajek mechanism that converts black-hole spin into the Poynting-flux jet."},{"cited_title":"Shibata, Numerical Relativity (World Scientific Pub- lishing Company, 2016)","cited_arxiv_id":null,"evidence_quote":"Observational sample of broad-lined type Ic supernovae used to compare ejecta mass, explosion energy, and nickel mass."}],"review_version":1}