Explodability matters: how realistic neutrino-driven explosions change explosive nucleosynthesis yields
Pith reviewed 2026-05-18 06:39 UTC · model grok-4.3
The pith
Realistic neutrino-driven supernova models produce lower iron-peak element yields than piston and bomb prescriptions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Piston and bomb models artificially increase the production of Fe-peak elements relative to neutrino-driven explosions, while differences in which progenitors explode create discrepancies primarily in the lighter elements.
What carries the argument
The GR1D+ spherically symmetric model with state-of-the-art neutrino transport and time-dependent mixing-length convection, which sets realistic explosion energy and mass cut before SkyNet yield calculation.
If this is right
- Galactic chemical evolution calculations should replace piston and bomb yields with neutrino-driven results to avoid overestimating iron-group elements.
- Explodability differences among progenitors shift production of oxygen, silicon, and other lighter elements more than simple models indicate.
- Yield tables for rotating, low-metallicity, and binary stars become available with energies and remnant masses matching type II-P, IIb, and Ib observations.
Where Pith is reading between the lines
- Incorporating these yields into chemical evolution codes could change predicted metallicity patterns in dwarf galaxies.
- Three-dimensional explosion simulations may refine the element-by-element differences found here.
- Binary evolution channels could amplify the lighter-element discrepancies if their explodability differs from single-star cases.
Load-bearing premise
The one-dimensional-plus model with mixing-length convection treatment captures the essential multi-dimensional dynamics and neutrino transport that fix the final explosion energy and mass cut.
What would settle it
Direct comparison of predicted nickel-56 or calcium yields against observed abundances in supernova remnants or metal-poor stars would show whether realistic models match better than piston and bomb results.
Figures
read the original abstract
Explosive nucleosynthesis is affected by many uncertainties, particularly regarding assumptions and prescriptions adopted during the evolution of the star. Moreover, simple explosion models are often used in the literature, which can introduce large errors in the assumed explosion energy and mass cut. In this paper, our goal is to analyze the explosion properties and nucleosynthesis of a large range of progenitors from three different stellar evolution codes: FRANEC, KEPLER, and MESA. In particular, we will show the differences between the neutrino-driven explosions simulated in this work with the much simpler bomb and piston models that are typically widely used in the literature. We will then focus on the impact of different explodabilities and different explosion dynamics on the nucleosynthetic yields. We adopt the neutrino-driven core-collapse supernova explosion code GR1D+, i.e. a spherically symmetric model with state-of-the-art microphysics and neutrino transport and a time-dependent mixing-length model for neutrino-driven convection. We carry out explosions up to several seconds after bounce, and then calculate the nucleosynthetic yields with the post-processing code SkyNet. We find that our 1D+ simulations yield explosion energies and remnant masses in agreement with observations of type II-P, IIb, and Ib supernovae, as well as with the most recent 3D simulations of the explosion. We provide a complete set of yields for all the stars simulated, including rotating, low-metallicity, and binary progenitors. Finally, we find that piston and bomb models, compared to more realistic neutrino-driven explosions, can artificially increase the production of Fe-peak elements, whereas the different explodability tends to cause discrepancies in the lighter elements.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper simulates neutrino-driven core-collapse supernovae for progenitors from FRANEC, KEPLER, and MESA using the GR1D+ code (1D with state-of-the-art microphysics, neutrino transport, and time-dependent mixing-length convection). Explosions are evolved to several seconds post-bounce, nucleosynthesis yields are post-processed with SkyNet, and results are compared to standard piston and bomb models. The central claims are that GR1D+ yields match observed Type II-P/IIb/Ib energies and remnant masses as well as recent 3D simulations, that piston/bomb models artificially boost Fe-peak production, and that explodability differences primarily affect lighter-element yields. A full yield table for rotating, low-metallicity, and binary progenitors is provided.
Significance. If the central comparisons hold, the work supplies a large, observationally anchored set of yields from more realistic neutrino-driven explosions and quantifies systematic offsets introduced by the piston/bomb approximations that dominate the literature. This directly improves inputs for galactic chemical evolution and supernova yield libraries, while the reported agreement with 3D results and observed energies strengthens the practical utility of the 1D+ dataset.
major comments (1)
- [Comparison of explosion models and nucleosynthesis yields] The section comparing piston/bomb models to GR1D+ runs does not state whether the simple models were executed with the identical final explosion energy and remnant mass (or mass cut) that each GR1D+ progenitor produces. Without per-progenitor matching of these boundary conditions, the reported Fe-peak overproduction in piston/bomb models cannot be unambiguously attributed to the artificial nature of those prescriptions rather than to inconsistent E_exp or M_rem values. This matching is load-bearing for the Fe-peak part of the central claim.
minor comments (3)
- [Introduction and results] The abstract and introduction refer to 'different explodabilities' driving lighter-element discrepancies, but the manuscript should add a short table or figure that isolates the explodability effect (e.g., by holding E_exp and mass cut fixed across models).
- [Methods and figures] Figure captions and text should explicitly note the post-bounce time at which explosion energy and remnant mass are measured for the GR1D+ runs, to facilitate direct comparison with piston/bomb setups.
- [Nucleosynthesis post-processing] A brief statement on the sensitivity of the SkyNet post-processing to the adopted mass cut and fallback prescription would help readers assess robustness of the lighter-element trends.
Simulated Author's Rebuttal
We thank the referee for their thorough review and for identifying this key point about the model comparisons. We have revised the manuscript to explicitly address the concern and strengthen the clarity of our central claims.
read point-by-point responses
-
Referee: The section comparing piston/bomb models to GR1D+ runs does not state whether the simple models were executed with the identical final explosion energy and remnant mass (or mass cut) that each GR1D+ progenitor produces. Without per-progenitor matching of these boundary conditions, the reported Fe-peak overproduction in piston/bomb models cannot be unambiguously attributed to the artificial nature of those prescriptions rather than to inconsistent E_exp or M_rem values. This matching is load-bearing for the Fe-peak part of the central claim.
Authors: We agree that explicit per-progenitor matching of explosion energy and remnant mass is essential for isolating the effects of the explosion mechanism. In our work, the piston and bomb models were run with the identical final explosion energy and mass cut (remnant mass) obtained from each GR1D+ simulation for the corresponding progenitor. This matching was performed to ensure a fair comparison and to attribute differences in Fe-peak yields specifically to the artificial nature of the piston/bomb prescriptions rather than to differences in E_exp or M_rem. However, we acknowledge that the original manuscript did not state this matching procedure with sufficient clarity. We have revised Section 3.2 (and the associated figure captions) to explicitly describe how the boundary conditions were matched on a per-progenitor basis and have added a sentence confirming that the GR1D+ values of E_exp and M_rem were adopted directly for the piston and bomb runs. With this clarification, the reported Fe-peak overproduction is now unambiguously linked to the simplified explosion prescriptions. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper performs forward 1D+ hydrodynamic simulations of neutrino-driven explosions with GR1D+ (including time-dependent mixing-length convection and state-of-the-art neutrino transport), followed by post-processing nucleosynthesis yields via SkyNet. Explosion energies and remnant masses are reported to agree with external observational benchmarks (Type II-P, IIb, Ib supernovae) and recent 3D simulations, rather than being fitted or defined in terms of the nucleosynthesis outputs. Direct comparisons to piston and bomb models are executed as separate runs; no equations, parameters, or self-citations reduce the reported yield differences (Fe-peak or lighter elements) to inputs by construction. The derivation chain remains self-contained against external benchmarks with no load-bearing self-referential steps.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The time-dependent mixing-length model for neutrino-driven convection in GR1D+ adequately represents the multi-dimensional convective engine that determines explodability.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
piston and bomb models, compared to more realistic neutrino-driven explosions, can artificially increase the production of Fe-peak elements
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.
Forward citations
Cited by 1 Pith paper
-
Production of heavy $\alpha$-elements and $^{44}$Ti in Cas A: comparison to abundances from 1D core-collapse supernova models and evidence for Carbon-Oxygen shell mergers
Core-collapse supernova models including C-O shell mergers best match observed elemental ratios in Cas A, indicating mergers occur and contribute up to 20-30% of 44Ti outside the reverse shock.
Reference graph
Works this paper leans on
-
[1]
Abdikamalov, E., Ott, C. D., Radice, D., et al. 2015, ApJ, 808, 70
work page 2015
-
[2]
Adams, S. M., Kochanek, C. S., Gerke, J. R., & Stanek, K. Z. 2017, MNRAS, 469, 1445
work page 2017
-
[3]
Andrassy, R., Herwig, F., Woodward, P., & Ritter, C. 2020, MNRAS, 491, 972
work page 2020
-
[4]
Andresen, H., O’Connor, E. P., Andersen, O. E., & Couch, S. M. 2024, A&A, 687, A55
work page 2024
- [5]
- [6]
- [7]
- [8]
-
[9]
Arnett, W. D. 1969, ApJ, 157, 1369
work page 1969
-
[10]
Arnett, W. D. 1982, ApJ, 253, 785
work page 1982
-
[11]
Arnett, W. D. 1987, ApJ, 319, 136
work page 1987
-
[12]
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
work page 2009
-
[13]
B., Baron, E., & Thielemann, F
Aufderheide, M. B., Baron, E., & Thielemann, F. K. 1991, ApJ, 370, 630
work page 1991
-
[14]
Battino, U., Roberti, L., Lawson, T. V ., Laird, A. M., & Todd, L. 2024, Universe, 10, 204
work page 2024
-
[15]
Blinnikov, S. I., Eastman, R., Bartunov, O. S., Popolitov, V . A., & Woosley, S. E. 1998, ApJ, 496, 454
work page 1998
- [16]
- [17]
-
[18]
Boccioli, L., Mathews, G. J., Suh, I.-S., & O’Connor, E. P. 2022, ApJ, 926, 147
work page 2022
- [19]
-
[20]
Boccioli, L., Roberti, L., Limongi, M., Mathews, G. J., & Chieffi, A. 2023, ApJ, 949, 17
work page 2023
-
[21]
Bohm-Vitense, E. 1958, Z. Astrophys., 46, 108
work page 1958
-
[22]
Bruenn, S. W. 1985, ApJS, 58, 771
work page 1985
- [23]
-
[24]
Buras, R., Rampp, M., Janka, H. T., & Kifonidis, K. 2006, A&A, 447, 1049
work page 2006
-
[25]
Burbidge, E. M., Burbidge, G. R., Fowler, W. A., & Hoyle, F. 1957, RvMP, 29, 547
work page 1957
-
[26]
Burrows, A., Reddy, S., & Thompson, T. A. 2006, NuPhA, 777, 356
work page 2006
- [27]
- [28]
-
[29]
Chan, C., Müller, B., Heger, A., Pakmor, R., & Springel, V . 2018, ApJL, 852, L19
work page 2018
-
[30]
Couch, S. M. & Ott, C. D. 2015, ApJ, 799, 5
work page 2015
- [31]
- [32]
-
[33]
Diehl, R. & Prantzos, N. 2020, in Handbook of nuclear physics, ed. I. Tanihata, H. Toki, & T. Kajino (Singapore: Springer Nature Singapore), 1–83
work page 2020
- [34]
-
[35]
Ebinger, K., Curtis, S., Ghosh, S., et al. 2020, ApJ, 888, 91 Eggenberger Andersen, O., O’Connor, E., Andresen, H., da Silva Schneider, A., & Couch, S. M. 2025, ApJ, 980, 53 Eggenberger Andersen, O., Zha, S., da Silva Schneider, A., et al. 2021, ApJ, 923, 201
work page 2020
-
[36]
Ertl, T., Janka, H. T., Woosley, S. E., Sukhbold, T., & Ugliano, M. 2016, ApJ, 818, 124
work page 2016
-
[37]
Falla, A., Roberti, L., Limongi, M., & Chieffi, A. 2025, ApJ, 991, 21
work page 2025
- [38]
-
[39]
Farmer, R., Laplace, E., de Mink, S. E., & Justham, S. 2021, ApJ, 923, 214
work page 2021
-
[40]
Farmer, R., Laplace, E., Ma, J.-z., de Mink, S. E., & Justham, S. 2023, ApJ, 948, 111
work page 2023
-
[41]
Fischer, T., Bastian, N.-U., Blaschke, D., et al. 2017, PASA, 34, e067
work page 2017
-
[42]
Fischer, T., Guo, G., Langanke, K., et al. 2024, PrPNP, 137, 104107
work page 2024
-
[43]
J., Lucente, G., Padilla-Gay, I., & Patwardhan, A
Friedland, A., Li, D. J., Lucente, G., Padilla-Gay, I., & Patwardhan, A. V . 2025, arXiv e-prints, arXiv:2508.02055 Fröhlich, C., Hauser, P., Liebendörfer, M., et al. 2006a, ApJ, 637, 415 Fröhlich, C., Martínez-Pinedo, G., Liebendörfer, M., et al. 2006b, PhRvL, 96, 142502
-
[44]
2017, in Handbook of supernovae, ed
Gal-Yam, A. 2017, in Handbook of supernovae, ed. A. W. Alsabti & P. Murdin (Cham: Springer International Publishing), 195–237
work page 2017
- [45]
- [46]
- [47]
-
[48]
Hix, W. R. & Thielemann, F. K. 1999, Journal of Computational and Applied Mathematics, 109, 321
work page 1999
-
[49]
Horowitz, C. J. 2002, PhRvD, 65, 043001
work page 2002
-
[50]
Horowitz, C. J., Caballero, O. L., Lin, Z., O’Connor, E., & Schwenk, A. 2017, PhRvC, 95, 025801
work page 2017
-
[51]
Imasheva, L., Janka, H.-T., & Weiss, A. 2023, MNRAS, 518, 1818
work page 2023
- [52]
-
[53]
2017, in Handbook of supernovae, ed
Janka, H.-T. 2017, in Handbook of supernovae, ed. A. W. Alsabti & P. Murdin (Cham: Springer International Publishing), 1095–1150
work page 2017
- [54]
- [55]
- [56]
-
[57]
V ., Pignatari, M., Stancliffe, R
Lawson, T. V ., Pignatari, M., Stancliffe, R. J., et al. 2022, MNRAS, 511, 886
work page 2022
- [58]
-
[59]
Li, X. 2024, Master’s thesis, Stockholm University, Department of Astronomy/ Stockholm University, Department of Astronomy Liebendörfer, M., Whitehouse, S. C., & Fischer, T. 2009, ApJ, 698, 1174
work page 2024
- [60]
- [61]
- [62]
- [63]
-
[64]
Litvinova, I. Y . & Nadezhin, D. K. 1985, Soviet Astronomy Letters, 11, 145
work page 1985
- [65]
- [66]
- [67]
- [68]
-
[69]
Meyer, B. S., Mathews, G. J., Howard, W. M., Woosley, S. E., & Hoffman, R. D. 1992, ApJ, 399, 656
work page 1992
-
[70]
Mezzacappa, A., Endeve, E., Messer, O. E. B., & Bruenn, S. W. 2020, LRCA, 6, 4
work page 2020
-
[71]
2024, ApJL, 965, L27 Müller, B
Milisavljevic, D., Temim, T., De Looze, I., et al. 2024, ApJL, 965, L27 Müller, B. 2020, LRCA, 6, 3 Müller, B., Janka, H.-T., & Marek, A. 2012, ApJ, 756, 84 Article number, page 16 of 23 Boccioli & Roberti: Neutrino-driven explosive nucleosynthesis Müller, B., Melson, T., Heger, A., & Janka, H.-T. 2017, MNRAS, 472, 491 Müller, B., Tauris, T. M., Heger, A....
work page 2024
-
[72]
2015, PASJ, 67, 107 Navó, G., Reichert, M., Obergaulinger, M., & Arcones, A
Nakamura, K., Takiwaki, T., Kuroda, T., & Kotake, K. 2015, PASJ, 67, 107 Navó, G., Reichert, M., Obergaulinger, M., & Arcones, A. 2023, ApJ, 951, 112
work page 2015
-
[73]
Neustadt, J. M. M., Kochanek, C. S., Stanek, K. Z., et al. 2021, MNRAS, 508, 516
work page 2021
- [74]
-
[75]
2013, ARA&A, 51, 457 O’Connor, E
Nomoto, K., Kobayashi, C., & Tominaga, N. 2013, ARA&A, 51, 457 O’Connor, E. 2015, ApJS, 219, 24 O’Connor, E. & Ott, C. D. 2010, CQGra, 27, 114103
work page 2013
-
[76]
Patton, R. A. & Sukhbold, T. 2020, MNRAS, 499, 2803
work page 2020
- [77]
- [78]
- [79]
-
[80]
Popov, D. V . 1993, ApJ, 414, 712
work page 1993
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.