We present results from the first generation of multi-dimensional general relativistic neutrino hydrodynamics simulations of core-collapse supernovae. A comparison with models computed using either the purely Newtonian approximation or the "effective gravitational potential" approach reveals appreciable quantitative differences in the heating conditions and the gravitational wave spectra. Our results underscore the important role of general relativity in the supernova problem (which appears to be on par with other important factors such as the dimensionality and the equation of state) both for our understanding of the explosion dynamics as well as for predictions of observable signatures.

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-resolution 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)

Presented in this paper is a technique that we propose for extracting the physical parameters of a rotating stellar core collapse from the observation of the associated gravitational wave signal from the collapse and core bounce. Data from interferometric gravitational wave detectors can be used to provide information on the mass of the progenitor model, precollapse rotation and the nuclear equation of state. We use waveform libraries provided by the latest numerical simulations of rotating stellar core collapse models in general relativity, and from them create an orthogonal set of eigenvectors using principal component analysis. Bayesian inference techniques are then used to reconstruct the associated gravitational wave signal that is assumed to be detected by an interferometric detector. Posterior probability distribution functions are derived for the amplitudes of the principal component analysis eigenvectors, and the pulse arrival time. We show how the reconstructed signal and the principal component analysis eigenvector amplitude estimates may provide information on the physical parameters associated with the core collapse event.

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

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.

Using numerical simulations based on solving the general relativistic hydrodynamic equations, we study the dynamics of a phase transition in the dense core of isolated rotating neutron stars, triggered by the back bending instability reached via angular momentum loss. In particular, we investigate the dynamics of a migration from an unstable configuration into a stable one, which leads to a mini-collapse of the neutron star and excites sizeable pulsations in its bulk until it acquires a new stable equilibrium state. We consider equations of state with softening at high densities, a simple analytic one with a mixed hadron-quark phase in an intermediate pressure interval and pure quark matter at very high densities, and a microphysical one that has a first-order phase transition, originating from kaon condensation. Although the marginally stable initial models are rigidly rotating, we observe that during the collapse (albeit little) differential rotation is created. We analyze the emission of gravitational radiation, which in some models is amplified by mode resonance effects, and assess its prospective detectability by interferometric detectors. We expect that the most favorable conditions for dynamic migration exist in very young magnetars. We find that the damping of the post-migration pulsations strongly depends on the character of the equation of state softening. The damping of pulsations in the models with the microphysical equation of state is caused by dissipation associated with matter flowing through the density jump at the edge of the dense core. If at work, this mechanism dominates over all other types of dissipation, like bulk viscosity in the exotic-phase core, gravitational radiation damping, or numerical viscosity.

We present the generalization of a recently introduced modified gravitational potential for self-gravitating fluids. The use of this potential allows for an accurate approximation of general relativistic effects in an otherwise Newtonian hydrodynamics code also in cases of rapid rotation. We test this approach in numerical simulations of astrophysical scenarios related to compact stars, like supernova core collapse with both a simplified and detailed microphysical description of matter, and rotating neutron stars in equilibrium. We assess the quality of the new potential, and demonstrate that it provides a significant improvement compared to previous formulations for such potentials. Newtonian simulations of compact objects employing such an effective relativistic potential predict inaccurate pulsation frequencies despite the excellent agreement of the collapse dynamics and structure of the compact objects with general relativistic results. We analyze and discuss the reason for this behavior.

We present results from detailed general relativistic simulations of stellar core collapse to a proto-neutron star, using two different microphysical nonzero-temperature nuclear equations of state as well as an approximate description of deleptonization during the collapse phase. Investigating a wide variety of rotation rates and profiles as well as masses of the progenitor stars and both equations of state, we confirm in this very general setup the recent finding that a generic gravitational wave burst signal is associated with core bounce, already known as type I in the literature. The previously suggested type II (or "multiple-bounce") waveform morphology does not occur. Despite this reduction to a single waveform type, we demonstrate that it is still possible to constrain the progenitor and postbounce rotation based on a combination of the maximum signal amplitude and the peak frequency of the emitted gravitational wave burst. Our models include to sufficient accuracy the currently known necessary physics for the collapse and bounce phase of core-collapse supernovae, yielding accurate and reliable gravitational wave signal templates for gravitational wave data analysis. In addition, we assess the possiblity of nonaxisymmetric instabilities in rotating nascent proto-neutron stars. We find strong evidence that in an iron core-collapse event the postbounce core cannot reach sufficiently rapid rotation to become subject to a classical bar-mode instability. However, many of our postbounce core models exhibit sufficiently rapid and differential rotation to become subject to the recently discovered dynamical instability at low rotation rates.

This paper presents results from axisymmetric simulations of magneto-rotational stellar core collapse to neutron stars in general relativity using the passive field approximation for the magnetic field. These simulations are performed using a new general relativistic numerical code specifically designed to study this astrophysical scenario. The code is based on the conformally-flat approximation of Einstein's field equations and conservative formulations of the magneto-hydrodynamics equations. The code has been recently upgraded to incorporate a tabulated, microphysical equation of state and an approximate deleptonization scheme. This allows us to perform the most realistic simulations of magneto-rotational core collapse to date, which are compared with simulations employing a simplified (hybrid) equation of state, widely used in the relativistic core collapse community. Furthermore, state-of-the-art (unmagnetized) initial models from stellar evolution are used. In general, stellar evolution models predict weak magnetic fields in the progenitors, which justifies our simplification of performing the computations under the approach that we call the passive field approximation for the magnetic field. Our results show that for the core collapse models with microphysics the saturation of the magnetic field cannot be reached within dynamical time scales by winding up the poloidal magnetic field into a toroidal one. We estimate the effect of other amplification mechanisms including the magneto-rotational instability (MRI) and several types of dynamos.

We perform general relativistic simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state as well as an approximate description of deleptonization. We show that for a wide variety of rotation rates and profiles the gravitational wave burst signals from the core bounce are of a generic type, known as Type I in the literature. In our systematic study, using both general relativity and Newtonian gravity, we identify and individually quantify the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the equation of state, and deleptonization. Such a generic type of signal templates will likely facilitate a more efficient search in current and future gravitational wave detectors of both interferometric and resonant type.

We present detailed results from performing general relativistic (GR) simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state (EoS) as well as an approximate description of deleptonization during the collapse phase. We show that for a wide variety of rotation rates and profiles the gravitational wave (GW) burst signals from the core bounce are of a generic type, already known as Type I in the literature. In addition, for most models the characteristic frequency of the GW burst signal lies in a narrow range around approximately 718 Hz. In our systematic study, using both GR and Newtonian gravity, we identify, individually quantify, and discuss in detail the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the EoS, and deleptonization. We also discuss the detectability prospects of such GW burst signals by GW detectors, and infer that such a generic type of signal templates will likely facilitate a more efficient search in current and future detectors of both interferometric and resonant type.

We present results from the first 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in one of our models that leads to prolonged narrow-band GW emission at ~930 Hz over several tens of milliseconds.

We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in two of our models that leads to prolonged narrow-band GW emission at ~ 920-930 Hz over several tens of milliseconds.

We continue our investigations of the magnetorotational collapse of stellar cores discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV potential used in our simulations catches several features of fully relativistic simulations quite well. It is able to correctly reproduce the behavior of models which show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. If this is not the case, the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential allowing for a more efficient amplification of the magnetic field by differential rotation. Sufficiently strong magnetic fields brake down the core's rotation and trigger a contraction phase to higher densities. Several models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce does not suffice to eject any mass, and the core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Instead it belongs to the signal type IV introduced by us in the first paper of this series. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and > 100 kHz denotes the generation of a jet-like outflow. [Abstract abbreviated]

We study non-linear axisymmetric pulsations of rotating relativistic stars using a general relativistic hydrodynamics code under the assumption of a conformal flatness. We compare our results to previous simulations where the spacetime dynamics was neglected. The pulsations are studied along various sequences of both uniformly and differentially rotating relativistic polytropes with index N = 1. We identify several modes, including the lowest-order l = 0, 2, and 4 axisymmetric modes, as well as several axisymmetric inertial modes. Differential rotation significantly lowers mode frequencies, increasing prospects for detection by current gravitational wave interferometers. We observe an extended avoided crossing between the l = 0 and l = 4 first overtones, which is important for correctly identifying mode frequencies in case of detection. For uniformly rotating stars near the mass-shedding limit, we confirm the existence of the mass-shedding-induced damping of pulsations, though the effect is not as strong as in the Cowling approximation. We also investigate non-linear harmonics of the linear modes and notice that rotation changes the pulsation frequencies in a way that would allow for various parametric instabilities between two or three modes to take place. We assess the detectability of each obtained mode by current gravitational wave detectors and outline how the empirical relations that have been constructed for gravitational wave asteroseismology could be extended to include the effects of rotation.

We review various approaches to approximating general relativistic effects in hydrodynamic simulations of stellar core collapse and post-bounce evolution. Different formulations of a modified Newtonian gravitational potential are presented. Such an effective relativistic potential can be used in an otherwise standard Newtonian hydrodynamic code. An alternative approximation of general relativity is the assumption of conformal flatness for the three-metric, and its extension by adding second post-Newtonian order terms. Using a code which evolves the coupled system of metric and fluid equations, we apply the various approximation methods to numerically simulate axisymmetric models for the collapse of rotating massive stellar cores. We compare the collapse dynamics and gravitational wave signals (which are extracted using the quadrupole formula), and thereby assess the quality of the individual approximation method. It is shown that while the use of an effective relativistic potential already poses a significant improvement compared to a genuinely Newtonian approach, the two conformal-flatness-based approximation methods yield even more accurate results, which are qualitatively and quantitatively very close to those of a fully general relativistic code even for rotating models which almost collapse to a black hole.

We investigate the possibility to approximate relativistic effects in hydrodynamical simulations of stellar core collapse and post-bounce evolution by using a modified gravitational potential in an otherwise standard Newtonian hydrodynamic code. Different modifications of the effective relativistic potential introduced by Rampp & Janka (2002) are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hydrodynamic simulations with two different codes are compared with fully relativistic results. One code is applied for one- and two-dimensional calculations with a simple equation of state, and employs either the modified effective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale supernova runs including a microphysical equation of state and neutrino transport based on the solution of the Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully relativistic solutions in spherical symmetry, leading to significant improvements compared to previously published approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation.

We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We demonstrate that the combination of a finite difference grid and a spectral grid can be successfully accomplished. This "Mariage des Maillages" (French for grid wedding) approach results in high accuracy of the metric solver and allows for fully 3D applications using computationally affordable resources, and ensures long term numerical stability of the evolution. We compare our new approach to two other, finite difference based, methods to solve the metric equations. A variety of tests in 2D and 3D is presented, involving highly perturbed neutron star spacetimes and (axisymmetric) stellar core collapse, demonstrating the ability to handle spacetimes with and without symmetries in strong gravity. These tests are also employed to assess gravitational waveform extraction, which is based on the quadrupole formula.

Core collapse supernovae are a promising source of detectable gravitational waves. Most of the existing (multidimensional) numerical simulations of core collapse in general relativity have been done using approximations of the Einstein field equations. As recently shown by Dimmelmeier et al (2002a,b), one of the most interesting such approximation is the so-called conformal flatness condition (CFC) of Isenberg, Wilson and Mathews. Building on this previous work we present here new results from numerical simulations of relativistic rotational core collapse in axisymmetry, aiming at improving the dynamics and the gravitational waveforms. The computer code used for these simulations evolves the coupled system of metric and fluid equations using the 3+1 formalism, specialized to a new framework for the gravitational field equations which we call CFC+. In this approach we add new degrees of freedom to the original CFC equations, which extend them by terms of second post-Newtonian order. The corrections for CFC+ are computed solving a system of elliptic linear equations. The new formalism is assessed with time evolutions of both rotating neutron stars in equilibrium and gravitational core collapse of rotating polytropes. Gravitational wave signals for a comprehensive sample of collapse models are extracted using either the quadrupole formula or directly from the metric. We discuss our results on the dynamics and the gravitational wave emission through a detailed comparison between CFC and CFC+ simulations. The main conclusion is that, for the neutron star spacetimes analyzed in the present work, no significant differences are found among CFC, CFC+, and full general relativity, which highlights the suitability of the former.

We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in axisymmetry and have computed the gravitational radiation emitted by such an event. 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. 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 and frequencies as the corresponding Newtonian ones. 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. 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.

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. 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 which 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.

Axisymmetric pulsations of rotating neutron stars can be excited in several scenarios, such as core-collapse, crust and core-quakes and binary mergers and could become detectable either in gravitational waves or high-energy radiation. Here, we present a comprehensive study of all low-order axisymmetric modes of uniformly and rapidly rotating relativistic stars. Initial stationary configurations are appropriately perturbed and are numerically evolved using an axisymmetric, nonlinear relativistic hydrodynamics code, assuming time-independence of the gravitational field (Cowling approximation). The simulations are performed using a high-resolution shock-capturing finite-difference scheme accurate enough to maintain the initial rotation law for a large number of rotational periods, even for stars at the mass-shedding limit. Through Fourier transforms of the time evolution of selected fluid variables, we compute the frequencies of quasi-radial and non-radial modes with spherical harmonic indices l=0,1,2 and 3, for a sequence of rotating stars from the non-rotating limit to the mass-shedding limit. The frequencies of the axisymmetric modes are affected significantly by rotation only when the rotation rate exceeds about 50% of the maximum allowed. As expected, at large rotation rates, apparent mode crossings between different modes appear. In addition to the above modes, several axisymmetric inertial modes are also excited in our numerical evolutions.

We present results from simulations of axisymmetric relativistic rotational core collapse. The general relativistic hydrodynamic equations are formulated in flux-conservative form and solved using a high-resolution shock-capturing scheme. The Einstein equations are approximated with a conformally flat 3-metric. We use the quadrupole formula to extract waveforms of the gravitational radiation emitted during the collapse. A comparison of our results with those of Newtonian simulations shows that the wave amplitudes agree within 30%. Surprisingly, in some cases, relativistic effects actually diminish the amplitude of the gravitational wave signal. We further find that the parameter range of models suffering multiple coherent bounces due to centrifugal forces is considerably smaller than in Newtonian simulations.