Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of Stellar Core Collapse, Bounce, and Postbounce Evolution with Boltzmann Neutrino Transport, and Implications for the Core Collapse Supernova Mechanism
82
0
0.0
(
0
)
Authors
A. Mezzacappa
Ask ChatGPT about the research
No Arabic abstract
In this paper, we present results from a simulation of stellar core collapse, bounce, and postbounce evolution with Boltzmann neutrino transport. We motivate the development of our Boltzmann solver in light of the sensitivity of the neutrino-heating core collapse supernova paradigm to details in the neutrino transport, particularly near the neutrinospheres, where the neutrinos are neither diffusing nor free streaming and a kinetic description is necessary, and in light of the mixed outcomes and transport approximations used in all prior supernova models in both one and two dimensions. We discuss the implications of our findings for the supernova mechanism and future supernova research. We also present the results of a Boltzmann transport prediction of the early neutrino light curves in the model included here.
rate research
Read More
We present self-consistent general relativistic simulations of stellar core collapse, bounce, and postbounce evolution for 13, 15, and 20 solar mass progenitors in spherical symmetry. Our simulations implement three-flavor Boltzmann neutrino transport and standard nuclear physics. The results are compared to our corresponding simulations with Newtonian hydrodynamics and O(v/c) Boltzmann transport.
General relativistic multi-group and multi-flavor Boltzmann neutrino transport in spherical symmetry adds a new level of detail to the numerical bridge between microscopic nuclear and weak interaction physics and the macroscopic evolution of the astrophysical object. Although no supernova explosions are obtained, we investigate the neutrino luminosities in various phases of the postbounce evolution for a wide range of progenitor stars between 13 and 40 solar masses. The signal probes the dynamics of material layered in and around the protoneutron star and is, within narrow limits, sensitive to improvements in the weak interaction physics. Only changes that dramatically exceed physical limitations allow experiments with exploding models. We discuss the differences in the neutrino signal and find the electron fraction in the innermost ejecta to exceed 0.5 as a consequence of thermal balance and weak equilibrium at the masscut.
We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashis nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there.
We investigate the criteria for successful core-collapse supernova explosions by the neutrino mechanism. We find that a critical-luminosity/mass-accretion-rate condition distinguishes non-exploding from exploding models in hydrodynamic one-dimensional (1D) and two-dimensional (2D) simulations. We present 95 such simulations that parametrically explore the dependence on neutrino luminosity, mass accretion rate, resolution, and dimensionality. While radial oscillations mediate the transition between 1D accretion (non-exploding) and exploding simulations, the non-radial standing accretion shock instability characterizes 2D simulations. We find that it is useful to compare the average dwell time of matter in the gain region with the corresponding heating timescale, but that tracking the residence time distribution function of tracer particles better describes the complex flows in multi-dimensional simulations. Integral quantities such as the net heating rate, heating efficiency, and mass in the gain region decrease with time in non-exploding models, but for 2D exploding models, increase before, during, and after explosion. At the onset of explosion in 2D, the heating efficiency is $sim$2% to $sim$5% and the mass in the gain region is $sim$0.005 M$_{sun}$ to $sim$0.01 M$_{sun}$. Importantly, we find that the critical luminosity for explosions in 2D is $sim$70% of the critical luminosity required in 1D. This result is not sensitive to resolution or whether the 2D computational domain is a quadrant or the full 180$^{circ}$. We suggest that the relaxation of the explosion condition in going from 1D to 2D (and to, perhaps, 3D) is of a general character and is not limited by the parametric nature of this study.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
