No Arabic abstract
We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino transport, including all relevant $mathcal{O}(v/c)$ terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate $mathcal{O}(v/c)$ terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 milliseconds after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ``ray-by-ray approach employed by all other groups may be compromising their results. We show that ``ray-by-ray calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion.
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 investigate the impact of rotation and magnetic fields on the dynamics and gravitational wave emission in 2D core-collapse supernova simulations with neutrino transport. We simulate 16 different models of $15,M_odot$ and $39,M_odot$ progenitor stars with various initial rotation profiles and initial magnetic fields strengths up to $10^{12}, mathrm{G}$, assuming a dipolar field geometry in the progenitor. Strong magnetic fields generally prove conducive to shock revival, though this trend is not without exceptions. The impact of rotation on the post-bounce dynamics is more variegated, in line with previous studies. A significant impact on the time-frequency structure of the gravitational wave signal is found only for rapid rotation or strong initial fields. For rapid rotation, the angular momentum gradient at the proto-neutron star surface can appreciably affect the frequency of the dominant mode, so that known analytic relations for the high-frequency emission band no longer hold. In case of two magnetorotational explosion models, the time-frequency structure of the post-bounce emission appears rather different from neutrino-driven explosions. In one of these two models, a new high-frequency emission component of significant amplitude emerges about $200, mathrm{ms}$ after the burst of gravitational wave emission around shock revival has subsided. This emission is characterised by broad-band power well into the kHz range. Its emission mechanism remains unclear and needs to be investigated further. We also estimate the maximum detection distances for our waveforms. The magnetorotational models do not stick out for higher detectability during the post-bounce and explosion phase.
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.
Most supernova explosions accompany the death of a massive star. These explosions give birth to neutron stars and black holes and eject solar masses of heavy elements. However, determining the mechanism of explosion has been a half-century journey of great complexity. In this paper, we present our perspective of the status of this theoretical quest and the physics and astrophysics upon which its resolution seems to depend. The delayed neutrino-heating mechanism is emerging as a robust solution, but there remain many issues to address, not the least of which involves the chaos of the dynamics, before victory can unambiguously be declared. It is impossible to review in detail all aspects of this multi-faceted, more-than-half-century-long theoretical quest. Rather, we here map out the major ingredients of explosion and the emerging systematics of the observables with progenitor mass, as we currently see them. Our discussion will of necessity be speculative in parts, and many of the ideas may not survive future scrutiny. Some statements may be viewed as informed predictions concerning the numerous observables that rightly exercise astronomers witnessing and diagnosing the supernova Universe. Importantly, the same explosion in the inside, by the same mechanism, can look very different in photons, depending upon the mass and radius of the star upon explosion. A 10$^{51}$-erg (one Bethe) explosion of a red supergiant with a massive hydrogen-rich envelope, a diminished hydrogen envelope, no hydrogen envelope, and, perhaps, no hydrogen envelope or helium shell all look very different, yet might have the same core and explosion evolution.
We investigate core-collapse supernova (CCSN) nucleosynthesis with self-consistent, axisymmetric (2D) simulations performed using the radiation-hydrodynamics code Chimera. Computational costs have traditionally constrained the evolution of the nuclear composition within multidimensional CCSN models to, at best, a 14-species $alpha$-network capable of tracking only $(alpha,gamma)$ reactions from $^{4}$He to $^{60}$Zn. Such a simplified network limits the ability to accurately evolve detailed composition and neutronization or calculate the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. However, limitations such as poor spatial resolution of the tracer particles, inconsistent thermodynamic evolution, including misestimation of expansion timescales, and uncertain determination of the multidimensional mass-cut at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of such uncertainties for four self-consistent axisymmetric CCSN models initiated from stellar metallicity, non-rotating progenitors of 12 $M_odot$, 15 $M_odot$, 20 $M_odot$, and 25 $M_odot$ and evolved with the smaller $alpha$-network to more than 1 s after the launch of an explosion.