pith. sign in

arxiv: 2605.19968 · v1 · pith:ZRSDB6KZnew · submitted 2026-05-19 · 🌌 astro-ph.SR · astro-ph.HE

Simulating the Convective Urca Process with Multiple Urca Pairs in a Simmering White Dwarf

Brendan Boyd , Ferran Poca-Amor\'os , Alan Calder , Dean M. Townsley This is my paper

Pith reviewed 2026-05-20 04:22 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.HE
keywords convective Urca processwhite dwarf convectionType Ia supernovaesimmering phasehydrodynamic simulationsUrca pairscarbon burningneutrino emission
0
0 comments X

The pith

Convective Urca reactions reduce mixing efficiency near the white dwarf convection boundary without shrinking the zone, and the A=23 pair dominates.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

During the simmering phase before a Type Ia supernova, carbon burning in a near-Chandrasekhar white dwarf drives convection while Urca weak reactions link neutrino losses to compositional changes. Three-dimensional low-Mach simulations were run with a carbon-burning network that includes the A=21, A=23, and A=25 Urca pairs, and runs with and without these reactions were compared to isolate their effects on the flow. The convective Urca process lowers mixing efficiency close to the convective boundary but leaves the overall size of the convection zone unchanged. The A=23 pair contributes the largest share of the effect among the three. These results clarify how the pre-explosion convection zone behaves when neutrino-emitting reactions operate inside the convective region.

Core claim

In 3D hydrodynamic simulations of the convection zone in a simmering white dwarf that include a comprehensive carbon-burning network along with the A=21, A=23, and A=25 Urca pairs, the convective Urca process reduces the efficiency of convective mixing near the convective boundary but does not restrict the size of the convection zone, and the A=23 Urca pair is the most important to the process.

What carries the argument

The convective Urca process, the coupling of weak nuclear reactions (beta decays and electron captures that emit neutrinos) with convective mixing in the white dwarf core.

If this is right

  • The overall extent of the convection zone in the simmering white dwarf stays the same when Urca reactions are active.
  • Convective mixing becomes less efficient in the layers immediately adjacent to the boundary.
  • The A=23 Urca pair produces the strongest influence on the process compared with the A=21 and A=25 pairs.
  • These boundary-layer changes can be used to refine one-dimensional models of white-dwarf evolution before explosion.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Altered mixing near the boundary could change the local carbon depletion rate and thereby shift the timing or location of ignition.
  • Stellar evolution codes that treat convection with mixing-length theory may need revised prescriptions at Urca-active boundaries.
  • Repeating the simulations at different central densities or with additional Urca pairs could reveal whether the dominance of A=23 persists under other conditions.

Load-bearing premise

The low-Mach hydrodynamic approximation and chosen nuclear network accurately capture the interaction between weak reactions and convective flows without dominant numerical artifacts at the boundary.

What would settle it

A simulation or observation showing identical mixing efficiency and velocity profiles right at the convective boundary whether or not the Urca pairs are included would falsify the central result.

Figures

Figures reproduced from arXiv: 2605.19968 by Alan Calder, Brendan Boyd, Dean M. Townsley, Ferran Poca-Amor\'os.

Figure 1
Figure 1. Figure 1: — A graphical visualization of the nuclear reaction network. Each node represents a different isotope in the network. Arrows represent reactions connecting two isotopes, with the direction of the arrow indicating the direction of the reaction. The horizontal axis indicates the proton number, Z. The vertical axis indicates the number of excess neutrons (positive) or protons (negative). Helium, protons, and … view at source ↗
Figure 2
Figure 2. Figure 2: — Temperature profiles vs. radius for a series of initial models. The blue curve represents the initial profile used in the presented simulations. The green curves represents the initial pro￾file used in PA2026, solid representing the FN1/NB1 model and dashed representing the FN2 model. The grey curve represents the initial profile used in B2025. The vertical black lines indicate the location of the A=21 (… view at source ↗
Figure 4
Figure 4. Figure 4: — The top plot shows the angle-averaged X( 20Ne) vs. radial bin. The bottom plot shows the rms velocity vs. radial bin. The light purple indicates the t=0 initial conditions. The dark purple curve represents the simulation at the end of the ini￾tial evolution, about t = 885 s. The vertical black lines indicate the location of the A=21 (dashed), A=23 (dot-dashed), and A=25 (solid) Urca shells. Overall, afte… view at source ↗
Figure 5
Figure 5. Figure 5: — The two plots tracks the mass enclosed by the convection zone (left) and rms velocity in the convection zone (right) over time. The grey curve shows the initial evolution, to t = 885 s, prior to evolving the base state. The blue curve represent the FN simulation and the orange curve represents the NB simulation. 1000 1500 2000 2500 3000 3500 Time (s) 2.8 3.0 3.2 3.4 3.6 12C Burn Rate (10 − 4 M hr −1) Ful… view at source ↗
Figure 6
Figure 6. Figure 6: — Rate of 12C mass burned in the simulations over time. The blue curves represents the FN simulation while the orange curves represents the NB simulation. The shaded regions corre￾spond to the time intervals with which we calculate time-averages of some parameters. using centered differences as described in Boyd et al. (2025a). In [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: — Slices through the center of the white dwarf, zoomed in on the convection zone. The colorbar represents the mass fraction of 20Ne. The left half of the slice shows the FN simulation and the right half shows the NB simulation. The three white circles indicate the location of the three Urca shells, A=21, A=23, and A=25 (from innermost to outermost). to beta-decay products. This transition does not oc￾cur r… view at source ↗
Figure 8
Figure 8. Figure 8: — Slices through the center of the white dwarf, zoomed in on the convection zone. The colorbar represents the magnitude of the tangential velocity. The left half of each slice shows the FN simulation and the right half shows the NB simulation. The three white circles indicate the location of the three Urca shells, A=21, A=23, and A=25 (from innermost to outermost). In comparison, the network presented in S… view at source ↗
Figure 11
Figure 11. Figure 11: — Density-weighted angle-averaged profiles of the ratio of convective gradients vs. binned radius. The dashed curves relate to the Schwarzschild criterion ∇/∇ad. The solid curves relate to the Ledoux criterion ∇/∇Led. The blue curves represent the FN simulation, while the orange represent the NB simulation. Each curve represents the time averaged value for the given simulation over the last ∼600 s, see [… view at source ↗
Figure 12
Figure 12. Figure 12: , further demonstrate the key differences in the structure of the two convection zones. The difference in the central entropy between the two simulations is due to differences in temperature and composition, which still largely reflects the initial composition seen in [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: — Density-weighted angle-averaged composition profiles for the A=21 (left), A=23 (center), and A=25 (right) Urca pairs vs. binned radius. Solid curves represent the relatively neutron-poor isotope (product of beta-decay) and the dashed curves represent the relatively neutron-rich isotope (product of electron capture). The blue curves represent the FN simulation, while the orange represent the NB simulatio… view at source ↗
Figure 15
Figure 15. Figure 15: — The total energy generation rate from nuclear re￾actions per radial bin vs. radial bin. The blue curves represent the FN simulation, while the orange represent the NB simulation. Each curve represents the time averaged value for the given simu￾lation over the last ∼600 s, see [PITH_FULL_IMAGE:figures/full_fig_p013_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: — The rms velocity profiles vs. binned radius. The blue curve represents the time averaged value for the FN simulation over the last ∼600 s, see [PITH_FULL_IMAGE:figures/full_fig_p014_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: — Density-weighted angle-averaged composition profiles for the A=23 Urca pair vs. binned radius. The blue curve repre￾sents the time averaged value for the FN simulation over the last ∼600 s, see [PITH_FULL_IMAGE:figures/full_fig_p014_17.png] view at source ↗
read the original abstract

Type Ia supernovae are bright thermonuclear explosions of one or more white dwarf stars. The exact origin and explosion mechanism for these supernovae is still poorly understood. In the near-Chandrasekhar mass progenitor model, a simmering phase precedes the explosion. During this simmering phase, central carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, a linking of weak nuclear reactions and convective mixing. Convective Urca has the potential to alter characteristics of the convection zone and thus alter the evolution of the white dwarf. To study the convective Urca process, we use the low Mach number hydrodynamic code MAESTROeX to run 3D simulations of the convection zone. We build off previous work to implement a more comprehensive carbon burning network and include the A=21, A=23, and A=25 Urca pairs in the simulations. We compare simulations with and without the convective Urca process to isolate the direct effects the process has on the convection zone. We find the convective Urca process reduces the efficiency of convective mixing near the the convective boundary, but does not restrict the size of the convection zone. We additionally find the A=23 Urca pair to be the most important Urca pair to the convective Urca process in these simulations. All together, our results better inform our understanding of this complex phenomena as well as demonstrates the range of potential convective structures, particularly at the convective boundary, of a simmering white dwarf.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper uses the low-Mach number code MAESTROeX to perform 3D simulations of convection in a simmering near-Chandrasekhar white dwarf, incorporating an extended carbon-burning network and the A=21, A=23, and A=25 Urca pairs. Direct comparisons between otherwise identical runs with and without the Urca process are used to isolate its effects; the central claims are that the convective Urca process reduces mixing efficiency near the convective boundary without shrinking the convection zone and that the A=23 pair dominates the process.

Significance. If the reported differences survive variations in boundary treatment and resolution, the work supplies concrete numerical evidence on how weak reactions couple to convective flows in the simmering phase, helping constrain the pre-explosion structure of Type Ia supernova progenitors. The inclusion of multiple Urca pairs and the with/without comparison approach are strengths that allow a direct assessment of relative importance.

major comments (1)
  1. [simulation setup and comparison methodology] Simulation setup and comparison methodology section: the headline result that Urca reduces mixing efficiency near the boundary rests on the assumption that MAESTROeX's low-Mach boundary treatment and any implicit filtering do not themselves suppress mixing at the interface. Without explicit tests (e.g., altered boundary formulations, grid stretching variations, or resolution studies that quantify the mixing-efficiency metric), it remains possible that the reported reduction is at least partly numerical rather than physical.
minor comments (1)
  1. [abstract] Abstract contains a repeated word: 'near the the convective boundary'.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment of its significance. We address the single major comment below.

read point-by-point responses

  1. Referee: [simulation setup and comparison methodology] Simulation setup and comparison methodology section: the headline result that Urca reduces mixing efficiency near the boundary rests on the assumption that MAESTROeX's low-Mach boundary treatment and any implicit filtering do not themselves suppress mixing at the interface. Without explicit tests (e.g., altered boundary formulations, grid stretching variations, or resolution studies that quantify the mixing-efficiency metric), it remains possible that the reported reduction is at least partly numerical rather than physical.

    Authors: We agree that explicit verification of numerical robustness is important for claims about mixing near the convective boundary. Our headline result, however, is obtained from otherwise identical runs that differ only by the presence or absence of the Urca reactions; the low-Mach boundary conditions, grid, and any implicit filtering are therefore common to both sets of simulations. Any systematic suppression of mixing at the interface would appear in both the with-Urca and without-Urca cases, so the differential reduction we measure is attributable to the Urca process itself. The MAESTROeX boundary treatment used here is the same as that validated in our earlier studies of convective white-dwarf interiors. In the revised manuscript we will expand the methods section with a concise description of the boundary formulation and its prior validation for convective-boundary problems, and we will add a short paragraph discussing why the differential comparison isolates the physical effect from shared numerical influences. revision: yes

Circularity Check

0 steps flagged

No significant circularity: results from controlled numerical comparisons

full rationale

The paper isolates the convective Urca effects through direct side-by-side 3D MAESTROeX runs that differ only in the inclusion of the Urca pairs and associated weak reactions. The reported reduction in mixing efficiency near the boundary, the lack of restriction on convection-zone size, and the dominance of the A=23 pair are outputs of the time-dependent hydrodynamic evolution under the low-Mach approximation and chosen network; they are not obtained by fitting parameters to the target observables, by self-referential definitions, or by load-bearing self-citations that presuppose the result. The methodology therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger reflects standard modeling choices in low-Mach stellar hydrodynamics rather than paper-specific fitted constants or new entities.

axioms (1)
  • domain assumption The low Mach number approximation remains valid throughout the convection zone and near its boundary.
    Invoked by the choice of MAESTROeX for these simmering-white-dwarf runs.

pith-pipeline@v0.9.0 · 5822 in / 1440 out tokens · 44811 ms · 2026-05-20T04:22:10.149397+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

50 extracted references · 50 canonical work pages · 1 internal anchor

  1. [1]

    1978, ApJ, 226, 1034, doi: 10.1086/156681

    Alastuey, A., & Jancovici, B. 1978, ApJ, 226, 1034, doi: 10.1086/156681

    work page doi:10.1086/156681 1978

  2. [2]

    Barkat, Z., & Wheeler, J. C. 1990, ApJ, 355, 602, doi: 10.1086/168793

    work page doi:10.1086/168793 1990

  3. [3]

    2025a, ApJ, 979, 216, doi: 10.3847/1538-4357/ad9bb0

    Boyd, B., Calder, A., Townsley, D., & Zingale, M. 2025a, ApJ, 979, 216, doi: 10.3847/1538-4357/ad9bb0

    work page doi:10.3847/1538-4357/ad9bb0

  4. [4]

    C., Townsley, D

    Boyd, B., Calder, A. C., Townsley, D. M., & Zingale, M. 2025b, in Journal of Physics Conference Series, Vol. 2997, Journal of Physics Conference Series, 012006, doi: 10.1088/1742-6596/2997/1/012006

    work page doi:10.1088/1742-6596/2997/1/012006

  5. [5]

    Bruenn, S. W. 1973, ApJ, 183, L125, doi: 10.1086/181269

    work page doi:10.1086/181269 1973

  6. [6]

    A., Brown, E

    Chamulak, D. A., Brown, E. F., Timmes, F. X., & Dupczak, K. 2008, ApJ, 677, 160, doi: 10.1086/528944

    work page doi:10.1086/528944 2008

  7. [7]

    , year = 2010, volume =

    Cyburt, R. H., Amthor, A. M., Ferguson, R., et al. 2010, ApJS, 189, 240, doi: 10.1088/0067-0049/189/1/240

    work page doi:10.1088/0067-0049/189/1/240 2010

  8. [8]

    A., Truran, J

    Denissenkov, P. A., Truran, J. W., Herwig, F., et al. 2015, MNRAS, 447, 2696, doi: 10.1093/mnras/stu2589

    work page doi:10.1093/mnras/stu2589 2015

  9. [9]

    2019a, The Journal of Open Source Software, 4, 1757, doi: 10.21105/joss.01757

    Fan, D., Nonaka, A., Almgren, A., et al. 2019a, The Journal of Open Source Software, 4, 1757, doi: 10.21105/joss.01757

    work page doi:10.21105/joss.01757

  10. [10]

    S., Harpole, A., & Zingale, M

    Fan, D., Nonaka, A., Almgren, A. S., Harpole, A., & Zingale, M. 2019b, ApJ, 887, 212, doi: 10.3847/1538-4357/ab4f75

    work page doi:10.3847/1538-4357/ab4f75

  11. [11]

    J., Challis, P

    Foley, R. J., Challis, P. J., Chornock, R., et al. 2013, ApJ, 767, 57, doi: 10.1088/0004-637X/767/1/57

    work page doi:10.1088/0004-637x/767/1/57 2013

  12. [12]

    FLASH: An Adaptive Mesh Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes

    Fryxell, B., Olson, K., Ricker, P., et al. 2000, ApJS, 131, 273, doi: 10.1086/317361

    work page doi:10.1086/317361 2000

  13. [13]

    C., Dewitt, H

    Graboske, H. C., Dewitt, H. E., Grossman, A. S., & Cooper, M. S. 1973, ApJ, 181, 457, doi: 10.1086/152062

    work page doi:10.1086/152062 1973

  14. [14]

    J., Kawaler, S

    Hansen, C. J., Kawaler, S. D., & Trimble, V. 2004, Stellar interiors : physical principles, structure, and evolution (Springer New York, NY), doi: 10.1007/978-1-4419-9110-2

    work page doi:10.1007/978-1-4419-9110-2 2004

  15. [15]

    Harris and K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

    work page doi:10.1038/s41586-020-2649-2 2020

  16. [16]

    Hunter, J. D. 2007, Computing in Science and Engg., 9, 90, doi: 10.1109/MCSE.2007.55

    work page doi:10.1109/mcse.2007.55 2007

  17. [17]

    1996, ApJ, 470, 1015, doi: 10.1086/177926

    Itoh, N., Nishikawa, A., & Kohyama, Y. 1996, ApJ, 470, 1015, doi: 10.1086/177926

    work page doi:10.1086/177926 1996

  18. [18]

    Itoh, N., Totsuji, H., Ichimaru, S., & Dewitt, H. E. 1979, ApJ, 234, 1079, doi: 10.1086/157590

    work page doi:10.1086/157590 1979

  19. [19]

    2012, Journal of the Atmospheric Sciences, 69, 961, doi: 10.1175/JAS-D-11-0110.1

    Klein, R., & Pauluis, O. 2012, Journal of the Atmospheric Sciences, 69, 961, doi: 10.1175/JAS-D-11-0110.1

    work page doi:10.1175/jas-d-11-0110.1 2012

  20. [20]

    Lesaffre, P., Podsiadlowski, P., & Tout, C. A. 2005, MNRAS, 356, 131, doi: 10.1111/j.1365-2966.2004.08428.x

    work page doi:10.1111/j.1365-2966.2004.08428.x 2005

  21. [21]

    K., & Han, Z

    Liu, Z.-W., R¨ opke, F. K., & Han, Z. 2023, Research in Astronomy and Astrophysics, 23, 082001, doi: 10.1088/1674-4527/acd89e

    work page doi:10.1088/1674-4527/acd89e 2023

  22. [22]

    M., Nonaka, A., Almgren, A

    Malone, C. M., Nonaka, A., Almgren, A. S., Bell, J. B., & Zingale, M. 2011, ApJ, 728, 118, doi: 10.1088/0004-637X/728/2/118

    work page doi:10.1088/0004-637x/728/2/118 2011

  23. [23]

    2014, ARA&A, 52, 107, doi: 10.1146/annurev-astro-082812-141031 Mart´ ınez-Rodr´ ıguez, H., Piro, A

    Maoz, D., Mannucci, F., & Nelemans, G. 2014, ARA&A, 52, 107, doi: 10.1146/annurev-astro-082812-141031 Mart´ ınez-Rodr´ ıguez, H., Piro, A. L., Schwab, J., & Badenes, C. 2016, ApJ, 825, 57, doi: 10.3847/0004-637X/825/1/57

    work page doi:10.1146/annurev-astro-082812-141031 2014

  24. [24]

    Simpy: A discrete-event simulation library,

    Meurer, A., Smith, C. P., Paprocki, M., et al. 2017, PeerJ Computer Science, 3, e103, doi: 10.7717/peerj-cs.103

    work page doi:10.7717/peerj-cs.103 2017

  25. [25]

    K., & Yokoi, K

    Nomoto, K., Thielemann, F. K., & Yokoi, K. 1984, ApJ, 286, 644, doi: 10.1086/162639 Paczy´ nski, B. 1972, Astrophys. Lett., 11, 53

    work page doi:10.1086/162639 1984

  26. [26]

    1999, ApJ, 517, 565

    Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565

    work page 1999

  27. [27]

    Phillips, M. M. 1993, ApJ, 413, L105

    work page 1993

  28. [28]

    2017, ApJ, 836, L9, doi: 10.3847/2041-8213/aa5c7e

    Piersanti, L., Bravo, E., Cristallo, S., et al. 2017, ApJ, 836, L9, doi: 10.3847/2041-8213/aa5c7e

    work page doi:10.3847/2041-8213/aa5c7e 2017

  29. [29]

    2022, ApJ, 926, 103, doi: 10.3847/1538-4357/ac403b

    Piersanti, L., Bravo, E., Straniero, O., Cristallo, S., & Dom´ ınguez, I. 2022, ApJ, 926, 103, doi: 10.3847/1538-4357/ac403b

    work page doi:10.3847/1538-4357/ac403b 2022

  30. [30]

    On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf

    Piro, A. L., & Bildsten, L. 2008, ApJ, 673, 1009, doi: 10.1086/524189 Poca-Amor´ os, F., Boyd, B., Townsley, D. M., & Calder, A. 2026, On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf. https://arxiv.org/abs/2602.21440 17

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/524189 2008

  31. [31]

    G., Filippenko, A

    Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009

    work page 1998

  32. [32]

    L., & Badenes, C

    Schwab, J., Mart´ ınez-Rodr´ ıguez, H., Piro, A. L., & Badenes, C. 2017, ApJ, 851, 105, doi: 10.3847/1538-4357/aa9a3c

    work page doi:10.3847/1538-4357/aa9a3c 2017

  33. [33]

    I., Johnson, E

    Smith, A. I., Johnson, E. T., Chen, Z., et al. 2023, The Astrophysical Journal, 947, 65, doi: 10.3847/1538-4357/acbaff

    work page doi:10.3847/1538-4357/acbaff 2023

  34. [34]

    Stein, J., Barkat, Z., & Wheeler, J. C. 1999, ApJ, 523, 381, doi: 10.1086/307715

    work page doi:10.1086/307715 1999

  35. [35]

    Stein, J., & Wheeler, J. C. 2006, ApJ, 643, 1190, doi: 10.1086/503246

    work page doi:10.1086/503246 2006

  36. [36]

    , year = 2016, volume =

    Suzuki, T., Toki, H., & Nomoto, K. 2016, ApJ, 817, 163, doi: 10.3847/0004-637X/817/2/163

    work page doi:10.3847/0004-637x/817/2/163 2016

  37. [37]

    2023, GNU Parallel 20231122 (’Grindav´ ık’), Zenodo, doi: 10.5281/zenodo.10199085

    Tange, O. 2023, GNU Parallel 20231122 (’Grindav´ ık’), Zenodo, doi: 10.5281/zenodo.10199085

    work page doi:10.5281/zenodo.10199085 2023

  38. [38]

    2017, in Handbook of Supernovae, ed

    Taubenberger, S. 2017, in Handbook of Supernovae, ed. A. W. Alsabti & P. Murdin (Springer International Publishing), 317, doi: 10.1007/978-3-319-21846-5 37

    work page doi:10.1007/978-3-319-21846-5 2017

  39. [39]

    Team, A.-A. M. D., Bishop, A., Fields, C. E., et al. 2024, AMReX-Astro/Microphysics: Release 24.08, 24.08, Zenodo, doi: 10.5281/zenodo.13150461

    work page doi:10.5281/zenodo.13150461 2024

  40. [40]

    X., Brown, E

    Timmes, F. X., Brown, E. F., & Truran, J. W. 2003, ApJ, 590, L83, doi: 10.1086/376721

    work page doi:10.1086/376721 2003

  41. [41]

    X., & Swesty, F

    Timmes, F. X., & Swesty, F. D. 2000, ApJS, 126, 501, doi: 10.1086/313304

    work page doi:10.1086/313304 2000

  42. [42]

    J., Smith, B

    Turk, M. J., Smith, B. D., Oishi, J. S., et al. 2011, The Astrophysical Journal Supplement Series, 192, 9, doi: 10.1088/0067-0049/192/1/9

    work page doi:10.1088/0067-0049/192/1/9 2011

  43. [43]

    1999, ApJ, 522, L43, doi: 10.1086/312213

    Umeda, H., Nomoto, K., Kobayashi, C., Hachisu, I., & Kato, M. 1999, ApJ, 522, L43, doi: 10.1086/312213

    work page doi:10.1086/312213 1999

  44. [44]

    M., Lecoanet, D., Brown, B

    Vasil, G. M., Lecoanet, D., Brown, B. P., Wood, T. S., & Zweibel, E. G. 2013, ApJ, 773, 169, doi: 10.1088/0004-637X/773/2/169 Wes McKinney. 2010, in Proceedings of the 9th Python in Science Conference, ed. St´ efan van der Walt & Jarrod Millman, 56 – 61, doi: 10.25080/Majora-92bf1922-00a

    work page doi:10.1088/0004-637x/773/2/169 2013

  45. [45]

    E., & Zingale, M

    Willcox, D. E., & Zingale, M. 2018, Journal of Open Source Software, 3, 588, doi: 10.21105/joss.00588

    work page doi:10.21105/joss.00588 2018

  46. [46]

    E., Taam, R

    Woosley, S. E., Taam, R. E., & Weaver, T. A. 1986, ApJ, 301, 601, doi: 10.1086/163926

    work page doi:10.1086/163926 1986

  47. [47]

    Zegers, R. G. T., Brown, E. F., Akimune, H., et al. 2008, Phys. Rev. C, 77, 024307, doi: 10.1103/PhysRevC.77.024307

    work page doi:10.1103/physrevc.77.024307 2008

  48. [48]

    2019, Journal of Open Source Software, 4, 1370, doi: 10.21105/joss.01370

    Zhang, W., Almgren, A., Beckner, V., et al. 2019, Journal of Open Source Software, 4, 1370, doi: 10.21105/joss.01370

    work page doi:10.21105/joss.01370 2019

  49. [49]

    S., Bell, J

    Zingale, M., Almgren, A. S., Bell, J. B., Nonaka, A., & Woosley, S. E. 2009, ApJ, 704, 196, doi: 10.1088/0004-637X/704/1/196

    work page doi:10.1088/0004-637x/704/1/196 2009

  50. [50]

    2024, ApJ, 977, 30, doi: 10.3847/1538-4357/ad8a66 This paper was built using the Open Journal of As- trophysics LATEX template

    Clark, A. 2024, ApJ, 977, 30, doi: 10.3847/1538-4357/ad8a66 This paper was built using the Open Journal of As- trophysics LATEX template. The OJA is a journal which provides fast and easy peer review for new papers in the astro-phsection of the arXiv, making the reviewing pro- cess simpler for authors and referees alike. Learn more at http://astro.theoj.org

    work page doi:10.3847/1538-4357/ad8a66 2024