Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 M sub( )blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times ...shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l = 2 and l = 1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987 A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to similar to 3300 km s super(-1), strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km s super(-1), and a final prolate anisotropy of the inner ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2 x 10 super(51) erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The first effect protects the (fastest) clumps from interacting with the strong reverse shock that forms below the He/H composition interface, by keeping their propagation timescale through the He-core shorter than the reverse shock formation time. This ensures that the outward motion of the clumps remains always subsonic, and that thus their energy dissipation is minimal (in contrast to the supersonic case). The second effect is responsible for the strong inward mixing of hydrogen: the aspherical shock deposits large amounts of vorticity into the He/H interface layer at early times (around t = 100 s). This triggers the growth of a strong Richtmyer-Meshkov instability that results in a global anisotropy of the inner ejecta at late times (i.e. around t = 10 000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987 A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event.
We study here the formation of heavy r-process nuclei in the high-entropy environment of rapidly expanding neutrino-driven winds from compact objects. In particular, we explore the sensitivity of the ...element creation in the $A \ga 130$ region to the low-temperature behavior of the outflows. For this purpose we employ a simplified model of the dynamics and of the thermodynamical evolution for radiation dominated, adiabatic outflows. It consists of a first stage of fast, exponential cooling with timescale $\tau_\mathrm{dyn}$, followed by a second phase of slower evolution, assuming either constant density and temperature or a power-law decay of these quantities. These cases describe a strong deceleration or decreasing acceleration of the transsonic outflows, respectively, and thus are supposed to capture the most relevant effects associated with a change in the wind expansion behavior at large radii, for example because of the collision with the slower, preceding supernova ejecta and the possible presence of a wind termination shock. We find that for given entropy, expansion timescale, and proton-to-baryon ratio not only the transition temperature between the two expansion phases can make a big difference in the formation of the platinum peak, but also the detailed cooling law during the later phase. Because the nuclear photodisintegration rates between about $2\times 10^8$ K and roughly 109 K are more sensitive to the temperature than the neutron-capture rates are to the free neutron density, a faster cooling but continuing high neutron density in this temperature regime allow the r-process path to move closer to the neutron-drip line. With low (γ,n)- but high β-decay rates, the r-processing does then not proceed through a (γ, n)-(n, γ) equilibrium but through a quasi-equilibrium of (n, γ)-reactions and β-decays, as recently also pointed out by Wanajo. Unless the transition temperature and corresponding (free neutron) density become too low ($T \la 2\times 10^8$ K), a lower temperature or faster temperature decline during the slow, late evolution phase therefore allow for a stronger appearance of the third abundance peak.
It has been recently proposed that the shocked surface layers of exploding O-Ne-Mg cores provide the conditions for r-process nucleosynthesis, because their rapid expansion and high entropies enable ...heavy r-process isotopes to form even in an environment with very low initial neutron excess of the matter. We show here that the most sophisticated available hydrodynamic simulations (in spherical and axial symmetry) do not support this new r-process scenario because they fail to provide the necessary conditions of temperature, entropy, and expansion timescale by significant factors. This suggests that, either the formation of r-process elements works differently than suggested by Ning et al. (ApJ, 667, L159, NQM07), or that some essential core properties with influence on the explosion dynamics might be different from those predicted by Nomoto's progenitor model.
Theory of core-collapse supernovae Janka, H.-Th; Langanke, K.; Marek, A. ...
Physics reports,
04/2007, Letnik:
442, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Advances in our understanding and the modeling of stellar core-collapse and supernova explosions over the past 15 years are reviewed, concentrating on the evolution of hydrodynamical simulations, the ...description of weak interactions and nuclear equation of state effects, and new insights into the nucleosynthesis occurring in the early phases of the explosion, in particular the neutrino-p process. The latter is enabled by the proton-richness of the early ejecta, which was discovered because significant progress has been made in the treatment of neutrino transport and weak interactions. This progress has led to a new generation of sophisticated Newtonian and relativistic hydrodynamics simulations in spherical symmetry. Based on these, it is now clear that the prompt bounce-shock mechanism is not the driver of supernova explosions, and that the delayed neutrino-heating mechanism can produce explosions without the aid of multi-dimensional processes only if the progenitor star has an ONeMg core inside a very dilute He-core, i.e., has a mass in the
8
–
10
M
⊙
range. Hydrodynamic instabilities of various kinds have indeed been recognized to occur in the supernova core and to be of potential importance for the explosion. Neutrino-driven explosions, however, have been seen in two-dimensional simulations with sophisticated neutrino transport so far only when the star has a small iron core and low density in the surrounding shells as being found in stars near
10
–
11
M
⊙
. The explosion mechanism of more massive progenitors is still a puzzle. It might involve effects of three-dimensional hydrodynamics or might point to the relevance of rapid rotation and magnetohydrodynamics, or to still incompletely explored properties of neutrinos and the high-density equation of state.
By means of three-dimensional hydrodynamic simulations with an Eulerian PPM code we investigate the time-dependent evolution and properties of accretion tori around nonrotating and rotating ...stellar-mass black holes, using a pseudo-Newtonian (Paczynski & Wiita or Artemova-Bjoernsson-Novikov) potential to approximate the effects of general relativity. The simulations are performed with three nested Cartesian grids to ensure sufficient resolution near the central black hole on the one hand and a large computational volume on the other. The black hole and torus are considered as the remnant of a binary neutron star or neutron-star black-hole merger. Referring to results from previous hydrodynamical simulations of such events, we assume the initial configurations to consist of a black hole with a mass of about 4\,M_\odot girded by a toroidal accretion disk with a mass in the range from about 0.01 M_\odot to 0.2 M_\odot. We simulate the torus evolution without and with physical shear viscosity, employing a simple \alpha-model for the gas viscosity. As in our previous work on merging neutron star binaries and neutron star/black hole binaries, we use the equation of state of Lattimer and Swesty. The energy loss and lepton number change due to neutrino emission from the hot torus are treated by a neutrino-trapping scheme. The energy deposition by neutrino-antineutrino annihilation around the disk is evaluated in a post-processing step. The time-dependent efficiency of converting gravitational energy to neutrinos, expressed by the ratio of neutrino luminosity to accretion rate of rest-mass energy, can reach maximum values of up to about 10%. The corresponding efficiency of converting neutrino energy into a pair-photon fireball by neutrino annihilation peaks at values of several percent. Interestingly, we find that the rate of neutrino- antineutrino annihilation decays with time much less steeply than the total neutrino luminosity does with the decreasing gas mass of the torus, because the ongoing protonization of the initially neutron-rich disk matter leads to a rather stable product of neutrino and antineutrino luminosities. The neutrino luminosity and total energy release of the torus increase steeply with higher viscosity, larger torus mass, and larger black hole spin in corotation with the disk, in particular when the spin parameter is a \ga 0.8. The latter dependence is moderated in case of a high disk viscosity. For rotation rates as expected for post-merger black holes (a\ga 0.5) and reasonable values of the alpha viscosity of the torus ( \alpha \sim 0.1), torus masses in the investigated range can release sufficient energy in neutrinos to account for the energetics of the well-localized short gamma-ray bursts recently detected by Hete and Swift, if collimation of the ultrarelativistic outflows into about 1% of the sky is invoked, as predicted by recent hydrodynamic jet simulations.
Context.We investigate the deposition of energy and momentum due to the annihilation of neutrinos (ν) and antineutrinos ($\bar\nu$) in the vicinity of steady, axisymmetric accretion tori around ...stellar-mass black holes (BHs). This process is widely considered as an energy source for driving ultrarelativistic outflows with the potential to produce gamma-ray bursts. Aims.We analyze the influence of general relativistic (GR) effects in combination with different neutrinosphere properties on the $\nu\bar\nu$-annihilation efficiency and spatial distribution of the energy deposition rate. Methods.Assuming axial symmetry, we numerically compute the annihilation rate 4-vector. For this purpose, we construct the local neutrino distribution by ray-tracing neutrino trajectories in a Kerr space-time using null geodesics. We vary the value of the dimensionless specific angular momentum a of the central BH, which provides the gravitational field in our models. We also study different shapes of the neutrinospheres, spheres, thin disks, and thick accretion tori, whose structure ranges from idealized tori to equilibrium non-selfgravitating matter distributions. Furthermore, we compute Newtonian models where the influence of the gravitational field on the annihilation process is neglected. Results.Compared to Newtonian calculations, GR effects increase the total annihilation rate measured by an observer at infinity by a factor of two when the neutrinosphere is a thin disk, but the increase is only ≈$25\%$ for toroidal and spherical neutrinospheres. Comparing cases with similar luminosities, thin disk models yield the highest energy deposition rates by $\nu\bar\nu$-annihilation, and spherical neutrinospheres the lowest ones, independently of whether GR effects are included. Increasing a from 0 to 1 enhances the energy deposition rate measured by an observer at infinity by roughly a factor of 2 due to the change of the inner radius of the neutrinosphere. General relativity and rotation cause important differences in the spatial distribution of the energy deposition rate by $\nu\bar\nu$-annihilation.
We investigate the possibility approximating relativistic effects in hydrodynamical simulations of stellar core collapse and post-bounce evolution by using a modified gravitational potential in an ...otherwise standard Newtonian hydrodynamic code. Different modifications of a previously introduced effective relativistic potential are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hydrodynamic simulations with two different codes are compared with fully relativistic results. One code is applied for one- and two-dimensional calculations with a simple equation of state, and employs either the modified effective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale supernova runs including a microphysical equation of state and neutrino transport based on the solution of the Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully relativistic solutions in spherical symmetry, leading to significant improvements compared to previously published approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation.
Performing two-dimensional hydrodynamic simulations including a detailed treatment of the equation of state of the stellar plasma and for the neutrino transport and interactions, we investigate here ...the interplay between different kinds of non-radial hydrodynamic instabilities that can play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, which is driven by the negative entropy gradients caused by neutrino heating or by variations in the shock strength in transient phases of shock expansion and contraction, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in those cases where the gas in the postshock region is rapidly advected towards the gain radius, the growth of such a buoyancy instability can be suppressed. In this situation the shock and postshock flow can nevertheless develop non- radial asymmetry with an oscillatory growth in the amplitude. This phenomenon has been termed "standing (or spherical) accretion shock instability" (SASI). It is shown here that the SASI oscillations can trigger convective instability, and like the latter, they lead to an increase in the average shock radius and in the mass of the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter accreted through the neutrino-heating layer and thus enhance the neutrino energy deposition in support of the neutrino-driven explosion mechanism. A rapidly contracting and more compact nascent neutron star turns out to be favorable for explosions, because the accretion luminosity and neutrino heating are greater and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI observed in our simulations agrees with the one estimated for the advective-acoustic cycle (AAC), in which perturbations are carried by the accretion flow from the shock to the neutron star and pressure waves close an amplifying global feedback loop. A variety of other features in our models, as well as differences in their behavior, can also be understood on the basis of the AAC hypothesis. The interpretation of the SASI in our simulations as a purely acoustic phenomenon, however, appears difficult.
Performing two-dimensional hydrodynamic simulations including a detailed treatment of the equation of state of the stellar plasma and for the neutrino transport and interactions, we investigate here ...the interplay between different kinds of non-radial hydrodynamic instabilities that can play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, which is driven by the negative entropy gradients caused by neutrino heating or by variations in the shock strength in transient phases of shock expansion and contraction, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in those cases where the gas in the postshock region is rapidly advected towards the gain radius, the growth of such a buoyancy instability can be suppressed. In this situation the shock and postshock flow can nevertheless develop non-radial asymmetry with an oscillatory growth in the amplitude. This phenomenon has been termed “standing (or spherical) accretion shock instability” (SASI). It is shown here that the SASI oscillations can trigger convective instability, and like the latter, they lead to an increase in the average shock radius and in the mass of the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter accreted through the neutrino-heating layer and thus enhance the neutrino energy deposition in support of the neutrino-driven explosion mechanism. A rapidly contracting and more compact nascent neutron star turns out to be favorable for explosions, because the accretion luminosity and neutrino heating are greater and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI observed in our simulations agrees with the one estimated for the advective-acoustic cycle (AAC), in which perturbations are carried by the accretion flow from the shock to the neutron star and pressure waves close an amplifying global feedback loop. A variety of other features in our models, as well as differences in their behavior, can also be understood on the basis of the AAC hypothesis. The interpretation of the SASI in our simulations as a purely acoustic phenomenon, however, appears difficult.
Two-dimensional hydrodynamic simulations of stellar core collapse are presented which for the first time were performed by solving the Boltzmann equation for the neutrino transport including a ...state-of-the-art description of neutrino interactions. Stellar rotation is also taken into account. Although convection develops below the neutrinosphere and in the neutrino-heated region behind the supernova shock, the models do not explode. This suggests missing physics, possibly with respect to the nuclear equation of state and weak interactions in the subnuclear regime. However, it might also indicate a fundamental problem with the neutrino-driven explosion mechanism.