Gravitational Waves from Gravitational Collapse

Kimberly C. B. New

Los Alamos National Laboratory MS B220, X-2 Los Alamos, New Mexico 87545 U.S.A.

2003-03-10

Abstract
Gravitational wave emission from the gravitational collapse of massive stars has been studied for more than three decades. Current state-of-the-art numerical investigations of collapse include those that use progenitors with more realistic angular momentum profiles, properly treat microphysics issues, account for general relativity, and examine non-axisymmetric effects in three dimensions. Such simulations predict that gravitational waves from various phenomena associated with gravitational collapse could be detectable with ground-based and space-based interferometric observatories.

1 Introduction

The field of gravitational wave (GW) astronomy will soon become a reality. The first generation of ground-based interferometric detectors (LIGO [37] , VIRGO [99] , GEO 600 [203] , TAMA 300 [150] ) are beginning their search for GWs. Toward the end of this decade, two of these detectors (LIGO, VIRGO) will undergo upgrades that should allow them to reach the sensitivities necessary to regularly detect emissions from astrophysical sources. A space-based interferometric detector, LISA [108] , could be launched in the early part of the next decade. One important class of sources for these observatories is stellar gravitational collapse. This class includes the accretion induced collapse (AIC) of white dwarf binary components and the core collapse of massive stars ( M > 8 M ⊙   ), very massive Population III stars ( M = 100 – 500 M ⊙   ), and supermassive stars (SMSs, M > 10 6 M ⊙   ). Some of these collapses result in explosions (Type II, Ib/c supernovae and hypernovae) and all leave behind neutron star or black hole remnants.
Strong GWs can be emitted during a gravitational collapse/explosion and, following the collapse, by the resulting compact remnant [197, 142, 143, 64, 174, 71, 72, 94] . GW emission during the collapse itself may result if the collapse or explosion involves aspherical bulk mass motion or convection.
Rotational or fragmentation instabilities encountered by the collapsing star will also produce GWs. Neutron star remnants of collapse may emit GWs due to the growth of rotational or r-mode instabilities. Black hole remnants will also be sources of GWs if they experience accretion induced ringing. All of these phenomena have the potential of being detected by gravitational wave observatories because they involve the rapid change of dense matter distributions.
Observation of gravitational collapse by gravitational wave detectors will provide unique information, complementary to that derived from electromagnetic and neutrino detectors. Gravitational radiation arises from the coherent superposition of mass motion, whereas electromagnetic emission is produced by the incoherent superposition of radiation from electrons, atoms, and molecules. Thus, GWs carry different kinds of information than other types of radiation. Furthermore, electromagnetic radiation interacts strongly with matter and thus gives a view of the collapse only from lower density regions near the surface of the star, and it is weakened by absorption as it travels to the detector. Neutrinos can escape from much further within the collapsing star, but even they are scattered by the highest density regions in the core. By contrast, gravitational waves can propagate from the innermost parts of the stellar core to detectors without attenuation by intervening matter.
The characteristics of the GW emission from gravitational collapse have been the subject of much study. Core collapse supernovae, in particular, have been investigated as sources of gravitational radiation for more than three decades (see, e.g., [165, 198, 166, 50, 148, 139, 188, 65, 137, 217, 160, 71, 72] ).
However, during this time research has produced estimates of GW strength that vary over orders of magnitude. This is due to the complex nature of core collapse. Important theoretical and numerical issues include
  • ∙ construction of accurate progenitor models, including realistic angular momentum distributions,
  • ∙ proper treatment of microphysics, including the use of realistic equations of state and neutrino transport,
  • ∙ simulation in three-dimensions to study non-axisymmetric effects,
  • ∙ inclusion of general relativistic effects,
  • ∙ inclusion of magnetic field effects, and
  • ∙ study of the effect of an envelope on core behavior.
To date, collapse simulations generally include state-of-the-art treatments of only one or two of the above physics issues (often because of numerical constraints). For example, those studies that include advanced microphysics have often been run with Newtonian gravity (and approximate evaluation of the GW emission; see Section  2.4 ). A 3D, general relativistic collapse simulation that includes all significant physics effects is not feasible at present. However, good progress has been made on the majority of the issues listed above; the more recent work will be reviewed in some detail here.
The remainder of this article is structured as follows. Each category of gravitational collapse will be discussed in a separate section (AIC in Section  2 , collapse of massive stars in Section  3 , collapse of Population III stars in Section  4 , and collapse of SMSs in Section  5 ). Each of these sections ( 2 ,  3 ,  4 ,  5 ) is divided into subsection topics: collapse scenario, formation rate, GW emission mechanisms, and numerical predictions of GW emission. In the subsections on numerical predictions, the detectability of the GW emission from various phenomena associated with collapse is examined.
In particular, the predicted characteristics of GW emission are compared to the sensitivities of LIGO (for sources with frequencies of 1   to 10 4 H z   ) and LISA (for sources with lower frequencies in the range of 10 − 4   to 1 H z   ).

2 Accretion Induced Collapse

2.1 Collapse scenario

As a white dwarf accretes matter from a binary companion, its density will increase. If the rate of the resulting compressional heating exceeds the white dwarf 's cooling rate, the star's temperature will increase as well. When the mass of the white dwarf exceeds the Chandrasekhar stability limit, it will begin to collapse.
If the temperature of the core is high enough, nuclear burning will begin and the stellar pressure will increase. However, electron capture will act to reduce both the temperature and pressure behind the burn front. The collapse outcome is determined by the relative strength of the nuclear burning and electron capture. If nuclear burning is strong enough, the white dwarf will explode as a Type Ia supernova (SN). If electron capture wins out, the collapse will continue and result in neutron star formation. Electron capture will dominate if the density at which nuclear ignition occurs exceeds a critical density ρ c r i t   . For C-O white dwarfs, ρ c r i t   is in the range 6 × 10 9 – 10 10 g c m − 3   [26] . For O-Ne-Mg white dwarfs, electron capture may be stronger than nuclear burning under most conditions (if a Rayleigh-Taylor instability does not produce a turbulent burn front) [155, 98] . The collapse of an O-Ne-Mg white dwarf begins when its central density reaches 4 × 10 9 g c m − 3   . (For more details about the conditions under which AIC occurs, see [155, 98, 26, 25, 125] .) The dynamics of the collapse itself are somewhat similar to the dynamics of core collapse SNe (without the resulting explosive disruption of the star). The collapse is halted when nuclear densities are reached in the core. The core bounces and sends a bounce-shock outward through the star. After several milliseconds, the shock stalls and the remainder of the star collapses through the shock front (see, e.g., [211] ). Less than 10 − 1 M ⊙   will likely be ejected by the star (due to the bounce itself or due to neutrino absorption/wind mechanisms) [69] .

2.2 Formation rate

The AIC occurrence rate is difficult to determine for a number of reasons. These include incomplete understanding of binary star evolution and accretion processes [40, 111] . Another uncertainty is whether the collapse of an accreting Chandrasekhar mass white dwarf results in a supernovae explosion or a complete AIC (with accompanying neutron star formation).
The AIC rate can be indirectly inferred from the observed amount of rare, neutron rich isotopes present in the Galaxy. These isotopes (formed via electron capture) are present in the portion ( ∼ 0.1 M ⊙   ) of the outer envelope ejected by the star during an AIC. If all of these isotopes present in the Galaxy are assumed to have originated in AICs, an upper limit of ∼ 10 − 5 y r − 1   can be set for the galactic AIC rate [69] .
Binary population synthesis analysis can be used to determine which accreting white dwarfs will undergo AIC. The results of Yungelson and Livio [214] predict that the galactic AIC rate is between 8 × 10 − 7   and 8 × 10 − 5 y r − 1   . Thus, a reasonable occurrence rate can be found for an observation distance of 100 M p c   .

2.3 GW emission mechanisms

During AIC, emission of GWs will occur if the infall of matter is aspherical. GWs will also be produced if the collapsing star or neutron star remnant encounters non-axisymmetric instabilities.
These include global rotational mode, r   -mode, and fragmentation instabilities.
Global rotational instabilities in fluids arise from non-axisymmetric modes e ± i m φ   , where m = 2   is known as the “bar-mode” [192, 6] . It is convenient to parameterize a system's susceptibility to these modes by the stability parameter β = T r o t / | W |   . Here, T r o t   is the rotational kinetic energy and W   is the gravitational potential energy. Dynamical rotational instabilities, driven by Newtonian hydrodynamics and gravity, develop on the order of the rotation period of the object. For the uniform-density, incompressible, uniformly rotating MacLaurin spheroids, the dynamical bar-mode instability sets in at β d ≈ 0.27   . For differentially rotating fluids with a polytropic equation of state, numerical simulations have determined that the stability limit β d ≈ 0.27   is valid for initial angular momentum distributions that are similar to those of MacLaurin spheroids [180, 54, 130, 97, 157, 92, 201] . If the object has an off-center density maximum, β d   could be as low as 0.10   [200, 210, 157, 41] . Furthermore, recent work by Shibata, Kariono, and Eriguchi [184] suggests that β d   could be as low as ∼ 0.01   for stars with a large degree of differential rotation. General relativity may enhance the dynamical bar-mode instability by slightly reducing β d   [182, 170] . Secular rotational instabilities are driven by dissipative processes such as gravitational radiation reaction and viscosity. When this type of instability arises, it develops on the timescale of the relevant dissipative mechanism, which can be much longer than the rotation period. The secular bar-mode instability limit for MacLaurin spheroids is β s ≈ 0.14   .
In rotating stars, gravitational radiation reaction drives the r   -modes toward unstable growth [5, 67] .
In hot, rapidly rotating neutron stars, this instability may not be suppressed by internal dissipative mechanisms (such as viscosity and magnetic fields) [121] . If not limited, the dimensionless amplitude α   of the dominant ( m = 2   ) r   -mode will grow to order unity within ten minutes of the formation of a neutron star rotating with a millisecond period. The emitted GWs carry away angular momentum, and will cause the newly formed neutron star to spin down over time. The spindown timescale and the strength of the GWs themselves are directly dependent on the maximum value α m a x   to which the amplitude is allowed to grow [122, 123] . Originally, it was thought that α m a x ∼ 1   . Later work indicated that α m a x   may be ≥ 3   [122, 189, 173, 123] . Some research suggests that magnetic fields, hyperon cooling, and hyperon bulk viscosity may limit the growth of the r   -mode instability, even in nascent neutron stars [110, 109, 163, 164, 123, 120, 80, 6] (significant uncertainties remain regarding the efficacy of these dissipative mechanisms). Furthermore, a recent study of a simple barotropic neutron star model by Arras et al. [8] suggests that multimode couplings could limit α m a x   to values ≪ 1   . If α m a x   is indeed ≪ 1   , GW emission from r   -modes in collapsed remnants is likely undetectable. For the sake of completeness, an analysis of GW emission from r   -modes (which assumes α m a x ∼ 1   ) is presented in the remainder of this paper. However, because it is quite doubtful that α m a x   is sizeable, r   -mode sources are omitted from figures comparing source strengths and detector sensitivities and from discussions of likely detectable sources in the concluding section.
There is some numerical evidence that a collapsing star may fragment into two or more orbiting clumps [75] . If this does indeed occur, the orbiting fragments would be a strong GW source.

2.4 Numerical predictions of GW emission

2.4.1 Full collapse simulations

The accretion-induced collapse of white dwarfs has been simulated by a number of groups [13, 133, 211, 69] .
The majority of these simulations have been Newtonian and have focused on mass ejection and neutrino and γ   -ray emission during the collapse and its aftermath (note that neglecting relativistic effects likely introduces an error of order ( v / c ) 2 ∼ 10 %   for the neutron star remnants of AIC; see below). The most sophisticated are those carried out by Fryer et al. [69] , as they include realistic equations of state, neutrino transport, and rotating progenitors.
As a part of their general evaluation of upper limits to GW emission from gravitational collapse, Fryer, Holz, and Hughes (hereafter, FHH) [71] examined an AIC simulation (model 3) of Fryer et al. [69] . FHH used both numerical and analytical techniques to estimate the peak amplitude h p k   , energy E G W   , and frequency f G W   of the gravitational radiation emitted during the collapse simulations they studied.
For direct numerical computation of the GWs emitted in these simulations, FHH used the quadrupole approximation, valid for nearly Newtonian sources [136] . This approximation is standardly used to compute the GW emission in Newtonian simulations. The reduced or traceless quadrupole moment of the source can be expressed as
I - i j = ∫ ρ ( x i x j − 1 3 δ i j r 2 ) d 3 r , (1)
where i , j = 1 , 2 , 3   are spatial indices and r = x 2 + y 2 + z 2   is the distance to the source. The two polarizations of the gravitational wave field, h +   and h ×   , can be computed in terms of I - . . . i j   , where an overdot indicates a time derivative d / d t   . The energy E G W   is a function of I - . . . i j   . Equation ( 1 ) can be used to calculate I - i j   from the results of a numerical simulation by direct summation over the computational grid. Numerical time derivatives of I - i j   can then be taken to compute h +   , h ×   , and E G W   . However, successive application of numerical derivatives generally introduces artificial noise. Methods for computing I - . . i j   and I - . . . i j   without taking numerical time derivatives have been developed [65, 20, 137] . These methods recast the time derivatives of I - i j   as spatial derivatives of hydrodynamic quantities computed in the collapse simulation (including the density, velocities, and gravitational potential). Thus, instantaneous values for I - . . i j   and I - . . . i j   can be computed on a single numerical time slice (see [137, 217] for details). Note that FHH define h p k   as the maximum value of the rms strain h = 〈 h + 2 + h × 2 〉   , where angular brackets indicate that averages have been taken over both wavelength and viewing angle on the sky.
Errors resulting from the neglect of general relativistic effects (in collapse evolutions as a whole and in GW emission estimations like the quadrupole approximation) are of order ( v / c ) 2   . These errors are typically ∼ 10 %   for neutron star remnants of AIC, < 30 %   for neutron star remnants of massive stellar collapse, and > 30 %   for black hole remnants. Neglect of general relativity in rotational collapse studies is of special concern because relativistic effects counteract the stabilizing effects of rotation (see Section  3.4 ).
Because the code used in the collapse simulations examined by FHH [69] was axisymmetric, their use of the numerical quadrupole approximation discussed above does not account for GW emission that may occur due to non-axisymmetric mass flow. The GWs computed directly from their simulations come only from polar oscillations (which are significant when the mass flow during collapse [or explosion] is largely aspherical).
In order to predict the GW emission produced by non-axisymmetric instabilities, FHH employed rough analytical estimates. The expressions they used to approximate the rms strain h   and the power P = d E G W / d t   of the GWs emitted by a star that has encountered the bar-mode instability are
h b a r = 32 45 G c 4 m r 2 ω 2 d (2)
and
P b a r = 32 45 G c 5 m 2 r 4 ω 6 . (3)
Here m   , 2 r   , and ω   are the mass, length, and angular frequency of the bar and d   is the distance to the source. FHH vary the mass m   assumed to be enclosed by the bar (which has a corresponding length 2 r   ) and compute the characteristics of the GW emission as a function of this enclosed mass.
For simplicity, FHH assumed that a fragmentation instability will cause a star to break into two clumps (although more clumps could certainly be produced). Their estimates for the rms strain and power radiated by the orbiting binary fragments are
h b i n = 128 5 G c 4 m r 2 ω 2 d (4)
and
P b i n = 128 5 G c 5 m 2 r 4 ω 6 . (5)
Here m   is the mass of a single fragment, 2 r   is the separation of the fragments, and ω   is their orbital frequency.
For their computation of the GWs radiated via r-modes, FHH used the method of Ho and Lai [89] (which assumes α m a x = 1   ) and calculated only the emission from the dominant m = 2   mode. This approach is detailed in FHH. If the neutron star mass and initial radius are taken to be 1.4 M ⊙   and 12.53 k m   , respectively, the resulting formula for the average GW strain is
h ( t ) = 1.8 × 10 − 24 α ( ν s 1 k H z ) ( 20 M p c d ) , (6)
where α   is the mode amplitude and ν s   is the spin frequency.
FHH's numerical quadrupole estimate of the GWs from polar oscillations in the AIC simulation of Fryer et al. [69] predicts a peak dimensionless amplitude h p k = 5.9 × 10 − 24   (for d = 100 M p c   ).
The energy E G W = 3 × 10 45 e r g   is emitted at a frequency of f G W ≈ 50 H z   . This amplitude is about an order of magnitude too small to be observed by the advanced LIGO-II detector. The sensitivity curve for the broadband configuration of the LIGO-II detector is shown in Figure  1 (see Appendix A of [71] for details on the computation of this curve). Note that the characteristic strain h   is plotted along the vertical axis in Figure  1 (and in LISA's sensitivity curve, shown in Figure  12 ). For burst sources, h = h p k   . For sources that persist for N   cycles, h = N h p k   .

Figure 1 : A comparison between the GW amplitude h ( f )   for various sources and the LIGO-II sensitivity curve. See the text for details regarding the computations of h   . The AIC sources are assumed to be located at a distance of 100 M p c   ; the SNe sources at 10 M p c   ; and the Population III sources at a luminosity distance of ∼ 50 G p c   . Secular bar-mode sources are identified with an (s), dynamical bar-modes with a (d).

The simulation of Fryer et al. [69] considered by FHH does produce an object with an off-center density maximum. However, because the maximum β   reached in the AIC simulation of Fryer et al.
was < 0.06   (and because the degree of differential rotation present in the remnant was low [68] ), the collapsing object is not likely to encounter dynamical rotational instabilities. The choice of initial angular momentum J   for the progenitor in an AIC simulation can affect this outcome. Fryer et al.'s choice of J = 10 49 g c m 2 s − 1   is such that the period of the white dwarf progenitor is 10 s   less than the shortest observed period 30 s   of a cataclysmic variable white dwarf [112] . If higher values of J   exist in accreting white dwarfs, bar-mode instabilities may be more likely to occur (see the discussion of work by Liu and Lindblom below).
According to FHH, the remnant of this AIC simulation will be susceptible to r   -mode growth.
Assuming α m a x ∼ 1   (which is likely not physical; see Section  2.3 ), they predict E G W > 10 52 e r g   .
FHH compute h ( f G W )   for coherent observation of the neutron star as it spins down over the course of a year. For a neutron star located at a distance of 100 M p c   , this track is always below the LIGO-II noise curve. The point on this track with the maximum h   , which corresponds to the beginning of r   -mode evolution, is shown in Figure  1 . The track moves down and to the left (i.e., h   and f   decrease) in this figure as the r   -mode evolution continues.

2.4.2 Analyses of equilibrium representations of collapse

In addition to full hydrodynamics collapse simulations, many studies of gravitational collapse have used hydrostatic equilibrium models to represent stars at various stages in the collapse process. Some investigators use sequences of equilibrium models to represent snapshots of the phases of collapse (e.g., [14, 153, 125] ). Others use individual equilibrium models as initial conditions for hydrodynamical simulations(e.g., [187, 157, 152, 41] ). Such simulations represent the approximate evolution of a model beginning at some intermediate phase during collapse or the evolution of a collapsed remnant. These studies do not typically follow the intricate details of the collapse itself. Instead, their goals include determining the stability of models against the development of non-axisymmetric modes and estimation of the characteristics of any resulting GW emission.
Liu and Lindblom [125, 124] have applied this equilibrium approach to AIC. Their investigation began with a study of equilibrium models built to represent neutron stars formed from AIC [125] .
These neutron star models were created via a two-step process, using a Newtonian version of Hachisu's self-consistent field method [76] . Hachisu's method ensures that the forces due to the centrifugal and gravitational potentials and the pressure are in balance in the equilibrium configuration.
Liu and Lindblom's process of building the nascent neutron stars began with the construction of rapidly rotating, pre-collapse white dwarf models. Their Models I and II are C-O white dwarfs with central densities ρ c = 10 10   and 6 × 10 9 g c m − 3   , respectively (recall this is the range of densities for which AIC is likely for C-O white dwarfs). Their Model III is an O-Ne-Mg white dwarf that has ρ c = 4 × 10 9 g c m − 3   (recall this is the density at which collapse is induced by electron capture). All three models are uniformly rotating, with the maximum allowed angular velocities. The models' values of total angular momentum are roughly 3 – 4 times that of Fryer et al.'s AIC progenitor Model 3 [69] . The realistic equation of state used to construct the white dwarfs is a Coulomb corrected, zero temperature, degenerate gas equation of state [171, 44] .
In the second step of their process, Liu and Lindblom [125] built equilibrium models of the collapsed neutron stars themselves. The mass, total angular momentum, and specific angular momentum distribution of each neutron star remnant is identical to that of its white dwarf progenitor (see Section 3 of [125] for justification of the specific angular momentum conservation assumption). These models were built with two different realistic neutron star equations of state.
Liu and Lindblom's cold neutron star remnants had values of the stability parameter β   ranging from 0.23 – 0.26   . It is interesting to compare these results with those of Zwerger and Müller [217] .
Zwerger and Müller performed axisymmetric hydrodynamics simulations of stars with polytropic equations of state ( P ∝ ρ Γ   ). Their initial models were Γ = 4 / 3   polytropes, representative of massive white dwarfs. All of their models started with ρ c = 10 10 g c m − 3   . Their model that was closest to being in uniform rotation (A1B3) had 22% less total angular momentum than Liu and Lindblom's Model I. The collapse simulations of Zwerger and Müller that started with model A1B3 all resulted in remnants with values of β < 0.07   . Comparison of the results of these two studies could indicate that the equation of state may play a significant role in determining the structure of collapsed remnants. Or it could suggest that the assumptions employed in the simplified investigation of Liu and Lindblom are not fully appropriate. Zwerger and Müller's work will be discussed in much more detail in Section  3 , as it was performed in the context of core collapse supernovae.
In a continuation of the work of Liu and Lindblom, Liu [124] used linearized hydrodynamics to perform a stability analysis of the cold neutron star AIC remnants of Liu and Lindblom [125] .
He found that only the remnant of the O-Ne-Mg white dwarf (Liu and Lindblom's Model III) developed the dynamical bar-mode ( m = 2   ) instability. This model had an initial β = 0.26   . Note that the m = 1   mode, observed by others to be the dominant mode in unstable models with values of β   much lower than 0.27   [200, 210, 157, 41] , did not grow in his simulation. Because Liu and Lindblom's Models I and II had lower values of β   , Liu identified the onset of instability for neutron stars formed via AIC as β d ≈ 0.25   . Liu estimated the peak amplitude of the GWs emitted by the Model III remnant to be h p k ≈ 1.4 × 10 − 24   and the LIGO-II signal-to-noise ratio (for a persistent signal like that seen in the work of [152] and [28] ) to be S / N ≤ 3   (for f G W ≈ 450 H z   ). These values are for a source located at 100 M p c   . He also predicted that the timescale for gravitational radiation to carry away enough angular momentum to eliminate the bar-mode is τ G W ∼ 7 s   ( ∼ 3 × 10 3   cycles). Thus, h ∼ 8 × 10 − 23   .
(Note that this value for h   is merely an upper limit as it assumes that the amplitude and frequency of the GWs do not change over the 7 s   during which they are emitted. Of course, they will change as angular momentum is carried away from the object via GW emission.) Such a signal may be marginally detectable with LIGO-II (see Figure  1 ). Details of the approximations on which these estimates are based can be found in [124] .
Liu cautions that his results hold if the magnetic field of the proto-neutron star is B ≤ 10 12 G   . If the magnetic field is larger, then it may have time to suppress some of the neutron star's differential rotation before it cools. This would make bar formation less likely. Such a large field could only result if the white dwarf progenitor's B   field was ≥ 10 8 G   . Observation-based estimates suggest that about 25 %   of white dwarfs in interacting close binaries (cataclysmic variables) are magnetic and that the field strengths for these stars are ∼ 10 7 – 3 × 10 8 G   [206] .

2.5 Going further

The AIC scenario is generally discussed in terms of the collapse of an accreting white dwarf to a neutron star. However, Shibata, Baumgarte, and Shapiro have examined the collapse of a rotating, supramassive neutron star to a black hole [183] . Such supramassive neutron stars (with masses greater than the maximum mass for a nonrotating neutron star) could be formed and pushed to collapse via accretion from a binary companion. They performed 3D, fully general relativistic hydrodynamics simulations of uniformly rotating neutron stars. Dynamical non-axisymmetric instabilities (such as the bar-mode) did not have time to grow in their simulations prior to black hole formation. Differentially rotating neutron star progenitors could have higher values of β   than the uniformly rotating models used in this study and may be susceptible to non-axisymmetric instabilities on a shorter timescale.

3 Collapse of Massive Stars

3.1 Collapse scenario

Stars with mass greater than ∼ 8 M ⊙   will undergo core collapse at the end of their thermonuclear burning life cycles. The collapse will trigger a supernova explosion if the star's mass is < 40 – 50 M ⊙   (for reviews, see [7, 18, 32] ). Core collapse SNe include Types II and Ib/Ic. SNe Ib/Ic are distinguished from Type II SNe by the absence of hydrogen in their spectra. SNe Ib/Ic are thought to result from the collapse of the cores of massive stars that have lost their hydrogen envelopes (and possibly part of their helium envelopes) by stellar winds or by mass transfer. The SN Ib/Ic progenitors that lose their outer envelopes via stellar winds are known as Wolf–Rayet stars and have initial masses ∼ > 30 M ⊙   ; those that undergo mass loss via mass transfer in binaries have progenitor masses in the range 12 M ⊙ – 18 M ⊙   [143] .
It is believed that the remnant of a core collapse SN is a neutron star if the mass of the progenitor is less than ∼ 20 – 25 M ⊙   [73, 71] . If the progenitor's mass is in the range 20 − 25 ∼ < M ∼ < 40 – 50 M ⊙   , not the entire star is ejected in the SN explosion. More than 2 M ⊙   will fall back onto the nascent neutron star and lead to black hole formation. If the progenitor's mass exceeds roughly 40 – 50 M ⊙   , then no explosion will occur and the star will collapse directly to a black hole. These objects are known as collapsars. However, it is unclear if high metallicity stars with M ∼ > 40 – 50 M ⊙   actually end their lives in collapse or are prevented from doing so by stellar wind driven mass loss [71] . Note that the limits on the progenitor masses quoted in this paragraph (especially the 40 – 50 M ⊙   lower limit for direct black hole formation) are uncertain because the progenitor mass dependence of the neutrino explosion mechanism (see below) is unknown [81, 147] .
The massive iron cores of SN II/Ib/Ic progenitors are supported by both thermal and electron degeneracy pressures. The density and temperature of such a core will eventually rise, due to the build up of matter consumed by thermonuclear burning, to the point where electron capture and photodissociation of nuclei begin. Dissociation lowers the photon and electron temperatures and thereby reduces the core's thermal support [62] . Electron capture reduces the electron degeneracy pressure. One or both of these processes will trigger the collapse of the core. The relative importance of dissociation and electron capture in instigating collapse is determined by the mass of the star [62] . The more massive the core, the bigger is the role played by dissociation.
Approximately 70 %   of the inner portion of the core collapses homologously and subsonically.
The outer core collapses at supersonic speeds [62, 137] . The maximum velocity of the outer regions of the core reaches ∼ 7 × 10 4 k m s − 1   . It takes just 1 s for an earth-sized core to collapse to a radius of 50 km [7] .
The inward collapse of the core is halted by nuclear forces when its central density ρ c   is 2 – 10 times the density of nuclear material [12, 11] . The core overshoots its equilibrium position and bounces. A shock wave is formed when the supersonically infalling outer layers hit the rebounding inner core. If the inner core pushes the shock outwards with energy E > 10 51 e r g   (supplied by the binding energy of the nascent neutron star), then the remainder of the star can be ejected in about 20 ms [7] . This so-called “prompt explosion” mechanism will only succeed if a very soft supra-nuclear equation of state is used in conjunction with a relatively small core ( M ∼ < 1.35 M ⊙   , derived from a very low mass progenitor) and a large portion of the collapse proceeds homologously [18, 143] .
Inclusion of general relativistic effects in collapse simulations can increase the success of the prompt mechanism in some cases [12, 190] .
Both dissociation of nuclei and electron capture can reduce the ejection energy, causing the prompt mechanism to fail. The shock will then stall at a radius in the range 100 – 200 km. Colgate and White [47] were the first to suggest that energy from neutrinos emitted by the collapsed core could revive the stalled shock. (See Burrows and Thompson [36] for a recent review of core collapse neutrino processes.) However, their simulations did not include the physics necessary to accurately model this “delayed explosion” mechanism. Wilson and collaborators were the first to perform collapse simulations with successful delayed ejections [24, 23, 207, 19, 209, 208] . However, their simulations and those of others had difficulty producing energies high enough to match observations [45, 30, 104] . Observations of SN 1987a show that significant mixing occurred during this supernovae (see Arnett et al. [7] for a review). Such mixing can be attributed to nonradial motion resulting from fluid instabilities. Convective instabilities play a significant role in the current picture of the delayed explosion mechanism. The outer regions of the nascent neutron star are convectively unstable after the shock stalls (for an interval of 10 – 100 ms after bounce) due to the presence of negative lepton and energy gradients [143] . This has been confirmed by both 2D and 3D simulations [33, 34, 104, 140, 141, 105, 143, 134, 100, 102, 159] . Convective motion is more effective at transporting neutrinos out of the proto-neutron star than is diffusion. Less than 10 %   of the neutrinos emitted by the neutron star need to be absorbed and converted to kinetic energy for the shock to be revived [143] . The “hot bubble” region above the surface of the neutron star also has been shown to be convectively unstable [45, 18, 46, 74] . Janka and Müller have demonstrated that convection in this region only aids the explosion if the neutrino luminosity is in a narrow region [106] . Some simulations that include advanced neutrino transport methods have cast doubt on the ability of convection to ensure the success of the delayed explosion mechanism [135] . However, Rampp and Janka [102, 159] have recently performed collapse simulations with sophisticated neutrino physics (including Boltzmann transport and state-of-the-art neutrino-matter interactions) that are very close to producing explosions. Their results indicate that successful explosions could be very sensitive to the details of the neutrino physics. Multi-dimensional simulations that include the full effects of general relativity, rotation, and an improved treatment of neutrino interactions will likely be necessary to properly model the delayed explosion mechanism [103, 102] .
In addition to the mixing seen in SN 1987a, observations of (i) polarization in the spectra of several core collapse SNe, (ii) jets in the Cas A remnant, and (iii) kicks in neutron stars suggest that supernovae are inherently aspherical (see [7, 3, 90] and references therein). Note that these asphericities could originate in the central explosion mechanism itself and/or the mechanism(s) for energy transfer between the core and ejecta [102] . These observations partly motivated the multi-dimensional studies of convection in the delayed explosion mechanism discussed above and have also spurred investigations of the role jets may play in SN explosions [128, 205, 212] . Recent 3D radiation hydrodynamics simulations performed by Höflich et al. [90] indicate that low velocity jets stalled inside SN envelopes can account for the observed asymmetries. Possible sources of such asymmetries include magnetic fields [15, 3] , anisotropic neutrino emission, and convection.

3.2 Formation rate

Type II/Ib/Ic supernovae are observed to occur in only spiral and irregular galaxies. The most thorough computation of SN rates is that of Cappellaro et al. [38] . Their sample includes 137 SNe from five different SN searches. They determined that the core collapse SN rate in the Galaxy is 0.6 × 10 − 3 – 1.6 × 10 − 2 y r − 1   . Thus, a reasonable occurrence rate can be found for an observation distance of 10 Mpc. (Note that a recent infrared survey estimates that the rate may actually be an order of magnitude higher [129] .) Approximately 5 – 40% of these SNe will leave behind black hole remnants [73] . The formation rate of collapsars, the massive cores that collapse to black holes without an accompanying SNe explosion, is unknown. This is because of the uncertainty in stellar wind-driven mass-loss rates [73] .

3.3 GW emission mechanisms

Gravitational radiation will be emitted during the collapse/explosion of a core collapse SN due to the star's changing quadrupole moment. A rough description of the possible evolution of the quadrupole moment is given in the remainder of this paragraph. During the first 100 – 250 ms of the collapse, as the core contracts and flattens, the magnitude of the quadrupole moment I - i j   will increase. The contraction speeds up over the next 20 ms and the density distribution becomes a centrally condensed torus [137] . In this phase the core's shrinking size dominates its increasing deformation and the magnitude of I - i j   decreases. As the core bounces, I - i j   changes rapidly due to the deceleration and rebound. If the bounce occurs because of nuclear pressure, its timescale will be < 1 m s   . If centrifugal forces play a role in halting the collapse, the bounce can last up to several m s   [137] . The magnitude of I - i j   will increase due to the core's expansion after bounce. As the resulting shock moves outwards, the unshocked portion of the core will undergo oscillations, causing I - i j   to oscillate as well. The shape of the core, the depth of the bounce, the bounce timescale, and the rotational energy of the core all strongly affect the GW emission. For further details see [62, 137] .
Convectively driven inhomogeneities in the density distribution of the outer regions of the nascent neutron star and anisotropic neutrino emission are other sources of GW emission during the collapse/explosion [35, 146] .
As discussed in the case of AIC, global rotational instabilities (such as the m = 2   bar-mode) may develop during the collapse itself or in a neutron star remnant. A neutron star remnant will likely also be susceptible to the radiation reaction driven r   -modes. Both of these types of instabilities will emit GWs, as will a fragmentation instability if one occurs. See Section  2.3 for further details regarding these instabilities.
If the collapsed remnant is a black hole, GWs will be radiated as the infall of the remaining stellar matter distorts the black hole's geometry. This “ringdown phase” will end when gravitational radiation has dissipated all of the black hole's accretion-induced distortion. Zanotti, Rezzolla, and Font [215] have recently suggested that the torus of matter surrounding the black hole may be an even stronger source of GWs than the collapse itself (see also [113, 204] ). When perturbed, such “toroidal neutron stars” may undergo regular oscillations and thus emit copious GWs.

3.4 Numerical predictions of GW emission

3.4.1 Historical investigations

The collapse of the progenitors of core collapse supernovae has been investigated as a source of gravitational radiation for more than three decades. In an early study published in 1971, Ruffini and Wheeler [165] identified mechanisms related to core collapse that could produce GWs and provided order-of-magnitude estimates of the characteristics of such emission.
Quantitative computations of GW emission during the infall phase of collapse were performed by Thuan and Ostriker [198] and Epstein and Wagoner [57, 56] , who simulated the collapse of oblate dust spheroids. Thuan and Ostriker used Newtonian gravity and computed the emitted radiation in the quadrupole approximation. Epstein and Wagoner discovered that post-Newtonian effects prolonged the collapse and thus lowered the GW luminosity. Subsequently, Novikov [156] and Shapiro and Saenz [177, 166] included internal pressure in their collapse simulations and were thus able to examine the GWs emitted as collapsing cores bounced at nuclear densities. The quadrupole GWs from the ringdown of the collapse remnant were initially investigated by the perturbation study of Turner and Wagoner [202] and later by Saenz and Shapiro [167, 168] .
Müller [139] calculated the quadrupole GW emission from 2D axisymmetric collapse based on the Newtonian simulations of Müller and Hillebrandt [144] (these simulations used a realistic equation of state and included differential rotation). He found that differential rotation enhanced the efficiency of the GW emission.
The first fully general relativistic investigations of stellar core collapse were Nakamura's 2D simulations of rotating collapse [148, 149] . However, because of the limits of his numerical formalism and computational resources, he was unable to compute the emitted gravitational radiation (the energy of this emission is quite small compared to the rest mass energy and thus was difficult to extract numerically). The results of this work indicate that collapse does not lead to black hole formation if the parameter q = J / M 2   exceeds unity (here J   and M   are the angular momentum and gravitational mass of the remnant).
Stark and Piran [188, 158] were the first to compute the GW emission from fully relativistic collapse simulations, using the ground-breaking formalism of Bardeen and Piran [10] . They followed the (pressure-cut induced) collapse of rotating polytropes in 2D. Their work focused in part on the conditions for black hole formation and the nature of the resulting ringdown waveform, which they found could be described by the quasi-normal modes of a rotating black hole. In each of their simulations, less than 1 %   of the gravitational mass was converted to GW energy.
Seidel and collaborators also studied the effects of general relativity on the GW emission during collapse and bounce [175, 176] . They employed a perturbative approach, valid only in the slowly rotating regime.
The gravitational radiation from non-axisymmetric collapse was investigated by Detweiler and Lindblom, who used a sequence of non-axisymmetric ellipsoids to represent the collapse evolution [50] . They found that the radiation from their analysis of non-axisymmetric collapse was emitted over a more narrow range of frequency than in previous studies of axisymmetric collapse.
For further discussion of the first two decades of study of the GW emission from stellar collapse see [63] . In the remainder of Section  3.4 , more recent investigations will be discussed.

3.4.2 Axisymmetric simulations

The core collapse simulations of Mönchmeyer et al. began with better iron core models and a more realistic microphysical treatment (including a realistic equation of state, electron capture, and a simple neutrino transport scheme) than any previous study of GW emission from axisymmetric stellar core collapse [137] . The shortcomings of their investigation included initial models that were not in rotational equilibrium, an equation of state that was somewhat stiff in the subnuclear regime, and the use of Newtonian gravity. Each of their four models had a different initial angular momentum profile. The rotational energies of the models ranged from 0.1 – 0.45 of the maximum possible rotational energy.
The collapses of three of the four models of Mönchmeyer et al. were halted by centrifugal forces at subnuclear densities. This type of low ρ c   bounce had been predicted by Shapiro and Lightman [178] and Tohline [199] (in the context of the “fizzler” scenario for failed supernovae; see also [82, 83, 96] ), and had been observed in earlier collapse simulations [145, 191] . Mönchmeyer and collaborators found that a bounce caused by centrifugal forces would last for several m s   , whereas a bounce at nuclear densities would occur in < 1 m s   . They also determined that a subnuclear bounce produced larger amplitude oscillations in density and radius, with larger oscillation periods, than a bounce initiated by nuclear forces alone. They pointed out that these differences in timescale and oscillatory behavior should affect the GW signal. Therefore, the GW emission could indicate whether the bounce was a result of centrifugal or nuclear forces.
Mönchmeyer et al. identified two different types of waveforms in their models (computed using the numerical quadrupole approximation discussed in Section  2.4 ). The waveforms they categorized as Type I (similar to those observed in previous collapse simulations [139, 65] ) are distinguished by a large amplitude peak at bounce and subsequent damped ringdown oscillations. They noted that Type I signals were produced by cores that bounced at nuclear densities (or bounced at subnuclear densities if the cores had small ratios of radial kinetic to rotational kinetic energies).
The quadrupole gravitational wave amplitude A 20 E 2   for a Type I waveform is shown in Figure  2 (see [194, 217] for expressions relating A 20 E 2   to h   ). The waveforms identified as Type II exhibit several maxima, which result from multiple bounces (see Figure  3 for an example of a Type II waveform).
Note that the waveforms displayed in Figures  2 and  3 are from the study of Zwerger and Müller [217] , discussed below.

Figure 2 : Type I waveform (quadrupole amplitude A 20 E 2   as a function of time) from one of Zwerger and Müller's [217] simulations of a collapsing polytrope. The vertical dotted line marks the time at which the first bounce occurred. (Figure 5d of [217] ; used with permission.)

Figure 3 : Type II waveform (quadrupole amplitude A 20 E 2   as a function of time) from one of Zwerger and Müller's [217] simulations of a collapsing polytrope. The vertical dotted line marks the time at which bounce occurred. (Figure 5a of [217] ; used with permission.)

The model of Mönchmeyer et al. that bounced due to nuclear forces had the highest GW amplitude of all of their models, h p k ∼ 10 − 23   for a source distance d = 10 M p c   , and the largest emitted energy E G W ∼ 10 47 e r g   . The accompanying power spectrum peaked in the frequency range 5 × 10 2 – 10 3 H z   .
The most extensive Newtonian survey of the parameter space of axisymmetric, rotational core collapse is that of Zwerger and Müller [217] . They simulated the collapse of 78 initial models with varying amounts of rotational kinetic energy (reflected in the initial value of the stability parameter β i   ), differential rotation, and equation of state stiffness. In order to make this large survey tractable, they used a simplified equation of state and did not explicitly account for electron capture or neutrino transport. Their initial models were constructed in rotational equilibrium via the method of Eriguchi and Müller [58] . The models had a polytropic equation of state, with initial adiabatic index Γ i = 4 / 3   . Collapse was induced by reducing the adiabatic index to a value Γ r   in the range 1.28 – 1.325. The equation of state used during the collapse evolution had both polytropic and thermal contributions.

Figure 4 : Type III waveform (quadrupole amplitude A 20 E 2   as a function of time) from one of Zwerger and Müller's [217] simulations of a collapsing polytrope. The vertical dotted line marks the time at which bounce occurred. (Figure 5e of [217] ; used with permission).

The major result of Zwerger and Müller's investigation was that the signal type of the emitted gravitational waveform in their runs was determined by the stiffness of the equation of state of the collapsing core (i.e., the value of Γ r   ). In their simulations, Type I signals (as labelled by Mönchmeyer et al. [137] ) were produced by models with relatively soft equations of state, Γ r ∼ < 1.31   . Type II signals were produced by the models with stiffer equations of state, Γ r ∼ > 1.32   . They found a smooth transition between these signal types if Γ r   was increased while all other parameters were held fixed. They also observed another class of signal, Type III, for their models with the lowest Γ r   ( = 1.28   ). Type III waveforms have a large positive peak just prior to bounce, a smaller negative peak just after bounce, and smaller subsequent oscillations with very short periods (see Figure  4 ). Type III signals were not observed in the evolution of strongly differentially rotating Γ r = 1.28   models. Their waveforms were computed with the same technique used in [137] .
In contrast to the results of Mönchmeyer et al. [137] , in Zwerger and Müller's investigation the value of ρ c   at bounce did not determine the signal type. Instead, the only effect on the waveform due to ρ c   was a decrease in h p k   in models that bounced at subnuclear densities. The effect of the initial value of β i   on h p k   was non-monotonic. For models with β i ∼ < 0.1   , h p k   increased with increasing β i   . This is because the deformation of the core is larger for faster rotators. However, for models with larger β i   , h p k   decreases as β i   increases. These models bounce at subnuclear densities. Thus, the resulting acceleration at bounce and the GW amplitude are smaller. Zwerger and Müller found that the maximum value of h p k   for a given sequence was reached when ρ c   at bounce was just less than ρ n u c   . The degree of differential rotation did not have a large effect on the emitted waveforms computed by Zwerger and Müller. However, they did find that models with soft equations of state emitted stronger signals as the degree of differential rotation increased.
The models of Zwerger and Müller that produced the largest GW signals fell into two categories:
those with stiff equations of state and β i < 0.01   ; and those with soft equations of state, β i ≥ 0.018   , and large degrees of differential rotation. The GW amplitudes emitted during their simulations were in the range 4 × 10 − 25 ∼ < h ∼ < 4 × 10 − 23   , for d = 10 M p c   (the model with the highest h   is identified in Figure  1 ). The corresponding energies ranged from 10 44 ∼ < E G W ∼ < 10 47 e r g   .
The peaks of their power spectra were between 500 Hz and 1 kHz. Such signals would fall just outside of the range of LIGO-II. Yamada and Sato [213] used techniques very similar to those of Zwerger and Müller [217] in their core collapse study. Their investigation revealed that the h p k   for Type I signals became saturated when the dimensionless angular momentum of the collapsing core, q = J / ( 2 G M / c )   , reached ∼ 0.5   .
They also found that h p k   was sensitive to the stiffness of the equation of state for densities just below ρ n u c   . The characteristics of the GW emission from their models were similar to those of Zwerger and Müller.

3.4.3 Non-axisymmetric simulations

The GW emission from non-axisymmetric hydrodynamics simulations of stellar collapse was first studied by Bonazzola and Marck [131, 22] . They used a Newtonian, pseudo-spectral hydrodynamics code to follow the collapse of polytropic models. Their simulations covered only the pre-bounce phase of the collapse. They found that the magnitudes of h p k   in their 3D simulations were within a factor of two of those from equivalent 2D simulations and that the gravitational radiation efficiency did not depend on the equation of state.
The first 3D hydrodynamics collapse simulations to study the GW emission well beyond the core bounce phase were performed by Rampp, Müller, and Ruffert [160] . These authors started their Newtonian simulations with the only model (A4B5G5) of Zwerger and Müller [217] that had a post-bounce value for the stability parameter β = 0.35   that significantly exceeded 0.27   (recall this is the value at which the dynamical bar instability sets in for MacLaurin spheroid-like models).
This model had the softest EOS ( Γ r = 1.28   ), highest β i = 0.04   , and largest degree of differential rotation of all of Zwerger and Müller's models. The model's initial density distribution had an off-center density maximum (and therefore a torus-like structure). Rampp, Müller, and Ruffert evolved this model with a 2D hydrodynamics code until its β   reached ≈ 0.1   . At that point, 2.5 ms prior to bounce, the configuration was mapped onto a 3D nested cubical grid structure and evolved with a 3D hydrodynamics code.
Before the 3D simulations started, non-axisymmetric density perturbations were imposed to seed the growth of any non-axisymmetric modes to which the configuration was unstable. When the imposed perturbation was random (5% in magnitude), the dominant mode that arose was m = 4   .
The growth of this particular mode was instigated by the cubical nature of the computational grid.
When an m = 3   perturbation was imposed (10% in magnitude), three clumps developed during the post-bounce evolution and produced three spiral arms. These arms carried mass and angular momentum away from the center of the core. The arms eventually merged into a bar-like structure (evidence of the presence of the m = 2   mode). Significant non-axisymmetric structure was visible only within the inner 40 km of the core. Their simulations were carried out to ∼ 14 m s   after bounce.
The amplitudes of the emitted gravitational radiation (computed in the quadrupole approximation) were only ∼ 2 %   different from those observed in the 2D simulation of Zwerger and Müller. Because of low angular resolution in the 3D runs, the energy emitted was only 65 %   of that emitted in the corresponding 2D simulation.
The findings of Centrella et al. [41] indicate it is possible that some of the post-bounce configurations of Zwerger and Müller, which have lower values of β   than the model studied by Rampp, Müller, and Ruffert [160] , may also be susceptible to non-axisymmetric instabilities. Centrella et al. have performed 3D hydrodynamics simulations of Γ = 1.3   polytropes to test the stability of configurations with off-center density maxima (as are present in many of the models of Zwerger and Müller [217] ).
The simulations carried out by Centrella and collaborators were not full collapse simulations, but rather began with differentially rotating equilibrium models. These simulations tracked the growth of any unstable non-axisymmetric modes that arose from the initial 1 %   random density perturbations that were imposed. Their results indicate that such models can become dynamically unstable at values of β ∼ > 0.14   . The observed instability had a dominant m = 1   mode.
Centrella et al. estimate that if a stellar core of mass M ∼ 1.4 M ⊙   and radius R ∼ 200 k m   encountered this instability, the values of h p k   from their models would be ∼ 2 × 10 − 24 – 2 × 10 − 23   , for d = 10 M p c   . The frequency at which h p k   occurred in their simulations was ∼ 200 H z   . This instability would have to persist for at least ∼ 15   cycles to be detected with LIGO-II. Brown [29] carried out an investigation of the growth of non-axisymmetric modes in post-bounce cores that was similar in many respects to that of Rampp, Müller, and Ruffert [160] . He performed 3D hydrodynamical simulations of the post-bounce configurations resulting from 2D simulations of core collapse. His pre-collapse initial models are Γ = 4 / 3   polytropes in rotational equilibrium. The differential rotation laws used to construct Brown's initial models were motivated by the stellar evolution study of Heger, Langer, and Woosley [86] . The angular velocity profiles of their pre-collapse progenitors were broad and Gaussian-like. Brown's initial models had peak angular velocities ranging from 0.8 – 2.4 times those of [86] . The model evolved by Rampp, Müller, and Ruffert [160] had much stronger differential rotation than any of Brown's models. To induce collapse, Brown reduced the adiabatic index of his models to Γ = 1.28   , the same value used by [160] .
Brown found that β   increased by a factor ∼ < 2   during his 2D collapse simulations. This is much less than the factor of ∼ 9   observed in the model studied by Rampp, Müller, and Ruffert [160] .
This is likely a result of the larger degree of differential rotation in the model of Rampp et al.
Brown performed 3D simulations of the two most rapidly rotating of his post-bounce models (models Ω 24   and Ω 20   , both of which had β > 0.27   after bounce) and of the model of Rampp et al.
(which, although it starts out with β = 0.35   , has a sustained β < 0.2   ). Brown refers to the Rampp et al. model as model RMR. Because Brown's models do not have off-center density maxima, they are not expected to be unstable to the m = 1   mode observed by Centrella et al. [41] . He imposed random 1 %   density perturbations at the start of all three of these 3D simulations (note that this perturbation was of a much smaller amplitude than those imposed by [160] ).
Brown's simulations determined that both his most rapidly rotating model Ω 24   (with post-bounce β > 0.35   ) and model RMR are unstable to growth of the m = 2   bar-mode. However, his model Ω 20   (with post-bounce β > 0.3   ) was stable. Brown observed no dominant m = 3   or m = 4   modes growing in model RMR at the times at which they were seen in the simulations of Rampp et al.
This suggests that the mode growth in their simulations was a result of the large perturbations they imposed. The m = 2   mode begins to grow in model RMR at about the same time as Rampp et al. stopped their evolutions. No substantial m = 1   growth was observed.
The results of Brown's study indicate that the overall β   of the post-bounce core may not be a good diagnostic for the onset of instability. He found, as did Rampp, Müller, and Ruffert [160] , that only the innermost portion of the core (with ρ > 10 10 g c m − 3   ) is susceptible to the bar-mode.
This is evident in the stability of his model Ω 20   . This model had an overall β > 0.3   , but an inner core with β i c = 0.15   . Brown also observed that the β   of the inner core does not have to exceed 0.27   for the model to encounter the bar-mode. Models Ω 24   and RMR had β i c ≈ 0.19   . He speculates that the inner cores of these later two models may be bar-unstable because interaction with their outer envelopes feeds the instability or because β d < 0.27   for such configurations.
Fryer and Warren [74] have recently performed the first 3D collapse simulations to follow the entire collapse through explosion. They used a smoothed particle hydrodynamics code, a realistic equation of state, the flux-limited diffusion approximation for neutrino transport, and Newtonian spherical gravity. Their initial model was nonrotating. Thus, no bar-mode instabilities could develop during their simulations. The only GW emitting mechanism present in their models was convection in the core. The maximum amplitude h   of this emission, computed in the quadrupole approximation, was ∼ 3 × 10 − 26   , for d = 10 M p c   [72] . Preliminary results from rotating simulations performed with this same 3D code do not show the development of bar or fragmentation instabilities [72] .
The GW emission from nonradial quasinormal mode oscillations in proto-neutron stars has recently been examined by Ferrari, Miniutti, and Pons [61] . They found that the frequencies of emission f G W   during the first second after formation (600 – 1100 Hz for the first fundamental and gravity modes) are significantly lower than the corresponding frequencies for cold neutron stars and thus reside in the bandwidths of terrestrial interferometers. However, for first generation interferometers to detect the GW emission from an oscillating proto-neutron star located at 10 Mpc, with a signal-to-noise ratio of 5, E G W   must be ∼ 10 − 3 – 10 − 2 M ⊙ c 2   . It is unlikely that this much energy is stored in these modes (the collapse itself may only emit ∼ 10 − 7 M ⊙ c 2   in gravitational waves [53] ).

3.4.4 General relativistic simulations

General relativistic effects oppose the stabilizing influence of rotation in pre-collapse cores.
Thus, stars that might be prevented from collapsing due to rotational support in the Newtonian limit may collapse when general relativistic effects are considered. Furthermore, general relativity will cause rotating stars undergoing collapse to bounce at higher densities than in the Newtonian case [192, 217, 160, 31] .
The full collapse simulations of Fryer and Heger [70] are the most sophisticated axisymmetric simulations from which the resultant GW emission has been studied [71, 72] . Fryer and Heger include the effects of general relativity, but assume (for the purposes of their gravity treatment only) that the mass distribution is spherical. The GW emission from these simulations was evaluated with either the quadrupole approximation or simpler estimates (see below).
The work of Fryer and Heger [70] is an improvement over past collapse investigations because it starts with rotating progenitors evolved to collapse with a stellar evolution code (which incorporates angular momentum transport via an approximate diffusion scheme) [84] , incorporates realistic equations of state and neutrino transport, and follows the collapse to late times. The values of total angular momentum of the inner cores of Fryer and Heger ( 0.95 – 1.9 × 10 49 g c m 2 s − 1   ) are lower than has often been assumed in studies of the GW emission from core collapse. Note that the total specific angular momentum of these core models may be lower by about a factor of 10 if magnetic fields were included in the evolution of the progenitors [3, 85, 87] .
FHH's [71] numerical quadrupole estimate of the GWs from polar oscillations in the collapse simulations of Fryer and Heger [70] predicts a peak dimensionless amplitude h p k = 4.1 × 10 − 23   (for d = 10 M p c   ), emitted at f G W ≈ 20 H z   . The radiated energy E G W ∼ 2 × 10 44 e r g   . This signal would be just out of the detectability range of the LIGO-II detector.
The cores in the simulations of Fryer and Heger [70] are not compact enough (or rotating rapidly enough) to develop bar instabilities during the collapse and initial bounce phases. However, the explosion phase ejects a good deal of low angular momentum material along the poles in their evolutions. Therefore, about 1 s after the collapse, β   becomes high enough in their models to exceed the secular bar instability limit. The β   of their model with the least angular momentum actually exceeds the dynamical bar instability limit as well (it contracts to a smaller radius and thus has a higher spin rate than the model with higher angular momentum). FHH (and [72] ) compute an upper limit (via Equation ( 2 )) to the emitted amplitude from their dynamically unstable model of h ∼ 3 × 10 − 22   (if coherent emission from a bar located at 10 Mpc persists for 100 cycles).
The corresponding frequency and maximum power are f G W ≈ 10 3 H z   and P G W = 10 53 e r g s − 1   .
LIGO-II should be able to detect such a signal (see Figure  1 , where FHH's upper limit to h   for this dynamical bar-mode is identified).
As mentioned above, the proto-neutron stars of Fryer and Heger are likely to be unstable to the development of secular bar instabilities. The GW emission from proto-neutron stars that are secularly unstable to the bar-mode has been examined by Lai and Shapiro [117, 115] . Because the timescale for secular evolution is so long, 3D hydrodynamics simulations of the nonlinear development of a secular bar can be impractical. To bypass this difficulty, Lai [115] considers only incompressible fluids, for which there are exact solutions for (Dedekind and Jacobi-like) bar development. He predicts that such a bar located at 10 Mpc would emit GWs with a peak characteristic amplitude h ∼ 10 − 21   , if the bar persists for 10 2 – 10 4   cycles. The maximum f G W   of the emitted radiation is in the range 10 2 – 10 3 H z   . This type of signal should be easily detected by LIGO-I (although detection may require a technique like the fast chirp transform method of Jenet and Prince [107] due to the complicated phase evolution of the emission).
FHH predict that a fragmentation instability is unlikely to develop during core collapse SNe because the cores have central density maxima (see also [72] ). However, they do give estimates (calculated via Equations ( 4 ,  5 )) for the amplitude, power, and frequency of the emission from such an instability: h p k ∼ 2 × 10 − 22   , P G W = 10 54 e r g s − 1   , f G W ≈ 2 × 10 3 H z   . Again, this signal would fall just beyond the upper limit of LIGO-II's frequency range. Drawing on analogy with fragmentation observed in star formation simulations, Davies et al. [49] have recently suggested that the fragmentation of collapsing massive stellar cores and the subsequent gravitational radiation-driven coalescence of the fragments may power gamma-ray bursts (see also [113] ). They say that due to kick velocity constraints, this would be possible only for supernovae that leave behind black hole remnants with M ∼ > 12 M ⊙   (formed via the merger of the fragments). Under the assumption that the collapsing core splits into two fragments (each with M = 0.7 M ⊙   ), they estimate that LIGO-II should detect the merger of ∼ 400   core collapse fragments per year. Note that the signal strength would be even greater if the total mass of the black hole remnant were indeed 12 M ⊙   (however the formation rate for these more massive objects is likely lower than that of lower mass cores).
The GW emission from r   -mode unstable neutron star remnants of core collapse SNe would be easily detectable if α m a x ∼ 1   (which is likely not physical; see Section  2.3 ). Multiple GW bursts will occur as material falls back onto the neutron star and results in repeat episodes of r   -mode growth (note that a single r   -mode episode can have multiple amplitude peaks [122] ). FHH calculate that the characteristic amplitude of the GW emission from this r   -mode evolution tracks from 6 – 1 × 10 − 22   , over a frequency range of 10 3 – 10 2 H z   (see Section  2.4 for details). They estimate the emitted energy to exceed 10 52 e r g   .
If the collapse remnant is a black hole, GWs will be emitted as the accreting black hole “rings down”. Even with very optimistic accretion scenarios, FHH conclude that such radiation will be of very low amplitude and beyond the upper frequency reach of LIGO-II (see [71] for details).
General relativity has been more fully accounted for in the core collapse studies of Dimmelmeier, Font, and Müller [51, 52, 53] , which build on the Newtonian, axisymmetric collapse simulations of Zwerger and Müller [217] . In all, they have followed the collapse evolution of 26 different models, with both Newtonian and general relativistic simulations. As in the work of Zwerger and Müller, the different models are characterized by varying degrees of differential rotation, initial rotation rates, and adiabatic indices. They use the conformally flat metric to approximate the space time geometry [48] in their relativistic hydrodynamics simulations. This approximation gives the exact solution to Einstein's equations in the case of spherical symmetry. Thus, as long as the collapse is not significantly aspherical, the approximation is relatively accurate. However, the conformally flat condition does eliminate GW emission from the spacetime. Because of this, Dimmelmeier, Font, and Müller used the quadrupole approximation to compute the characteristics of the emitted GW signal (see [217] for details).
The general relativistic simulations of Dimmelmeier et al. showed the three different types of collapse evolution (and corresponding gravitational radiation signal) seen in the Newtonian simulations of Zwerger and Müller (regular collapse – Type I signal; multiple bounce collapse – Type II signal; and rapid collapse – Type III signal). However, relativistic effects sometimes led to a different collapse type than in the Newtonian case. This is because general relativity did indeed counteract the stabilizing effects of rotation and led to much higher bounce densities (up to 700% higher). They found that multiple bounce collapse is much rarer in general relativistic simulations (occurring in only two of their models). When multiple bounce does occur, relativistic effects shorten the time interval between bounces by up to a factor of four. Movies of the simulations of four models from Dimmelmeier et al. [53] are shown in Figures  5 ,  6 ,  7 , and  8 . The four evolutions shown include a regular collapse (Movie  5 ), a rapid collapse (Movie  6 ), a multiple bounce collapse (Movie  7 ), and a very rapidly and differentially rotating collapse (Movie  8 ). The left frames of each movie contain the 2D evolution of the logarithmic density. The upper and lower right frames display the evolutions of the gravitational wave amplitude and the maximum density, respectively.
These movies can also be viewed at [132] .

Figure 5 : Movie of the evolution of the regular collapse model A3B2G4 of Dimmelmeier et al. [53] . The left frame contains the 2D evolution of the logarithmic density. The upper and lower right frames display the evolutions of the gravitational wave amplitude and the maximum density, respectively.

Figure 6 : Same as Movie  5 , but for rapid collapse model A3B2G5 of Dimmelmeier et al. [53] .

Figure 7 : Same as Movie  5 , but for multiple collapse model A2B4G1 of Dimmelmeier et al. [53] .

Figure 8 : Same as Movie  5 , but for rapid, differentially rotating collapse model A4B5G5 of Dimmelmeier et al. [53] .

Dimmelmeier et al. found that models for which the collapse type was the same in both Newtonian and relativistic simulations had lower GW amplitudes h p k   in the relativistic case. This is because the Newtonian models were less compact at bounce and thus had material with higher densities and velocities at larger radii. Both higher and lower values of h p k   were observed in models for which the collapse type changed. Overall, the range of h p k   ( 4 × 10 − 24 – 3 × 10 − 23   , for a source located at 10 Mpc) seen in the relativistic simulations was quite close to the corresponding Newtonian range. The average E G W   was somewhat higher in the relativistic case ( 1.5 × 10 47 e r g   compared to the Newtonian value of 6.4 × 10 46 e r g   ). The overall range of GW frequencies observed in their relativistic simulations (60 – 1000 Hz) was close to the Newtonian range. They did note that relativistic effects always caused the characteristic frequency of emission, f G W   , to increase (up to five-fold). For most of their models, this increase in f G W   was not accompanied by an increase in h p k   . This means that relativistic effects could decrease the detectability of GW signals from some core collapses. However, the GW emission from the models of Dimmelmeier et al. could be detected by the first generation of ground-based interferometric detectors if the sources were fortuitously located in the Local Group of galaxies. A catalog containing the signals and spectra of the GW emission from all of their models can be found at [138] .
Fully general relativistic collapse simulations (i.e., without the conformally flat approximation) have been performed by Shibata [181] . He used an axisymmetric code that solves the Einstein equations in Cartesian coordinates and the hydrodynamics equations in cylindrical coordinates.
The use of the Cartesian grid eliminates the presence of singularities and allows for stable, long-duration axisymmetric simulations [4] . The focus of this work was the effect of rotation on the criteria for prompt black hole formation. Shibata found that if the parameter q = J / M 2   is less than 0.5   , black hole formation occurred for rest masses slightly greater than the maximum mass of spherical stars. However, for 0.5 < q < 1   , the maximum stable rest mass is increased by ∼ 70 – 80 %   . The results are only weakly dependent on the initial rotation profile. Shibata did not compute the GW emission in his collapse simulations.
The new general relativistic hydrodynamics simulations of Zanotti, Rezzolla, and Font [215] suggest that a torus of neutron star matter surrounding a black hole remnant may be a stronger source of GWs than the collapse itself. They used a high resolution shock-capturing hydrodynamics method in conjunction with a static (Schwarzschild) spacetime to follow the evolution of “toroidal neutron stars”. Their results indicate that if a toroidal neutron star (with constant specific angular momentum) is perturbed, it could undergo regular oscillations. They estimate that the resulting GW emission would have a characteristic amplitude h c   ranging from 6 × 10 − 24 – 5 × 10 − 23   , for ratios of torus mass to black hole mass in the range 0.1 – 0.5. (These amplitude values are likely underestimated because the simulations of Zanotti et al. are axisymmetric.) The corresponding frequency of emission is f G W ∼ 200 H z   . The values of h c   and f G W   quoted here are for a source located at 10 Mpc. This emission would be just outside the range of LIGO-II (see Figure  1 ).
Further numerical investigations, which study tori with non-constant angular momenta and include the effects of self-gravity and black hole rotation, are needed to confirm these predictions. Movies from the simulations of Zanotti et al. can be viewed at [162] .
Magnetized tori around rapidly spinning black holes (formed via either core collapse or neutron star-black hole coalescence) have recently been examined in the theoretical study of van Putten and Levinson [204] . They find that such a torus–black hole system can exist in a suspended state of accretion if the ratio of poloidal magnetic field energy to kinetic energy E B / E k   is less than 0.1   .
They estimate that ∼ 10 %   of the spin energy of the black hole will be converted to gravitational radiation energy through multipole mass moment instabilities that develop in the torus. If a magnetized torus–black hole system located at 10 Mpc is observed for 2 × 10 4   rotation periods, the characteristic amplitude of the GW emission is ∼ 6 × 10 − 20   . It is possible that this emission could take place at several frequencies. Observations of x-ray lines from gamma-ray bursts (which are possibly produced by these types of systems) could constrain these frequencies by providing information regarding the angular velocities of the tori (preliminary estimates from observations suggest f G W ∼ 500 H z   , placing the radiation into a range detectable by LIGO-I).

3.4.5 Simulations of convective instabilities

Convectively driven inhomogeneities in the density distribution of the outer regions of the nascent neutron star and anisotropic neutrino emission are other sources of GW emission during the collapse/explosion [35, 143] . GW emission from these processes results from small-scale asphericities, unlike the large-scale motions responsible for GW emission from aspherical collapse and non-axisymmetric global instabilities. Note that Rayleigh–Taylor instabilities also induce time-dependent quadrupole moments at composition interfaces in the stellar envelope. However, the resultant GW emission is too weak to be detected because the Rayleigh–Taylor instabilities occur at very large radii [143] .
Burrows and Hayes [35] have suggested that the large proper motions observed in pulsars may result from asymmetrical core collapse SNe. They have computed the resultant GW emission under the assumption that such asymmetries were caused by convection and anisotropic neutrino emission during the collapse and explosion. Hydrodynamic simulations and theoretical investigations suggest that asymmetries present in the star prior to collapse (in part due to convection during silicon and oxygen burning) will be amplified during the collapse [16, 116] . Based on these arguments, Burrows and Hayes induced a mass dipole anisotropy of < 0.1 %   in the pre-collapse core used in their “exploratory” collapse simulation. The initial quadrupole anisotropy present in their model was 0.012   . They imposed no aspherical perturbations on the initial velocities. The collapse was followed in 1D inside a radius of 15 km and in 2D (with azimuthal symmetry) to larger radii. To reduce the computational burden of this exploratory simulation, they artificially hardened the neutrino spectrum to drive the explosion at a faster rate (note that this practice was observed to minimize the GW emission in the simulations of Müller and Janka [146] ; see below).
The initial asymmetries were observed to grow in Burrows and Hayes' calculation. The nascent neutron star developed a recoil speed of ∼ 530 k m s − 1   . Aspherical mass motion was responsible for the vast majority of the recoil. Only 16% of the recoil velocity could be attributed to anisotropic neutrino emission. The gravitational waveform from this simulation (including separate matter and neutrino contributions) is shown in Figure  9 . The contributions of matter motion and neutrino emission to the GW amplitude h   were of opposite sign. Anisotropic neutrino radiation was basically the sole contributor to h   for the first 20 ms, after which the contribution from mass motion became significant. The peak amplitude calculated was h p k ∼ 3 × 10 − 24   , for a source located at 10 Mpc. Most ( ∼ 70 %   ) of the emitted E G W ∼ 2.0 × 10 45 e r g   came from mass motion (because the mass quadrupole was changing more rapidly as the bounce occurred). The frequency of the GWs emitted was in the range 10 < f G W < 100 H z   . Burrows and Hayes predicted LIGO-II should be able to detect this emission, with a signal-to-noise ratio of 10, from a core collapse that fortuitously happened to occur in our Galaxy (with d = 10 k p c   ).

Figure 9 : The gravitational waveform (including separate matter and neutrino contributions) from the collapse simulations of Burrows and Hayes [35] . The curves plot the product of the gravitational wave amplitude and distance to the source as a function of time. (Figure 3 of [35] ; used with permission.)

The study by Nazin and Postnov [151] predicts a lower limit for E G W   emitted during an asymmetric core collapse SN (where such asymmetries could be induced by both aspherical mass motion and neutrino emission). They assume that observed pulsar kicks are solely due to asymmetric collapse. They suggest that the energy associated with the kick ( M v 2 / 2   , where M   and v   are the mass and velocity of the neutron star) can be set as a lower limit for E G W   (which can be computed without having to know the mechanism behind the asymmetric collapse). From observed pulsar proper motions, they estimate the degree of asymmetry ε   present in the collapse and the corresponding characteristic GW amplitude ( h ∝ ε   ). This amplitude is 3 × 10 − 25   for a source located at 10 Mpc and emitting at f G W = 1 k H z   . Müller and Janka performed both 2D and 3D simulations of convective instabilities in the proto-neutron star and hot bubble regions during the first second of the explosion phase of a Type II SN [146] . They numerically computed the GW emission from the convection-induced aspherical mass motion and neutrino emission in the quadrupole approximation (for details, see Section 3 of their paper).

Figure 10 : Convective instabilities inside the proto-neutron star in the 2D simulation of Müller and Janka [146] . The evolutions of the temperature (left panels) and logarithmic density (right panels) distributions are shown for the radial region 15 – 95 k m   . The upper and lower panels correspond to times 12   and 21 m s   , respectively, after the start of the simulation. The temperature values range from 2.5 × 10 10   to 1.8 × 10 11 K   . The values of the logarithm of the density range from 10.5   to 13.3 g c m − 3   . The temperature and density both increase as the colors change from blue to green, yellow, and red. (Figure 7 of [146] ; used with permission.)

For typical iron core masses, the convectively unstable region in the proto-neutron star extends over the inner 0.7 – 1.20 M ⊙   of the core mass (this corresponds to a radial range of ∼ 10 – 50 k m   ).
The convection in this region, which begins approximately 10 – 20 ms after the shock forms and may last for ≈ 20 m s – 1 s   , is caused by unstable gradients in entropy and/or lepton number resulting from the stalling of the prompt shock and deleptonization outside the neutrino sphere.
Müller and Janka's simulations of convection in this region began with the 1D, non-rotating, 12 ms post-bounce model of Hillebrandt [88] . This model included general relativistic corrections that had to be relaxed away prior to the start of the Newtonian simulations. Neutrino transport was neglected in these runs (see Section 2.1 of [146] for justification); however, a sophisticated equation of state was utilized. Figure  10 shows the evolution of the temperature and density distributions in the 2D simulation of Müller and Janka. The peak GW amplitude resulting from convective mass motions in these simulations of the proto-neutron star was ≈ 3 × 10 − 24   in 2D and ≈ 2 × 10 − 25   in 3D, for d = 10 M p c   . The emitted energy was 9.8 × 10 44 e r g   in 2D and 1.3 × 10 42 e r g   in 3D. The power spectrum peaked at frequencies of 200 – 600 Hz in 2D and 100 – 200 Hz in 3D. Such signals would not be detectable with LIGO-II. The reasons for the differences between the 2D and 3D results include smaller convective elements and less underand overshooting in 3D. The relatively low angular resolution of the 3D simulations may have also played a role. The quadrupole gravitational wave amplitude A 20 E 2   from the 2D simulation is shown in the upper left panel of Figure  11 (see [217, 194] for expressions relating A 20 E 2   to h   ).

Figure 11 : Quadrupole amplitudes A 20 E 2 [ c m ]   from convective instabilities in various models of [146] . The upper left panel is the amplitude from a 2D simulation of proto-neutron star convection. The other three panels are amplitudes from 2D simulations of hot bubble convection. The imposed neutrino flux in the hot bubble simulations increases from the top right model through the bottom right model. (Figure 18 of [146] ; used with permission.)

Convection in the hot bubble region between the shock and neutrino sphere arises because of an unstable entropy gradient resulting from neutrino heating. This unstable region extends over the inner mass range 1.25 – 1.40 M ⊙   (corresponding to a radial range of ≈ 100 – 1000 k m   ). Convection in the hot bubble begins ≈ 50 – 80 m s   after shock formation and lasts for ≈ 100 – 500 m s   .
Only 2D simulations were performed in this case. These runs started with a 25 ms post-bounce model provided by Bruenn. A simple neutrino transport scheme was used in the runs and an imposed neutrino flux was located inside the neutrino sphere. Due to computational constraints, the computational domain did not include the entire convectively unstable region inside the proto-neutron star (thus this set of simulations only accurately models the convection in the hot bubble region, not in the proto-neutron star).
The peak GW amplitude resulting from these 2D simulations of convective mass motions in the hot bubble region was h p k ≈ 5 × 10 − 25   , for d = 10 M p c   . The emitted energy was ∼ < 2 × 10 42 e r g   . The energy spectrum peaked at frequencies of 50 – 200 Hz. As the explosion energy was increased (by increasing the imposed neutrino flux), the violent convective motions turn into simple rapid expansion. The resultant frequencies drop to f G W ∼ 10 H z   . The amplitude of such a signal would be too low to be detectable with LIGO-II. Müller and Janka estimated the GW emission from the convection induced anisotropic neutrino radiation in their simulations (see [146] for details). They could not directly compute the characteristics of the GW emission from neutrinos because their proto-neutron star simulations did not include neutrino transport and the neutrino transport in their hot bubble simulations was only 1D. They found that the amplitude of the GWs emitted can be a factor of 5 – 10 higher than the GW amplitudes resulting from convective mass motion. However, their results indicate that the E G W   coming from anisotropic neutrino emission is only a small fraction of that emitted by convective mass motion.
Note that the qualitative results (including h   and E G W   magnitudes) of Burrows and Hayes' [35] “exploratory” 2D calculation are in rough agreement with the more detailed 2D simulations of Müller and Janka. However, both investigations have shortcomings. Necessary improvements would include more realistic initial models (which include rotation and are produced with multi-dimensional collapse simulations to better estimate initial asphericities), use of sophisticated neutrino transport schemes, inclusion of general relativistic effects, and higher-resolution 3D runs.

3.5 Going further

The background of GW emission from a population of core collapse SNe at cosmological distances may be detectable by LIGO-II, according to Ferrari, Matarrese, and Schneider [60] . They determined the SN rate as a function of redshift using observations to determine the evolution of the star formation rate. Only collapses that lead to black hole formation were considered. This simplified the study because the GW emission from such collapses is generally a function of just the black hole mass and angular momentum. They found that the stochastic background from these sources is not continuous and suggest that this could be used to optimize detection strategies. The maximum GW spectral strain amplitude they computed was in the range 10 − 28 – 10 − 27 H z   , at frequencies of a few times 10 2 H z   . Such a signal may be detected by a pair of LIGO-II detectors.

4 Collapse of Population III Stars

4.1 Collapse scenario

The first generation of stars to form in the early universe are known as Population III stars (formed at redshifts z ∼ > 5   ). Theoretical and computational evidence suggests that Population III stars may have had masses ∼ > 100 M ⊙   [2, 1, 75] . Since these massive stars contained no metals, it was possible for them to form directly and then evolve with very low stellar winds and thus very little mass loss. If the mass of a nonrotating Population III star is ∼ > 260 M ⊙   , its fate is to collapse directly to a black hole at the end of its life [75] . If rotational support prevents the star from direct collapse to a black hole, explosive thermonuclear burning will cause the star to undergo a giant hypernova explosion. Prior to black hole formation, the rotating, collapsed core will have a mass of 50 – 70 M ⊙   and a radius of 1000 – 2000 km. Note that because these massive stars evolve so quickly (in a few million years [9] ), the events associated with their deaths will take place at roughly the redshifts of their births.

4.2 Formation rate

The formation rate of Population III stars can be indirectly estimated from the re-ionization fraction of the early universe, which was re-ionized by light from these stars [91, 39] . Using estimates of the ultraviolet light emission from Population III stars, their ionization efficiency, and the re-ionization fraction of the early universe, one can determine that about 0.01% – 1% of the universe's baryonic matter was found in these very massive stars. This corresponds to ∼ 10 4 – 10 7   Population III stars in a 10 11 M ⊙   galaxy and thus a collapse rate that is ∼ < 10 − 3 y r − 1   . Thus, a reasonable occurrence rate can be found for an observation (luminosity) distance of ∼ 50 G p c   (which corresponds to a redshift of z = 5   , in the cosmology used by [71] ). However, uncertainties in the assumptions make this formation rate uncertain by a few orders of magnitude [1, 75, 2, 71] . It is unknown how many Population III stars will collapse directly to black holes.

4.3 GW emission mechanisms

The GW emission mechanisms related to the collapse of Population III stars are a subset of those discussed in the sections on AIC and SNe/collapsars. These include aspherical collapse, global rotational and fragmentation instabilities that may arise during the collapse/explosion and in the collapse remnant (prior to black hole formation), and the “ring-down” of the remnant black hole.

4.4 Numerical predictions of GW emission

The GW emission from the collapse of Population III stars has recently been investigated by Fryer and collaborators (Fryer, Woosley, and Heger [75] , FHH [71] , and Fryer, Holz, Hughes, and Warren [72] ). The collapse simulations of Fryer, Woosley, and Heger again started with rotating collapse progenitors that had been evolved with a stellar evolution code [84] . The initial models used by the evolution code were in rigid rotation with a surface ratio of centrifugal to gravitational forces of 20% (this ratio is seen in current observations of O stars).
The results of Fryer, Woosley, and Heger suggest that the collapse remnant (prior to black hole formation) is susceptible to the development of a secular bar-mode instability. However, at z > 5   , the GW emission would be redshifted out of LIGO-II's frequency range. At z = 5   , h p k = 8 × 10 − 23   , with a corresponding frequency of 10 Hz [71, 72] . Even if such a signal persists for a hundred cycles, it probably would be undetectable by LIGO-II. Note that these signal strengths are orders of magnitude lower than the qualitative estimates of signal strength given in Carr, Bond, and Arnett [39] .
LIGO-II may be able to detect the GW emission from binary clumps formed via a fragmentation instability. If such a signal is emitted at z = 5   and persists for 10 cycles, h   would be ∼ 10 − 22   , over a frequency range of 10 – 100 Hz [71, 72] . The likelihood of the development of a fragmentation instability is diminished by the fact that the off-center density maxima present in the simulations of Fryer, Woosley, and Heger are not very pronounced [68] .
The “ring-down” of the black hole remnant will likely be strong because Fryer, Woosley, and Heger observe a high accretion rate after collapse. FHH estimate that for a source located at z = 20   , the GWs would be redshifted out of LIGO-II's bandwidth. However, for a source at z = 5   , h p k ∼ 6 × 10 − 23   and the frequency range is 20 – 70 Hz. This signal may be marginally detectable with LIGO-II (see Figure  1 ).

5 Collapse of Supermassive Stars

5.1 Collapse scenario

There is a large body of observational evidence that supermassive black holes (SMBHs, M ∼ > 10 6 M ⊙   ) exist in the centers of many, if not most galaxies (see, e.g., the reviews of Rees [161] and Macchetto [127] ).
The masses of SMBHs in the centers of more than 45 galaxies have been estimated from observations [59] and there are more than 30 galaxies in which the presence of a SMBH has been confirmed [114] .
One of the possible formation mechanisms for SMBHs involves the gravitational collapse of supermassive stars (SMSs). The timescale for this formation channel is short enough to account for the presence of SMBHs at redshifts z > 6   [101] . Supermassive stars may contract directly out of the primordial gas, if radiation and/or magnetic field pressure prevent fragmentation [79, 55, 78, 126, 27, 1] .
Alternatively, they may build up from fragments of stellar collisions in clusters [172, 17] . Supermassive stars are radiation dominated, isentropic and convective [180, 216, 126] . Thus, they are well represented by an n = 3   polytrope. If the star's mass exceeds 10 6 M ⊙   , nuclear burning and electron/positron annihilation are not important.
After formation, an SMS will evolve through a phase of quasistationary cooling and contraction.
If the SMS is rotating when it forms, conservation of angular momentum requires that it spins up as it contracts. There are two possible evolutionary regimes for a cooling SMS. The path taken by an SMS depends on the strength of its viscosity and magnetic fields and on the nature of its angular momentum distribution.
In the first regime, viscosity or magnetic fields are strong enough to enforce uniform rotation throughout the star as it contracts. Baumgarte and Shapiro [14] have studied the evolution of a uniformly rotating SMS up to the onset of relativistic instability. They demonstrated that a uniformly rotating, cooling SMS will eventually spin up to its mass shedding limit. The mass shedding limit is encountered when matter at the star's equator rotates with the Keplerian velocity.
The limit can be represented as β s h e d = ( T / | W | ) s h e d   . In this case, β s h e d = 9 × 10 − 3   . The star will then evolve along a mass shedding sequence, losing both mass and angular momentum. It will eventually contract to the onset of relativistic instability [95, 42, 43, 180, 101] .
Baumgarte and Shapiro used both a second-order, post-Newtonian approximation and a fully general relativistic numerical code to determine that the onset of relativistic instability occurs at a ratio of R / M ∼ 450   , where R   is the star's radius and G = c = 1   in the remainder of this section.
Note that a second-order, post-Newtonian approximation was needed because rotation stabilizes the destabilizing role of nonlinear gravity at the first post-Newtonian level. If the mass of the star exceeds 10 6 M ⊙   , the star will then collapse and possibly form a SMBH. If the star is less massive, nuclear reactions may lead to explosion instead of collapse.
The major result of Baumgarte and Shapiro's work is that the universal values of the following ratios exist for the critical configuration at the onset of relativistic instability: T / | W |   , R / M   , and J / M 2   . These ratios are completely independent of the mass of the star or its prior evolution.
Because uniformly rotating SMSs will begin to collapse from a universal configuration, the subsequent collapse and the resulting gravitational waveform will be unique.
In the opposite evolutionary regime, neither viscosity nor magnetic fields are strong enough to enforce uniform rotation throughout the cooling SMS as it contracts. In this case, it has been shown that the angular momentum distribution is conserved on cylinders during contraction [21] . Because viscosity and magnetic fields are weak, there is no means of redistributing angular momentum in the star. So, even if the star starts out rotating uniformly, it cannot remain so.
The star will then rotate differentially as it cools and contracts. In this case, the subsequent evolution depends on the star's initial angular momentum distribution, which is largely unknown.
One possible outcome is that the star will spin up to mass-shedding (at a different value of β s h e d   than a uniformly rotating star) and then follow an evolutionary path that may be similar to that described by Baumgarte and Shapiro [14] . The alternative outcome is that the star will encounter the dynamical bar instability prior to reaching the mass-shedding limit. New and Shapiro [153, 154] have demonstrated that a bar-mode phase is likely to be encountered by differentially rotating SMSs with a wide range of initial angular momentum distributions. This mode will transport mass and angular momentum outward and thus may hasten the onset of collapse.

5.2 Formation rate

An estimate of the rate of the collapse of SMSs can be derived from the quasar luminosity function. Haehnelt [77] has used the quasar luminosity function to compute the rate of GW bursts from supermassive black holes, assuming that each quasar emits one such burst during its lifetime (and that each quasar is a supermassive black hole). If it is assumed that each of these bursts is due to the formation of a supermassive black hole via the collapse of a SMS, then Haehnelt's rate estimates can be used as estimates of the rate of SMS collapse. This rate is likely an overestimate of the SMS collapse rate because many SMBHs may have been formed via merger. Haehnelt predicts that the integrated event rate through redshift z = 4.5   ranges from ∼ 10 − 6 y r − 1   for M = 10 8 M ⊙   objects to ∼ 1 y r − 1   for M = 10 6 M ⊙   objects. Thus, as in the case of Population III stars, a reasonable occurrence rate can be found for an observation (luminosity) distance of 50 Gpc.

5.3 GW emission mechanisms

The GW emission mechanisms related to the collapse of SMSs are a subset of those discussed in the sections on AIC, SNe/collapsars, and Population III stellar collapse. These include (i) aspherical collapse, (ii) global rotational and fragmentation instabilities that may arise during the collapse/explosion and in the collapsed remnant (prior to black hole formation), and (iii) the “ring-down” of the remnant black hole.

5.4 Numerical predictions of GW emission

The outcome of SMS collapse can be determined only with numerical, relativistic 3D hydrodynamics simulations.
Until recently, such simulations had been published only for nearly spherical collapse. The spherical simulations of Shapiro and Teukolsky [179] produced collapse evolutions that were nearly homologous. In this case, the collapse time τ c o l l   is roughly the free-fall time at the horizon
τ c o l l = ( R 3 4 π M ) 1 / 2 = 14 s ( M 10 6 M ⊙ ) − 1 . (7)
The peak GW frequency f G W = τ c o l l − 1   is then 10 − 2 H z   , if the mass of the star is 10 6 M ⊙   . This is in the middle of LISA's frequency band of 10 − 4 – 1 H z   [196, 66] .
The amplitude h   of this burst signal can be roughly estimated in terms of the star's quadrupole moment
h ≤ ε 2 M 2 R d
≤ ε ⋅ 1 × 10 − 18 ( M 10 6 M ⊙ ) ( d 50 G p c ) − 1 . (8)
Here d   is the distance to the star and ε ∼ T / | W |   is a measure of the star's deviation from spherical symmetry. In this case, ε   will be much less than one near the horizon, since the collapse is nearly spherical.
There are two possible aspherical collapse outcomes that can be discussed. The first outcome is direct collapse to a SMBH. In this case, ε   will be on the order of one near the horizon. Thus, according to Equation  8 , the peak amplitude of the GW burst signal will be
h p k ∼ 1 × 10 − 18 ( M 10 6 M ⊙ ) ( d 50 G p c ) − 1 . (9)
Alternatively, the star may encounter the dynamical bar mode instability prior to complete collapse. Baumgarte and Shapiro [14] have estimated that a uniformly rotating SMS will reach β ∼ 0.27   when R / M = 15   . The frequency of the quasiperiodic gravitational radiation emitted by the bar can be estimated in terms of its rotation frequency to be
f G W = 2 f b a r ∼ 2 ( G M R 3 ) 1 / 2
= 2 × 10 − 3 H z ( M 10 6 M ⊙ ) − 1 , (10)
when R / M = 15   . The corresponding h p k   , again estimated in terms of the star's quadrupole moment, is
h p k ≤ 2 M 2 R d
≤ 1 × 10 − 19 ( M 10 6 M ⊙ ) ( d 50 G p c ) − 1 . (11)
The LISA sensitivity curve is shown in Figure  12 (see [93] for details on the computation of this curve; a mission time of 3 y e a r s   has been assumed). The GW signal from this dynamical bar-mode could be detected with LISA.

Figure 12 : A comparison between the GW amplitude h ( f )   for various sources and the LISA noise curve. See the text for details regarding the computations of h   . The SMS sources are assumed to be located at a luminosity distance of 50 G p c   . The bar-mode source is a dynamical bar-mode.

Shibata and Shapiro [185] have recently published a fully general relativistic, axisymmetric simulation of the collapse of a rapidly, rigidly rotating SMS. They found that the collapse remained homologous during the early part of the evolution. An apparent horizon does appear in their simulation, indicating the formation of a black hole. Because of the symmetry condition used in their run, non-axisymmetric instabilities were unable to develop.
The collapse of a uniformly rotating SMS has been investigated with post-Newtonian hydrodynamics, in 3+1 dimensions, by Saijo, Baumgarte, Shapiro, and Shibata [169] . Their numerical scheme used a post-Newtonian approximation to the Einstein equations, but solved the fully relativistic hydrodynamics equations. Their initial model was an n = 3   polytrope.
The results of Saijo et al. indicate that the collapse of a uniformly rotating SMS is coherent (i.e., no fragmentation instability develops). The collapse evolution of density contours from their model is shown in Figure  13 . Although the work of Baumgarte and Shapiro [14] suggests that a bar instability should develop prior to BH formation, no bar development was observed by Saijo et al.
They use the quadrupole approximation to estimate a mean GW amplitude from the collapse itself:
h = 4 × 10 − 21   , for a 10 6 M ⊙   star located at a distance of 50 Gpc. Their estimate for f G W   at the time of BH formation is 3 × 10 − 3 H z   . This signal would be detectable with LISA (see Figure  12 ).

Figure 13 : Meridional plane density contours from the SMS collapse simulation of Saijo, Baumgarte, Shapiro, and Shibata [169] . The contour lines denote densities ρ = ρ c × d ( 1 − i / 16 )   , where ρ c   is the central density. The frames are plotted at ( t / t D   , ρ c   , d   )=(a)( 5.0628 × 10 − 4   , 8.254 × 10 − 9   , 10 − 7   ), (b)( 2.50259   , 1.225 × 10 − 4   , 10 − 5   ), (c)( 2.05360   , 8.328 × 10 − 3   , 5.585 × 10 − 7   ), (d)( 2.50405   , 3.425 × 10 − 2   , 1.357 × 10 − 7   ), respectively. Here t   , t D   , and M 0   are the time, dynamical time ( = R e 3 / M   , where R e   is the initial equatorial radius and M   is the mass), and rest mass. (Figure 15 of [169] ; used with permission.)

Saijo et al. also consider the GW emission from the ringdown of the black hole remnant.
For the l = m = 2   quasi-normal mode of a Kerr black hole with a / M = 0.9   , they estimate the characteristic frequency and amplitude of the emission to be f G W ∼ 2 × 10 − 2 H z   and h ∼ 1 × 10 − 20 [ ( △ E G W / M ) / 10 − 4 ] 1 / 2   , for an M = 10 6 M ⊙   source located at a luminosity distance of 50 Gpc (see [118, 195, 186] for details). Here, △ E G W / M   is the radiated energy efficiency and may be ∼ < 7 × 10 − 4   [188] . This GW signal is within LISA's range of sensitivity (see Figure  12 ).

6 Summary

It is hoped that as gravitational collapse simulations become more sophisticated, the historically widely varying estimates of the magnitude of GW emission from collapse may start to converge.
Steady progress in this field has been made in the last decade. Some researchers have begun to use progenitor models produced with stellar evolution codes, which thus have more realistic angular momentum profiles, as starting points for collapse simulations [71] . This reduces the need for collapse studies that include large surveys of the angular momentum parameter space. Other progress made in the numerical study of collapse includes the use of realistic equations of state [71] , advanced neutrino transport and interaction schemes [102, 119, 159, 193] , and the performance of 3D Newtonian [160, 29, 74] and improved 2D general relativistic simulations [181] .
There is still much work to be done toward the goal of self-consistent, 3D general relativistic collapse simulations. Accurate progenitor modelling and collapse simulations must include the effects of magnetic fields, as they can significantly alter the amount of angular momentum and differential rotation present in collapsing stars. Many of the more advanced studies, which include proper microphysics treatment and/or general relativistic effects, have been limited to axisymmetry.
Full 3D simulations are necessary to compute the characteristics of the GW emission from non-axisymmetric collapse phenomena. Furthermore, simulations that follow both the collapse and the evolution of the collapsed remnant are necessary to consistently predict GW emission. One benefit of long duration simulations is that they will facilitate the investigation of the effects of the envelope on any instabilities that develop in the collapsing core or remnant. Of course, lengthy 3D simulations are computationally intensive. This burden may be reduced by the use of advanced numerical techniques, including adaptive mesh refinement and parallel algorithms.
The current numerical simulations of gravitational collapse indicate that interferometric observatories could detect GWs emitted by some collapse phenomena. LIGO-I may be able to detect GWs from secular bar-mode instabilities in core collapse SNe [115] and magnetized tori surrounding black hole collapse remnants [204] . LIGO-II could observe GWs from dynamical bar-mode instabilities in AIC [124] and core collapse SNe [71] , and possibly from the fragmentation of very massive SN cores that merge to form BHs [49] . LISA should be able to detect the collapse (and any bar-mode instabilities that develop during the collapse) of SMSs [14] and the ringdown of black hole remnants of collapsed Population III stars [71] and SMSs [169] . These observations will provide unique information about gravitational collapse and its associated progenitors and remnants.

7 Acknowledgements

It is a pleasure to thank Paul Bradley, Adam Burrows, Harald Dimmelmeier, Chris Fryer, Alex Heger, Scott Hughes, Hui Li, Ewald Müller, and Stuart Shapiro for helpful conversations and/or permission to reprint figures/movies from their published works. I also gratefully acknowledge Ewald Müller and a second referee for their beneficial reviews of this article. This work was performed under the auspices of the US Department of Energy by the Los Alamos National Laboratory under contract W-7405-ENG-36. References

  1. Abel, T., Bryan, G.L., and Norman, M.L., “The Formation and Fragmentation of Primordial Molecular Clouds”, Astrophys. J., 540, 39–44, (2000). Related online version (cited on 7 February 2000): . ☻ open access ✓
  2. Abel, T., Bryan, G.L., and Norman, M.L., “The formation of the first star in the universe”, Science, 295, 93–98, (2002). Related online version (cited on 4 December 2001): . ☻ open access ✓
  3. Akiyama, S., and Wheeler, J.C., “Magnetic Fields in Supernovae”, in Fryer, C.L., ed., Stellar Collapse, Astrophysics and Space Science Library, (Kluwer Academic Publishers, Dordrecht, Netherlands; Boston, U.S.A., 2004). Related online version (cited on 20 November 2002): . published after original LRR review. ☻ open access ✓
  4. Alcubierre, M., Brügmann, B., Holz, D.E., Takahashi, R., Brandt, S.R., Seidel, E., and Thornburg, J., “Symmetry without symmetry: Numerical simulations of axisymmetric systems using Cartesian grids”, Int. J. Mod. Phys. D, 10, 273–289, (2001). Related online version (cited on 4 August 1999): . ☻ open access ✓
  5. Andersson, N., “A New Class of Unstable Modes of Rotating Relativistic Stars”, Astrophys. J., 502, 708–713, (1998). Related online version (cited on 24 June 1997): . ☻ open access ✓
  6. Andersson, N., “Gravitational waves from instabilities in relativistic stars”, (November, 2002). URL (cited on 4 November 2002): . ☻ open access ✓
  7. Arnett, W.D., Bahcall, J.N., Kirshner, R.P., and Woosley, S.E., “Supernova 1987a”, Annu. Rev. Astron. Astrophys., 27, 629–700, (1989).
  8. Arras, P., Flanagan, É.É., Morsink, S.M., Schenk, A.K., Teukolsky, S.A., and Wasserman, I., “Saturation of the r   -mode instability”, Astrophys. J., 591, 1129–1151, (2003). published after original LRR review.
  9. Baraffe, I., Heger, A., and Woosley, S.E., “On the stability of very massive supermassive stars”, Astrophys. J., 550, 890–896, (2001). Related online version (cited on 26 September 2000): . ☻ open access ✓
  10. Bardeen, J.M., and Piran, T., “General relativistic axisymmetric rotating systems: Coordinates and equations”, Phys. Rep., 96, 205–250, (1983).
  11. Baron, E., Cooperstein, J., and Kahana, S., “Supernovae and the nuclear equation of state at high densities”, Nucl. Phys. A, 440, 744–754, (1985).
  12. Baron, E., Cooperstein, J., and Kahana, S., “Type-II Supernovae in 12 M ⊙   and 15 M ⊙   Stars: The Equation of State and General Relativity”, Phys. Rev. Lett., 55, 126–129, (1985).
  13. Baron, E., Cooperstein, J., Kahana, S., and Nomoto, K., “Collapsing white dwarfs”, Astrophys. J., 320, 304–307, (1987).
  14. Baumgarte, T.W., and Shapiro, S.L., “Evolution of rotating supermassive stars to the onset of collapse”, Astrophys. J., 526, 941–952, (1999). Related online version (cited on 14 September 1999): . ☻ open access ✓
  15. Baumgarte, T.W., and Shapiro, S.L., “Collapse of a Magnetized Star to a Black Hole”, Astrophys. J., 585, 930–947, (2003). URL (cited on 14 November 2002): . ☻ open access ✓
  16. Bazan, G., and Arnett, D., “Convection, nucleosynthesis, and core collapse”, Astrophys. J. Lett., 433, L41–L43, (1994).
  17. Begelman, M.C., and Rees, M.J., “The fate of dense stellar systems”, Mon. Not. R. Astron. Soc., 185, 847–860, (1978).
  18. Bethe, H.A., “Supernova mechanisms”, Rev. Mod. Phys., 62, 801–866, (1990).
  19. Bethe, H.A., and Wilson, J.R., “Revival of a stalled supernova shock by neutrino heating”, Astrophys. J., 295, 14–23, (1985).
  20. Blanchet, L., Damour, T., and Schäfer, G., “Post-Newtonian hydrodynamics and post-Newtonian gravitational wave generation for numerical relativity”, Mon. Not. R. Astron. Soc., 242, 289–305, (1990).
  21. Bodenheimer, P., and Ostriker, J.P., “Rapidly rotating stars. VIII. Zero-viscosity polytropic sequences”, Astrophys. J., 180, 159–170, (1973).
  22. Bonazzola, S., and Marck, J.-A., “Efficiency of gravitational radiation from axisymmetric and 3D stellar collapse I. Polytropic case”, Astron. Astrophys., 267, 623–633, (1993).
  23. Bowers, R.L., and Wilson, J.R., “Collapse of iron stellar cores”, Astrophys. J., 263, 366–376, (1982).
  24. Bowers, R.L., and Wilson, J.R., “A numerical model for stellar core collapse”, Astrophys. J. Suppl. Ser., 50, 115–159, (1982).
  25. Brachwitz, F., Dean, D.J., Hix, W.R., Iwamoto, K., Langanke, K., Martínez-Pinedo, G., Nomoto, K., Strayer, M.R., Thielemann, F.-K., and Umeda, H., “The role of electron captures in Chandrasekhar-mass models for Type Ia supernovae”, Astrophys. J., 536, 934–947, (2000). Related online version (cited on 26 January 2000): . ☻ open access ✓
  26. Bravo, E., and García-Senz, D., “Coulomb corrections to the equation of state of nuclear statistical equilibrium matter: Implications for SNIa nucleosynthesis and the accretion-induced collapse of white dwarfs”, Mon. Not. R. Astron. Soc., 307, 984–992, (1999).
  27. Bromm, V., Coppi, P.S., and Larson, R.B., “Forming the First Stars in the Universe: The Fragmentation of Primordial Gas”, Astrophys. J. Lett., 527, L5–L8, (1999). Related online version (cited on 13 October 1999): . ☻ open access ✓
  28. Brown, J.D., “Gravitational waves from the dynamical bar instability in a rapidly rotating star”, Phys. Rev. D, 62, 084024–1–11, (2000). Related online version (cited on 1 April 2000): . ☻ open access ✓
  29. Brown, J.D., “Rotational instabilities in post-collapse stellar cores”, in Centrella, J.M., ed., Astrophysical Sources for Ground-Based Gravitational Wave Detectors, Philadelphia, Pennsylvania, 30 October 1 November 2000, vol. 575 of AIP Conference Proceedings, 234–245, (American Institute of Physics, Melville, U.S.A., 2001).
  30. Bruenn, S.W., “Numerical Simulations of Core Collapse Supernovae”, in Guidry, M.W., and Strayer, M.R., eds., Nuclear Physics in the Universe, Proceedings of the First Symposium on Nuclear Physics in the Universe held in Oak Ridge, Tennessee, USA 24–26 September 1992, 31–50, (Institute of Physics, Bristol, U.K.; Philadelphia, U.S.A., 1993).
  31. Bruenn, S.W., De Nisco, K.R., and Mezzacappa, A., “General relativistic effects in the core collapse supernova mechanism”, Astrophys. J., 560, 326–338, (2001). Related online version (cited on 23 January 2001): . ☻ open access ✓
  32. Burrows, A., “Supernova explosions in the Universe”, Nature, 403, 727–733, (2000).
  33. Burrows, A., and Fryxell, B.A., “An instability in neutron stars at birth”, Science, 258, 430–434, (1992).
  34. Burrows, A., and Fryxell, B.A., “A convective trigger for supernova explosions”, Astrophys. J. Lett., 418, L33–L35, (1993).
  35. Burrows, A., and Hayes, J., “Pulsar recoil and gravitational radiation due to asymmetrical stellar collapse and explosion”, Phys. Rev. Lett., 76, 352–355, (1996).
  36. Burrows, A., and Thompson, T.A., “Neutrino-Matter Interaction Rates in Supernovae”, in Fryer, C.L., ed., Stellar Collapse, Astrophysics and Space Science Library, (Kluwer Academic Publishers, Dordrecht, Netherlands; Boston, U.S.A., 2004). Related online version (cited on 18 November 2002): . published after original LRR review. ☻ open access ✓
  37. California Institute of Technology, “LIGO Laboratory Home Page”, project homepage. URL (cited on 4 October 2002): . ☻ open access ✓
  38. Cappellaro, E., Evans, R., and Turatto, M., “A new determination of supernova rates and a comparison with indicators for galactic star formation”, Astron. Astrophys., 351, 459–466, (1999).
  39. Carr, B.J., Bond, J.R., and Arnett, W.D., “Cosmological consequences of population III stars”, Astrophys. J., 277, 445–469, (1984).
  40. Cassisi, S., Iben Jr, I., and Tornambé, A., “Hydrogen-accreting carbon-oxygen white dwarfs”, Astrophys. J., 496, 376–385, (1998).
  41. Centrella, J.M., New, K.C.B., Lowe, L., and Brown, J.D., “Dynamical rotational instability at low T / W   ”, Astrophys. J. Lett., 550, L193–L196, (2001). Related online version (cited on 27 October 2000): . ☻ open access ✓
  42. Chandrasekhar, S., “Dynamical Instability of Gaseous Masses Approaching the Schwarzschild Limit in General Relativity”, Phys. Rev. Lett., 12, 114–116, (1964). Erratum 437–438.
  43. Chandrasekhar, S., “The dynamical instability of gaseous masses approaching the Schwarzschild limit in general relativity”, Astrophys. J., 140, 417–433, (1964).
  44. Chandrasekhar, S., An Introduction to the Study of Stellar Structure, (Dover, New York, U.S.A., 1967).
  45. Colgate, S.A., “Supernovae: Hot bubbles drive explosions”, Nature, 341, 489–490, (1989).
  46. Colgate, S.A., Herant, M., and Benz, W., “Neutron star accretion and the neutrino fireball”, Phys. Rep., 227, 157–174, (1993).
  47. Colgate, S.A., and White, R.H., “The hydrodynamic behavior of supernovae explosions”, Astrophys. J., 143, 626–681, (1966).
  48. Cook, G.B., Shapiro, S.L., and Teukolsky, S.A., “Testing a simplified version of Einstein's equations for numerical relativity”, Phys. Rev. D, 53(10), 5533–5540, (May, 1996). Related online version (cited on 5 December 1995): . ☻ open access ✓
  49. Davies, M.B., King, A.R., Rosswog, S., and Wynn, G., “Gamma-ray bursts, supernova kicks, and gravitational radiation”, Astrophys. J. Lett., 579, L63–L66, (2002). Related online version (cited on 22 April 2022): . ☻ open access ✓
  50. Detweiler, S.L., and Lindblom, L., “On the evolution of the homogeneous ellipsoidal figures. II – Gravitational collapse and gravitational radiation”, Astrophys. J., 250, 739–749, (1981).
  51. Dimmelmeier, H., Font, J.A., and Müller, E., “Gravitational waves from relativistic rotational core collapse”, Astrophys. J. Lett., 560, L163–L166, (2001). Related online version (cited on 21 April 2001): . ☻ open access ✓
  52. Dimmelmeier, H., Font, J.A., and Müller, E., “Relativistic simulations of rotational core collapse. I. Methods, initial models and code tests”, Astron. Astrophys., 388, 917–935, (2002). Related online version (cited on 17 April 2002): . ☻ open access ✓
  53. Dimmelmeier, H., Font, J.A., and Müller, E., “Relativistic simulations of rotational core collapse II. Collapse dynamics and gravitational radiation”, Astron. Astrophys., 393, 523–542, (2002). Related online version (cited on 17 April 2002): . ☻ open access ✓
  54. Durisen, R.H., and Tohline, J.E., “Fission of rapidly rotating fluid systems”, in Black, D., and Matthews, M., eds., Protostars and Planets II, 534–575, (University of Arizona Press, Tucson, U.S.A., 1985).
  55. Eisenstein, D.J., and Loeb, A., “Origin of quasar progenitors from the collapse of low-spin cosmological perturbations”, Astrophys. J., 443, 11–17, (1995). Related online version (cited on 11 January 1994): . ☻ open access ✓
  56. Epstein, R., The post-Newtonian theory of the generation of gravitational radiation and its application to stellar collapse, Ph.D. Thesis, (Stanford University, Stanford, U.S.A., 1976).
  57. Epstein, R., and Wagoner, R.V., “Post-Newtonian generation of gravitational waves”, Astrophys. J., 197, 717–723, (1975).
  58. Eriguchi, Y., and Müller, E., “Equilibrium models of differentially rotating polytropes and the collapse of rotating stellar cores”, Astron. Astrophys., 147, 161–168, (1985).
  59. Ferrarese, L., and Merritt, D., “A fundamental relation between supermassive black holes and their host galaxies”, Astrophys. J., 539, L9–L12, (2000). Related online version (cited on 4 June 2000): . ☻ open access ✓
  60. Ferrari, V., Matarrese, S., and Schneider, R., “Gravitational wave background from a cosmological population of core-collapse supernovae”, Mon. Not. R. Astron. Soc., 303, 247–257, (1999). Related online version (cited on 15 October 1998): . ☻ open access ✓
  61. Ferrari, V., Miniutti, G., and Pons, J.A., “Gravitational waves from newly born, hot neutron stars”, (October, 2002). URL (cited on 26 October 2002): . ☻ open access ✓
  62. Finn, L.S., “Supernovae, gravitational radiation, and the quadrupole formula”, in Evans, C.R., Finn, L.S., and Hobill, D.W., eds., Frontiers in Numerical Relativity, International workshop devoted to research in numerical relativity, held in Urbana-Champaign in May 1988, 126–145, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1989).
  63. Finn, L.S., “Detectability of gravitational radiation from stellar-core collapse”, in Buchler, J.R., Detweiler, S.L., and Ipser, J.R., eds., Nonlinear problems in relativity and cosmology, 6th Florida Workshop in Nonlinear Astronomy, held on October 2-4, 1990 in Gainesville, Florida, 156–172, (New York Academy of Sciences, New York, U.S.A., 1991).
  64. Finn, L.S., “Gravitional Radiation Sources and Signatures”, in Dixon, L.J., ed., Gravity: From the Hubble Length to the Planck Length, Proceedings of the 26th SLAC Summer Institute on Particle Physics (SSI 98), Stanford, USA, 3-14 August 1998, vol. 538 of SLAC-R, (SLAC, Springfield, U.S.A., 2001). URL (cited on 30 March 1999): . also at: http://arXiv.org/abs/gr-qc/9903107. ☻ open access ✓
  65. Finn, L.S., and Evans, C.R., “Determining gravitational radiation from Newtonian self-gravitating systems”, Astrophys. J., 351, 588–600, (1990).
  66. Folkner, W.M., ed., Laser Interferometer Space Antenna: Second International LISA Symposium on the Detection and Observation of Gravitational Waves in Space, Pasadena 1998, vol. 456 of AIP Conference Proceedings, (Springer, New York, U.S.A., 1993).
  67. Friedman, J.L., and Morsink, S.M., “Axial Instability of Rotating Relativistic Stars”, Astrophys. J., 502, 714–720, (1998). Related online version (cited on 23 June 1997): . ☻ open access ✓
  68. Fryer, C.L., personal communication, (2002).
  69. Fryer, C.L., Benz, W., Herant, M., and Colgate, S.A., “What can the accretion-induced collapse of white dwarfs really explain?”, Astrophys. J., 516, 892–899, (1999). Related online version (cited on 2 December 1998): . ☻ open access ✓
  70. Fryer, C.L., and Heger, A., “Core-collapse simulations of rotating stars”, Astrophys. J., 541, 1033–1050, (2000). Related online version (cited on 30 July 1999): . ☻ open access ✓
  71. Fryer, C.L., Holz, D.E., and Hughes, S.A., “Gravitational wave emission from core collapse of massive stars”, Astrophys. J., 565, 430–446, (2002). Related online version (cited on 25 September 2001): . ☻ open access ✓
  72. Fryer, C.L., Holz, D.E., Hughes, S.A., and Warren, M.S., “Stellar Collapse and Gravitational Waves”, in Fryer, C.L., ed., Stellar Collapse, Astrophysics and Space Science Library, (Kluwer Academic Publishers, Dordrecht, Netherlands; Boston, U.S.A., 2004). URL (cited on 27 November 2002): . published after original LRR review. ☻ open access ✓
  73. Fryer, C.L., and Kalogera, V., “Theoretical black hole mass distributions”, Astrophys. J., 554, 548–560, (2001).
  74. Fryer, C.L., and Warren, M.S., “Modeling Core-Collapse Supernovae in Three Dimensions”, Astrophys. J., 574, L65–L68, (2002). URL (cited on 3 June 2002): . ☻ open access ✓
  75. Fryer, C.L., Woosley, S.E., and Heger, A., “Pair-instability supernovae, gravity waves, and gamma-ray transients”, Astrophys. J., 550, 372–382, (2001). Related online version (cited on 13 July 2000): . ☻ open access ✓
  76. Hachisu, I., “A versatile method for obtaining structures of rapidly rotating stars”, Astrophys. J. Suppl. Ser., 61, 479–507, (1986).
  77. Haehnelt, M.G., “Low-frequency gravitational waves from supermassive black holes”, Mon. Not. R. Astron. Soc., 269, 199–208, (1994).
  78. Haehnelt, M.G., Natarajan, P., and Rees, M.J., “High-redshift galaxies, their active nuclei and central black holes”, Mon. Not. R. Astron. Soc., 300, 817–827, (1998). Related online version (cited on 18 December 1997): . ☻ open access ✓
  79. Haehnelt, M.G., and Rees, M.J., “The formation of nuclei in newly formed galaxies and the evolution of the quasar population”, Mon. Not. R. Astron. Soc., 263, 168–178, (1993).
  80. Haensel, P., Levenfish, K.P., and Yakovlev, D.G., “Bulk viscosity in superfluid neutron star cores. III. Effects of Σ −   hyperons”, Astron. Astrophys., 381, 1080–1089, (2002). Related online version (cited on 26 October 2001): . ☻ open access ✓
  81. Hamuy, M., “Observed and physical properties of core-collapse supernovae”, (September, 2002). URL (cited on 10 September 2002): . ☻ open access ✓
  82. Hayashi, A., Eriguchi, Y., and Hashimoto, M., “On the possibility of the nonexplosive core contraction of massive stars: New evolutionary paths from rotating white dwarfs to rotating neutron stars”, Astrophys. J., 492, 286–297, (1998).
  83. Hayashi, A., Eriguchi, Y., and Hashimoto, M., “On the possibility of the nonexplosive core contraction of massive stars. II. General relativistic analysis”, Astrophys. J., 521, 376–381, (1999).
  84. Heger, A., The presupernova evolution of rotating massive stars, Ph.D. Thesis, (Technische Universität München, Munich, Germany, 1998).
  85. Heger, A., personal communication, (2002).
  86. Heger, A., Langer, N., and Woosley, S.E., “Presupernova Evolution of Rotating Massive Stars. I. Numerical Method and Evolution of the Internal Stellar Structure”, Astrophys. J., 528, 368–396, (2000). Related online version (cited on 12 April 1999): . ☻ open access ✓
  87. Heger, A., and Woosley, S.E., “On the Progenitors of Collapsars”, in Ricker, G.R., and Vanderspek, R.K., eds., Gamma-Ray Burst and Afterglow Astronomy 2001: A Workshop Celebrating the First Year of the HETE Mission, Woods Hole, Massachusetts, USA, 5-9 November 2001, vol. 662 of AIP Conference Proceedings, 214–216, (American Institute of Physics, Melville, U.S.A., 2003). URL (cited on 1 June 2002): . published after original LRR review. ☻ open access ✓
  88. Hillebrandt, W., “Stellar Collapse and Supernova Explosions”, in Pacini, F., ed., High Energy Phenomena around Collapsed Stars, Proceedings of the NATO Advanced Study Institute, Cargèse, Corsica, France, September 2-13, 1985, vol. 195 of NATO Science Series, 73–104, (Reidel, Dordrecht, Netherlands; Boston, U.S.A., 1987).
  89. Ho, W.C.G., and Lai, D., “ r   -Mode oscillations and spin-down of young rotating magnetic neutron stars”, Astrophys. J., 543, 386–394, (2000). Related online version (cited on 15 December 1999): . ☻ open access ✓
  90. Höflich, P., Khokhlov, A., Wang, L., Wheeler, J.C., and Baade, D., “Aspherical Supernovae Explosions”, in van der Hucht, K.A., Herrero, A., and Esteban, C., eds., A Massive Star Odyssey, from Main Sequence to Supernova, Lanzarote, Canary Islands, Spain, June 24-28, 2002, vol. 212 of IAU Symposia, (Astronomical Society of the Pacific, San Francisco, U.S.A., 2003). Related online version (cited on 12 July 2002): . ☻ open access ✓
  91. Hogan, C.J., “A model of pregalactic evolution”, Mon. Not. R. Astron. Soc., 188, 781–790, (1979).
  92. Houser, J.L., “The effect of rotation on the gravitational radiation and dynamical instability of stiff stellar cores”, Mon. Not. R. Astron. Soc., 299, 1069–1086, (1998).
  93. Hughes, S.A., “Untangling the merger history of massive black holes with LISA”, Mon. Not. R. Astron. Soc., 331, 805–816, (2002). Related online version (cited on 30 August 2001): . ☻ open access ✓
  94. Hughes, S.A., Márka, S., Bender, P.L., and Hogan, C.J., “New physics and astronomy with the new gravitational-wave observatories”, in Graf, N., ed., Proceedings of Snowmass 2001, The Future of Particle Physics, 30 June 21 July 2001, Snowmass Village, Colorado, vol. C010630, P402, (SLAC eConf, Stanford, U.S.A., 2001). URL (cited on 15 October 2001): . ☻ open access ✓
  95. Iben, I., “Massive stars in quasi-static equilibrium”, Astrophys. J., 138, 1090–1096, (1963).
  96. Imamura, J.N., and Durisen, R.H., “The Dominance of Dynamic Barlike Instabilities in the Evolution of a Massive Stellar Core Collapse That “Fizzles””, Astrophys. J., 549, 1062–1075, (2001).
  97. Imamura, J.N., Friedman, J.L., and Durisen, R.H., “Secular stability limits for rotating polytropic stars”, Astrophys. J., 294, 474–478, (1985).
  98. Isern, J., Canal, R., and Labay, J., “The outcome of explosive ignition of ONeMg cores: supernovae, neutron stars, or “iron” white dwarfs?”, Astrophys. J. Lett., 372, L83–L86, (1991).
  99. Istituto Nazionale di Fisica Nucleare, “The Virgo Project”, project homepage. URL (cited on 4 October 2002): . ☻ open access ✓
  100. Janka, H.-T., “Conditions for shock revival by neutrino heating in core-collapse supernovae”, Astron. Astrophys., 368, 527–560, (2001). Related online version (cited on 28 August 2000): . ☻ open access ✓
  101. Janka, H.-T., “Supermassive Stars: Fact or Fiction?”, in Chui, C.K., Siuniaev, R.A., and Churazov, E., eds., Lighthouses of the Universe: The Most Luminous Celestial Objects and Their Use for Cosmology, Proceedings of the MPA/ESO/MPE/USM Joint Astronomy Conference, held in Garching, Germany, 6-10 August 2001, ESO Astrophysics Symposia, 357–368, (Springer, Berlin, Germany; New York, U.S.A.; et al., 2002). Related online version (cited on 1 February 2002): . ☻ open access ✓
  102. Janka, H.-T., Buras, R., Kifonidis, K., Rampp, M., and Plewa, T., “Explosion Mechanisms of Massive Stars”, in Fryer, C.L., ed., Stellar Collapse, Astrophysics and Space Science Library, (Kluwer Academic Publishers, Dordrecht, Netherlands; Boston, U.S.A., 2004). Related online version (cited on 13 December 2002): . published after original LRR review. ☻ open access ✓
  103. Janka, H.-T., Kifonidis, K., and Rampp, M., “Supernova explosions and neutron star formation”, in Blaschke, D., Glendenning, N. K., and Sedrakian, A., eds., Physics of Neutron Star Interiors, Proceedings of a Workshop on Physics of Neutron Star Interiors (ECT, Trento, June 19 July 7, 2000), vol. 578 of Lecture Notes in Physics, 333–363, (Springer, Berlin, Germany; New York, U.S.A., 2001). Related online version (cited on 1 March 2001): . ☻ open access ✓
  104. Janka, H.-T., and Müller, E., “Neutrino-driven type-II supernovae: Neutrino heating and post bounce dynamics”, in Suzuki, Y., and Nakamura, K., eds., Frontiers of Neutrino Astrophysics, Proceedings of the International Symposium on Neutrino Astrophysics held on October 19-22, 1992, Takayama / Kamioka, Japan, vol. 5 of Frontiers Science Series, 203–217, (Universal Academy Press, Tokyo, Japan, 1993).
  105. Janka, H.-T., and Müller, E., “Dynamics of Type-II supernovae”, in McCray, R., and Wang, Z., eds., Supernovae and Supernovae Remnants, Proceedings of the IAU Colloquium 145, held in Xian, China, May 24-29, 1993, IAU Colloquium 145, 109–118, (Cambridge University Press, Cambridge, U.K., 1996).
  106. Janka, H.-T., and Müller, E., “Neutrino heating, convection, and the mechanism of Type-II supernova explosions”, Astron. Astrophys., 306, 167–198, (1996).
  107. Jenet, F.A., and Prince, T.A., “Detection of variable frequency signals using a fast chirp transform”, Phys. Rev. D, 62, 122001–1–10, (2000). Related online version (cited on 7 December 2000): . ☻ open access ✓
  108. Jet Propulsion Laboratory, “Laser Interferometer Space Antenna Home Page”, project homepage, (October, 2002). URL (cited on 4 October 2002): . ☻ open access ✓
  109. Jones, P.B., “Bulk viscosity of neutron-star matter”, Phys. Rev. D, 64, 084003–1–7, (2001).
  110. Jones, P.B., “Comment on “Gravitational Radiation Instability in Hot Young Neutron Stars””, Phys. Rev. Lett., 86, 1384, (2001).
  111. Kato, M., and Hachisu, I., “A new estimation of mass accumulation efficiency in helium shell flashes toward Type Ia supernova explosions”, Astrophys. J. Lett., 513, L41–L44, (1999). Related online version (cited on 8 January 1999): . ☻ open access ✓
  112. King, A.R., and Lasota, J.-P., “Spin evolution and magnetic fields in cataclysmic variables”, Astrophys. J., 378, 674–681, (1991).
  113. Kobayashi, S., and Mészáros, P., “Gravitational radiation from gamma-ray burst progenitors”, (October, 2002). URL (cited on 9 October 2002): . ☻ open access ✓
  114. Kormendy, J., “Supermassive Black Holes in Disk Galaxies”, in Funes, J.G., and Corsini, E.M., eds., Galaxy Disks and Disk Galaxies, Proceedings of a conference sponsored by the Vatican Observatory, held at the Pontifical Gregorian University in Rome, Italy, 12-16 June 2000, vol. 230 of ASP Conference Series, 247–256, (Astronomical Society of the Pacific, San Francisco, U.S.A., 2001).
  115. Lai, D., “Secular bar-mode evolution and gravitational waves from neutron stars”, in Centrella, J.M., ed., Astrophysical Sources for Ground-based Gravitational Wave Detectors, Philadelphia, Pennsylvania, USA, 30 October 1 November 2000, vol. 575 of AIP Conference Proceedings, 246–257, (American Institute of Physics, Melville, U.S.A., 2001).
  116. Lai, D., and Goldreich, P., “Growth of perturbations in gravitational collapse and accretion”, Astrophys. J., 535, 402–411, (2000). Related online version (cited on 25 June 1999): . ☻ open access ✓
  117. Lai, D., and Shapiro, S.L., “Gravitational radiation from rapidly rotating nascent neutron stars”, Astrophys. J., 442, 259–272, (1995). Related online version (cited on 17 August 1994): . ☻ open access ✓
  118. Leaver, E.W., “An analytic representation for the quasi-normal modes of Kerr black holes”, Proc. R. Soc. London, Ser. A, 402, 285–298, (1985).
  119. Liebendörfer, M., Messer, O.E.B., Mezzacappa, A., Bruenn, S.W., Cardall, C.Y., and Thielemann, F.-K., “A finite difference representation of neutrino radiation hydrodynamics for spherically symmetric general relativistic supernova simulations”, (July, 2002). URL (cited on 1 July 2002): . ☻ open access ✓
  120. Lindblom, L., and Owen, B.J., “Effect of hyperon bulk viscosity on neutron-star r-modes”, Phys. Rev. D, 65, 063006–1–15, (2002). Related online version (cited on 25 October 2001): . ☻ open access ✓
  121. Lindblom, L., Owen, B.J., and Morsink, S.M., “Gravitational Radiation Instability in Hot Young Neutron Stars”, Phys. Rev. Lett., 80, 4843–4846, (1998). Related online version (cited on 13 March 1998): . ☻ open access ✓
  122. Lindblom, L., Tohline, J.E., and Vallisneri, M., “Nonlinear Evolution of the r-Modes in Neutron Stars”, Phys. Rev. Lett., 86, 1152–1155, (2001). Related online version (cited on 31 October 2000): . ☻ open access ✓
  123. Lindblom, L., Tohline, J.E., and Vallisneri, M., “Numerical evolutions of nonlinear r-modes in neutron stars”, Phys. Rev. D, 65, 084039–1–15, (2002). Related online version (cited on 20 September 2001): . ☻ open access ✓
  124. Liu, Y.T., “Dynamical instability of new-born neutron stars as sources of gravitational radiation”, Phys. Rev. D, 65, 124003–1–14, (2002). Related online version (cited on 21 September 2001): . ☻ open access ✓
  125. Liu, Y.T., and Lindblom, L., “Models of rapidly rotating neutron stars: remnants of accretion-induced collapse”, Mon. Not. R. Astron. Soc., 324, 1063–1073, (2001). Related online version (cited on 9 December 2000): . ☻ open access ✓
  126. Loeb, A., and Rasio, F.A., “Collapse of primordial gas clouds and the formation of quasar black holes”, Astrophys. J., 432, 52–61, (1994). Related online version (cited on 16 January 1994): . ☻ open access ✓
  127. Macchetto, F.D., “Supermassive black holes and galaxy morphology”, Astrophys. Space Sci., 269, 269–291, (1999). Related online version (cited on 5 October 1999): . ☻ open access ✓
  128. MacFadyen, A.I., and Woosley, S.E., “Collapsars: Gamma-Ray Bursts and Explosions in “Failed Supernovae””, Astrophys. J., 524, 262–289, (1999).
  129. Maiolino, R., Vanzi, L., Mannucci, F., Cresci, G., Ghinassi, F., and Della Valle, M., “Discovery of two infrared supernovae: a new window on the SN search”, (April, 2002). URL (cited on 6 April 2002): . ☻ open access ✓
  130. Managan, R.A., “On the secular instability of axisymmetric rotating stars to gravitational radiation reaction”, Astrophys. J., 294, 463–473, (1985).
  131. Marck, J.-A., and Bonazzola, S., “Gravitational radiation from three-dimensional gravitational stellar core collapse”, in D'Inverno, R., ed., Approaches to Numerical Relativity, Proceedings of the International Workshop on Numerical Relativity, Southampton, England, 16-20 December 1991, 247, (Cambridge University Press, Cambridge, U.K., 1992).
  132. Max-Planck-Institut für Astrophysik, “General relativistic simulations of rotational supernova collapse”, project homepage, (July, 2002). URL (cited on 3 July 2002): . ☻ open access ✓
  133. Mayle, R., and Wilson, J.R., “Supernovae from collapse of oxygen-magnesium-neon cores”, Astrophys. J., 334, 909–926, (1988).
  134. Mezzacappa, A., Calder, A.C., Bruenn, S.W., Blondin, J.M., Guidry, M.W., Strayer, M.R., and Umar, A.S., “The interplay between proto-neutron star convection and neutrino transport in core collapse supernovae”, Astrophys. J., 493, 848–862, (1998). Related online version (cited on 18 September 1997): . ☻ open access ✓
  135. Mezzacappa, A., Calder, A.C., Bruenn, S.W., Blondin, J.M., Guidry, M.W., Strayer, M.R., and Umar, A.S., “An investigation of neutrino-driven convection and the core collapse supernova mechanism using multigroup neutrino transport”, Astrophys. J., 495, 911–926, (1998). Related online version (cited on 18 September 1997): . ☻ open access ✓
  136. Misner, C.W., Thorne, K.S., and Wheeler, J.A., Gravitation, (W.H. Freeman, San Francisco, U.S.A., 1973).
  137. Mönchmeyer, R., Schäfer, G., Müller, E., and Kates, R.E., “Gravitational waves from the collapse of rotating stellar cores”, Astron. Astrophys., 246, 417–440, (1991).
  138. Müller, E., “MPA Hydro Gang Homepage”, project homepage, Max Planck Institute for Astrophysics. URL (cited on 17 April 2002): . old url at cited date: http://www.mpa-garching.mpg.de/Hydro/hydro.html. ☻ open access ✓
  139. Müller, E., “Gravitational radiation from collapsing rotating stellar cores”, Astron. Astrophys., 114, 53–59, (1982).
  140. Müller, E., in Hillebrandt, W., and Müller, E., eds., Proceedings of the 7th Workshop on Nuclear Astrophysics (Ringberg Castle, Tegernsee, Germany, March 22-27, 1993), 27, (Max-Planck-Institut für Astrophysik, Garching, Germany, 1993).
  141. Müller, E., “Multidimensional hydrodynamical simulations of supernova explosions”, in Bludman, S.A., Mochkovitch, R., and Zinn-Justin, J., eds., Supernovae, Proceedings of the Les Houches Summer School, Session LIV, 31 July 1 September 1990, vol. 54 of Les Houches Summer School Proceedings, 395, (Elsevier, Amsterdam, Netherlands, 1994).
  142. Müller, E., “Gravitational waves from core collapse supernovae”, in Marck, J.-A., and Lasota, J.-P., eds., Relativistic Gravitation and Gravitational Radiation, Proceedings of the Les Houches School of Physics, held in Les Houches, Haute Savoie, 26 September 6 October, 1995, Cambridge Contemporary Astrophysics, 273–308, (Cambridge University Press, Cambridge, U.K., 1997).
  143. Müller, E., “Simulation of astrophysical fluid flow”, in Steiner, O., and Gautschy, A., eds., Computational Methods for Astrophysical Fluid Flow, Lecture Notes 1997 of the Swiss Society for Astronomy and Astrophysics (SSAA), held March 3–8, 1997 in Les Diablerets, Switzerland, vol. 27 of Saas-Fee Advanced Course, 343–494, (Springer, Berlin, Germany; New York, U.S.A., 1998).
  144. Müller, E., and Hillebrandt, W., “The collapse of rotating stellar cores”, Astron. Astrophys., 103, 358–366, (1981).
  145. Müller, E., Hillebrandt, W., and Rozyczka, M., “Stellar Collapse – Adiabatic hydrodynamics and shock wave propagation”, Astron. Astrophys., 81, 288–292, (1980).
  146. Müller, E., and Janka, H.-T., “Gravitational radiation from convective instabilities in Type II supernova explosions”, Astron. Astrophys., 317, 140–163, (1997).
  147. Nadyozhin, D.K., “Physical Properties of SNe IIP Derived from a Comparison of Theoretical Models with Observations”, Astron. Astrophys., submitted, (2002).
  148. Nakamura, T., “General Relativistic Collapse of Axially Symmetric Stars Leading to the Formation of Rotating Black Holes”, Prog. Theor. Phys., 65, 1876–1890, (1981).
  149. Nakamura, T., “General Relativistic Collapse of Accreting Neutron Stars with Rotation”, Prog. Theor. Phys., 70, 1144–1147, (1983).
  150. National Astronomical Observatory of Japan, “TAMA Project”, project homepage. URL (cited on 4 October 2002): . ☻ open access ✓
  151. Nazin, S.N., and Postnov, K.A., “High neutron star birth velocities and gravitational radiation during supernova explosions”, Astron. Astrophys., 317, L79–L81, (1997). Related online version (cited on 15 January 1997): . ☻ open access ✓
  152. New, K.C.B., Centrella, J.M., and Tohline, J.E., “Gravitational waves from long-duration simulations of the dynamical bar instability”, Phys. Rev. D, 62, 064019–1–16, (2000). Related online version (cited on 30 November 1999): . ☻ open access ✓
  153. New, K.C.B., and Shapiro, S.L., “Evolution of differentially rotating supermassive stars to the onset of bar instability”, Astrophys. J., 548, 439–446, (2001). Related online version (cited on 9 October 2000): . ☻ open access ✓
  154. New, K.C.B., and Shapiro, S.L., “The formation of supermassive black holes and the evolution of supermassive stars”, Class. Quantum Grav., 18, 3965–3976, (2001). Related online version (cited on 6 September 2000): . ☻ open access ✓
  155. Nomoto, K., and Kondo, Y., “Conditions for accretion-induced collapse of white dwarfs”, Astrophys. J. Lett., 367, L19–L22, (1991).
  156. Novikov, I.D., “Gravitational radiation from a star collapsing into a disk”, Sov. Astron., 19, 398, (1975).
  157. Pickett, B.K., Durisen, R.H., and Davis, G., “The dynamic stability of rotating protostars and protostellar disks. I. The effects of the angular momentum distribution”, Astrophys. J., 458, 714–738, (1996).
  158. Piran, T., and Stark, R.F., “Numerical relativity, rotating gravitational collapse and gravitational radiation”, in Centrella, J.M., ed., Dynamical Spacetimes and Numerical Relativity, Proceedings of the workshop held at Drexel University on October 7–11, 1985, 40–73, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1986).
  159. Rampp, M., and Janka, H.-T., “Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations”, Astron. Astrophys., 396, 361–392, (2002). Related online version (cited on 7 March 2002): . ☻ open access ✓
  160. Rampp, M., Müller, E., and Ruffert, M., “Simulations of non-axisymmetric rotational core collapse”, Astron. Astrophys., 332, 969–983, (1998). Related online version (cited on 11 November 1997): . ☻ open access ✓
  161. Rees, M.J., “Astrophysical Evidence for Black Holes”, in Wald, R.M., ed., Black Holes and Relativistic Stars, Proceedings of the Symposium dedicated to the memory of Subrahmanyan Chandrasekhar, held in Chicago, December 14–15, 1996, 79–101, (University of Chicago Press, Chicago, U.S.A.; London, U.K., 1998).
  162. Rezzolla, L., “Relativistic Astrophysics movies at SISSA”, personal homepage, SISSA / ISAS. URL (cited on 4 October 2002): . ☻ open access ✓
  163. Rezzolla, L., Lamb, F.K., Marković, D., and Shapiro, S.L., “Properties of r   modes in rotating magnetic neutron stars. I. Kinematic secular effects and magnetic evolution equations”, Phys. Rev. D, 64, 104013–1–12, (2001). Related online version (cited on 17 July 2001): . ☻ open access ✓
  164. Rezzolla, L., Lamb, F.K., Marković, D., and Shapiro, S.L., “Properties of r modes in rotating magnetic neutron stars. II. Evolution of the r modes and stellar magnetic field”, Phys. Rev. D, 64, 104014–1–13, (2001). Related online version (cited on 17 July 2001): . ☻ open access ✓
  165. Ruffini, R., and Wheeler, J.A., “Relativistic Cosmology from Space Platforms”, in Hardy, V., and H. Moore, H., eds., Proceedings of the Conference on Space Physics, 45–174, (ESRO, Paris, France, 1971).
  166. Saenz, R.A., and Shapiro, S.L., “Gravitational radiation from stellar collapse – Ellipsoidal models”, Astrophys. J., 221, 286–303, (1978).
  167. Saenz, R.A., and Shapiro, S.L., “Gravitational and neutrion radiation from stellar core collapse – Improved ellipsoidal model calculations”, Astrophys. J., 229, 1107–1125, (1979).
  168. Saenz, R.A., and Shapiro, S.L., “Gravitational radiation from stellar core collapse. III – Damped ellipsoidal oscillations”, Astrophys. J., 244, 1033–1038, (1981).
  169. Saijo, M., Baumgarte, T.W., Shapiro, S.L., and Shibata, M., “Collapse of a rotating supermassive star to a supermassive black hole: Post-Newtonian simulations”, Astrophys. J., 569, 349–361, (2002). Related online version (cited on 6 February 2002): . ☻ open access ✓
  170. Saijo, M., Shibata, M., Baumgarte, T.W., and Shapiro, S.L., “Dynamical Bar Instability in Rotating Stars: Effect of General Relativity”, Astrophys. J., 548, 919–931, (2001). Related online version (cited on 10 October 2000): . ☻ open access ✓
  171. Salpeter, E.E., “Energy and pressure of a zero-temperature plasma”, Astrophys. J., 134, 669–682, (1961).
  172. Sanders, R.H., “The effects of stellar collisions in dense stellar systems”, Astrophys. J., 162, 791–809, (1970).
  173. Schenk, A.K., Arras, P., Flanagan, É.É., Teukolsky, S.A., and Wasserman, I., “Nonlinear mode coupling in rotating stars and the r   -mode instability in neutron stars”, Phys. Rev. D, 65, 024001–1–43, (2002). Related online version (cited on 23 January 2001): . ☻ open access ✓
  174. Schutz, B.F., “Gravitational wave astronomy”, Class. Quantum Grav., 16, A131–A156, (1999). Related online version (cited on 9 November 1999): . ☻ open access ✓
  175. Seidel, E., and Moore, T., “Gravitational radiation from realistic relativistic stars: Odd-parity fluid perturbations”, Phys. Rev. D, 35, 2287–2296, (1987).
  176. Seidel, E., and Moore, T., “Gravitational radiation from perturbations of stellar core collapse models”, in Evans, C.R., Finn, L.S., and Hobill, D.W., eds., Frontiers in Numerical Relativity, 146–162, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1988).
  177. Shapiro, S.L., “Gravitational radiation from stellar collapse – The initial burst”, Astrophys. J., 214, 566–575, (1977).
  178. Shapiro, S.L., and Lightman, A.P., “Rapidly rotating, post-Newtonian neutron stars”, Astrophys. J., 207, 263–278, (1976).
  179. Shapiro, S.L., and Teukolsky, S.A., “Gravitational collapse of supermassive stars to black holes – Numerical solution of the Einstein equations”, Astrophys. J., 234, L177–L181, (1979).
  180. Shapiro, S.L., and Teukolsky, S.A., Black Holes, White Dwarfs and Neutron Stars, (John Wiley & Sons, New York, U.S.A., 1983).
  181. Shibata, M., “Axisymmetric Simulations of Rotating Stellar Collapse in Full General Relativity – Criteria for Prompt Collapse to Black Holes –”, Prog. Theor. Phys., 104, 325–358, (2000). Related online version (cited on 19 July 2000): . ☻ open access ✓
  182. Shibata, M., Baumgarte, T.W., and Shapiro, S.L., “The bar-mode instability in differentially rotating neutron stars: Simulations in full general relativity”, Astrophys. J., 542, 453–463, (2000). Related online version (cited on 18 May 2000): . ☻ open access ✓
  183. Shibata, M., Baumgarte, T.W., and Shapiro, S.L., “Stability and collapse of rapidly rotating, supramassive neutron stars: 3D simulations in general relativity”, Phys. Rev. D, 61, 044012–1–11, (2000). Related online version (cited on 16 November 1999): . ☻ open access ✓
  184. Shibata, M., Karino, S., and Eriguchi, Y., “Dynamical instability of differentially rotating stars”, Mon. Not. R. Astron. Soc., 334, L27–L32, (2002). URL (cited on 6 June 2002): . ☻ open access ✓
  185. Shibata, M., and Shapiro, S.L., “Collapse of a rotating supermassive star to a supermassive black hole: Fully relativistic simulations”, Astrophys. J. Lett., 572, L39–L43, (2002).
  186. Shibata, M., Shapiro, S.L., and Uryu, K., “Equilibrium and stability of supermassive stars in binary systems”, Phys. Rev. D, 64, 24004–1–14, (2001). Related online version (cited on 25 April 2001): . ☻ open access ✓
  187. Smith, S., Houser, J.L., and Centrella, J.M., “Simulations of nonaxisymmetric instability in a rotating star: A comparison between Eulerian and Smooth Particle Hydrodynamics”, Astrophys. J., 458, 236–256, (1996). Related online version (cited on 9 October 1995): . ☻ open access ✓
  188. Stark, R.F., and Piran, T., “Gravitational-Wave Emission from Rotating Gravitational Collapse”, Phys. Rev. Lett., 55, 891–894, (1985).
  189. Stergioulas, N., and Font, J.A., “Nonlinear r   -modes in rapidly rotating relativistic stars”, Phys. Rev. Lett., 86, 1148–1151, (2001). Related online version (cited on 31 July 2000): . ☻ open access ✓
  190. Swesty, F.D., Lattimer, J.M., and Myra, E.S., “The role of the equation of state in the `prompt' phase of type II supernovae”, Astrophys. J., 425, 195–204, (1994).
  191. Symbalisty, E.M.D., “Magnetorotational iron core collapse”, Astrophys. J., 285, 729–746, (1984).
  192. Tassoul, J.-L., Theory of Rotating Stars, (Princeton University Press, Princeton, U.S.A., 1978).
  193. Thompson, T.A., Burrows, A., and Pinto, P.A., “Shock breakout in core-collapse supernovae and its neutrino signature”, (November, 2002). URL (cited on 10 November 2002): . ☻ open access ✓
  194. Thorne, K.S., “Multipole expansions of gravitational radiation”, Rev. Mod. Phys., 52, 299–340, (1980).
  195. Thorne, K.S., “Gravitational radiation”, in Hawking, S.W., and Israel, W., eds., Three Hundred Years of Gravitation, 330–458, (Cambridge University Press, Cambridge, U.K.; New York, U.S.A., 1987).
  196. Thorne, K.S., “Gravitational radiation”, in Böhringer, H., Morfill, G.E., and Trümper, J.E., eds., 17th Texas Symposium on Relativistic Astrophysics and Cosmology, vol. 759 of Annals of the New York Academy of Sciences, 127–152, (New York Academy of Sciences, New York, U.S.A., 1995).
  197. Thorne, K.S., “Gravitational Waves from Compact Bodies”, in van Paradijs, J., van den Heuvel, E.P.J., and Kuulkers, E., eds., Compact Stars in Binaries, Proceedings of the 165th Symposium of the International Astronomical Union, held in The Hague, the Netherlands, August 15-19, 1994, vol. 165 of IAU Symposia, 153–183, (Kluwer, Dordrecht, Netherlands; Boston, U.S.A., 1996). Related online version (cited on 30 June 1995): . ☻ open access ✓
  198. Thuan, T.X., and Ostriker, J.P., “Gravitational radiation from stellar collapse”, Astrophys. J. Lett., 191, L105–L107, (1974).
  199. Tohline, J.E., “The collapse of rotating stellar cores – Equilbria between white dwarf and neutron star densities”, Astrophys. J., 285, 721–728, (1984).
  200. Tohline, J.E., and Hachisu, I., “The breakup of self-gravitating rings, tori, and thick accretion disks”, Astrophys. J., 361, 394–407, (1990).
  201. Toman, J., Imamura, J.N., Pickett, B.K., and Durisen, R.H., “Nonaxisymmetric dynamic instabilities of rotating polytropes. I. The Kelvin modes”, Astrophys. J., 497, 370–387, (1998).
  202. Turner, M.S., and Wagoner, R.V., “Gravitational radiation from slowly-rotating supernovae – Preliminary results”, in Smarr, L., ed., Sources of Gravitational Radiation, Proceedings of the Battelle Seattle Workshop, July 24-August 4, 1978, 383–407, (Cambridge University Press, Cambridge, U.K., 1979).
  203. University of Hannover, “GEO 600 Home Page”, project homepage. URL (cited on 4 October 2002): . ☻ open access ✓
  204. van Putten, M.H.P.M., and Levinson, A., “Theory and astrophysical consequences of a magnetized torus around a rapidly rotating black hole”, Astrophys. J., 584, 937–953, (2003). URL (cited on 12 December 2002): . ☻ open access ✓
  205. Wheeler, J.C., Yi, I., Höflich, P., and Wang, L., “Asymmetric Supernovae, Pulsars, Magnetars, and Gamma-Ray Bursts”, Astrophys. J., 537, 810–823, (2000).
  206. Wickramasinghe, D.T., and Ferrario, L., “Magnetism in Isolated and Binary White Dwarfs”, Publ. Astron. Soc. Pac., 112, 873–924, (2000).
  207. Wilson, J.R., in Centrella, J.M., LeBlanc, J.M., and Bowers, J.L., eds., Numerical Astrophysics, Proceedings of a symposium in honor of James R. Wilson, held at the University of Illinois in October, 1982, 422–434, (Jones and Barlett, Boston, U.S.A., 1985).
  208. Wilson, J.R., and Mayle, R., “Convection in core collapse supernovae”, Phys. Rep., 163, 63–78, (1988).
  209. Wilson, J.R., Mayle, R., Woosley, S.E., and Weaver, T.A., “Stellar Core Collapse and Supernovae”, Ann. N.Y. Acad. Sci., 470, 267–293, (1986).
  210. Woodward, J., Tohline, J.E., and Hachisu, I., “The stability of thick, self-gravitating disks in protostellar systems”, Astrophys. J., 420, 247–267, (1994).
  211. Woosley, S.E., and Baron, E., “The collapse of white dwarfs to neutron stars”, Astrophys. J., 391, 228–235, (1992).
  212. Woosley, S.E., Zhang, W., and Heger, A., “Gamma-Ray Bursts and Jet-Powered Supernovae”, in Hillebrandt, W., and Leibundgut, B., eds., From Twilight to Highlight: The Physics of Supernovae, Proceedings of the ESO/MPA/MPE Workshop Held at Garching, Germany, 29-31 July 2002, ESO Astrophysics Symposia, 87–96, (Springer, Berlin, Germany; New York, U.S.A.; et al., 2003). Related online version (cited on 4 November 2002): . ☻ open access ✓
  213. Yamada, S., and Sato, K., “Gravitational Radiation from Rotational Collapse of a Supernova Core”, Astrophys. J., 450, 245–252, (1995).
  214. Yungelson, L.R., and Livio, M., “Type Ia Supernovae: An Examination of Potential Progenitors and the Redshift Distribution”, Astrophys. J., 497, 168–177, (1998). Related online version (cited on 18 November 1997): . ☻ open access ✓
  215. Zanotti, O., Rezzolla, L., and Font, J.A., “Quasi-periodic accretion and gravitational waves from oscillating “toroidal neutron stars” around a Schwarzschild black hole”, Mon. Not. R. Astron. Soc., 341, 832–848, (2003).
  216. Zel'dovich, Y.B., and Novikov, I.D., Relativistic Astrophysics, vol. 1, (University of Chicago Press, Chicago, U.S.A., 1971).
  217. Zwerger, T., and Müller, E., “Dynamics and gravitational wave signature of axisymmetric rotational core collapse”, Astron. Astrophys., 320, 209–227, (1997).

Note: The reference version of this article is published by Living Reviews in Relativity