No Arabic abstract
We investigate the explosion of stars with zero-age main-sequence masses between 20 and 35 solar masses and varying degrees of rotation and magnetic fields including ones commonly considered progenitors of gamma-ray bursts (GRBs). The simulations, combining special relativistic magnetohydrodynamics, a general relativistic approximate gravitational potential, and two-moment neutrino transport, demonstrate the viability of different scenarios for the post-bounce evolution. Having formed a highly massive proto-neutron star (PNS), several models launch successful explosions, either by the standard supernova mechanism based on neutrino heating and hydrodynamic instabilities or by magnetorotational processes. It is, however, quite common for the PNS to collapse to a black hole (BH) within a few seconds. Others might produce proto-magnetar-driven explosions. We explore several ways to describe the different explosion mechanisms. The competition between the timescales for advection of gas through the gain layer and heating by neutrinos provides an approximate explanation for models with insignificant magnetic fields. The fidelity of this explosion criterion in the case of rapid rotation can be improved by accounting for the strong deviations from spherical symmetry and mixing between pole and equator. We furthermore study an alternative description including the ram pressure of the gas falling through the shock. Magnetically driven explosions tend to arise from a strongly magnetised region around the polar axis. In these cases, the onset of the explosion corresponds to the equality between the advection timescale and the timescale for the propagation of Alfven waves through the gain layer.
We explore the influence of non-axisymmetric modes on the dynamics of the collapsed core of rotating, magnetized high-mass stars in three-dimensional simulations of a rapidly rotating star with an initial mass of $M_{ZAMS}$ = 35 solar masses endowed with four different pre-collapse configurations of the magnetic field, ranging from moderate to very strong field strength and including the field predicted by the stellar evolution model. The model with the weakest magnetic field achieves shock revival due to neutrino heating in a gain layer characterized by a large-scale, hydrodynamic m = 1 spiral mode. Later on, the growing magnetic field of the proto-neutron star launches weak outflows into the early ejecta. Their orientation follows the evolution of the rotational axis of the proto-neutron star, which starts to tilt from the original orientation due to the asymmetric accretion flows impinging on its surface. The models with stronger magnetization generate mildly relativistic, magnetically driven polar outflows propagating over a distance of $10^4$ km within a few 100 ms. These jets are stabilized against disruptive non-axisymmetric instabilities by their fast propagation and by the shear of their toroidal magnetic field. Within the simulation times of around 1 s, the explosions reach moderate energies and the growth of the proto-neutron star masses ceases at values substantially below the threshold for black hole formation, which, in combination with the high rotational energies, might suggest a possible later proto-magnetar activity.
We assess the variance of the post-collapse evolution remnants of compact, massive, low-metallicity stars, under small changes in the degrees of rotation and magnetic field of selected pre-supernova cores. These stellar models are commonly considered progenitors of long gamma-ray bursts. The fate of the proto-neutron star (PNS) formed after collapse, whose mass may continuously grow due to accretion, critically depends on the poloidal magnetic field strength at bounce. Should the poloidal magnetic field be sufficiently weak, the PNS collapses to a black hole (BH) within a few seconds. Models on this evolutionary track contain promising collapsar engines. Poloidal magnetic fields smooth over large radial scales (e.g. dipolar fields) or slightly augmented with respect to the original pre-supernova core yield long-lasting PNSs. In these models, BH formation is avoided or staved off for a long time, hence, they may produce proto-magnetars (PMs). Some of our PM candidates have been run for $lesssim 10,$s after core bounce, but they have not entered the Kelvin-Helmholtz phase yet. Among these models, some display episodic events of spin-down during which we find properties broadly compatible with the theoretical expectations for PMs ($M_{PNS} approx 1.85,M_odot - 2.5,M_odot$, $bar{P}_{PNS} approx 1.5 - 4,$ms, and $b^{rm surf}_{PNS} lesssim 10^{15},$G) and their very collimated supernova ejecta has nearly reached the stellar surface with (still growing) explosion energies $gtrsim 2times 10^{51},$erg.
I summarize what we have learned about the nature of stars that ultimately explode as core-collapse supernovae from the examination of images taken prior to the explosion. By registering pre-supernova and post-supernova images, usually taken at high resolution using either space-based optical detectors, or ground-based infrared detectors equipped with laser guide star adaptive optics systems, nearly three dozen core-collapse supernovae have now had the properties of their progenitor stars either directly measured or (more commonly) constrained by establishing upper limits on their luminosities. These studies enable direct comparison with stellar evolution models that, in turn, permit estimates of the progenitor stars physical characteristics to be made. I review progenitor characteristics (or constraints) inferred from this work for each of the major core-collapse supernova types (II-Plateau, II-Linear, IIb, IIn, Ib/c), with a particular focus on the analytical techniques used and the processes through which conclusions have been drawn. Brief discussion of a few individual events is also provided, including SN 2005gl, a type IIn supernova that is shown to have had an extremely luminous -- and thus very massive -- progenitor that exploded shortly after a violent, luminous blue variable-like eruption phase, contrary to standard theoretical predictions.
Recent multi-dimensional simulations of core-collapse supernovae are producing successful explosions and explosion-energy predictions. In general, the explosion-energy evolution is monotonic and relatively smooth, suggesting a possible analytic solution. We derive analytic solutions for the expansion of the gain region under the following assumptions: spherical symmetry, one-zone shell, and powered by neutrinos and $alpha$ particle recombination. We consider two hypotheses: I) explosion energy is powered by neutrinos and $alpha$ recombination, II) explosion energy is powered by neutrinos alone. Under these assumptions, we derive the fundamental dimensionless parameters and analytic scalings. For the neutrino-only hypothesis (II), the asymptotic explosion energy scales as $E_{infty} approx 1.5 M_g v_0^2 eta^{2/3}$, where $M_g$ is the gain mass, $v_0$ is the free-fall velocity at the shock, and $eta$ is a ratio of the heating and dynamical time scales. Including both neutrinos and recombination (hypothesis I), the asymptotic explosion energy is $E_{infty} approx M_g v_0^2 (1.5eta^{2/3} + beta f(rho_0))$, where $beta$ is the dimensionless recombination parameter. We use Bayesian inference to fit these analytic models to simulations. Both hypotheses fit the simulations of the lowest progenitor masses that tend to explode spherically. The fits do not prefer hypothesis I or II; however, prior investigations suggest that $alpha$ recombination is important. As expected, neither hypothesis fits the higher-mass simulations that exhibit aspherical explosions. In summary, this explosion-energy theory is consistent with the spherical explosions of low progenitor masses; the inconsistency with higher progenitor-mass simulations suggests that a theory for them must include aspherical dynamics.
We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_odot$, 15 $M_odot$, and 27 $M_odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.