We introduce novel high order well-balanced finite volume methods for the full compressible Euler system with gravity source term. They require no \`a priori knowledge of the hydrostatic solution which is to be well-balanced and are not restricted to certain classes of hydrostatic solutions. In one spatial dimension and under the assumption of an ideal gas law we construct a method that exactly balances a high order discretization of any hydrostatic state. The method is extended to general equations of state and two spatial dimensions using a local high order approximation of any hydrostatic state in each cell. The proposed methods are simple, flexible, and robust. Numerical tests verify that the well-balanced method improves the capability to accurately resolve small perturbations on hydrostatic states.

We present an Advanced Spectral Leakage (ASL) scheme to model neutrinos in the context of core-collapse supernovae and compact binary mergers. Based on previous gray leakage schemes, the ASL scheme computes the neutrino cooling rates by interpolating local production and diffusion rates (relevant in optically thin and thick regimes, respectively), separately for discretized values of the neutrino energy. Neutrino trapped components are also modeled, based on equilibrium and timescale arguments. The better accuracy achieved by the spectral treatment allows a more reliable computation of neutrino heating rates in optically thin conditions. The scheme has been calibrated and tested against Boltzmann transport in the context of Newtonian spherically symmetric models of core-collapse supernovae. ASL shows a very good qualitative and a partial quantitative agreement, for key quantities from collapse to a few hundreds of milliseconds after core bounce. We have proved the adaptability and flexibility of our ASL scheme coupling it to an axisymmetric Eulerian and to a three-dimensional SPH code to simulate core-collapse. Therefore, the neutrino treatment presented here is ideal for large parameter-space explorations, parametric studies, high-resolution tests, code developments, and long-term modeling of asymmetric configurations, where more detailed neutrino treatments are not available or currently computationally too expensive.

We explore heavy element nucleosynthesis in the explosion of massive stars which are triggered by a quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. The present study is based on general relativistic radiation hydrodynamics simulations with three-flavor Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation of state based on the MIT bag model for strange quark matter. The quark-hadron phase transition inside the stellar core forms a shock wave propagating towards the surface of the proto-neutron star. The shock wave results in an explosion and ejects neutron-rich matter which is piled up or accreting on the proto-neutron star. Later, during the cooling phase, the proto-neutron star develops a proton-rich neutrino-driven wind. We present a detailed analysis of the nucleosynthesis outcome in both neutron-rich and proton-rich ejecta and compare our integrated nucleosynthesis with observations of metal poor stars.

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.