The Impact of Different Neutrino Transport Methods on Multidimensional Core-collapse Supernova Simulations

Kuo-Chuan Pan,Carlos Mattes,Evan P. O'Connor,Sean M. Couch,Albino Perego,Almudena Arcones
DOI: https://doi.org/10.1088/1361-6471/aaed51
2018-10-15
Abstract:Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations.
High Energy Astrophysical Phenomena
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is the influence of different neutrino transport methods on multi - dimensional core - collapse supernova simulations. Specifically, by comparing three approximate neutrino transport schemes (the M1 closure method, the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method) in 1 - D and 2 - D core - collapse supernova simulations, the author explores the performance differences of these methods in the supernova explosion mechanism. The key points of concern in the paper include: 1. **Selection of neutrino transport methods**: The paper selects three different neutrino transport methods, namely: - **M1 closure method**: It is solved based on the first two moments of the Boltzmann equation and is the most accurate but computationally most expensive method. - **Isotropic Diffusion Source Approximation (IDSA)**: The neutrino distribution function is decomposed into free - streaming and trapped parts, which are solved separately and coupled through a diffusion source term. - **Advanced Spectral Leakage (ASL) method**: Based on the grey leakage scheme, it carries discretized neutrino energy spectrum information and models the trapped neutrino components. 2. **Consistency of simulation settings**: To ensure the fairness of the comparison, the author uses the same grid settings, hydrodynamic and gravitational solvers in all three transport methods, thereby eliminating the interference of other factors. 3. **Comparison of key physical quantities**: The paper focuses on analyzing the influence of different transport methods on the following physical quantities: - **Proto - neutron star (PNS) radius and mass**: It is found that different transport methods have a relatively small influence on the PNS radius and mass. - **Neutrino luminosity and average energy**: It is found that there are quantitative differences of about 10 - 20% in neutrino luminosity and average energy caused by different transport methods, but the overall trend is relatively consistent. - **Properties of the gain region**: Including the gain region mass and net heating rate, as well as the PNS core convection intensity, these properties have large differences among different transport methods. - **Computational cost**: The computational costs of different transport methods are compared to evaluate their feasibility in practical applications. 4. **Other influencing factors**: The paper also explores the influence of the initial model (progenitor star), nuclear equation of state (EoS) and random perturbations on the simulation results, and finds that the influence of these factors on key physical quantities is comparable to that of different transport methods. In summary, this paper aims to reveal the advantages and disadvantages of different neutrino transport methods in core - collapse supernova simulations through systematic comparative research, providing a reference for future research.