No Arabic abstract
Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, Agile-Boltztran of the Oak Ridge-Basel group and Vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 solar mass star with Newtonian gravity and a 15 solar mass star with relativistic gravity.
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.
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 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.
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our codes capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code.