An erroneous interpretation of the hydrodynamical results led to an incorrect determination of the fallback masses in Ugliano et al. (2012), which also (on a smaller level) affects the neutron star masses provided in that paper. This problem was already addressed and corrected in the follow-up works by Ertl et al. (2015) and Sukhbold et al. (2015). Therefore, the reader is advised to use the new data of the latter two publications. In the remaining text of this Erratum we present the differences of the old and new fallback results in detail and explain the origin of the mistake in the original analysis by Ugliano et al. (2012).

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.

We perform hydrodynamic supernova simulations in spherical symmetry for over 100 single stars of solar metallicity to explore the progenitor-explosion and progenitor-remnant connections established by the neutrino-driven mechanism. We use an approximative treatment of neutrino transport and replace the high-density interior of the neutron star (NS) by an inner boundary condition based on an analytic proto-NS core-cooling model, whose free parameters are chosen such that explosion energy, nickel production, and energy release by the compact remnant of progenitors around 20 solar masses are compatible with Supernova 1987A. Thus we are able to simulate the accretion phase, initiation of the explosion, subsequent neutrino-driven wind phase for 15-20 s, and the further evolution of the blast wave for hours to days until fallback is completed. Our results challenge long-standing paradigms. We find that remnant mass, launch time, and properties of the explosion depend strongly on the stellar structure and exhibit large variability even in narrow intervals of the progenitors' zero-age-main-sequence mass. While all progenitors with masses below about 15 solar masses yield NSs, black hole (BH) as well as NS formation is possible for more massive stars, where partial loss of the hydrogen envelope leads to weak reverse shocks and weak fallback. Our NS baryonic masses of ~1.2-2.0 solar masses and BH masses >6 solar masses are compatible with a possible lack of low-mass BHs in the empirical distribution. Neutrino heating accounts for SN energies between some 10^{50} erg and about 2*10^{51} erg, but can hardly explain more energetic explosions and nickel masses higher than 0.1-0.2 solar masses. These seem to require an alternative SN mechanism.

We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for the equations of ideal MHD. The solution to the Riemann problem is approximated by a five wave pattern, comprised of two outermost fast shocks, two rotational discontinuities and a contact surface in the middle. The proposed scheme is considerably more elaborate than in the classical case since the normal velocity is no longer constant across the rotational modes. Still, proper closure to the Rankine-Hugoniot jump conditions can be attained by solving a nonlinear scalar equation in the total pressure variable which, for the chosen configuration, has to be constant over the whole Riemann fan. The accuracy of the new Riemann solver is validated against one dimensional tests and multidimensional applications. It is shown that our new solver considerably improves over the popular HLL solver or the recently proposed HLLC schemes.