No Arabic abstract
We apply our recently developed code for spherically symmetric, fully general relativistic (GR) Lagrangian hydrodynamics and multigroup flux-limited diffusion neutrino transport to examine the effects of GR on the hydrodynamics and transport during collapse, bounce, and the critical shock reheating phase of core collapse supernovae. Comparisons of models computed with GR versus Newtonian hydrodynamics show that collapse to bounce takes slightly less time in the GR limit, and that the shock propagates slightly farther out in radius before receding. After a secondary quasistatic rise in the shock radius, the shock radius declines considerably more rapidly in the GR simulations than in the corresponding Newtonian simulations. During the shock reheating phase, core collapse computed with GR hydrodynamics results in a substantially more compact structure from the center out to the stagnated shock. The inflow speed of material behind the shock is also increased. Comparisons also show that the luminosity and rms energy of any neutrino flavor during the shock reheating phase increases when switching from Newtonian to GR hydrodynamics, and decreases when switching from Newtonian to GR transport. This latter decrease in neutrino luminosities and rms energies is less in magnitude than the increase that arise when switching from Newtonian to GR hydrodynamics, with the result that a fully GR simulation gives higher neutrino luminosities and harder neutrino spectra than a fully Newtonian simulation of the same precollapse model.
An attempt is made to assess the significance of rotation in the core-collapse supernova phenomenon, from both observational and theoretical point of view. The data on supernovae particularly indicative of the role of rotation in the collapse-triggered explosion is emphasized. The problem of including the rotation of presupernova core into the supernova theory is considered. A two-dimensional classification scheme of core-collapse supernovae is proposed which unifies classical supernovae of type Ib/c and type II, hypernovae and some GRB events.
We continue our investigations of the magnetorotational collapse of stellar cores discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV potential used in our simulations catches several features of fully relativistic simulations quite well. It is able to correctly reproduce the behavior of models which show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. If this is not the case, the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential allowing for a more efficient amplification of the magnetic field by differential rotation. Sufficiently strong magnetic fields brake down the cores rotation and trigger a contraction phase to higher densities. Several models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce does not suffice to eject any mass, and the core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Instead it belongs to the signal type IV introduced by us in the first paper of this series. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and > 100 kHz denotes the generation of a jet-like outflow. [Abstract abbreviated]
Here we present the results from two sets of simulations, in two and three spatial dimensions. In two dimensions, the simulations include multifrequency flux-limited diffusion neutrino transport in the ray-by-ray-plus approximation, two-dimensional self gravity in the Newtonian limit, and nuclear burning through a 14-isotope alpha network. The three-dimensional simulations are model simulations constructed to reflect the post stellar core bounce conditions during neutrino shock reheating at the onset of explosion. They are hydrodynamics-only models that focus on critical aspects of the shock stability and dynamics and their impact on the supernova mechanism and explosion. In two dimensions, we obtain explosions (although in one case weak) for two progenitors (11 and 15 Solar mass models). Moreover, in both cases the explosion is initiated when the inner edge of the oxygen layer accretes through the shock. Thus, the shock is not revived while in the iron core, as previously discussed in the literature. The three-dimensional studies of the development of the stationary accretion shock instability (SASI) demonstrate the fundamentally new dynamics allowed when simulations are performed in three spatial dimensions. The predominant l=1 SASI mode gives way to a stable m=1 mode, which in turn has significant ramifications for the distribution of angular momentum in the region between the shock and proto-neutron star and, ultimately, for the spin of the remnant neutron star. Moreover, the three-dimensional simulations make clear, given the increased number of degrees of freedom, that two-dimensional models are severely limited by artificially imposed symmetries.
We explore with self-consistent 2D F{sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.
We present a new mechanism for core-collapse supernova explosions that relies upon acoustic power generated in the inner core as the driver. In our simulation using an 11-solar-mass progenitor, a strong advective-acoustic oscillation a la Foglizzo with a period of ~25-30 milliseconds (ms) arises ~200 ms after bounce. Its growth saturates due to the generation of secondary shocks, and kinks in the resulting shock structure funnel and regulate subsequent accretion onto the inner core. However, this instability is not the primary agent of explosion. Rather, it is the acoustic power generated in the inner turbulent region and most importantly by the excitation and sonic damping of core g-mode oscillations. An l=1 mode with a period of ~3 ms grows to be prominent around ~500 ms after bounce. The accreting protoneutron star is a self-excited oscillator. The associated acoustic power seen in our 11-solar-mass simulation is sufficient to drive the explosion. The angular distribution of the emitted sound is fundamentally aspherical. The sound pulses radiated from the core steepen into shock waves that merge as they propagate into the outer mantle and deposit their energy and momentum with high efficiency. The core oscillation acts like a transducer to convert accretion energy into sound. An advantage of the acoustic mechanism is that acoustic power does not abate until accretion subsides, so that it is available as long as it may be needed to explode the star. [abridged]