No Arabic abstract
We compare Newtonian three-flavor multigroup Boltzmann (MGBT) and (Bruenns) multigroup flux-limited diffusion (MGFLD) neutrino transport in postbounce core collapse supernova environments. We focus our study on quantities central to the postbounce neutrino heating mechanism for reviving the stalled shock. Stationary-state three-flavor neutrino distributions are developed in thermally and hydrodynamically frozen time slices obtained from core collapse and bounce simulations that implement Lagrangian hydrodynamics and MGFLD neutrino transport. Most important, we find, for a region above the gain radius, net heating rates for MGBT that are as much as ~2 times the corresponding MGFLD rates, and net cooling rates below the gain radius that are typically ~0.8 times the MGFLD rates. These differences stem from differences in the neutrino luminosities and mean inverse flux factors, which can be as much as 11% and 24%, respectively. They are greatest at earlier postbounce times for a given progenitor mass and, for a given postbounce time, greater for greater progenitor mass. We discuss the ramifications these new results have for the supernova mechanism.
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.
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.
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 investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with an ``optimistic 15 solar mass precollapse model, which is representative of the class of stars with compact iron cores. This model is evolved through core collapse and bounce in one dimension using multigroup (neutrino-energy--dependent) flux-limited diffusion (MGFLD) neutrino transport and Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 106 ms after bounce at a time when there is a well-developed gain region, and proceeds for 400 ms. We couple two-dimensional (PPM) hydrodynamics to one-dimensional MGFLD neutrino transport. At 225 ms after bounce we see large-scale convection behind the shock, characterized by high-entropy, mushroom-like, expanding upflows and dense, low-entropy, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial convection velocities become supersonic just below the shock, reaching magnitudes in excess of 10^9 cm/sec. Eventually, however, the shock recedes to smaller radii, and at about 500 ms after bounce there is no evidence in our simulation of an explosion or of a developing explosion. Failure in our ``optimistic 15 solar mass Newtonian model leads us to conclude that it is unlikely, at least in our approximation, that neutrino-driven convection will lead to explosions for more massive stars with fatter iron cores or in cases in which general relativity is included.
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.