pith. sign in

arxiv: 2302.12089 · v5 · submitted 2023-02-23 · 🌌 astro-ph.HE · astro-ph.IM

Time integration for neutrino radiation transport using minimally implicit Runge-Kutta methods

Samuel Santos-P\'erez , Martin Obergaulinger , Isabel Cordero-Carri\'on This is my paper

Pith reviewed 2026-05-24 09:23 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IM
keywords neutrino radiation transportminimally implicit Runge-Kuttatime integrationstiff equationsradiation hydrodynamicscore-collapse supernovaeneutrino-matter interactions
0
0 comments X

The pith

Minimally implicit Runge-Kutta methods integrate neutrino radiation transport equations stably at explicit computational cost.

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

This paper develops minimally implicit Runge-Kutta methods for the time integration of equations describing the interaction between neutrinos and matter in astrophysical systems. The methods are constructed so the implicit operator can be inverted analytically rather than numerically, keeping the cost comparable to explicit schemes. Derivation relies on the physical behavior of variables when reactions become stiff, such as in optically thick environments. This approach is tested on neutrino-matter reactions and applied to core-collapse supernova simulations. A reader would care because many astrophysical evolutions involve stiff radiation hydrodynamics that are computationally expensive to simulate accurately.

Core claim

The authors present minimally implicit Runge-Kutta methods for integrating the coupled hydrodynamic and radiative transfer equations. By accounting for the physical behavior of the evolved variables in the stiff regime, these methods allow the implicit operator to be inverted analytically. As a result, the computational cost equals that of an explicit method while ensuring stability and accuracy in the treatment of neutrino reactions with matter.

What carries the argument

Minimally implicit Runge-Kutta methods derived from stiff-regime physical behavior that enable analytical inversion of the implicit operator.

If this is right

  • The methods maintain stability and accuracy in stiff regimes without additional approximations.
  • They can be applied directly to realistic core-collapse supernovae simulations.
  • The approach reduces the computational burden of radiation hydrodynamics calculations to explicit-method levels.
  • Tests on neutrino-matter reactions validate the method's performance.

Where Pith is reading between the lines

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

  • Similar minimally implicit techniques could be adapted for other radiation transport problems involving photons or other particles.
  • The analytical inversion might allow higher-resolution simulations of supernova dynamics that were previously limited by computational cost.
  • Extensions to multi-dimensional or relativistic cases could follow from the same physical-behavior principle.

Load-bearing premise

The physical behavior of the evolved variables in the stiff regime permits derivation of methods where the implicit operator inversion is analytical and preserves stability and accuracy without additional approximations or post-hoc adjustments.

What would settle it

Running a stiff neutrino transport test case with the new methods and comparing the results to those from a standard implicit solver; disagreement in stability or accuracy beyond expected truncation error would falsify the claims.

Figures

Figures reproduced from arXiv: 2302.12089 by Isabel Cordero-Carri\'on, Martin Obergaulinger, Samuel Santos-P\'erez.

Figure 1
Figure 1. Figure 1: Results for E using the MIRK1 method with κ = 102 and CFL=1/2 5 . Several spatial resolutions at t = 5 are displayed. The whole spatial domain is shown at the left panel and a zoom near the origin is displayed at the right panel. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 r 0 0.2 0.4 0.6 0.8 1 1.2 F n=50 n=100 n=200 n=400 exact solution 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 r 0.8 0.85 0.9 0.95 1… view at source ↗
Figure 2
Figure 2. Figure 2: Results for F using the MIRK1 method with κ = 102 and CFL=1/2 5 . Several spatial resolutions at t = 5 are displayed. The whole spatial domain is shown at the left panel and a zoom near the maximum value for F is displayed at the right panel. Let us compute estimates of the order of convergence. We denote ε(∆t, ∆r) the L2 norm of the error of the numerical solution for E at a temporal res￾olution ∆t and a … view at source ↗
Figure 3
Figure 3. Figure 3: L2 norm errors of E between the exact and numerical solutions divided by E(t, r = 0) as function of ∆t, using the MIRK1 method and κ = 102 . CFL=1/2 k with k = 0, 1, . . . , 6 (values from right to left for each curve) have been used. We find two regimes with a specific dependence of the error on ∆t: for a larger time step, ε(∆t, ∆r) is described by a power law, ε(∆t, ∆r) ∝ (∆t) pt , while the error become… view at source ↗
Figure 4
Figure 4. Figure 4: Results for E using the MIRK2 method with κ = 102 and CFL=1/2 5 . Several spatial resolutions at t = 5 are displayed. The whole spatial domain is shown at the left panel and a zoom near the origin is displayed at the right panel. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 r 0 0.2 0.4 0.6 0.8 1 1.2 F n=50 n=100 n=200 n=400 exact solution 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 r 0.8 0.85 0.9 0.95 1… view at source ↗
Figure 5
Figure 5. Figure 5: Results for F using the MIRK2 method with κ = 102 and CFL=1/2 5 . Several spatial resolutions at t = 5 are displayed. The whole spatial domain is shown at the left panel and a zoom near the maximum value for F is displayed at the right panel. The errors ε(∆t, ∆r) for the same spatial and temporal resolutions as for the MIRK1 method are shown in [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: L2 norm errors of E between the exact and numerical solutions divided by E(t, r = 0) as function of ∆t, using the MIRK2 method and κ = 102 . CFL=1/2 k with k = 0, 1, . . . , 6 (values from right to left for each curve) have been used. for the spatial error at fixed CFL factor (see the second row of Tab. 1) is observed, before entering the regime in which the temporal error dominates over the spatial one. A… view at source ↗
Figure 7
Figure 7. Figure 7: Results for E using the MIRK1 method with κ = 105 and CFL=1/2 5 . Several spatial resolutions at t = 10 are displayed. t ( t, r) 2 -20.0 2 -19.0 2 -18.0 2 -17.0 2 -16.0 2 -15.0 2 -14.0 2 -13.0 2 -12.0 2 -11.0 2 -10.0 2-9.0 2-8.0 2-7.0 2-6.0 2-5.0 2-4.0 2-3.0 n=3200 n=6400 n=12800 n=25600 [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: L2 norm errors of E between the exact and numerical solutions divided by E(t, r = 0) as function of ∆t, using the MIRK1 method and κ = 105 . CFL=1/2 k with k = 0, 1, . . . , 5 (values from right to left for each curve) have been used. The case k = 6 for n = 3200, 6400 has been also included. 20 [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Results for E using the MIRK2 method with κ = 105 and CFL=1/2 5 . Several spatial resolutions at t = 10 are displayed. An important difference between the results obtained for the MIRK1 and the MIRK2 methods is that in the case of the MIRK2 method all the considered resolutions lie within the convergence regime, while in the case of the MIRK1 method this only happens when higher resolutions are considered.… view at source ↗
Figure 10
Figure 10. Figure 10: L2 norm errors of E between the exact and numerical solutions divided by E(t, r = 0) as function of ∆t, using the MIRK2 method and κ = 105 . CFL=1/2 k with k = 0, 1, . . . , 5 (values from right to left for each curve) have been used. The case k = 6 for n = 3200, 6400 has been also included. the order of convergence (for both MIRK1 and MIRK2), but resulting in similar values. 4.2. Toy model for a proto-ne… view at source ↗
Figure 11
Figure 11. Figure 11: Initial data of the toy model for PNS cooling. Top panel: d [PITH_FULL_IMAGE:figures/full_fig_p025_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Radial profiles of the RK2 model of the PNS toy model at fo [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Time evolution of the neutrino luminosity of the PNS toy mod [PITH_FULL_IMAGE:figures/full_fig_p028_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Same as Fig. 12. The lines show the results of simulation MIR [PITH_FULL_IMAGE:figures/full_fig_p029_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Same as Fig. 13, but comparing run MIRK2-1 run (lines) to [PITH_FULL_IMAGE:figures/full_fig_p030_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Same as Fig. 12. The lines show the results of simulation MIR [PITH_FULL_IMAGE:figures/full_fig_p031_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Evolution of important variables of our simulations. Models [PITH_FULL_IMAGE:figures/full_fig_p033_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Evolution of important variables of our simulations. Models [PITH_FULL_IMAGE:figures/full_fig_p034_18.png] view at source ↗
Figure 20
Figure 20. Figure 20: The density profiles at representative epochs during the e [PITH_FULL_IMAGE:figures/full_fig_p035_20.png] view at source ↗
Figure 19
Figure 19. Figure 19: Radial profiles of selected models at a few times after bou [PITH_FULL_IMAGE:figures/full_fig_p036_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Radial profiles of selected models at a few times after bou [PITH_FULL_IMAGE:figures/full_fig_p037_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Comparison of unstable simulations, as indicated in the lege [PITH_FULL_IMAGE:figures/full_fig_p039_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: Same as Fig. 17, but comparison of simulations that evolve [PITH_FULL_IMAGE:figures/full_fig_p040_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: Same as Fig. 18, but comparison of simulations that evolve [PITH_FULL_IMAGE:figures/full_fig_p041_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Same as Fig. 19, but comparison of simulations that evolve [PITH_FULL_IMAGE:figures/full_fig_p042_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: Same as Fig. 20, but comparison of simulations that evolve [PITH_FULL_IMAGE:figures/full_fig_p043_25.png] view at source ↗
read the original abstract

The evolution of many astrophysical systems is dominated by the interaction between matter and radiation such as photons or neutrinos. The dynamics can be described by the evolution equations of radiation hydrodynamics in which reactions between matter particles and radiation quanta couples the hydrodynamic equations to those of radiative transfer (see Munier & Weaver (1986a) and Munier & Weaver (1986b)). The numerical treatment has to account for their potential stiffness (e.g., in optically thick environments). In this article, we will present a new method to numerically integrate these equations in a stable way by using minimally implicit Runge-Kutta methods. With these methods, the inversion of the implicit operator can be done analytically, so the computational cost is equivalent to that of an explicit method. We strongly take into account the physical behavior of the evolved variables in the limit of the stiff regime in the derivation of the methods. We will show the results of applying these methods to the reactions between neutrinos and matter in some tests and also in realistic core-collapse supernovae simulations.

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

2 major / 2 minor

Summary. The manuscript introduces minimally implicit Runge-Kutta methods for time integration of the coupled radiation-hydrodynamics equations governing neutrino-matter interactions. By incorporating the physical stiff-limit behavior of the evolved variables, the methods are constructed so that the implicit operator inverts analytically, yielding computational cost comparable to explicit schemes while aiming to preserve stability and accuracy. The approach is demonstrated on neutrino reaction tests and applied to realistic core-collapse supernova simulations.

Significance. If the analytical inversion holds without hidden approximations that degrade accuracy or stability, the method would provide an efficient alternative for handling stiff source terms in radiation transport, potentially enabling higher-resolution astrophysical simulations at reduced cost. The grounding in physical stiff-regime limits is a constructive feature for relaxation-type systems, though the absence of detailed error analysis limits assessment of its robustness.

major comments (2)
  1. [Abstract and derivation section] Abstract and derivation section: the central claim that the implicit operator inverts analytically due to stiff-regime physical behavior is asserted without an explicit derivation, reduced system equations, or stability proof; this is load-bearing for the efficiency and accuracy assertions.
  2. [Results section (tests and CCSN simulations)] Results section (tests and CCSN simulations): no quantitative error analysis, convergence rates, or direct comparisons to standard IMEX or fully implicit methods are provided to confirm that stability and accuracy are retained without post-hoc adjustments.
minor comments (2)
  1. [Abstract] The abstract would benefit from a brief statement of the specific test problems used and the number of neutrino species/groups considered.
  2. [Introduction] Notation for the radiation moments and coupling terms should be defined consistently with standard references such as Munier & Weaver (1986a,b) to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful review and constructive comments on our manuscript. We address each major comment below and outline revisions that will strengthen the presentation while preserving the core contributions of the minimally implicit Runge-Kutta approach.

read point-by-point responses

  1. Referee: [Abstract and derivation section] Abstract and derivation section: the central claim that the implicit operator inverts analytically due to stiff-regime physical behavior is asserted without an explicit derivation, reduced system equations, or stability proof; this is load-bearing for the efficiency and accuracy assertions.

    Authors: We agree that the manuscript would benefit from a more explicit step-by-step derivation of the analytical inversion. In the revised version we will add the reduced system equations for the stiff limit, detail how the physical behavior of the evolved variables (e.g., equilibrium values of neutrino number densities and energies) leads to an analytically invertible implicit operator, and include a short stability argument showing that the scheme remains stable by construction in the stiff regime. These additions will be placed in a new subsection of the methods section without altering the reported results. revision: yes

  2. Referee: [Results section (tests and CCSN simulations)] Results section (tests and CCSN simulations): no quantitative error analysis, convergence rates, or direct comparisons to standard IMEX or fully implicit methods are provided to confirm that stability and accuracy are retained without post-hoc adjustments.

    Authors: The current tests demonstrate qualitative stability and physical consistency on both idealized reaction problems and full CCSN runs, but we acknowledge the absence of quantitative error metrics and direct method comparisons. In the revision we will add L1 and L2 error norms versus reference solutions for the test suite, report observed convergence rates under grid refinement, and include timing and accuracy comparisons against a standard IMEX Runge-Kutta scheme on a subset of the CCSN models. These additions will be confined to an expanded results section and will not require new simulations. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives minimally implicit Runge-Kutta methods by incorporating the physical stiff-limit behavior of neutrino-matter coupling into the method construction, allowing analytic inversion of the implicit operator. No steps reduce to self-definition, fitted inputs renamed as predictions, or load-bearing self-citations. The central claim rests on standard IMEX construction for relaxation systems applied to the radiation hydrodynamics equations, with the stiff-regime analysis providing independent physical input rather than tautological closure. The abstract and described approach show no reduction of the result to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the domain assumption that stiff-regime physics allows analytical inversion; no free parameters or invented entities are mentioned in the abstract.

axioms (2)
  • standard math Standard Runge-Kutta stability and consistency properties hold for the derived schemes
    The methods are extensions of the RK framework applied to stiff ODEs from radiation hydrodynamics.
  • domain assumption Physical behavior of variables in the stiff limit permits analytical inversion of the implicit operator
    Explicitly stated in the abstract as the basis for method derivation.

pith-pipeline@v0.9.0 · 5722 in / 1229 out tokens · 25228 ms · 2026-05-24T09:23:49.607803+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    With these methods, the inversion of the implicit operator can be done analytically, so the computational cost is equivalent to that of an explicit method. We strongly take into account the physical behavior of the evolved variables in the limit of the stiff regime in the derivation of the methods.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    Conditions a < 1/2, b < 1/2 must be fulfilled for the spectral radius of the updated matrix to be strictly bounded by 1... With this choice, it is satisfied that E^{n+1} = E_eq^n = E_eq^{n+1} + O(Δt).

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

35 extracted references · 35 canonical work pages · 13 internal anchors

  1. [1]

    Munier, R

    A. Munier, R. Weaver, Radiation transfer in the fluid frame: A cov ariant formulation: Part i: Radiation hydrodynamics, Comput. Phys. Rep. 3 (3) (1986) 127–164

    work page 1986

  2. [2]

    Munier, R

    A. Munier, R. Weaver, Radiation transfer in the fluid frame: A cov ariant formulation: Part ii: The radiation transfer equation, Comput. Phy s. Rep. 3 (3) (1986) 165–208

    work page 1986

  3. [3]

    Cordero-Carri´ on, S

    I. Cordero-Carri´ on, S. Santos-P´ erez, C. Mart ´ ınez-Vidallach, Numerical evolution of the resistive relativistic magnetohydrodynamic equations: A minimally implicit runge-kutta scheme, Appl Math Comput 443 (2023) 127774. doi:10.1016/j.amc.2022.127774

    work page doi:10.1016/j.amc.2022.127774 2023

  4. [4]

    J. F. Mahlmann, M. A. Aloy, V. Mewes, P. Cerda-Duran, Computa - tional general relativistic force-free electrodynamics-ii. charac terization of numerical diffusivity, Astronomy & Astrophysics 647 (2021) A58

    work page 2021

  5. [5]

    Koellermeier, G

    J. Koellermeier, G. Samaey, Spatially adaptive projective integra tion schemes for stiff hyperbolic balance laws with spectral gaps, J. Com put. Math. 8 (2022) 295–325. 48

    work page 2022

  6. [6]

    Koellermeier, E

    J. Koellermeier, E. Pimentel-Garc ´ ıa, Steady states and well-bala nced schemes for shallow water moment equations with topography, App l Math Comput 427 (2022) 127166

    work page 2022

  7. [7]

    Pareschi, G

    L. Pareschi, G. Russo, Implicit–explicit runge–kutta schemes an d appli- cations to hyperbolic systems with relaxation, J Sci Comp 25 (1) (20 05) 129–155

    work page

  8. [8]

    M. R. Izquierdo, L. Pareschi, B. Mi˜ nano, J. Mass´ o, C. Palenzuela, Global high-order numerical schemes for the time evolution of the genera l rel- ativistic radiation magneto-hydrodynamics equations, ArXiv Gener al Relativity and Quantum Cosmology e-prints, arXiv:gr-qc/00027 (20 22). arXiv:gr-qc/00027

    work page

  9. [9]

    Foucart, E

    F. Foucart, E. O’Connor, L. Roberts, L. E. Kidder, H. P. Pfeiffe r, M. A. Scheel, Impact of an improved neutrino energy estimate on outflow s in neutron star merger simulations, Physical Review D 94 (2016) 1230 16. doi:10.1103/PhysRevD.94.123016

    work page doi:10.1103/physrevd.94.123016 2016

  10. [10]

    Radice, S

    D. Radice, S. Bernuzzi, A. Perego, R. Haas, A new moment-bas ed general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers, Monthly Notices of the RAS 5 12 (2022) 1499–1521. doi:10.1093/mnras/stac589

    work page doi:10.1093/mnras/stac589 2022

  11. [11]

    O. Just, M. Obergaulinger, H.-T. Janka, A new multidimen- sional, energy-dependent two-moment transport code for neut rino- hydrodynamics, Monthly Notices of the RAS 453 (4) (2015) 3386–3 413

    work page 2015

  12. [12]

    Foucart, Monte carlo closure for moment-based transpor t schemes in general relativistic radiation hydrodynamic simulations, Monthly No- tices of the RAS 475 (2018) 4186–4207

    F. Foucart, Monte carlo closure for moment-based transpor t schemes in general relativistic radiation hydrodynamic simulations, Monthly No- tices of the RAS 475 (2018) 4186–4207. doi:10.1093/mnras/sty108

    work page doi:10.1093/mnras/sty108 2018

  13. [13]

    Foucart, M

    F. Foucart, M. D. Duez, F. H´ ebert, L. E. Kidder, P. Kovarik, H. P. Pfeiffer, M. A. Scheel, Implementation of monte carlo transport in t he general relativistic spec code, Astrophysical Journal 920 (2) (2 021) 82. doi:10.3847/1538-4357/ac1737

    work page doi:10.3847/1538-4357/ac1737

  14. [14]

    S. Jin, C. Levermore, Numerical Schemes for Hyperbolic Conse rvation Laws with Stiff Realaxation Terms, JCoPh 126 (1996) 446–467. 49

    work page 1996

  15. [15]

    J. A. Pons, J. M. Ib´ a˜ nez, J. A. Miralles, Hyperbolic character of the angular moment equations of radia- tive transfer and numerical methods, Monthly Notices of the RAS 317 (2000) 550–562. arXiv:astro-ph/0005310, doi:10.1046/j.1365-8711.2000.03679.x

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1046/j.1365-8711.2000.03679.x 2000

  16. [16]

    S. Jin, L. Pareschi, G. Toscani, Uniformly accurate diffusive rela xation schemes for multiscale transport equations, SIAM J. Num. Anal. 38 (2000) 913

    work page 2000

  17. [17]

    A radiation-hydrodynamics scheme valid from the transport to the diffusion limit

    E. Audit, P. Charrier, J. Chi` eze, B. Dubroca, A radiation- hydrodynamics scheme valid from the transport to the diffusion limit, ArXiv Astrophysics e-prints, arXiv:astro-ph/0206281 (#jun# 20 02). arXiv:astro-ph/0206281

    work page internal anchor Pith review Pith/arXiv arXiv

  18. [18]

    R. J. LeVeque, Numerical methods for conservation laws, Vol. 214, Springer, 1992

    work page 1992

  19. [19]

    On the measurements of numerical viscosity and resistivity in Eulerian MHD codes

    T. Rembiasz, M. Obergaulinger, P. Cerd´ a-Dur´ an, M.-´A. Aloy, E. M¨ uller, On the Measurements of Numerical Viscosity and Resistivity in Eule- rian MHD Codes, Astrophysical Journal, Supplement 230 (2) (2017 ) 18. arXiv:1611.05858, doi:10.3847/1538-4365/aa6254

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4365/aa6254 2017

  20. [20]

    J. Verwer, Convergence and order reduction of diagonally implic it Runge-Kutta schemes in the method of lines, Numerical Analysis: Pr o- ceedings of the Dundee Conference on Numerical Analysis (1975) 2 20– 237doi:https://ir.cwi.nl/pub/2470

    work page 1975

  21. [21]

    H. T. Janka, Conditions for shock revival by neutrino heating in core- collapse supernovae, Astronomy and Astrophysics 368 (2001) 52 7–560. arXiv:astro-ph/0008432, doi:10.1051/0004-6361:20010012

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361:20010012 2001

  22. [22]

    Rampp, H.-T

    M. Rampp, H.-T. Janka, Radiation hydrodynamics with neutri- nos. Variable Eddington factor method for core-collapse super- nova simulations, Astronomy and Astrophysics 396 (2002) 361–39 2. doi:10.1051/0004-6361:20021398

    work page doi:10.1051/0004-6361:20021398 2002

  23. [23]

    A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time

    M. Liebend¨ orfer, O. E. B. Messer, A. Mezzacappa, S. W. Bru enn, C. Y. Cardall, F. Thielemann, A Finite Difference Representation of Neutrin o 50 Radiation Hydrodynamics in Spherically Symmetric General Relativis- tic Spacetime, Astrophysical Journal, Supplement 150 (2004) 263 –316. arXiv:arXiv:astro-ph/0207036, doi:10.1086/380191

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/380191 2004

  24. [24]

    Stellar core collapse in full general relativity with microphysics - Formulation and Spherical collapse test -

    Y. Sekiguchi, Stellar Core Collapse in Full General Relativity with Mi- crophysics – Formulation and Spherical Collapse Test, Prog. Theor . Phys. 124 (2010) 331–379. arXiv:1009.3320

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  25. [25]

    A New Open-Source Code for Spherically-Symmetric Stellar Collapse to Neutron Stars and Black Holes

    E. O’Connor, C. D. Ott, A new open-source code for spherically symmetric stellar collapse to neutron stars and black holes, Class Quantum Gravity 27 (11) (2010) 114103–+. arXiv:0912.2393, doi:10.1088/0264-9381/27/11/114103

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/27/11/114103 2010

  26. [26]

    A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry

    B. M¨ uller, H. Janka, H. Dimmelmeier, A New Multi-dimensional Gen- eral Relativistic Neutrino Hydrodynamic Code for Core-collapse Su- pernovae. I. Method and Code Tests in Spherical Symmetry, Astr o- physical Journal, Supplement 189 (2010) 104–133. arXiv:1001.4841, doi:10.1088/0067-0049/189/1/104

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0067-0049/189/1/104 2010

  27. [27]

    O’Connor, An Open-source Neutrino Radiation Hydrodynamic s Code for Core-collapse Supernovae, Astrophysical Journal, Suppleme nt 219 (2015) 24

    E. O’Connor, An Open-source Neutrino Radiation Hydrodynamic s Code for Core-collapse Supernovae, Astrophysical Journal, Suppleme nt 219 (2015) 24. doi:10.1088/0067-0049/219/2/24

    work page doi:10.1088/0067-0049/219/2/24 2015

  28. [28]

    A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars

    T. Kuroda, T. Takiwaki, K. Kotake, A New Multi-energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to the Gravitational Collapse of Massive Stars, Astro- physical Journal, Supplement 222 (2016) 20. arXiv:1501.06330, doi:10.3847/0067-0049/222/2/20

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/0067-0049/222/2/20 2016

  29. [29]

    An advanced leakage scheme for neutrino treatment in astrophysical simulations

    A. Perego, R. M. Cabez´ on, R. K¨ appeli, An Advanced Leakage Scheme for Neutrino Treatment in Astrophysical Simulations, Astr o- physical Journal, Supplement 223 (2016) 22. arXiv:1511.08519, doi:10.3847/0067-0049/223/2/22

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/0067-0049/223/2/22 2016

  30. [30]

    O’Connor, R

    E. O’Connor, R. Bollig, A. Burrows, S. Couch, T. Fischer, H.-T. Janka, K. Kotake, E. J. Lentz, M. Liebend¨ orfer, O. E. B. Messe r, A. Mezzacappa, T. Takiwaki, D. Vartanyan, Global comparison of core-collapse supernova simulations in spherical symmetry, J Phys G Nucl Part Phys 45 (2018) 104001, aDS Bibcode: 2018JPhG...45j400 1O. doi:10.1088/1361-6471/aadeae. 51

    work page doi:10.1088/1361-6471/aadeae 2018

  31. [31]

    O. Just, R. Bollig, H.-T. Janka, M. Obergaulinger, R. Glas, S. Nag ataki, Core-collapse supernova simulations in one and two dimensions: com- parison of codes and approximations, Monthly Notices of the RAS 48 1 (2018) 4786–4814. doi:10.1093/mnras/sty2578

    work page doi:10.1093/mnras/sty2578 2018

  32. [32]

    M. P. Laiu, E. Endeve, R. Chu, J. A. Harris, O. E. B. Messer, A DG-IMEX Method for Two-moment Neutrino Transport: Non- linear Solvers for Neutrino-Matter Coupling, Astrophysical Jour- nal, Supplement 253 (2021) 52, aDS Bibcode: 2021ApJS..253...52L. doi:10.3847/1538-4365/abe2a8

    work page doi:10.3847/1538-4365/abe2a8 2021

  33. [33]

    Exploring the relativistic regime with Newtonian hydrodynamics: An improved effective gravitational potential for supernova simulations

    A. Marek, H. Dimmelmeier, H.-T. Janka, E. M¨ uller, R. Buras, Exploring the relativistic regime with Newtonian hydrodynamics: an improved ef - fective gravitational potential for supernova simulations, Astro nomy and Astrophysics 445 (2006) 273–289. arXiv:arXiv:astro-ph/0502161, doi:10.1051/0004-6361:20052840

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361:20052840 2006

  34. [34]

    A. W. Steiner, M. Hempel, T. Fischer, Core-collapse Super- nova Equations of State Based on Neutron Star Observa- tions, Astrophysical Journal 774 (2013) 17. arXiv:1207.2184, doi:10.1088/0004-637X/774/1/17

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/774/1/17 2013

  35. [35]

    Time integration for neutrino radiation transport using minimally implicit Runge-Kutta methods

    M. J. Castro, C. Par´ es, Well-balanced high-order finite volume methods for systems of balance laws, J Sci Comp 82 (2) (2020) 48. 52 This figure "pnstoy-RK2-profiles.png" is available in "png" format from: http://arxiv.org/ps/2302.12089v3

    work page internal anchor Pith review Pith/arXiv arXiv 2020