We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.

We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged]

Core-collapse supernovae are among the most energetic cosmic cataclysms. They are prodigious emitters of neutrinos and quite likely strong galactic sources of gravitational waves. Observation of both neutrinos and gravitational waves from the next galactic or near extragalactic core-collapse supernova will yield a wealth of information on the explosion mechanism, but also on the structure and angular momentum of the progenitor star, and on aspects of fundamental physics such as the equation of state of nuclear matter at high densities and low entropies. In this contribution to the proceedings of the Neutrino 2012 conference, we summarize recent progress made in the theoretical understanding and modeling of core-collapse supernovae. In this, our emphasis is on multi-dimensional processes involved in the explosion mechanism such as neutrino-driven convection and the standing accretion shock instability. As an example of how supernova neutrinos can be used to probe fundamental physics, we discuss how the rise time of the electron antineutrino flux observed in detectors can be used to probe the neutrino mass hierarchy. Finally, we lay out aspects of the neutrino and gravitational-wave signature of core-collapse supernovae and discuss the power of combined analysis of neutrino and gravitational wave data from the next galactic core-collapse supernova.

We present two methods for integrating forced geodesic equations in the Kerr spacetime, which can accommodate arbitrary forces. As a test case, we compute inspirals under a simple drag force, mimicking the presence of gas. We verify that both methods give the same results for this simple force. We find that drag generally causes eccentricity to increase throughout the inspiral. This is a relativistic effect qualitatively opposite to what is seen in gravitational-radiation-driven inspirals, and similar to what is observed in hydrodynamic simulations of gaseous binaries. We provide an analytic explanation by deriving the leading order relativistic correction to the Newtonian dynamics. If observed, an increasing eccentricity would provide clear evidence that the inspiral was occurring in a non-vacuum environment. Our two methods are especially useful for evolving orbits in the adiabatic regime. Both use the method of osculating orbits, in which each point on the orbit is characterized by the parameters of the geodesic with the same instantaneous position and velocity. Both methods describe the orbit in terms of the geodesic energy, axial angular momentum, Carter constant, azimuthal phase, and two angular variables that increase monotonically and are relativistic generalizations of the eccentric anomaly. The two methods differ in their treatment of the orbital phases and the representation of the force. In one method the geodesic phase and phase constant are evolved together as a single orbital phase parameter, and the force is expressed in terms of its components on the Kinnersley orthonormal tetrad. In the second method, the phase constants of the geodesic motion are evolved separately and the force is expressed in terms of its Boyer-Lindquist components. This second approach is a generalization of earlier work by Pound and Poisson for planar forces in a Schwarzschild background.

Gravitational waves from test masses bound to geodesic orbits of rotating black holes are simulated, using Teukolsky's black hole perturbation formalism, for about ten thousand generic orbital configurations. Each binary radiates power exclusively in modes with frequencies that are integer-linear-combinations of the orbit's three fundamental frequencies. The following general spectral properties are found with a survey of orbits: (i) 99% of the radiated power is typically carried by a few hundred modes, and at most by about a thousand modes, (ii) the dominant frequencies can be grouped into a small number of families defined by fixing two of the three integer frequency multipliers, and (iii) the specifics of these trends can be qualitatively inferred from the geometry of the orbit under consideration. Detections using triperiodic analytic templates modeled on these general properties would constitute a verification of radiation from an adiabatic sequence of black hole orbits and would recover the evolution of the fundamental orbital frequencies. In an analogy with ordinary spectroscopy, this would compare to observing the Bohr model's atomic hydrogen spectrum without being able to rule out alternative atomic theories or nuclei. The suitability of such a detection technique is demonstrated using snapshots computed at 12-hour intervals throughout the last three years before merger of a kludged inspiral. Because of circularization, the number of excited modes decreases as the binary evolves. A hypothetical detection algorithm that tracks mode families dominating the first 12 hours of the inspiral would capture 98% of the total power over the remaining three years.

This is the second in a series of papers whose aim is to generate ``adiabatic'' gravitational waveforms from the inspiral of stellar-mass compact objects into massive black holes. In earlier work, we presented an accurate (2+1)D finite-difference time-domain code to solve the Teukolsky equation, which evolves curvature perturbations near rotating (Kerr) black holes. The key new ingredient there was a simple but accurate model of the singular source term based on a discrete representation of the Dirac-delta function and its derivatives. Our earlier work was intended as a proof of concept, using simple circular, equatorial geodesic orbits as a testbed. Such a source is effectively static, in that the smaller body remains at the same coordinate radius and orbital inclination over an orbit. (It of course moves through axial angle, but we separate that degree of freedom from the problem. Our numerical grid has only radial, polar, and time coordinates.) We now extend the time-domain code so that it can accommodate dynamic sources that move on a variety of physically interesting world lines. We validate the code with extensive comparison to frequency-domain waveforms for cases in which the source moves along generic (inclined and eccentric) bound geodesic orbits. We also demonstrate the ability of the time-domain code to accommodate sources moving on interesting non-geodesic worldlines. We do this by computing the waveform produced by a test mass following a ``kludged'' inspiral trajectory, made of bound geodesic segments driven toward merger by an approximate radiation loss formula.

I review the status of research, conducted by a variety of independent groups, aimed at the eventual observation of Extreme Mass Ratio Inspirals (EMRIs) with gravitational wave detectors. EMRIs are binary systems in which one of the objects is much more massive than the other, and which are in a state of dynamical evolution that is dominated by the effects of gravitational radiation. Although these systems are highly relativistic, with the smaller object moving relative to the larger at nearly light-speed, they are well described by perturbative calculations which exploit the mass ratio as a natural small parameter. I review the use of such approximations to generate waveforms needed by data analysis algorithms for observation. I also briefly review the status of developing the data analysis algorithms themselves. Although this article is almost entirely a review of previous work, it includes (as an appendix) a new analytical estimate for the time over which the influence of radiation on the binary itself is observationally negligible.

Using black hole perturbation theory, we calculate the gravitational waves produced by test particles moving on bound geodesic orbits about rotating black holes. The orbits we consider are generic - simultaneously eccentric and inclined. The waves can be described as having radial, polar, and azimuthal "voices", each of which can be made to dominate by varying eccentricity and inclination. Although each voice is generally apparent in the waveform, the radial voice is prone to overpowering the others. We also compute the radiative fluxes of energy and axial angular momentum at infinity and through the event horizon. These fluxes, coupled to a prescription for the radiative evolution of the Carter constant, will be used in future work to adiabatically evolve through a sequence of generic orbits. This will enable the calculation of inspiral waveforms that, while lacking certain important features, will approximate those expected from astrophysical extreme mass ratio captures sufficiently well to aid development of measurement algorithms on a relatively short timescale.

In many astrophysical problems, it is important to understand the behavior of functions that come from rotating (Kerr) black hole orbits. It can be particularly useful to work with the frequency domain representation of those functions, in order to bring out their harmonic dependence upon the fundamental orbital frequencies of Kerr black holes. Although, as has recently been shown by W. Schmidt, such a frequency domain representation must exist, the coupled nature of a black hole orbit's $r$ and $\theta$ motions makes it difficult to construct such a representation in practice. Combining Schmidt's description with a clever choice of timelike coordinate suggested by Y. Mino, we have developed a simple procedure that sidesteps this difficulty. One first Fourier expands all quantities using Mino's time coordinate $\lambda$. In particular, the observer's time $t$ is decomposed with $\lambda$. The frequency domain description is then built from the $\lambda$-Fourier expansion and the expansion of $t$. We have found this procedure to be quite simple to implement, and to be applicable to a wide class of functionals. We test the procedure using a simple test function, and then apply it in a particularly interesting case, the Weyl curvature scalar $\psi_4$ used in black hole perturbation theory.

A key source for LISA will be the inspiral of compact objects into supermassive black holes. Recently Mino has shown that in the adiabatic limit, gravitational waveforms for these sources can be computed using for the radiation reaction force the gradient of one half the difference between the retarded and advanced metric perturbations. Using post-Newtonian expansions, we argue that the resulting waveforms should be sufficiently accurate for signal detection with LISA. Data-analysis templates will require higher accuracy, going beyond adiabaticity; this remains a significant challenge. We describe an explicit computational procedure for obtaining waveforms based on Mino's result, for the case of a point particle coupled to a scalar field. We derive an expression for the time-averaged time derivative of the Carter constant, and verify that the expression correctly predicts that circular orbits remain circular while evolving via radiation reaction. The derivation uses detailed properties of mode expansions, Green's functions and bound geodesic orbits in the Kerr spacetime, which we review in detail. This paper is about three quarters review and one quarter new material. The intent is to give a complete and self-contained treatment of scalar radiation reaction in the Kerr spacetime, in a single unified notation, starting with the Kerr metric and ending with formulae for the time evolution of all three constants of the motion that are sufficiently explicit to be used immediately in a numerical code.

We describe progress evolving an important limit of binary orbits in general relativity, that of a stellar mass compact object gradually spiraling into a much larger, massive black hole. These systems are of great interest for gravitational wave observations. We have developed tools to compute for the first time the radiated fluxes of energy and angular momentum, as well as instantaneous snapshot waveforms, for generic geodesic orbits. For special classes of orbits, we compute the orbital evolution and waveforms for the complete inspiral by imposing global conservation of energy and angular momentum. For fully generic orbits, inspirals and waveforms can be obtained by augmenting our approach with a prescription for the self force in the adiabatic limit derived by Mino. The resulting waveforms should be sufficiently accurate to be used in future gravitational-wave searches.

We address the issue of finding an optimal detection method for a discontinuous or intermittent gravitational wave stochastic background. Such a signal might sound something like popcorn popping. We derive an appropriate version of the maximum likelihood detection statistic, and compare its performance to that of the standard cross-correlation statistic both analytically and with Monte Carlo simulations. The maximum likelihood statistic performs better than the cross-correlation statistic when the background is sufficiently non-Gaussian. For both ground and space based detectors, this results in a gain factor, ranging roughly from 1 to 3, in the minimum gravitational-wave energy density necessary for detection, depending on the duty cycle of the background. Our analysis is exploratory, as we assume that the time structure of the events cannot be resolved, and we assume white, Gaussian noise in two collocated, aligned detectors. Before this detection method can be used in practice with real detector data, further work is required to generalize our analysis to accommodate separated, misaligned detectors with realistic, colored, non-Gaussian noise.

We derive a detection method for a stochastic background of gravitational waves produced by events where the ratio of the average time between events to the average duration of an event is large. Such a signal would sound something like popcorn popping. Our derivation is based on the somewhat unrealistic assumption that the duration of an event is smaller than the detector time resolution.