-
On the characteristics of fast neutrino flavor instabilities in three-dimensional core-collapse supernova models
Authors:
Sajad Abbar,
Francesco Capozzi,
Robert Glas,
H. -Thomas Janka,
Irene Tamborra
Abstract:
We assess the occurrence of fast neutrino flavor instabilities in two three-dimensional state-of-the-art core-collapse supernova simulations performed using a two-moment three-species neutrino transport scheme: one with an exploding 9$\mathrm{M_{\odot}}$ and one with a non-exploding 20$\mathrm{M_{\odot}}$ model. Apart from confirming the presence of fast instabilities occurring within the neutrino…
▽ More
We assess the occurrence of fast neutrino flavor instabilities in two three-dimensional state-of-the-art core-collapse supernova simulations performed using a two-moment three-species neutrino transport scheme: one with an exploding 9$\mathrm{M_{\odot}}$ and one with a non-exploding 20$\mathrm{M_{\odot}}$ model. Apart from confirming the presence of fast instabilities occurring within the neutrino decoupling and the supernova pre-shock regions, we detect flavor instabilities in the post-shock region for the exploding model. These instabilities are likely to be scattering-induced. In addition, the failure in achieving a successful explosion in the heavier supernova model seems to seriously hinder the occurrence of fast instabilities in the post-shock region. This is a consequence of the large matter densities behind the stalled or retreating shock, which implies high neutrino scattering rates and thus more isotropic distributions of neutrinos and antineutrinos. Our findings suggest that the supernova model properties and the fate of the explosion can remarkably affect the occurrence of fast instabilities. Hence, a larger set of realistic hydrodynamical simulations of the stellar collapse is needed in order to make reliable predictions on the flavor conversion physics.
△ Less
Submitted 11 March, 2021; v1 submitted 11 December, 2020;
originally announced December 2020.
-
Gravitational-wave Signals From Three-dimensional Supernova Simulations With Different Neutrino-Transport Methods
Authors:
Haakon Andresen,
Robert Glas,
H-Thomas Janka
Abstract:
We compare gravitational-wave (GW) signals from eight three-dimensional simulations of core-collapse supernovae, using two different progenitors with zero-age main sequence masses of 9 and 20 solar masses. The collapse of each progenitor was simulated four times, at two different grid resolutions and with two different neutrino transport methods, using the Aenus-Alcar code. The main goal of this s…
▽ More
We compare gravitational-wave (GW) signals from eight three-dimensional simulations of core-collapse supernovae, using two different progenitors with zero-age main sequence masses of 9 and 20 solar masses. The collapse of each progenitor was simulated four times, at two different grid resolutions and with two different neutrino transport methods, using the Aenus-Alcar code. The main goal of this study is to assess the validity of recent concerns that the so-called "Ray-by-Ray+" (RbR+) approximation is problematic in core-collapse simulations and can adversely affect theoretical GW predictions. Therefore, signals from simulations using RbR+ are compared to signals from corresponding simulations using a fully multidimensional (FMD) transport scheme. The 9 solar-mass progenitor successfully explodes, whereas the 20 solar-mass model does not. Both the standing accretion shock instability and hot-bubble convection develop in the postshock layer of the non-exploding models. In the exploding models, neutrino-driven convection in the postshock flow is established around 100 ms after core bounce and lasts until the onset of shock revival. We can, therefore, judge the impact of the numerical resolution and neutrino transport under all conditions typically seen in non-rotating core-collapse simulations. We find excellent qualitative agreement in all GW features. We find minor quantitative differences between simulations, but find no systematic differences between simulations using different transport schemes. Resolution-dependent differences in the hydrodynamic behaviour of low-resolution and high-resolution models have a greater impact on the GW signals than consequences of the different transport methods. Furthermore, increasing the resolution decreases the discrepancies between models with different neutrino transport.
△ Less
Submitted 12 April, 2021; v1 submitted 20 November, 2020;
originally announced November 2020.
-
Fast Neutrino Flavor Instability in the Neutron-star Convection Layer of Three-dimensional Supernova Models
Authors:
Robert Glas,
H. -Thomas Janka,
Francesco Capozzi,
Manibrata Sen,
Basudeb Dasgupta,
Alessandro Mirizzi,
Guenter Sigl
Abstract:
Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission sel…
▽ More
Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged)
△ Less
Submitted 17 January, 2020; v1 submitted 30 November, 2019;
originally announced December 2019.
-
Effects of LESA in Three-Dimensional Supernova Simulations with Multi-Dimensional and Ray-by-Ray-plus Neutrino Transport
Authors:
Robert Glas,
H. -Thomas Janka,
Tobias Melson,
Georg Stockinger,
Oliver Just
Abstract:
A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar ne…
▽ More
A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar neutrino/hydrodynamics code, which treats the energy- and velocity-dependent transport of neutrinos of all flavors by a two-moment scheme with algebraic M1 closure. For each of the progenitors, results with fully multi-dimensional (FMD) neutrino transport and with ray-by-ray-plus (RbR+) approximation are considered for two different grid resolutions. While the 9 solar-mass models develop explosions, the 20 solar-mass progenitor does not explode with the employed version of simplified neutrino opacities. In all 3D models we observe the growth of substantial dipole amplitudes of the lepton-number (electron neutrino minus antineutrino) flux with stable or slowly time-evolving direction and overall properties fully consistent with the LESA phenomenon. Models with RbR+ transport develop LESA dipoles somewhat faster and with temporarily higher amplitudes, but the FMD calculations exhibit cleaner hemispheric asymmetries with a far more dominant dipole. In contrast, the RbR+ results display much wider multipole spectra of the neutrino-emission anisotropies with significant power also in the quadrupole and higher-order modes. Our results disprove speculations that LESA is a numerical artifact of RbR+ transport. We also discuss LESA as consequence of a dipolar convection flow inside of the nascent neutron star and establish, tentatively, a connection to Chandrasekhar's linear theory of thermal instability in spherical shells.
△ Less
Submitted 7 June, 2019; v1 submitted 26 September, 2018;
originally announced September 2018.
-
Three-Dimensional Core-Collapse Supernova Simulations with Multi-Dimensional Neutrino Transport Compared to the Ray-by-Ray-plus Approximation
Authors:
Robert Glas,
Oliver Just,
H. -Thomas Janka,
Martin Obergaulinger
Abstract:
Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ tr…
▽ More
Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ transport for both tested grid resolutions in the cases of a 20 solar-mass progenitor, which does not explode with the employed simplified set of neutrino opacities, and of an exploding 9 solar-mass model. This is in stark contrast to corresponding axisymmetric (2D) simulations, which confirm previous claims that the RbR+ approximation can foster explosions in 2D in particular in models with powerful axial sloshing of the stalled shock due to the standing accretion shock instability (SASI). However, while local and instantaneous variations of neutrino fluxes and heating rates can still be considerably higher with RbR+ transport in 3D, the time-averaged quantities are very similar to FMD results because of the absence of a fixed, artificial symmetry axis that channels the flow. Therefore, except for stochastic fluctuations, the neutrino signals and the post-bounce evolution of 3D simulations with FMD and RbR+ transport are also very similar, in particular for our calculations with the better grid resolution. Higher spatial resolution has clearly a more important impact than the differences by the two transport treatments. Our results back up the use of the RbR+ approximation for neutrino transport in 3D SN modeling.
△ Less
Submitted 5 February, 2019; v1 submitted 26 September, 2018;
originally announced September 2018.
-
Core-collapse supernova simulations in one and two dimensions: comparison of codes and approximations
Authors:
Oliver Just,
Robert Bollig,
Hans-Thomas Janka,
Martin Obergaulinger,
Robert Glas,
Shigehiro Nagataki
Abstract:
We present spherically symmetric (1D) and axisymmetric (2D) supernova simulations for a convection-dominated 9 Msun and a 20 Msun progenitor that develops violent activity by the standing-accretion-shock instability (SASI). We compare in detail the Aenus-Alcar code, which uses fully multidimensional two-moment neutrino transport with an M1 closure, with a ray-by-ray-plus (RbR+) version of this cod…
▽ More
We present spherically symmetric (1D) and axisymmetric (2D) supernova simulations for a convection-dominated 9 Msun and a 20 Msun progenitor that develops violent activity by the standing-accretion-shock instability (SASI). We compare in detail the Aenus-Alcar code, which uses fully multidimensional two-moment neutrino transport with an M1 closure, with a ray-by-ray-plus (RbR+) version of this code and with the Prometheus-Vertex code that employs RbR+ two-moment transport with a Boltzmann closure. Besides testing consequences of ignored non-radial neutrino-flux components in the RbR+ approximation, we also discuss the influence of various transport ingredients applied or not applied in recent literature, namely simplified neutrino-pair processes, neutrino-electron scattering, velocity-dependent and gravitational-redshift terms, and strangeness and many-body corrections for neutrino-nucleon scattering. Alcar and Vertex show excellent agreement in 1D and 2D despite a slightly but systematically smaller radius (~1km) and stronger convection of the proto-neutron star with Alcar. As found previously, the RbR+ approximation is conducive to explosions, but much less severely in the convection-dominated 9 Msun case than in the marginally exploding 20 Msun model, where the onset time of explosion also exhibits big stochastic variations, and the RbR+ approximation has no distinctly stronger supportive effect than simplified pair processes or strangeness and many-body corrections. Neglecting neutrino-electron scattering has clearly unfavorable effects for explosions, while ignoring velocity and gravitational-redshift effects can both promote or delay the explosion. The ratio of advection timescale to neutrino-heating timescale in 1D simulations is a sensitive indicator of the influence of physics ingredients on explosions also in multidimensional simulations.
△ Less
Submitted 18 September, 2018; v1 submitted 4 May, 2018;
originally announced May 2018.