We perform spherically symmetric simulations of the core collapse of a single progenitor star of zero age main sequence mass $M_{\rm ZAMS} = 15 \, M_{\odot}$ with two models of heavy sterile neutrinos in the mass range of hundred MeV$/c^2$. According to both models, these hypothetical particles are copiously produced in the center, stream outwards a subsequently decay releasing energy into final states (including neutrinos) of the Standard Model. We find that they can lead to a successful explosion in otherwise non-exploding progenitors. Depending on their unknown parameters (e.g., mass and coupling constants with matter), we obtain either no explosion or an explosion of one of two types, i.e., through heating of gas downstream of the stalled shock wave, similarly to the standard scenario for supernova explosions or through heating of gas at higher radii that ejects matter from the outer core or the envelope while the center continues to accrete matter. In both cases, the explosion energies can be very high. We presume that this new type of explosion would produce an electromagnetic signal that significantly differs from common events because of the relative absence of heavy elements in the ejecta. The combination of core-collapse simulations and astrophysical observations may further constrain the parameters of the sterile neutrinos.

We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through one- and two-dimensional test problems.

We discuss the relevance of the magnetorotational instability (MRI) in core-collapse supernovae (CCSNe). Our recent numerical studies show that in CCSNe, the MRI is terminated by parasitic instabilities of the Kelvin-Helmholtz type. To determine whether the MRI can amplify initially weak magnetic fields to dynamically relevant strengths in CCSNe, we performed three-dimensional simulations of a region close to the surface of a differentially rotating proto-neutron star in non-ideal magnetohydrodynamics with two different numerical codes. We find that under the conditions prevailing in proto-neutron stars, the MRI can amplify the magnetic field by (only) one order of magnitude. This severely limits the role of MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields.

Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfven waves, and the tearing mode instability using the MHD code Aenus. By comparing the simu- lation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of tearing modes we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast-magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) systems. We also determine the dependance of the numerical viscosity and resistivity on the time integration method, the spatial reconstruction scheme and (to a lesser extent) the Riemann solver employed in the simulations. From the measured results we infer the numerical resolution (as a function of the spatial reconstruction method) required to properly resolve the growth and saturation level of the magnetic field amplified by the magnetorotational instability in the post-collapsed core of massive stars. Our results show that it is to the best advantage to resort to ultra-high order methods (e.g., 9th-order Monotonicity Preserving method) to tackle this problem properly, in particular in three dimensional simulations.

Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of the amplification factor were proposed within the parasitic model. To determine the amplification factor and test which criterion is a better predictor of the MRI termination, we perform three-dimensional shearing-disc and shearing-box simula- tions of a region close to the surface of a differentially rotating proto-neutron star in non-ideal MHD with two different numerical codes. We find that independently of the initial magnetic field strength, the MRI channel modes can amplify the magnetic field by, at most, a factor of 100. Under the conditions found in proto-neutron stars a more realistic value for the mag- netic field amplification is of the order of 10. This severely limits the role of the MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields. A further amplification should therefore rely on other physical processes, such as for example an MRI-driven turbulent dynamo.

We study the influence of numerical methods and grid resolution on the termination of the magnetorotational instability (MRI) by means of parasitic instabilities in three-dimensional shearing-disc simulations reproducing typical conditions found in core-collapse supernovae. Whether or not the MRI is able to amplify weak magnetic fields in this context strongly depends, among other factors, on the amplitude at which its growth terminates. The qualitative results of our study do not depend on the numerical scheme. In all our models, MRI termination is caused by Kelvin-Helmholtz instabilities, consistent with theoretical predictions. Quantitatively, however, there are differences, but numerical convergence can be achieved even at relatively low grid resolutions if high-order reconstruction methods are used.

The magnetorotational instability (MRI) can be a powerful mechanism amplifying the magnetic field in core-collapse supernovae. Whether initially weak magnetic fields can be amplified by this instability to dynamically relevant strengths is still a matter of debate. One of the main uncertainties concerns the process that terminates the growth of the instability. Parasitic instabilities of both Kelvin-Helmholtz and tearing-mode type have been suggested to play a crucial role in this process, disrupting MRI channel flows and quenching magnetic field amplification. We perform two-dimensional and three-dimensional sheering-disc simulations of a differentially rotating protoneutron star layer in non-ideal magnetohydrodynamics with unprecedented high numerical accuracy, finding that Kelvin-Helmholtz parasitic modes dominate tearing modes in the regime of large hydrodynamic and magnetic Reynolds numbers, as encountered close to the surface of protoneutron stars. They also determine the maximum magnetic field stress achievable during the exponential growth of the MRI. Our results are consistent with the theory of parasitic instabilities based on a local stability analysis. To simulate the Kelvin-Helmholtz instabilities properly, a very high numerical resolution is necessary. Using ninth-order spatial reconstruction schemes, we find that at least eight grid zones per MRI channel are necessary to simulate the growth phase of the MRI and reach an accuracy of ~10 per cent in the growth rate, while more than ~60 zones per channel are required to achieve convergent results for the value of the magnetic stress at MRI termination.

We compare Newtonian and general relativistic descriptions of the stationary accretion of self-gravitating fluids onto compact bodies. Spherical symmetry and thin gas approximation are assumed. Luminosity depends, amongst other factors, on the temperature and the contribution of gas to the total mass, in both -- general relativistic ($L_{GR}$) and Newtonian ($L_N$) -- models. We discover a remarkable universal behaviour for transonic flows: the ratio of respective luminosities $L_{GR}/L_N$ is independent of the fractional mass of the gas and depends on asymptotic temperature. It is close to 1 in the regime of low asymptotic temperatures and can grow by one order of magnitude for high temperatures. These conclusions are valid for a wide range of polytropic equations of state.