Gravitational Waves from Gravitational Collapse

Gravitational wave emission from stellar collapse 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. This review covers the entire range of stellar collapse sources of gravitational waves: from the accretion induced collapse of a white dwarf through the collapse down to neutron stars or black holes of massive stars to the collapse of supermassive stars. Electronic Supplementary Material Supplementary material is available for this article at 10.12942/lrr-2003-2.


Introduction
The field of gravitational wave (GW) astronomy will soon become a reality. The first generation of ground-based interferometric detectors (LIGO [45], VIRGO [126], GEO 600 [104], TAMA 300 [182]) are beginning their search for GWs. Toward the end of this decade, two of these detectors (LIGO, VIRGO) will begin upgrades that should allow them to reach the sensitivities necessary to regularly detect emissions from astrophysical sources. A space-based interferometric detector, LISA [181], 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 covers an entire spectrum of stellar masses, from the accretion induced collapse (AIC) of a white dwarf through the collapse of massive stars (M > 8 M ) including the "collapsar" engine believed to power long-duration gamma-ray bursts [261], 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 [244,172,173,73,214,91,88,120]. 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. Asymmetric neutrino emission can also produce a strong gravitational wave signature. 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 or if the disks around the black hole develop instabilities. 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. In contrast, gravitational waves can propagate from the innermost parts of the stellar core to detectors without attenuation by intervening matter. With their weak interaction cross-sections, neutrinos can probe the same region probed by GWs. But whereas neutrinos are extremely sensitive to details in the microphysics (equation of state and cross-sections), GWs are most sensitive to physics driving the mass motions (e.g., rotation). Combined, the neutrino and the GW signals can teach us much about the conditions in the collapsing core and ultimately the physics that governs stellar collapse (e.g., [7,87]).
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., [203,245,204,57,179,171,233,74,170,271,198,86,88]). 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, Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 • 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, collapsar models in Section 4, collapse of Population III stars in Section 5, and collapse of SMSs in Section 6). Each of these sections (2,3,4,5,6) 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 Hz) and LISA (for sources with lower frequencies in the range of 10 −4 to 1 Hz).

Collapse scenario
Stars with masses below 8 M end their lives ejecting their envelopes in a possible planetary nebula, leaving behind a white dwarf that gradually cools and fades away. For those white dwarfs in binaries, binary accretion can reheat the white dwarf. If the accreted material ignites degenerately, the resultant nuclear explosion produces a nova and ejects all of the accreted material. But if this material burns non-degenerately, the white dwarf will gain mass. When the mass of the white dwarf exceeds the Chandrasekhar stability limit, it will begin to collapse.
Two possible fates await this collapsing white dwarf. If the collapsing core can achieve high enough temperatures, nuclear burning will begin and the stellar pressure will increase. This nuclear burning can drive nature's largest nuclear bomb, producing type Ia supernovae. However, neutrino emission from electron capture (e.g., URCA processes) can damp this burning until the core has collapsed too deeply into its gravitational well for nuclear burning to turn around the collapse. The result of this collapse is the formation of a neutron star. Electron capture will dominate if the density at which nuclear ignition occurs exceeds a critical density ρ crit . For C-O white dwarfs, ρ crit is in the range 6 × 10 9 -10 10 g cm −3 [32]. 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) [187,127]. The collapse of an O-Ne-Mg white dwarf begins when its central density reaches 4 × 10 9 g cm −3 . (For more details about the conditions under which AIC occurs, see [187,127,32,31,156].) The dynamics of the collapse itself are somewhat similar to the dynamics of core collapse SNe: the collapse proceeds until the core reaches nuclear densities, the core then bounces and sends out a bounce shock that stalls when it becomes optically thin to neutrinos and loses its thermal energy. However, there is very little envelope around this core to prevent an explosion and the shock can easily be revived to drive a low-mass explosion. Exactly how much matter is ejected depends upon the details of the collapse calculation (compare [114,263,82]). Less than 10 −1 M will likely be ejected by the star (due to the bounce itself or due to neutrino absorption/wind mechanisms) [82].

Formation rate
The AIC occurrence rate is difficult to determine for a number of reasons. These include incomplete understanding of binary star evolution and determining how much matter is truly accreted onto the white dwarf versus the amount that is ejected through novae [48,137,234]. 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). Figure 1 shows one estimate of the region in the space of initial white dwarf mass and accretion rate that produces AICs. New results continue to alter the dividing lines between these fates [234].
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. The exact yield is sensitive to the neutron fraction in this ejecta, which depends both on the neutrino transport and the electron capture rates, but if all of these isotopes present in the Galaxy are assumed to have originated in AICs, an upper limit of ∼ 10 −5 yr −1 can be set for the Galactic AIC rate [82].
Binary population synthesis analysis can be used to determine which accreting white dwarfs will undergo AIC. The results of Yungelson and Livio [268] predict that the galactic AIC rate is between 8 × 10 −7 and 8 × 10 −5 yr −1 . Thus, a reasonable occurrence rate can be found for an observation distance of 100 Mpc.

GW emission mechanisms
During AIC, emission of GWs will occur if the infall of matter is aspherical. The convective time is likely to be brief, so it is unlikely that post-bounce convection will produce GWs. However, GWs will also be produced if the collapsing star or neutron star remnant develops rotational or pulsational instabilities [226,227,235,252,254]. These include global rotational mode, r-mode, f -mode, and fragmentation instabilities.
Global rotational instabilities in fluids arise from non-axisymmetric modes e ±imφ , where m = 2 is known as the "bar-mode" [239,6]. It is convenient to parameterize a system's susceptibility to these modes by the stability parameter β = T rot /|W |. Here, T rot 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 [221,63,162,125,192,118,248]. If the object has an off-center density maximum, β d could be as low as 0.10 [247,260,192,49]. General relativity may enhance the dynamical bar-mode instability by slightly reducing β d [223,209]. 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 (e.g., [225]. The secular bar-mode instability limit for MacLaurin spheroids is β s ≈ 0.14.
In an attempt to reduce these high rotation requirements, there has been increasing work studying bar-mode instabilities driven by dynamical sheer in differentially rotating neutron stars. Sheer instabilities excite the co-rotating f -mode. If viscous forces don't damp this instability altogether, it is possible that this instability can occur for β d -values as low as ∼ 0.01 for stars with a large degree of differential rotation.
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 In rotating stars, gravitational radiation reaction drives the r-modes toward unstable growth [5,77]. In hot, rapidly rotating neutron stars, this instability may not be suppressed by internal dissipative mechanisms (such as viscosity and magnetic fields) [152]. 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 α max to which the amplitude is allowed to grow [153,154]. Originally, it was thought that α max ∼ 1. Later work indicated that α max may be ≥ 3 [153,236,213,154]. 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 [136,135,201,202,154,151,102,6] (significant uncertainties remain regarding the efficacy of these dissipative mechanisms). Furthermore, a study of a simple barotropic neutron star model by Arras et al. [9] suggests that multimode couplings could limit α max to values 1. If α max is indeed 1 (see also [97]), 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 α max ∼ 1) is presented in the remainder of this paper. However, because it is quite doubtful that α max 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 [96]. If this does indeed occur, the orbiting fragments would be a strong GW source.

Full collapse simulations
The accretion-induced collapse of white dwarfs has been simulated by a number of groups [14,166,263,82]. 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. [82], 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) [86] examined an AIC simulation (model 3) of Fryer et al. [82]. FHH used both numerical and analytical techniques to estimate the peak amplitude h pk , energy E GW , and frequency f GW 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 [169]. 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 where i, j = 1, 2, 3 are spatial indices and r = x 2 + y 2 + z 2 is the distance to the source.  [170,271] for details). Note that FHH define h pk 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 [82] 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 = dE GW /dt of the GWs emitted by a star that has encountered the bar-mode instability are and Here m, 2r, 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 2r) 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 and Here m is the mass of a single fragment, 2r 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 [115] (which assumes α max = 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 km, respectively, the resulting formula for the average GW strain is where α is the mode amplitude and ν s is the spin frequency.
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 FHH's numerical quadrupole estimate of the GWs from polar oscillations in the AIC simulation of Fryer et al. [82] predicts a peak dimensionless amplitude h pk = 5.9 × 10 −24 (for d = 100 Mpc). The energy E GW = 3 × 10 45 erg is emitted at a frequency of f GW ≈ 50 Hz. 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 2 (see Appendix A of [86] for details on the computation of this curve). Note that the characteristic strain h is plotted along the vertical axis in Figure 2 (and in LISA's sensitivity curve, shown in Figure 18). For burst sources, h = h pk . For sources that persist for N cycles, h = √ N h pk . The simulation of Fryer et al. [82] 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), 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 cm 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 [139]. Higher angular momenta are certainly possible. 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 α max ∼ 1 (which is likely not physical; see Section 2.3), they predict E GW > 10 52 erg. FHH compute h(f GW ) 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 Mpc, this track is always below the LIGO-II noise curve. The point on this track with the maximum h, which corresponds to the beginning Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 of r-mode evolution, is shown in Figure 2. The track moves down and to the left ( i.e., h and f decrease) in this figure as the r-mode evolution continues.

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., [16,185,156]). Others use individual equilibrium models as initial conditions for hydrodynamical simulations (e.g., [231,192,184,49]). 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 nonaxisymmetric modes and estimation of the characteristics of any resulting GW emission.
Liu and Lindblom [156,155] have applied this equilibrium approach to AIC. Their investigation began with a study of equilibrium models built to represent neutron stars formed from AIC [156]. These neutron star models were created via a two-step process, using a Newtonian version of Hachisu's self-consistent field method [98]. 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 cm −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 cm −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 [82]. The realistic equation of state used to construct the white dwarfs is a Coulomb corrected, zero temperature, degenerate gas equation of state [210,52].
In the second step of their process, Liu and Lindblom [156] 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 [156] 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 [271]. 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 cm −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 [155] used linearized hydrodynamics to perform a stability analysis of the cold neutron star AIC remnants of Liu and Lindblom [156]. He found that only the remnant of the O-Ne-Mg white dwarf (Liu and Lindblom's Model III) Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 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 [247,260,192,49], 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 pk ≈ 1.4 × 10 −24 and the LIGO-II signal-to-noise ratio (for a persistent signal like that seen in the work of [184] and [34]) to be S/N ≤ 3 (for f GW ≈ 450 Hz). These values are for a source located at 100 Mpc. He also predicted that the timescale for gravitational radiation to carry away enough angular momentum to eliminate the bar-mode is τ GW ∼ 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 2). Details of the approximations on which these estimates are based can be found in [155].
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 [256].

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 [224]. 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.

Collapse scenario
Stars with mass greater than ∼ 8 -10 M will undergo core collapse at the end of their thermonuclear burning life cycles. The gravitational energy released during the collapse is believed to be the power source behind a large subset of supernovae (for reviews, see [8,23,79,81,143]). 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 [108,266]; those that undergo mass loss via mass transfer in binaries can have progenitor masses as low as 12 M [194].
If the stellar mass is less than ∼ 20 -25 M [78,89,86,108], it is believed that the star will produce a strong supernova, leaving behind a neutron star remnant. This section studies the gravitational waves produced during the stellar collapse and supernova explosion mechanism. Without mass loss, stars above 40 -50 M are thought to collapse directly to black holes without producing a supernova explosion [78,108]. If rotating, the collapsing stars may form an accretion disk around their black hole core. One of the leading models for long-duration gamma-ray bursts (the "collapsar" engine) argues that the energy extracted from the disk or the black hole spin can drive a relativistic jet [262]. We will discuss these collapsars and their resultant gravitational waves further in Section 4. However, if we include the effects of winds, these stars lose much of their mass through winds during their nuclear-burning lifetimes. Their fate will be closer to that of a 20 M and are likely to produce supernovae and, at least initially, neutron stars, not black holes [89,108]. If the progenitor's mass is in the range 20 − 25 M 40 -50 M , the entire star is not ejected in the SN explosion. More than 2 M will fall back onto the nascent neutron star and lead to black hole formation. These objects also produce a variant of the collapsar engine for gamma-ray bursts (GRBs). Although the discussion of GW production from the collapse and supernova explosion phase will be discussed in this section, the GWs produced during the fallback and black hole accretion disk phase will be discussed in Section 4. 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 [103,178].
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 [71]. 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 [71]. 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 [71,170]. The maximum velocity of the outer regions of the core reaches ∼ 7 × 10 4 km s −1 . It takes just 1 s for an earth-sized core to collapse to a radius of 50 km [8].
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 [13,12]. 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 erg (supplied by the binding energy of the nascent neutron star), then the remainder of the star can be ejected Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 in about 20 ms [8]. This so-called "prompt explosion" mechanism has been succesful numerically only when 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 [23,173]. Inclusion of general relativistic effects in collapse simulations can increase the success of the prompt mechanism in some cases [13,237].
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 [55] suggested that energy from neutrinos emitted by the collapsed core could revive the stalled shock. (See Burrows and Thompson [44] for a 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 [30,29,257,24,259,258]. However, their 1-dimensional simulations and those of others had difficulty producing energies high enough to match observations [53,36,131]. It was not until Herant and collaborators modelled the collapse and bounce phase with neutrinotransport in 2-dimensions that convection began to be accepted as a necessary puzzle piece in understanding the supernova explosion mechanism [112].
Observations of SN 1987a show that significant mixing occurred during this supernovae (see Arnett et al. [8] 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 [173]. This has been confirmed by both 2D and 3D simulations [112,42,132,173,167,78,83,128,130,197,91,85,26,39,92,80,212,253]. 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 [173]. The "hot bubble" region above the surface of the neutron star also has been shown to be convectively unstable [53,23,54,91]. Janka and Müller have demonstrated that convection in this region only aids the explosion if the neutrino luminosity is in a narrow region [133]. 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 [168,130,197,39] and this problem is far from solved. But progress not only in neutrino transport, but in understanding new features in the convection [26,43,75] some including the effect of magnetic fields is leading to new ideas about the supernova mechanism [3]. We will discuss the effects of this new physics on the GW signal at the end of this section.
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 [8,3,116,122,121] 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 [130]. If due to the mechanism itself, these asymmetries may provide clues into the true engine behind supernova explosions. Already, the observations partly motivated the multi-dimensional studies of convection in the delayed explosion mechanism as well as work on magnetic field engines [17,255,3]. Observations have also driven the work on jets and collapsars 4. Höflich et al. [116] have argued that low velocity jets stalled inside SN envelopes can account for the observed asymmetries. Hungerford and collaborators have argued that the asymmetries required are not so extreme [122,121], arguments that have now been confirmed [138]. This debate is crucial to our understanding of the supernova mechanism. If jets are required, magnetic field mechanisms are the likely source of the asymmetry. If jets are not required, the convective engine can produce asymmetries via a number of channels from low mode convection [26,212,43] to rotation(e.g., [83]) to asymmetric collapse [41,80]. Most GW emission calculations of stellar collapse have focused on the collapse of rotating, massive stars with rotation Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 periods at least as high as those assumed by Fryer & Heger [83], implicitly assuming that the asymmetries are driving by rotation.

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. [46]. Their sample includes 137 SNe from five different SN searches. They determined that the core collapse SN rate in the Galaxy is 6 × 10 −3 -1.6 × 10 −2 yr −1 . Thus, a reasonable occurrence rate can be found for an observation distance of 10 Mpc. (Note that a infrared survey estimates that the rate may actually be an order of magnitude higher [161].) What this rate might exclude are those stars that collapse to form neutron stars, but have sufficiently weak explosions that they eject very little nickel and hence are not very bright. Including these systems as well will increase the total number of stellar collapses initially into neutron stars by less than 5 -40% [89]. However, it does include those stars that collapse to black holes directly and form GRBs (producing supernova-like outbursts). We will discuss the GWs from such collapses in Section 4. Excluding them from this sample will not change the above numbers.

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 ij will increase. The contraction speeds up over the next 20 ms and the density distribution becomes a centrally condensed torus [170]. In this phase the core's shrinking size dominates its increasing deformation and the magnitude of I ij decreases. As the core bounces, I ij changes rapidly due to the deceleration and rebound. If the bounce occurs because of nuclear pressure, its timescale will be < 1 ms. If centrifugal forces play a role in halting the collapse, the bounce can last up to several ms [170]. The magnitude of I ij 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 ij 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 [71,170]. 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 (see [41,176,80,158] for reviews).
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.

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 [203] identified mechanisms related to core collapse that could produce GWs and provided order-of-magnitude estimates of the characteristics of such emission.
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 Quantitative computations of GW emission during the infall phase of collapse were performed by Thuan and Ostriker [245] and Epstein and Wagoner [66,65], 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 [188] and Shapiro and Saenz [218,204] 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 [249] and later by Saenz and Shapiro [205,206].
Müller [171] calculated the quadrupole GW emission from 2D axisymmetric collapse based on the Newtonian simulations of Müller and Hillebrandt [174] (these simulations used a realistic equation of state and included differential rotation). He found that differential rotation enhanced the efficiency of the GW emission.
Stark and Piran [233,193] were the first to compute the GW emission from fully relativistic collapse simulations, using the ground-breaking formalism of Bardeen and Piran [11]. 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 [215,216]. 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 [57]. 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 [72]. In the remainder of Section 3.4, more recent investigations will be discussed.

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 [170]. 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 [219] and Tohline [246] (in the context of the "fizzler" scenario for failed supernovae; see also [105,106,124]), and had been observed in earlier collapse simulations [175,238]. Mönchmeyer and collaborators found that a bounce caused by centrifugal forces would last for several ms, whereas a bounce at nuclear densities would occur in < 1 ms. 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 Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 the numerical quadrupole approximation discussed in Section 2.4). The waveforms they categorized as Type I (similar to those observed in previous collapse simulations [171,74]) 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 E2 20 for a Type I waveform is shown in Figure 3 (see [241,271] for expressions relating A E2 20 to h). The waveforms identified as Type II exhibit several maxima, which result from multiple bounces (see Figure 4 for an example of a Type II waveform). Note that the waveforms displayed in Figures 3 and 4 are from the study of Zwerger and Müller [271], discussed below. The model of Mönchmeyer et al. that bounced due to nuclear forces had the highest GW amplitude of all of their models, h pk ∼ 10 −23 for a source distance d = 10 Mpc, and the largest emitted energy E GW ∼ 10 47 erg. The accompanying power spectrum peaked in the frequency range 5 × 10 2 -10 3 Hz.
The most extensive Newtonian survey of the parameter space of axisymmetric, rotational core collapse is that of Zwerger and Müller [271]. 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 [67]. 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 Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2  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. [170]) 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 5). Type III signals were not observed in the evolution of strongly differentially rotating Γ r = 1.28 models and were also not seen in subsequent investigations [144,190]. Their waveforms were computed with the same technique used in [170].
In contrast to the results of Mönchmeyer et al. [170], 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 pk in models that bounced at subnuclear densities. The effect of the initial value of β i on h pk was non-monotonic. For models with β i 0.1, h pk increased with increasing β i . This is because the deformation of the core is larger for faster rotators. However, for models with larger β i , h pk decreases as β i increases. These models bounce at subnuclear Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 densities. Thus, the resulting acceleration at bounce and the GW amplitude are smaller. Zwerger and Müller found that the maximum value of h pk for a given sequence was reached when ρ c at bounce was just less than ρ nuc . 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 Mpc (the model with the highest h is identified in Figure 2). The corresponding energies ranged from 10 44 E GW 10 47 erg. 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. Magnetic fields lower these amplitudes by ∼ 10% [145], but realistic stellar profiles can lower the amplitudes by a factor of 4 − 10 [190,177], restricting the detectability of supernovae to within our Galaxy ( 10 kpc).
Yamada and Sato [264] used techniques very similar to those of Zwerger and Müller [271] in their core collapse study. Their investigation revealed that the h pk for Type I signals became saturated when the dimensionless angular momentum of the collapsing core, q = J/(2GM/c), reached ∼ 0.5. They also found that h pk was sensitive to the stiffness of the equation of state for densities just below ρ nuc . The characteristics of the GW emission from their models were similar to those of Zwerger and Müller.

Non-axisymmetric simulations
The GW emission from non-axisymmetric hydrodynamics simulations of stellar collapse was first studied by Bonazzola and Marck [163,28]. 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 pk 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 [198]. These authors started their Newtonian simulations with the only model (A4B5G5) of Zwerger and Müller [271] 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 ms after bounce.
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 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. [49] 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 [198], 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 [271]). 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 km encountered this instability, the values of h pk from their models would be ∼ 2 × 10 −24 -2 × 10 −23 , for d = 10 Mpc. The frequency at which h pk occurred in their simulations was ∼ 200 Hz. This instability would have to persist for at least ∼ 15 cycles to be detected with LIGO-II.
Brown [35] carried out an investigation of the growth of non-axisymmetric modes in postbounce cores that was similar in many respects to that of Rampp, Müller, and Ruffert [198]. 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 [109]. 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 [109]. The model evolved by Rampp, Müller, and Ruffert [198] 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 [198].
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 [198]. 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 [49]. 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 [198]).
Brown's simulations determined that both his most rapidly rotating model Ω24 (with postbounce β > 0.35) and model RMR are unstable to growth of the m = 2 bar-mode. However, his model Ω20 (with post-bounce β > 0. 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 [198], that only the innermost portion of the core (with ρ > 10 10 g cm −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 β ic = 0.15. Brown also observed that the β of the inner core does not have to exceed 0.27 Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 for the model to encounter the bar-mode. Models Ω24 and RMR had β ic ≈ 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 [91] 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 Mpc [88]. In later work, Fryer & Warren [92] included full Newtonian gravity through a tree algorithm and studied the rotating progenitors from Fryer & Heger [83]. By the launch of the explosion, no bar instabilities had developed. This was because of several effects: they used slowly rotating, but presumably realistic, progenitors [83], the explosion occured quickly for their models ( 100 ms) and, finally, because much of the high angular momentum material did not make it into the inner core. These models have been further studied for the GW signals [87]. The fastest rotating models achieved a signal of h ∼ 2 × 10 −24 for d = 10 Mpc and characteristic frequencies of f GW ∼ 1000 Hz. For supernovae occuring within the Galaxy, such a signal is detectable by LIGO-II.
Fryer and collaborators have also modeled asymmetric collapse and asymmetric explosion calculations in 3 dimensions [80,90]. These calculations will be discussed in Section 3.4.5.
The GW emission from nonradial quasinormal mode oscillations in proto-neutron stars has been examined by Ferrari, Miniutti, and Pons [70]. They found that the frequencies of emission f GW 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-tonoise ratio of 5, E GW 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 [60]).

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 [239,271,198,37].
The full collapse simulations of Fryer and Heger [83] are the most sophisticated axisymmetric simulations from which the resultant GW emission has been studied [86,88]. 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 [83] 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) [107], 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 cm 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,232,110].
FHH's [86] numerical quadrupole estimate of the GWs from polar oscillations in the collapse Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 simulations of Fryer and Heger [83] predicts a peak dimensionless amplitude h pk = 4.1 × 10 −23 (for d = 10 Mpc), emitted at f GW ≈ 20 Hz. The radiated energy E GW ∼ 2 × 10 44 erg. This signal would be just out of the detectability range of the LIGO-II detector.
The cores in the simulations of Fryer and Heger [83] 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 [88]) 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 GW ≈ 10 3 Hz and P GW = 10 53 erg s −1 . LIGO-II should be able to detect such a signal (see Figure 2, 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 [148,146]. 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 [146] 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 GW of the emitted radiation is in the range 10 2 -10 3 Hz. 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 [134] due to the complicated phase evolution of the emission).
Alternatively, Ou et al. [191] bypassed the long secular timescale by increasing the driving force of the instability. They found that a bar instability was maintained for several orbits before sheer flows, producing GW emission that would have a signal-to-noise ratio greater than 8 for LIGO-II out to 32 Mpc. A movie of this simulation is shown in Figure 6. FHH predict that a fragmentation instability is unlikely to develop during core collapse SNe because the cores have central density maxima (see also [88]). However, they do give estimates [calculated via Equations (4, 5)] for the amplitude, power, and frequency of the emission from such an instability: h pk ∼ 2 × 10 −22 , P GW = 10 54 erg s −1 , f GW ≈ 2 × 10 3 Hz. Again, this signal would fall just beyond the upper limit of LIGO-II's frequency range.
The GW emission from r-mode unstable neutron star remnants of core collapse SNe would be easily detectable if α max ∼ 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 [153]). 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 Hz (see Section 2.4 for details). They estimate the emitted energy to exceed 10 52 erg.
General relativity has been more fully accounted for in the core collapse studies of Dimmelmeier, Font, and Müller [58,59,60] and Shibata and collaborators [230], which build on the Newtonian, axisymmetric collapse simulations of Zwerger and Müller [271]. 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 [56] 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 [271] 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. [60] are shown in Figures 7, 8, 9, and 10. The four evolutions shown include a regular collapse (Movie 7), a rapid collapse (Movie 8), a multiple bounce collapse (Movie 9), and a very rapidly and differentially rotating collapse (Movie 10). 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 [165].
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 pk 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 pk were observed in models for which the collapse type changed. Overall, the range of h pk (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 GW was somewhat higher in the relativistic case (1.5 × 10 47 erg compared to the Newtonian value of 6.4 × 10 46 erg). 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 GW , to increase (up to five-fold). Studies of non-linear pulsations in neutron stars expect high frequencies between Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 1.8 -3.6 kHz [235]. For most of their models, this increase in f GW was not accompanied by an increase in h pk . 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 [164].

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 [41,173]. 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 timedependent 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 [173].
Since convection was suggested as a key ingredient into the explosion, it has been postulated that asymmetries in the convection can produce the large proper motions observed in the pulsar population [111]. Convection asymmetries can either be produced by asymmetries in the progenitor star that grow during collapse or by instabilities in the convection itself. Burrows and Hayes [41] proposed that asymmetries in the collapse could produce the pulsar velocities. The idea behind this work was 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 [20,147]. These asymmetries will then drive asymmetries in the convection and ultimately, the supernova explosion. Burrows and Hayes [41] found that not only could they produce strong motions in the nascent neutron star, but detectable gravitational wave signals. The peak amplitude calculated was h pk ∼ 3 × 10 −24 , for a source located at 10 Mpc.
Fryer [80] was unable to produce the large neutron star velocities seen by Burrows and Hayes [41] even after significantly increasing the level of asymmetry in the initial star in excess of 25%. This discrepancy is now known to be due to the crude 2-dimensional model and gravity scheme used by Burrows and Hayes [40]. However, the gravitational wave signal produced by both simulations is comparable. Figure 11 shows the gravitational waveform from the Burrows & Hayes simulation (including separate matter and neutrino contributions). Figures 12, 13 show the matter and neutrino contributions respectively to the gravitational wave forms for the Fryer results. The gravitational wave amplitude is dominated by the neutrino component and can exceed h pk ∼ 6 × 10 −24 , for a source located at 10 Mpc in Fryer's most extreme example.
The study by Nazin and Postnov [183] predicts a lower limit for E GW 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 GW (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 GW = 1 kHz.
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 [176]. They numerically computed the GW emission from the convection-induced aspherical Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2  mass motion and neutrino emission in the quadrupole approximation (for details, see Section 3 of their paper).
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 km). The convection in this region, which begins approximately 10 -20 ms after the shock forms and may last for ≈ 20 ms -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 postbounce model of Hillebrandt [113]. 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 [176] for justification); however, a sophisticated equation of state was utilized. Figure 14 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 Mpc. More recent calculations get amplitudes of ≈ 10 −26 in 2D [177] and ≈ 3 − 5 × 10 −26 in 3D [87]. The emitted energy was 9.8 × 10 44 erg in 2D and 1.3 × 10 42 erg 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 under-and overshooting in 3D. The relatively low angular resolution of the 3D simulations may have also played a role. The quadrupole gravitational wave amplitude A E2 20 from the 2D simulation is shown in the upper left panel of Figure 15 (see [271,241] for expressions relating A E2 20 to h). Convection in the hot bubble region between the shock and neutrino sphere arises because of an unstable entropy gradient resulting from the outward moving shock and subsequent neutrino heating. Figure 16 shows a movie of the development of this entropy-driven convection. This unstable region extends over the inner mass range 1.25 -1.40 M (corresponding to a radial range of ≈ 100 -1000 km). Convection in the hot bubble begins ≈ 50 -80 ms after shock formation Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 Figure 14: Convective instabilities inside the proto-neutron star in the 2D simulation of Müller and Janka [176]. The evolutions of the temperature (left panels) and logarithmic density (right panels) distributions are shown for the radial region 15 -95 km. The upper and lower panels correspond to times 12 and 21 ms, 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 cm −3 . The temperature and density both increase as the colors change from blue to green, yellow, and red. (Figure 7  and lasts for ≈ 100 -500 ms. 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 pk ≈ 5 × 10 −25 , for d = 10 Mpc. The emitted energy was 2 × 10 42 erg. 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 GW ∼ 10 Hz. The amplitude of such a signal would be too low to be detectable with LIGO-II.
The case for GWs from convection induced asymmetric neutrino emission has also varied with time. Müller and Janka estimated the GW emission from the convection induced anisotropic neutrino radiation in their simulations (see [176] for details). They found that the amplitude of the GWs emitted can be a factor of 5 -10 higher than the GW amplitudes resulting from convective Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 mass motion. Müller et al. (2004) [177] argue now that the GWs produced by asymmetric neutrino emission is less than that of the convective motions.
Our understanding of the convective engine is evolving with time. Scheck et al. [212] found that the convective cells could merge with time, producing a single lobe convective instability, fulfilling the prediction by Herant [111]. Blondin et al. [26] argue that a standing shock instability could develop to drive low-mode convection and Burrows et al. [43] argue that it is the shocks produced in this convection that truly drives the supernova explosion. This convection can drive oscillations in the neutron star which may also be a source for GWs (see Fig. 17).

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 [69]. 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 Hz, at frequencies of a few times 10 2 Hz. Such a signal may be detected by a pair of LIGO-II detectors. Buonanno et al. [38] have recently redone such a study with all the current results on supernovae, arguing that the GW background would be detectable by second-generation (e.g., Big Bang Observatory) space-based detectors, noting that Pop III (black-hole forming stars) could well dominate the total background.

Collapse scenario
Without winds, stars above ∼ 40 -50 M will not produce supernova explosions, but collapse directly to black holes [78]. Stars above ∼ 20 -25 M are believed to produce only weak supernova explosions. These explosions are unable to eject all of the star's mantle, and will ultimately also collapse to black holes as this mantle falls back onto the core [78]. If these stars are rotating, they will develop accretion disks around the newly formed black hole, opening up new ways to produce both GWs and baryonic explosions.
These black hole accretion disks formed by stellar collapse are currently the favored engine behind long-duration gamma-ray bursts [262,160]. It is believed that either neutrinos rising from the disk annihilate and drive an explosion or magnetic fields wound up in the disk can produce jets that then evolve into the highly relativistic outflows observed in the GRBs [195]. Relativistic and magnetic field effects are far more important in these calculations than those for normal stellar collapse models. And yet, most simulations of GRBs to date do not include the detailed microphysics, magnetic fields (although see Proga and collaborators [196]), and relativity at the same level of sophistication as the simulations of core-collapse supernovae. Even so, this topic has garnered some interest in the GW community and this interest will only grow as the simulations become more accurate.

Formation rate
The fraction of massive stars collapsing to black holes can be estimated theoretically and lies between 5 -40% [89]. But it is much more difficult to determine the fraction of these black-hole forming stars that also have enough rotation to produce a GRB. Likewise, the wide-field monitoring of the night sky of GRBs produces an accurate observed rate of GRBs, but this rate must be corrected by beaming factors to get a true long-duration GRB rate. Observations of supernova remnants have allowed GRB observers to place some constraints on the GRB rate relative to the Type Ib/c supernova rate: Rate GRB /Rate SNIb/c < 3% [22]. This predicts a Galactic rate that is below ∼ 10 −4 yr −1 and the true value is likely to lie between ∼ 10 −5 -10 −4 yr −1 . Correlating the observed GRB rate to the rate of black hole accretion disks formed in stellar collapse is not straightforward. First, it may be that long-duration GRBs are not formed from stellar collapse, but from stellar mergers [93,94]. If this is the case, then we do not expect any GW signal from the collapse itself, but any disk instabilities may still produce GW signals. Second, there may be a number of black hole accretion disk systems produced that do not produce gamma-ray bursts because they can not achieve the high relativistic flows, increasing the total rate of GW sources.

GW emission mechanisms
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 [269] have suggested that the torus of matter surrounding the black hole may be an even stronger source of GWs than the collapse itself (see also [140,141,251,250]).  have further found that the GW signal will be polarized and this polarization may be detectable.

Numerical Predictions of GW emission Mechanisms
The first fully general relativistic investigations of stellar core collapse were Nakamura's 2D simulations of rotating collapse [179,180]. 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).
Fully general relativistic collapse simulations (i.e., without the conformally flat approximation) have also been performed by Shibata [222]. 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, longduration 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. More recent results suggest that this limit can be eased for differentially rotating massive stellar cores, and systems with spin parameter q as high as 2.5 may collapse to form black holes [217]. Shibata did not compute the GW emission in his collapse simulations, but in a recent study using axisymmetric calculations, GW signals have been calculataed focusing on this collapse phase [267].
Duez et al. [62] found that if a black hole does form, but the disk is spinning rapidly, that the disk will fragment and its subsequent accretion will be in spurts, causing a "splash" onto the black hole, producing ringing and GW emission. Their result implies very strong gravitational wave amplitudes 10 −21 at distances of 10 Mpc. Black hole ringing was also estimated by FHH, where they too assumed discrete accretion events. They found that even with very optimistic accretion scenarios, that such radiation will be of very low amplitude and beyond the upper frequency reach of LIGO-II (see [86] for details).
The new general relativistic hydrodynamics simulations of Zanotti, Rezzolla, and Font [269] 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.  Figure 2). 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 [200].
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 [251]. 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, Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 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 GW ∼ 500 Hz, placing the radiation into a range detectable by LIGO-I [251].

Going further
As we noted in Section 3.5, the background of GW emission is likely dominated by the population of stars that collapse to black holes. We refer the reader back to this section on progress in this arena.
Baumgarte and collaborators have studied the neutrino signal from the collapse of a hot neutron star down to a black hole [18,19,15]. This collapse phase has a specific neutrino signal, that when compared with the gravitational wave observations, can tell us a lot about the behavior of matter at nuclear densities.

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,96]. 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 [96]. 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 [10]), the events associated with their deaths will take place at roughly the redshifts of their births.

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 [117,47]. 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 yr −1 . Thus, a reasonable occurrence rate can be found for an observation (luminosity) distance of ∼ 50 Gpc (which corresponds to a redshift of z = 5, in the cosmology used by [86]). However, uncertainties in the assumptions make this formation rate uncertain by a few orders of magnitude [1,96,2,86]. Indeed, the latest results suggest that the initial mass function begins to fall off dramatically above 100 M and the number of stars fomred with masses above 25 M may be many orders of magnitude below this rate [189]. Instead, Population III stars may be dominated by the collapse of stars in the 20 -100 M . These stars may also collapse ultimately to black holes. For discussion of these objects, we refer the reader back to Section 4.

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.

Numerical predictions of GW emission
The GW emission from the collapse of Population III stars has been investigated by Fryer and collaborators (Fryer, Woosley, and Heger [96], FHH [86], and Fryer, Holz, Hughes, and Warren [88]). The collapse simulations of Fryer, Woosley, and Heger again started with rotating collapse progenitors that had been evolved with a stellar evolution code [107]. 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 pk = 8×10 −23 , Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 with a corresponding frequency of 10 Hz [86,88]. 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 [47]. 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 [86,88]. 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.
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 pk ∼ 6 × 10 −23 and the frequency range is 20 -70 Hz. This signal may be marginally detectable with LIGO-II (see Figure 2).

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 [199] and Macchetto [159]). The masses of SMBHs in the centers of more than 45 galaxies have been estimated from observations [68] and there are more than 30 galaxies in which the presence of a SMBH has been confirmed [142].
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 [129]. Supermassive stars may contract directly out of the primordial gas, if radiation and/or magnetic field pressure prevent fragmentation [101,64,100,157,33,1]. Alternatively, they may build up from fragments of stellar collisions in clusters [211,21]. Supermassive stars are radiation dominated, isentropic and convective [221,270,157]. 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 [16] 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 β shed = (T /|W |) shed . In this case, β shed = 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 [123,50,51,221,129].
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 [27]. 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 β shed Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 than a uniformly rotating star) and then follow an evolutionary path that may be similar to that described by Baumgarte and Shapiro [16]. The alternative outcome is that the star will encounter the dynamical bar instability prior to reaching the mass-shedding limit. New and Shapiro [185,186] 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.

Formation rate
An estimate of the rate of the collapse of SMSs can be derived from the quasar luminosity function. Haehnelt [99] 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 yr −1 for M = 10 8 M objects to ∼ 1 yr −1 for M = 10 6 M objects. Thus, as in the case of Population III stars, a reasonable occurrence rate can be determined for an observation (luminosity) distance of 50 Gpc.

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.

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 [220] produced collapse evolutions that were nearly homologous. In this case, the collapse time τ coll is roughly the free-fall time at the horizon The peak GW frequency f GW = τ −1 coll is then 10 −2 Hz, if the mass of the star is 10 6 M . This is in the middle of LISA's frequency band of 10 −4 -1 Hz [243,76].
The amplitude h of this burst signal can be roughly estimated in terms of the star's quadrupole moment 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.
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 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 pk ∼ 1 × 10 −18 M 10 6 M d 50 Gpc Alternatively, the star may encounter the dynamical bar mode instability prior to complete collapse. Baumgarte and Shapiro [16] 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 when R/M = 15. The corresponding h pk , again estimated in terms of the star's quadrupole moment, is The LISA sensitivity curve is shown in Figure 18 (see [119] for details on the computation of this curve; a mission time of 3 years has been assumed). The GW signal from this dynamical bar-mode could be detected with LISA. Shibata and Shapiro [228] have 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 [208]. 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. (confirmed in conformally flat simulations [207]) 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 19. Although the work of Baumgarte and Shapiro [16] 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 GW at the time of BH formation is 3 × 10 −3 Hz. This signal would be detectable with LISA (see Figure 18).
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 amplitude of emission to be h ∼ 1.2 × 10 −20 [( E GW /M )/10 −4 ] 1/2 at f GW ∼ 2 × 10 −2 Hz for an M = 10 6 M source located at a luminosity distance of 50 Gpc (see [149,242,229] for details).
Here, E GW /M is the radiated energy efficiency and may be 7 × 10 −4 [233]. This GW signal is within LISA's range of sensitivity (see Figure 18).
Living Reviews in Relativity http://www.livingreviews.org/lrr-2003-2 Figure 18: 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 Gpc. The bar-mode source is a dynamical bar-mode.

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 [86]. 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 [86], advanced neutrino transport and interaction schemes [130,150,197,240], and the performance of 3D Newtonian [198,35,91,92] and improved methods for general relativistic simulations [222,61].
There is still much work to be done toward the goal of self-consistent, 3D general relativistic collapse simulations. The most rapidly rotating progenitors (with possibly the strongest GW signals) may be produced only in specific binary systems [95,84,265]. If the GW signal is dominated by convection, then it is critical that scientists actually understand the true supernova mechanism. 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 [146] and magnetized tori surrounding black hole collapse remnants [251]. LIGO-II could observe GWs from dynamical bar-mode instabilities in AIC [155] and core collapse SNe [86]. LISA should be able to detect the collapse (and any bar-mode instabilities that develop during the collapse) of SMSs [16] and the ringdown of black hole remnants of collapsed Population III stars [86] and SMSs [208]. These observations will provide unique information about gravitational collapse, the supernova and gamma-ray burst explosion mechanisms, and their associated progenitors and remnants.

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