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.

Abstract.We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in axisymmetry and have computed the gravitational radiation emitted by such an event. The Einstein equations are formulated using the conformally flat metric approximation, and the corresponding hydrodynamic equations are written as a first-order flux-conservative hyperbolic system. Details of the methodology and of the numerical code have been given in an accompanying paper. We have simulated the evolution of 26 models in both Newtonian and relativistic gravity. The initial configurations are differentially rotating relativistic 4/3-polytropes in equilibrium which have a central density of 10 10 g cm −3 . Collapse is initiated by decreasing the adiabatic index to some prescribed fixed value. The equation of state consists of a polytropic and a thermal part for a more realistic treatment of shock waves. Any microphysics like electron capture and neutrino transport is neglected. Our simulations show that the three different types of rotational supernova core collapse and gravitational waveforms identified in previous Newtonian simulations (regular collapse, multiple bounce collapse, and rapid collapse) are also present in relativistic gravity. However, rotational core collapse with multiple bounces is only possible in a much narrower parameter range in relativistic gravity. The relativistic models cover almost the same range of gravitational wave amplitudes (4 × 10 −21 ≤ h TT ≤ 3 × 10 −20 for a source at a distance of 10 kpc) and frequencies (60 Hz ≤ ν ≤ 1000 Hz) as the corresponding Newtonian ones. Averaged over all models, the total energy radiated in the form of gravitational waves is 8.2 × 10 −8 M c 2 in the relativistic case, and 3.6 × 10 −8 M c 2 in the Newtonian case. For all collapse models that are of the same type in both Newtonian and relativistic gravity, the gravitational wave signal is of lower amplitude. If the collapse type changes, either weaker or stronger signals are found in the relativistic case. For a given model, relativistic gravity can cause a large increase of the characteristic signal frequency of up to a factor of five, which may have important consequences for the signal detection. Our study implies that the prospects for detection of gravitational wave signals from axisymmetric supernova rotational core collapse do not improve when taking into account relativistic gravity. The gravitational wave signals obtained in our study are within the sensitivity range of the first generation laser interferometer detectors if the source is located within the Local Group. An online catalogue containing the gravitational wave signal amplitudes and spectra of all our models is available at the URL

Abstract. We describe an axisymmetric general relativistic code for rotational core collapse. The code evolves the coupled system of metric and fluid equations using the ADM 3 + 1 formalism and a conformally flat metric approximation of the Einstein equations. Within this approximation the ADM 3 + 1 equations reduce to a set of five coupled non-linear elliptic equations for the metric components. The equations are discretized on a 2D grid in spherical polar coordinates and are solved by means of a Newton-Raphson iteration using a block elimination scheme to solve the diagonally dominant, sparse linear system arising within each iteration step. The relativistic hydrodynamics equations are formulated as a first-order flux-conservative hyperbolic system and are integrated using high-resolution shock-capturing schemes based on Riemann solvers. We assess the quality of the conformally flat metric approximation for relativistic core collapse and present a comprehensive set of tests that the code successfully passed. The tests include relativistic shock tubes, the preservation of the rotation profile and of the equilibrium of rapidly and differentially rotating neutron stars (approximated as rotating polytropes), spherical relativistic core collapse, and the conservation of rest-mass and angular momentum in dynamic spacetimes. The application of the code to relativistic rotational core collapse, with emphasis on the gravitational waveform signature, is presented in an accompanying paper.

We present three-dimensional (3D) simulations of supernova explosions of nonrotating stars, triggered by the delayed neutrinoheating mechanism with a suitable choice of the core-neutrino luminosity. Our results show that asymmetric mass ejection caused by hydrodynamic instabilities can accelerate the neutron star (NS) up to recoil velocities of more than 700 km s −1 by the "gravitational tug-boat mechanism", which is sufficient to explain most observed pulsar space velocities. The associated NS spin periods for our nonrotating progenitors are about 100 ms to 8000 ms without any obvious correlation between spin and kick magnitudes or directions. This suggests that faster spins and a possible spin-kick alignment might require angular momentum in the progenitor core prior to collapse. Our simulations for the first time demonstrate a clear correlation between the size of the NS kick and anisotropic production and distribution of heavy elements created by explosive burning behind the shock. In the case of large pulsar kicks, the explosion is significantly stronger opposite to the kick vector. Therefore the bulk of the explosively fused iron-group elements, in particular nickel, are ejected mostly in large clumps against the kick direction. This contrasts with the case of low recoil velocity, where the nickelrich lumps are more isotropically distributed. Explosively produced intermediate-mass nuclei heavier than 28 Si (like 40 Ca and 44 Ti) also exhibit significant enhancement in the hemisphere opposite to the direction of fast NS motion, while the distribution of 12 C, 16 O, and 20 Ne is not affected, and that of 24 Mg only marginally. Mapping the spatial distribution of the heavy elements in supernova remnants with identified pulsar motion may offer an important diagnostic test of the kick mechanism. Unlike kick scenarios based on anisotropic neutrino emission, our hydrodynamical acceleration model predicts enhanced ejection of iron-group elements and of their nuclear precursors in the opposite direction to the NS recoil.

Using a collapsar progenitor model of MacFadyen & Woosley, we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multidimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range of 1050-1051 ergs s-1 within a 30 degrees cone around the rotation axis. The jet flow is strongly beamed (approximately less than a few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor (R*=2.98x1010 cm) intact. At breakout, the maximum Lorentz factor of the jet flow is 33. After breakout, the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then become constant, rhoext=10-5 g cm-3. Outside the star, the flow begins to expand laterally also (v approximately c), but the beam remains very well collimated. At a distance of 2.54 R*, where the simulation ends, the Lorentz factor has increased to 44.