We present a new general relativistic (GR) code for hydrodynamic supernova simulations with neutrino transport in spherical and azimuthal symmetry (1D/2D). The code is a combination of the CoCoNuT hydro module, which is a Riemann-solver based, high-r
esolution shock-capturing method, and the three-flavor, energy-dependent neutrino transport scheme VERTEX. VERTEX integrates the neutrino moment equations with a variable Eddington factor closure computed from a model Boltzmann equation and uses the ray-by-ray plus approximation in 2D, assuming the neutrino distribution to be axially symmetric around the radial direction, and thus the neutrino flux to be radial. Our spacetime treatment employs the ADM 3+1 formalism with the conformal flatness condition for the spatial three-metric. This approach is exact in 1D and has been shown to yield very accurate results also for rotational stellar collapse. We introduce new formulations of the energy equation to improve total energy conservation in relativistic and Newtonian hydro simulations with Eulerian finite-volume codes. Moreover, a modified version of the VERTEX scheme is developed that simultaneously conserves energy and lepton number with better accuracy and higher numerical stability. To verify our code, we conduct a series of tests, including a detailed comparison with published 1D results for stellar core collapse. Long-time simulations of proto-neutron star cooling over several seconds both demonstrate the robustness of the new CoCoNuT-VERTEX code and show the approximate treatment of GR effects by means of an effective gravitational potential as in PROMETHEUS-VERTEX to be remarkably accurate in 1D. (abridged)
An increase in the central density of a neutron star may trigger a phase transition from hadronic matter to deconfined quark matter in the core, causing it to collapse to a more compact hybrid-star configuration. We present a study of this, building
on previous work by Lin et al. (2006). We follow them in considering a supersonic phase transition and using a simplified equation of state, but our calculations are general relativistic (using 2D simulations in the conformally flat approximation) as compared with their 3D Newtonian treatment. We also improved the treatment of the initial phase transformation, avoiding the introduction of artificial convection. As before, we find that the emitted gravitational-wave spectrum is dominated by the fundamental quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much smaller than suggested previously, which is disappointing for the detection prospects. However, we see significantly smaller damping and observe a nonlinear mode resonance which substantially enhances the emission in some cases. We explain the damping mechanisms operating, giving a different view from the previous work. Finally, we discuss the detectability of the gravitational waves, showing that the signal-to-noise ratio for current or second generation interferometers could be high enough to detect such events in our Galaxy, although third generation detectors would be needed to observe them out to the Virgo cluster, which would be necessary for having a reasonable event rate.
