Showing papers by "Hans-Thomas Janka published in 2006"

Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport. I. Numerical method and results for a 15 M star

Abstract: Supernova models with a full spectral treatment of the neutrino transport are presented, employing the PROMETHEUS/VERTEX neutrinohydrodynamics code with a variable Eddington factor closure of the O(v/c) moments equations of neutrino number, energy, and momentum. Our ray-by-ray plus approximation developed for two- (or three-) dimensional problems assumes that the local neutrino distribution function is azimuthally symmetric around the radial direction, which implies that the nonradial flux components disappear. Other terms containing the angular velocity components are retained in the moments equations and establish a coupling of the transport at different latitudes by lateral derivatives. Also lateral components of the neutrino pressure gradients are included in the hydrodynamics equations. This approximative approach for neutrino transport in multi-dimensional environments is motivated and critically assessed with respect to its capabilities, limitations, and inaccuracies in the context of supernova simulations. In this first paper of a series, one- (1D) and two-dimensional (2D) core-collapse calculations for a (nonrotating) 15 M ○. star are discussed, uncertainties in the treatment of the equation of state - numerical and physical - are tested, Newtonian results are compared with simulations using a general relativistic potential, bremsstrahlung and interactions of neutrinos of different flavors are investigated, and the standard approximation in neutrino-nucleon interactions with zero energy transfer is replaced by rates that include corrections due to nucleon recoil, thermal motions, weak magnetism, and nucleon correlations. Models with the full implementation of the ray-by-ray plus spectral transport were found not to explode, neither in spherical symmetry nor in 2D when the computational grid is constrained to a lateral wedge (<±45°) around the equator. The success of previous two-dimensional simulations with grey, flux-limited neutrino diffusion can therefore not be confirmed. An explosion is obtained in 2D for the considered 15 M ○. progenitor, when the radial velocity terms in the neutrino momentum equation are omitted. This manipulation increases the neutrino energy density in the convective gain layer by about 20-30% and thus the integral neutrino energy deposition in this region by about a factor of two compared to the non-exploding 2D model with the full transport. The spectral treatment of the transport and detailed description of charged-current processes leads to proton-rich neutrino-heated ejecta, removing the problem that previous explosion models with approximate neutrino treatment overproduced N = 50 closed neutron shell nuclei by large factors.

Multidimensional supernova simulations with approximative neutrino transport. I. Neutron star kicks and the anisotropy of neutrino-driven explosions in two spatial dimensions

Abstract: We study hydrodynamic instabilities during the first seconds of core-collapse supernovae by means of 2D simulations with approximative neutrino transport and boundary conditions that parameterize the effects of the contracting neutron star and allow us to obtain sufficiently strong neutrino heating and, hence, neutrino-driven explosions. Confirming more idealised studies, as well as supernova simulations with spectral transport, we find that random seed perturbations can grow by hydrodynamic instabilities to a globally asymmetric mass distribution in the region between the nascent neutron star and the accretion shock, leading to a dominance of dipole ($l=1$) and quadrupole ($l=2$) modes in the explosion ejecta, provided the onset of the supernova explosion is sufficiently slower than the growth time scale of the low-mode instability. By gravitational and hydrodynamic forces, the anisotropic mass distribution causes an acceleration of the nascent neutron star, which lasts for several seconds and can propel the neutron star to velocities of more than 1000 km s -1 . Because the explosion anisotropies develop chaotically and change by small differences in the fluid flow, the magnitude of the kick varies stochastically. No systematic dependence of the average neutron star velocity on the explosion energy or the properties of the considered progenitors is found. Instead, the anisotropy of the mass ejection, and hence of the kick, seems to increase when the nascent neutron star contracts more quickly, and thus low-mode instabilities can grow more rapidly. Our more than 70 models separate into two groups, one with high and the other with low neutron star velocities and accelerations after one second of post-bounce evolution, depending on whether the $l=1$ mode is dominant in the ejecta or not. This leads to a bimodality of the distribution when the neutron star velocities are extrapolated to their terminal values. Establishing a link to the measured distribution of pulsar velocities, however, requires a much larger set of calculations and ultimately 3D modelling.

Nucleosynthesis in Early Supernova Winds. II. The Role of Neutrinos

Abstract: One of the outstanding unsolved riddles of nuclear astrophysics is the origin of the so-called p-process nuclei from A = 92 to 126. Both the lighter and heavier p-process nuclei are adequately produced in the neon and oxygen shells of ordinary Type II supernovae, but the origin of these intermediate isotopes, especially 92,94Mo and 96,98Ru, has long been mysterious. Here we explore the production of these nuclei in the neutrino-driven wind from a young neutron star. We consider such early times that the wind still contains a proton excess because the rates for νe and positron captures on neutrons are faster than those for the inverse captures on protons. Following a suggestion by Frohlich and coworkers, we also include the possibility that—in addition to the protons, α-particles, and heavy seed—a small flux of neutrons is maintained by the reaction p(e, e+)n. This flux of neutrons is critical in bridging the long waiting points along the path of the rp-process by (n, p) reactions. Using the unmodified ejecta histories from a recent two-dimensional supernova model by Janka and coworkers, we find synthesis of p-rich nuclei up to 102Pd, although our calculations do not show efficient production of 92Mo. If the entropy of these ejecta is increased by a factor of 2, the synthesis extends to 120Te. Still larger increases in entropy, which might reflect the role of magnetic fields or vibrational energy input neglected in the hydrodynamical model, result in the production of nuclei up to A ≈ 170. Elements synthesized in these more extreme outflows include numerous s- and p-process nuclei, and even some r-process nuclei can be synthesized in these proton-rich conditions.

Torus formation in neutron star mergers and well-localized short gamma-ray bursts

Abstract: Merging neutron stars (NSs) are hot candidates for the still enigmatic sources of short gamma-ray bursts (GRBs). If the central engines of the huge energy release are accreting relic black holes (BHs) of such mergers, it is important to understand how the properties of the BH-torus systems, in particular disc masses and mass and rotation rate of the compact remnant, are linked to the characterizing parameters of the NS binaries. For this purpose, we present relativistic smoothed particle hydrodynamic simulations with conformally flat approximation of the Einstein field equations and a physical, non-zero temperature equation of state. Thick disc formation is highlighted as a dynamical process caused by angular momentum transfer through tidal torques during the merging process of asymmetric systems or in the rapidly spinning triaxial post-merger object. Our simulations support the possibility that the first well-localized short and hard GRBs 050509b, 050709, 050724, 050813 have originated from NS merger events and are powered by neutrino-antineutrino annihilation around a relic BH-torus system. Using model parameters based on this assumption, we show that the measured GRB energies and durations lead to estimates for the accreted masses and BH mass accretion rates which are compatible with theoretical expectations. In particular, the low-energy output and short duration of GRB 050509b set a very strict upper limit of less than 100 ms for the time interval after the merging until the merger remnant has collapsed to a BH, leaving an accretion torus with a small mass of only ∼0.01 M ⊙ . This favours a (nearly) symmetric NS+NS binary with a typical mass as progenitor system.

Relativistic neutron star merger simulations with non-zero temperature equations of state I. Variation of binary parameters and equation of state

Abstract: An extended set of binary neutron star (NS) merger simulations is performed with an approximative conformally flat treatment of general relativity to systematically investigate the influence of the nuclear equation of state (EoS), the neutron star masses, and the NS spin states prior to merging. We employ the two non-zero temperature EoSs of Shen et al. (1998a,b) and Lattimer & Swesty (1991). In addition, we use the cold EoS of Akmal et al. (1998) with a simple ideal-gas-like extension according to Shibata & Taniguchi (2006), and an ideal-gas EoS with parameters fitted to the supernuclear part of the Shen-EoS. We estimate the mass sitting in a dilute high-angular momentum ``torus'' around the future black hole (BH). The dynamics and outcome of the models is found to depend strongly on the EoS and on the binary parameters. Larger torus masses are found for asymmetric systems (up to ~0.3 M_sun for a mass ratio of 0.55), for large initial NSs, and for a NS spin state which corresponds to a larger total angular momentum. We find that the postmerger remnant collapses either immediately or after a short time when employing the soft EoS of Lattimer& Swesty, whereas no sign of post-merging collapse is found within tens of dynamical timescales for all other EoSs used. The typical temperatures in the torus are found to be about 3-10 MeV depending on the strength of the shear motion at the collision interface between the NSs and thus depending on the initial NS spins. About 10^{-3}-10^{-2} M_sun of NS matter become gravitationally unbound during or right after the merging process. This matter consists of a hot/high-entropy component from the collision interface and (only in case of asymmetric systems) of a cool/low-entropy component from the spiral arm tips. (abridged)

Rotating Collapse of Stellar Iron Cores in General Relativity

Abstract: We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in two of our models that leads to prolonged narrow-band GW emission at ~ 920-930 Hz over several tens of milliseconds.

Neutrino-driven wind and wind termination shock in supernova cores

Abstract: The neutrino-driven wind from a nascent neutron star at the center of a supernova expands into the earlier ejecta of the explosion. Upon collision with this slower matter the wind material is decelerated in a wind termination shock. By means of hydrodynamic simulations in spherical symmetry we demonstrate that this can lead to a large increase of the wind entropy, density, and temperature, and to a strong deceleration of the wind expansion. The consequences of this phenomenon for the possible r-process nucleosynthesis in the late wind still need to be explored in detail. Two-dimensional models show that the wind-ejecta collision is highly anisotropic and could lead to a directional dependence of the nucleosynthesis even if the neutrino-driven wind itself is spherically symmetric.