No Arabic abstract
We study the effects of the magnetic field on the dynamics of non-rotating stellar cores by performing two-dimensional (2D), magnetohydrodynamics (MHD) simulations. To this end, we have updated our neutrino-radiation-hydrodynamics supernova code to include MHD employing a divergence cleaning method with both careful treatments of finite volume and area reconstructions. By changing the initial strength of the magnetic field, the evolution of $15.0$, $18.4$ and $27.0$ $M_odot$ presupernova progenitors is investigated. An intriguing finding in our study is that the neutrino-driven explosion occurs regardless of the strength of the initial magnetic field. For the 2D models presented in this work, the neutrino heating is the main driver for the explosion, whereas the magnetic field secondary contributes to the pre-explosion dynamics. Our results show that the strong magnetic field weakens the growth of the neutrino-driven turbulence in the small scale compared to the weak magnetic field. This results in the slower increase of the turbulent kinetic energy in the postshock region, leading to the slightly delayed onset of the shock revival for models with the stronger initial magnetic field.
Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($gtrsim 50 %$), but that the time-integrated signal has a relatively modest asymmetry ($lesssim 20 %$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $sim$20% if the distance to the CCSN is $lesssim 6$ kpc.
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.
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 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.