[Paper Review] Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations

Markus Rampp, Hans‐Thomas Janka|ArXiv.org|Mar 7, 2002

Neutrino Physics ResearchPhysics and Astronomy86 references213 citations

TL;DR

This paper presents a novel numerical code for simulating neutrino transport in core-collapse supernovae using a variable Eddington factor method to solve the time- and energy-dependent Boltzmann equation. By iteratively solving moment equations and a model Boltzmann equation on a ray grid with time-implicit integration, the method enables accurate modeling of neutrino energy and lepton number exchange, with modular design supporting multi-dimensional and relativistic applications.

ABSTRACT

Neutrino transport and neutrino interactions in dense matter play a crucial role in stellar core collapse, supernova explosions and neutron star formation. Here we present a detailed description of a new numerical code for treating the time and energy dependent neutrino transport in hydrodynamical simulations of such events. The code is based on a variable Eddington factor method to deal with the integro-differential character of the Boltzmann equation. The moments of the neutrino distribution function and the energy and lepton number exchange with the stellar medium are determined by iteratively solving the zeroth and first order moment equations in combination with a model Boltzmann equation. The latter is discretized on a grid of tangent rays. The integration of the transport equations and the neutrino source terms is performed in a time-implicit way. In the present version of the program, the transport part is coupled to an explicit hydrodynamics code which follows the evolution of the stellar plasma by a finite-volume method with piecewise parabolic interpolation, using a Riemann solver for calculating the hydrodynamic states. The neutrino source terms are implemented in an operator-split step. Neutrino transport and hydrodynamics can be calculated with different spatial grids and different time steppings. The structure of the described code is modular and offers a high degree of flexibility for an application to relativistic and multi-dimensional problems at different levels of refinement and accuracy. We critically evaluate results for a number of test cases, including neutrino transport in rapidly moving stellar media and approximate relativistic core collapse, and suggest a path for generalizing the code to be used in multi-dimensional simulations of convection in neutron stars and supernovae.

Motivation & Objective

• To develop a robust, flexible numerical framework for simulating neutrino transport in core-collapse supernovae with high accuracy.
• To address the computational challenge of solving the integro-differential Boltzmann equation in multi-dimensional, time-dependent hydrodynamical environments.
• To enable accurate modeling of neutrino energy and lepton number exchange with dense stellar matter during core collapse and shock revival.
• To provide a modular, extensible code structure suitable for future application to relativistic and multi-dimensional supernova simulations.
• To validate the method through test cases involving rapid flows and approximate relativistic collapse, demonstrating its reliability for complex astrophysical scenarios.

Proposed method

• Uses a variable Eddington factor method to approximate the angular dependence of the neutrino distribution function, reducing the full Boltzmann equation to a set of moment equations.
• Solves the zeroth- and first-order moment equations iteratively, coupled with a model Boltzmann equation discretized on a grid of tangent rays.
• Employs time-implicit integration for the transport and source term equations to ensure stability in stiff, rapidly evolving systems.
• Couples the neutrino transport to an explicit finite-volume hydrodynamics code using piecewise parabolic reconstruction and a Riemann solver for hydrodynamic states.
• Implements neutrino source terms via operator splitting, allowing independent spatial grids and time steps for hydrodynamics and neutrino transport.
• Models nuclear composition transitions via three distinct density-temperature regimes, with continuous transitions governed by dissociation/recombination factors $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$.

Experimental results

Research questions

• RQ1 How can the integro-differential Boltzmann equation for neutrino transport be efficiently and accurately solved in time- and energy-dependent core-collapse supernova simulations?
• RQ2 Can a variable Eddington factor method with iterative moment solving and ray-based discretization achieve sufficient accuracy for modeling neutrino heating and energy deposition?
• RQ3 How does the modular design of the code support future extensions to multi-dimensional and relativistic simulations?
• RQ4 What is the impact of nuclear composition transitions (e.g., $^{56}$Ni formation, $\alpha$-particle recombination, heavy nucleus dissociation) on energy and lepton transport during collapse?
• RQ5 How do the numerical results compare with known test cases, such as neutrino transport in rapidly moving media or approximate relativistic collapse?

Key findings

• The variable Eddington factor method successfully captures the essential physics of neutrino transport in core-collapse supernovae, including energy and lepton number exchange with the stellar medium.
• The time-implicit integration scheme ensures numerical stability even in highly dynamic, stiff environments such as those encountered during core bounce and shock formation.
• The code's modular architecture allows independent spatial and temporal discretization for hydrodynamics and neutrino transport, enhancing computational flexibility and efficiency.
• Test simulations of neutrino transport in rapidly moving media and approximate relativistic collapse show good agreement with expected physical behavior, validating the method’s accuracy.
• The treatment of nuclear composition transitions via $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$ factors ensures consistent energy conservation and thermodynamic consistency during phase changes in the plasma.
• The method is well-suited for future application to multi-dimensional simulations of convection and neutrino-driven explosion mechanisms in neutron stars and supernovae.