We report here on recent progress in understanding the birth conditions of neutron stars and the way how supernovae explode. More sophisticated numerical models have led to the discovery of new phenomena in the supernova core, for example a generic hydrodynamic instability of the stagnant supernova shock against low-mode nonradial deformation and the excitation of gravity-wave activity in the surface and core of the nascent neutron star. Both can have supportive or decisive influence on the inauguration of the explosion, the former by improving the conditions for energy deposition by neutrino heating in the postshock gas, the latter by supplying the developing blast with a flux of acoustic power that adds to the energy transfer by neutrinos. While recent two-dimensional models suggest that the neutrino-driven mechanism may be viable for stars from about 8 solar masses to at least 15 solar masses, acoustic energy input has been advocated as an alternative if neutrino heating fails. Magnetohydrodynamic effects constitute another way to trigger explosions in connection with the collapse of sufficiently rapidly rotating stellar cores, perhaps linked to the birth of magnetars. The global explosion asymmetries seen in the recent simulations offer an explanation of even the highest measured kick velocities of young neutron stars.

By 2D hydrodynamic simulations including a detailed equation of state and neutrino transport, we investigate the interplay between different non-radial hydrodynamic instabilities that play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, which is driven by negative entropy gradients caused by neutrino heating or by time variations of the shock strength, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in 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 such a situation the shocked flow nevertheless can develop non-radial asymmetry with an oscillatory growth of the amplitude. This phenomenon has been termed ``standing 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 of the average shock radius and of the mass in the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter 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 NS turns out to be favorable for explosions, because the accretion luminosity and neutrino heating are larger and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI and a variety of other features in our simulations agree with estimates 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. (abridged)

We investigate the behavior and consequences of the reverse shock that terminates the supersonic expansion of the baryonic wind which is driven by neutrino heating off the surface of (non-magnetized) new-born neutron stars in supernova cores. To this end we perform long-time hydrodynamic simulations in spherical symmetry. In agreement with previous relativistic wind studies, we find that the neutrino-driven outflow accelerates to supersonic velocities and in case of a compact, about 1.4 solar mass (gravitational mass) neutron star with a radius of about 10 km, the wind reaches entropies of about 100 k_B per nucleon. The wind, however, is strongly influenced by the environment of the supernova core. It is decelerated and shock-heated abruptly by a termination shock that forms when the supersonic outflow collides with the slower preceding supernova ejecta. The radial position of this reverse shock varies with time and depends on the strength of the neutrino wind and the different conditions in progenitor stars with different masses and structure. Its basic properties and behavior can be understood by simple analytic considerations. We demonstrate that the entropy of matter going through the reverse shock can increase to a multiple of the asymptotic wind value. Seconds after the onset of the explosion it therefore can exceed 400 k_B per nucleon. The temperature of the shocked wind has typically dropped to about or less than 10^9 K, and density and temperature in the shock-decelerated matter continue to decrease only very slowly. Such conditions might strongly affect the important phases of supernova nucleosynthesis in a time and progenitor dependent way. (abridged)

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.

We analyze the linear stability of a stalled accretion shock in a perfect gas with a parametrized cooling function L ~ rho^{beta-alpha} P^alpha. The instability is dominated by the l=1 mode if the shock radius exceeds 2-3 times the accretor radius, depending on the parameters of the cooling function. The growth rate and oscillation period are comparable to those observed in the numerical simulations of Blondin & Mezzacappa (2006). The instability mechanism is analyzed by separately measuring the efficiencies of the purely acoustic cycle and the advective-acoustic cycle. These efficiencies are estimated directly from the eigenspectrum, and also through a WKB analysis in the high frequency limit. Both methods prove that the advective-acoustic cycle is unstable, and that the purely acoustic cycle is stable. Extrapolating these results to low frequency leads us to interpret the dominant mode as an advective-acoustic instability, different from the purely acoustic interpretation of Blondin & Mezzacappa (2006). A simplified characterization of the instability is proposed, based on an advective-acoustic cycle between the shock and the radius r_nabla where the velocity gradients of the stationary flow are strongest. The importance of the coupling region in this mechanism calls for a better understanding of the conditions for an efficient advective-acoustic coupling in a decelerated, nonadiabatic flow, in order to extend these results to core-collapse supernovae.

A toy model is analyzed in order to evaluate the linear stability of the gain region immediately behind a stalled accretion shock, after core bounce. This model demonstrates that a negative entropy gradient is not sufficient to warrant linear instability. The stability criterion is governed by the ratio \chi of the advection time through the gain region divided by the local timescale of buoyancy. The gain region is linearly stable if \chi< 3. The classical convective instability is recovered in the limit \chi\gg 3. For \chi>3, perturbations are unstable in a limited range of horizontal wavelengths centered around twice the vertical size H of the gain region. The threshold horizontal wavenumbers k_{min} and k_{max} follow simple scaling laws such that Hk_{min}\propto 1/{\chi} and Hk_{max}\propto \chi. The convective stability of the l=1 mode in spherical accretion is discussed, in relation with the asymmetric explosion of core collapse supernovae. The advective stabilization of long wavelength perturbations weakens the possible influence of convection alone on a global l=1 mode.

By means of two-dimensional (2D) simulations we study hydrodynamic instabilities during the first seconds of neutrino-driven supernova explosions, using a PPM hydrodynamics code, supplemented with a gray, non-equilibrium approximation of radial neutrino transport. We consider three 15 solar mass progenitors with different structures and one rotating model, in which we replace the dense core of the newly formed neutron star (NS) by a contracting inner grid boundary, and trigger neutrino-driven explosions by systematically varying the neutrino fluxes emitted at this boundary. Confirming more idealized 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, leading to a dominance of dipole (l=1) and quadrupole (l=2) modes in the explosion ejecta. Anisotropic gravitational and hydrodynamic forces are found to accelerate the NS on a timescale of 2-3 seconds. Since the explosion anisotropies develop chaotically, the magnitude of the corresponding kick varies stochastically in response to small differences in the fluid flow. Our more than 70 models separate into two groups, one with high and the other with low NS velocities and accelerations after 1s 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 NS velocities are extrapolated to their terminal values. The fast group has an average velocity of about 500 km/s and peak values in excess of 1000 km/s. Establishing a link to the measured distribution of pulsar velocities, however, requires a much larger set of calculations and ultimately 3D modeling. (abridged)

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass 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 1987A, 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 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**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 latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A 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. (abbreviated)

New simulations demonstrate that low-mode, nonradial hydrodynamic instabilities of the accretion shock help starting hot-bubble convection in supernovae and thus support explosions by the neutrino-heating mechanism. The prevailing conditions depend on the high-density equation of state which governs stellar core collapse, core bounce, and neutron star formation. Tests of this sensitivity to nuclear physics variations are shown for spherically symmetric models. Implications of current explosion models for r-process nucleosynthesis are addressed.

A few seconds after bounce in a core-collapse supernova, the shock wave passes the density region corresponding to resonant neutrino oscillations with the ``atmospheric'' neutrino mass difference. The transient violation of the adiabaticity condition manifests itself in an observable modulation of the neutrino signal from a future galactic supernova. In addition to the shock wave propagation effects that were previously studied, a reverse shock forms when the supersonically expanding neutrino-driven wind collides with the slower earlier supernova ejecta. This implies that for some period the neutrinos pass two subsequent density discontinuities, giving rise to a ``double dip'' feature in the average neutrino energy as a function of time. We study this effect both analytically and numerically and find that it allows one to trace the positions of the forward and reverse shocks. We show that the energy dependent neutrino conversion probabilities allow one to detect oscillations even if the energy spectra of different neutrino flavors are the same as long as the fluxes differ. These features are observable in the \bar\nu_e signal for an inverted and in the \nu_e signal for a normal neutrino mass hierarchy, provided the 13-mixing angle is ``large'' (sin^2\theta_{13}\gg 10^{-5}).

Two- and three-dimensional simulations demonstrate that hydrodynamic instabilities can lead to low-mode (l=1,2) asymmetries of the fluid flow in the neutrino-heated layer behind the supernova shock. This provides a natural explanation for aspherical mass ejection and for pulsar recoil velocities even in excess of 1000 km/s. We propose that the bimodality of the pulsar velocity distribution might be a consequence of a dominant l=1 mode in case of the fast component, while higher-mode anisotropy characterizes the postshock flow and SN ejecta during the birth of the slow neutron stars. We argue that the observed large asymmetries of supernovae and the measured high velocities of young pulsars therefore do not imply rapid rotation of the iron core of the progenitor star, nor do they require strong magnetic fields to play a crucial role in the explosion. Anisotropic neutrino emission from accretion contributes to the neutron star acceleration on a minor level, and pulsar kicks do not make a good case for non-standard neutrino physics in the nascent neutron star.

We show by two-dimensional and first three-dimensional simulations of neutrino-driven supernova explosions that low (l=1,2) modes can dominate the flow pattern in the convective postshock region on timescales of hundreds of milliseconds after core bounce. This can lead to large global anisotropy of the supernova explosion and pulsar kicks in excess of 500 km/s.

Assuming that the neutrino luminosity from the neutron star core is sufficiently high to drive supernova explosions by the neutrino-heating mechanism, we show that low-mode (l = 1, 2) convection can develop from random seed perturbations behind the shock. A slow onset of the explosion is crucial, requiring the core luminosity to vary slowly with time, in contrast to the burst-like exponential decay assumed in previous work. Gravitational and hydrodynamic forces by the globally asymmetric supernova ejecta were found to accelerate the remnant neutron star on a timescale of more than a second to velocities above 500 km/s, in agreement with observed pulsar proper motions.

We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogeneous environment. The equation of state used in the simulations accounts for an arbitrary mixture of electrons, protons and electron-positron pairs. Using the hydrodynamic models we have also computed synthetic radio maps and the thermal Bremsstrahlung X-ray emission from their cavities. Although there is a difference of about three orders of magnitude in the temperatures of the cavities inflated by the simulated jets, we find that both the morphology and the dynamic behaviour are almost independent on the assumed composition of the jets. Their evolution proceeds in two distinct epochs. During the first one multidimensional effects are unimportant and the jets propagate ballistically. The second epoch starts when the first larger vortices are produced near the jet head causing the beam cross section to increase and the jet to decelerate. The evolution of the cocoon and cavity is in agreement with a simple theoretical model. The beam velocities are relativistic ($\Gamma \simeq 4$) at kiloparsec scales supporting the idea that the X-ray emission of several extragalactic jets may be due to relativistically boosted CMB photons. The radio emission of all models is dominated by the contribution of the hot spots. All models exhibit a depression in the X-rays surface brightness of the cavity interior in agreement with recent observations.