The defining feature of a classical black hole is being a perfect absorber. Any evidence showing otherwise would indicate a departure from the standard black-hole picture. Energy and angular momentum absorption by the horizon of a black hole is responsible for tidal heating in a binary. This effect is particularly important in the latest stages of an extreme mass ratio inspiral around a spinning supermassive object, one of the main targets of the future LISA mission. We study how this effect can be used to probe the nature of supermassive objects in a model independent way. We compute the orbital dephasing and the gravitational-wave signal emitted by a point particle in circular, equatorial motion around a spinning supermassive object to the leading order in the mass ratio. Absence of absorption by the central object can affect the gravitational-wave signal dramatically, especially at high spin. This effect will make it possible to put an unparalleled upper bound on the reflectivity of exotic compact objects, at the level of ${\cal O}(0.01)\%$. This stringent bound would exclude the possibility of observing echoes in the ringdown of a supermassive binary merger.

The coalescence of two black holes generates gravitational waves that carry detailed information about the properties of those black holes and their binary configuration. The final coalescence cycles are in the form of a {\it ringdown}: a superposition of quasi-normal modes of the merged remnant black hole. Each mode has an oscillation frequency and decay time that in general relativity is determined by the remnant's mass and spin. Measuring the frequency and decay time of multiple modes makes it possible to measure the remnant's mass and spin, and to test the waves against the predictions of gravity theories. In this {\it Letter}, we show that the relative amplitudes of these modes encodes information about a binary's {\it geometry}. Focusing on the large mass-ratio limit, which provides a simple-to-use tool for effectively exploring parameter space, we demonstrate how a binary's geometry is encoded in the relative amplitudes of these modes, and how to parameterize the modes in this limit. Although more work is needed to assess how well this carries over to less extreme mass ratios, our results indicate that measuring multiple ringdown modes from coalescence may aid in measuring important source properties, such as the misalignment of its members' spins and orbit.

The last gravitational waves emitted in the coalescence of two black holes are quasi-normal ringing modes of the merged remnant. In general relativity, the mass and the spin of the remnant black hole uniquely determine the frequency and damping time of each radiated mode. The amplitudes of these modes are determined by the mass ratio of the system and the geometry of the coalescence. This paper is part I of an analysis that aims to compute the "excitation factors" associated with misaligned binary black hole coalescence. To simplify the analysis, we consider a large mass ratio system consisting of a non-spinning body of mass $\mu$ that inspirals on a quasi-circular trajectory into a Kerr black hole of mass $M$ and spin parameter $a$, with $\mu/M \ll 1$. Our goal is to understand how different modes are excited as a function of the black hole spin $a$ and an angle $I$ which characterizes the misalignment of the orbit with the black hole's spin axis. Though the large mass ratio limit does not describe the binaries that are being observed by gravitational-wave detectors today, this limit makes it possible to quickly and easily explore the binary parameter space, and to develop insight into how the system's late ringing waves depend on the binary's geometry. In this first analysis, we develop the worldline which the small body follows as it inspirals and then plunges into the large black hole. Our analysis generalizes earlier work by Ori and Thorne to describe how a non-equatorial circular inspiral transitions into a plunging trajectory that falls into the black hole. The worldlines which we develop here are used in part II as input to a time-domain black hole perturbation solver. This solver computes the gravitational waves generated by such inspirals and plunges, making it possible to characterize the modes which the coalescence excites.

Using inspiral and plunge trajectories we construct with a generalized Ori-Thorne algorithm, we use a time-domain black hole perturbation theory code to compute the corresponding gravitational waves. The last cycles of these waveforms are a superposition of Kerr quasinormal modes. In this paper, we examine how the modes' excitations vary as a function of source parameters, such as the larger black hole's spin and the geometry of the smaller body's inspiral and plunge. We find that the mixture of quasinormal modes that characterize the final gravitational waves from a coalescence is entirely determined by the spin $a$ of the larger black hole, an angle $I$ which characterizes the misalignment of the orbital plane from the black hole's spin axis, a second angle $\theta_{\rm fin}$ which describes the location at which the small body crosses the black hole's event horizon, and the direction sgn$(\dot\theta_{\rm fin})$ of the body's final motion. If these large-mass-ratio results hold at less extreme mass ratios, then measuring multiple ringdown modes of binary black hole coalescence gravitational waves may provide important information about the source's binary properties, such as the misalignment of the orbit's angular momentum with black hole spin. This may be particularly useful for large mass binaries, for which the early inspiral waves are out of the detectors' most sensitive band.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

Bound orbits of black holes are very well understood. Given a Kerr black hole of mass $M$ and spin $S = aM^2$, it is simple to characterize its orbits as functions of the orbit's geometry. How do the orbits change if the black hole is itself evolving? How do the orbits change if the orbiting body evolves? In this paper, we consider a process that changes a black hole's mass and spin, acting such that the spacetime is described by the Kerr solution at any moment, or that changes the orbiting body's mass. Provided this change happens slowly, the orbit's actions ($J_r, J_\theta, J_\phi$) are {\it adiabatic invariants}, and thus are constant during this process. By enforcing adiabatic invariance of the actions, we deduce how an orbit evolves due to changes in the black hole's mass and spin and in the orbiting body's mass. We demonstrate the impact of these results with several examples: how an orbit responds if accretion changes a black hole's mass and spin; how it responds if the orbiting body's mass changes due to accretion; and how the inspiral of a small body into a black hole is affected by change to the hole's mass and spin due to the gravitational radiation absorbed by the event horizon. In all cases, the effect is very small, but can be an order of magnitude or more larger than what was found in previous work which did not take into account how the orbit responds due to these effects.

The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

In previous work, we developed tools for quantifying the tidal distortion of a black hole's event horizon due to an orbiting companion. These tools use techniques which require large mass ratios (companion mass $\mu$ much smaller than black hole mass $M$), but can be used for arbitrary bound orbits, and for any black hole spin. We also showed how to visualize these distorted black holes by embedding their horizons in a global Euclidean 3-space, ${\mathbb{E}}^3$. Such visualizations illustrate interesting and important information about horizon dynamics. Unfortunately, we could not visualize black holes with spin parameter $a_* > \sqrt{3}/2 \approx 0.866$: such holes cannot be globally embedded into ${\mathbb{E}}^3$. In this paper, we overcome this difficulty by showing how to embed the horizons of tidally distorted Kerr black holes in a hyperbolic 3-space, ${\mathbb{H}}^3$. We use black hole perturbation theory to compute the Gaussian curvatures of tidally distorted event horizons, from which we build a two-dimensional metric of their distorted horizons. We develop a numerical method for embedding the tidally distorted horizons in ${\mathbb{H}}^3$. As an application, we give a sequence of embeddings into ${\mathbb{H}}^3$ of a tidally interacting black hole with spin $a_*=0.9999$. A small amplitude, high frequency oscillation seen in previous work shows up particularly clearly in these embeddings.

The Doppler effect is a shift in the frequency of waves emitted from an object moving relative to the observer. By observing and analysing the Doppler shift in electromagnetic waves from astronomical objects, astronomers gain greater insight into the structure and operation of our universe. In this paper, a simple technique is described for teaching the basics of the Doppler effect to undergraduate astrophysics students using acoustic waves. An advantage of the technique is that it produces a visual representation of the acoustic Doppler shift. The equipment comprises a 40 kHz acoustic transmitter and a microphone. The sound is bounced off a computer fan and the signal collected by a DrDAQ ADC and processed by a spectrum analyser. Widening of the spectrum is observed as the fan power supply potential is increased from 4 to 12 V.

This paper describes a novel activity to model the dynamics of a Jupiter-mass, trans-Neptunian planet of a highly eccentric orbit. Despite a history rooted in modern astronomy, "Planet X", a hypothesised hidden planet lurking in our outer Solar System, has often been touted by conspiracy theorists as the cause of past mass extinction events on Earth, as well as other modern-day doomsday scenarios. Frequently dismissed as pseudoscience by astronomers, these stories continue to draw the attention of the public by provoking mass media coverage. Targeted at junior undergraduate levels, this activity allows students to debunk some of the myths surrounding Planet X by using simulation software to demonstrate that such a large-mass planet with extreme eccentricity would be unable to enter our Solar System unnoticed, let alone maintain a stable orbit.

We model the inspiral of a compact object into a more massive black hole rotating very near the theoretical maximum. We find that once the body enters the near-horizon regime the gravitational radiation is characterized by a constant frequency, equal to (twice) the horizon frequency, with an exponentially damped profile. This contrasts with the usual "chirping" behavior and, if detected, would constitute a "smoking gun" for a near-extremal black hole in nature.

In a previous paper, we developed tools for studying the horizon geometry of a Kerr black hole that is tidally distorted by a binary companion using techniques that require large mass ratios but can be applied to any bound orbit and allow for arbitrary black hole spin. We now apply these tools to generic Kerr black hole orbits. This allows us to investigate horizon dynamics: the tidal field perturbing the horizon's geometry varies over a generic orbit, with significant variations for eccentric orbits. Many of the features of the horizon's behavior found previously carry over to the dynamical case in a natural way. In particular, we find significant offsets between the applied tide and the horizon's response. This leads to bulging in the horizon's geometry which can lag or lead the orbit, depending upon the hole's rotation and the orbit's geometry. An interesting and apparently new feature we find are small-amplitude, high-frequency oscillations in the horizon's response. We have not been able to identify a mechanism for producing these oscillations, but find that they appear most clearly when rapidly rotating black holes are distorted by very strong-field orbits.

A small body orbiting a black hole follows a trajectory that, at leading order, is a geodesic of the black hole spacetime. Much effort has gone into computing "self force" corrections to this motion, arising from the small body's own contributions to the system's spacetime. Another correction to the motion arises from coupling of the small body's spin to the black hole's spacetime curvature. Spin-curvature coupling drives a precession of the small body, and introduces a "force" (relative to the geodesic) which shifts the small body's worldline. These effects scale with the small body's spin at leading order. In this paper, we show that the equations which govern spin-curvature coupling can be analyzed with a frequency-domain decomposition, at least to leading order in the small body's spin. We show how to compute the frequency of precession along generic orbits, and how to describe the small body's precession and motion in the frequency domain. We illustrate this approach with a number of examples. This approach is likely to be useful for understanding spin coupling effects in the extreme mass ratio limit, and may provide insight into modeling spin effects in the strong field for non-extreme mass ratios.

Extreme mass ratio inspirals (EMRIs) show a strong separation of timescales, with the time characterizing inspiral, $T_{\rm i}$, much longer than any time $T_{\rm o}$ characterizing orbital motions. The ratio of these timescales (which is essentially an EMRI's mass ratio) can be regarded as a parameter that controls a perturbative expansion. Here we describe the value and limitations of an "adiabatic" description of these binaries, which uses only the leading terms arising from such a two-timescale expansion. An adiabatic approach breaks down when orbits evolve through resonances, with important dynamical and observational consequences. We describe the shortfalls of an approach that only includes the adiabatic contributions to EMRI evolution, and outline what must be done to evolve these systems through resonance and to improve our ability to model EMRI systems more generally.

We report periods and JHKL observations for 648 oxygen-rich Mira variables found in two outer bulge fields at b=-7 degrees and l=+/-8 degrees and combine these with data on 8057 inner bulge Miras from the OGLE, Macho and 2MASS surveys, which are concentrated closer to the Galactic centre. Distance moduli are estimated for all these stars. Evidence is given showing that the bulge structure is a function of age. The longer period Miras (log P > 2.6, age about 5 Gyr and younger) show clear evidence of a bar structure inclined to the line of sight in both the inner and outer regions. The distribution of the shorter period (metal-rich globular cluster age) Miras, appears spheroidal in the outer bulge. In the inner region these old stars are also distributed differently from the younger ones and possibly suggest a more complex structure. These data suggest a distance to the Galactic centre, R0, of 8.9 kpc with an estimated uncertainty of 0.4 kpc. The possible effect of helium enrichment on our conclusions is discussed.

This paper describes a simple activity for plotting and characterizing the light curve from an exoplanet transit event by way of differential photometry analysis. Using free digital imaging software, participants analyse a series of telescope images with the goal of calculating various exoplanet parameters, including its size, orbital radius and habitability. The activity has been designed for a high school or undergraduate university level and introduces fundamental concepts in astrophysics and an understanding of the basis for exoplanetary science, the transit method and digital photometry.

Tidal coupling between members of a compact binary system can have an interesting and important influence on that binary's dynamical inspiral. Tidal coupling also distorts the binary's members, changing them (at lowest order) from spheres to ellipsoids. At least in the limit of fluid bodies and Newtonian gravity, there are simple connections between the geometry of the distorted ellipsoid and the impact of tides on the orbit's evolution. In this paper, we develop tools for investigating tidal distortions of rapidly rotating black holes using techniques that are good for strong-field, fast-motion binary orbits. We use black hole perturbation theory, so our results assume extreme mass ratios. We develop tools to compute the distortion to a black hole's curvature for any spin parameter, and for tidal fields arising from any bound orbit, in the frequency domain. We also develop tools to visualize the horizon's distortion for black hole spin $a/M \le \sqrt{3}/2$ (leaving the more complicated $a/M > \sqrt{3}/2$ case to a future analysis). We then study how a Kerr black hole's event horizon is distorted by a small body in a circular, equatorial orbit. We find that the connection between the geometry of tidal distortion and the orbit's evolution is not as simple as in the Newtonian limit.

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

The first direct observation of gravitational waves' action upon matter has recently been reported by the BICEP2 experiment. Advanced ground-based gravitational-wave detectors are being installed. They will soon be commissioned, and then begin searches for high-frequency gravitational waves at a sensitivity level that is widely expected to reach events involving compact objects like stellar mass black holes and neutron stars. Pulsar timing arrays continue to improve the bounds on gravitational waves at nanohertz frequencies, and may detect a signal on roughly the same timescale as ground-based detectors. The science case for space-based interferometers targeting millihertz sources is very strong. The decade of gravitational-wave discovery is poised to begin. In this writeup of a talk given at the 2013 TAUP conference, we will briefly review the physics of gravitational waves and gravitational-wave detectors, and then discuss the promise of these measurements for making cosmological measurements in the near future.

The inspiral of a stellar mass ($1 - 100\,M_\odot$) compact body into a massive ($10^5 - 10^7\,M_\odot$) black hole has been a focus of much effort, both for the promise of such systems as astrophysical sources of gravitational waves, and because they are a clean limit of the general relativistic two-body problem. Our understanding of this problem has advanced significantly in recent years, with much progress in modeling the "self force" arising from the small body's interaction with its own spacetime deformation. Recent work has shown that this self interaction is especially interesting when the frequencies associated with the orbit's $\theta$ and $r$ motions are in an integer ratio: $\Omega_\theta/\Omega_r = \beta_\theta/\beta_r$, with $\beta_\theta$ and $\beta_r$ both integers. In this paper, we show that key aspects of the self interaction for such "resonant" orbits can be understood with a relatively simple Teukolsky-equation-based calculation of gravitational-wave fluxes. We show that fluxes from resonant orbits depend on the relative phase of radial and angular motions. The purpose of this paper is to illustrate in simple terms how this phase dependence arises using tools that are good for strong-field orbits, and to present a first study of how strongly the fluxes vary as a function of this phase and other orbital parameters. Future work will use the full dissipative self force to examine resonant and near resonant strong-field effects in greater depth, which will be needed to characterize how a binary evolves through orbital resonances.

Transient orbital resonances have recently been identified as potentially important to the inspiral of small bodies into large black holes. These resonances occur as the inspiral evolves through moments in which two fundamental orbital frequencies, $\Omega_\theta$ and $\Omega_r$, are in a small integer ratio to one another. Previous work has demonstrated that a binary's parameters are "kicked" each time the inspiral passes through a resonance, changing the orbit's characteristics relative to a model that neglects resonant effects. In this paper, we use exact Kerr geodesics coupled to an accurate but approximate model of inspiral to survey orbital parameter space and estimate how commonly one encounters long-lived orbital resonances. We find that the most important resonances last for a few hundred orbital cycles at mass ratio $10^{-6}$, and that resonances are almost certain to occur during the time that a large mass ratio binary would be a target of gravitational-wave observations. Resonances appear to be ubiquitous in large mass ratio inspiral, and to last long enough that they are likely to affect binary evolution in observationally important ways.

Recent observations have accumulated compelling evidence that some short gamma-ray bursts (SGRBs) are associated with the mergers of neutron star (NS) binaries. This would indicate that the SGRB event is associated with a gravitational-wave (GW) signal corresponding to the final inspiral of the compact binary. In addition, the radioactive decay of elements produced in NS binary mergers may result in transients visible in the optical and infrared with peak luminosities on hours-days timescales. Simultaneous observations of the inspiral GWs and signatures in the electromagnetic band may allow us to directly and independently determine both the luminosity distance and redshift to a binary. These standard sirens (the GW analog of standard candles) have the potential to provide an accurate measurement of the low-redshift Hubble flow. In addition, these systems are absolutely calibrated by general relativity, and therefore do not experience the same set of astrophysical systematics found in traditional standard candles, nor do the measurements rely on a distance ladder. We show that 15 observable GW and EM events should allow the Hubble constant to be measured with 5% precision using a network of detectors that includes advanced LIGO and Virgo. Measuring 30 beamed GW-SGRB events could constrain H_0 to better than 1%. When comparing to standard Gaussian likelihood analysis, we find that each event's non-Gaussian posterior in H_0 helps reduce the overall measurement errors in H_0 for an ensemble of NS binary mergers.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries.

The future space-based gravitational wave detector LISA will be able to measure parameters of coalescing massive black hole binaries, often to extremely high accuracy. Previous work has demonstrated that the black hole spins can have a strong impact on the accuracy of parameter measurement. Relativistic spin-induced precession modulates the waveform in a manner which can break degeneracies between parameters, in principle significantly improving how well they are measured. Recent studies have indicated, however, that spin precession may be weak for an important subset of astrophysical binary black holes: those in which the spins are aligned due to interactions with gas. In this paper, we examine how well a binary's parameters can be measured when its spins are partially aligned and compare results using waveforms that include higher post-Newtonian harmonics to those that are truncated at leading quadrupole order. We find that the weakened precession can substantially degrade parameter estimation. This degradation is particularly devastating for the extrinsic parameters sky position and distance. Absent higher harmonics, LISA typically localizes the sky position of a nearly aligned binary a factor of $\sim 6$ less accurately than for one in which the spin orientations are random. Our knowledge of a source's sky position will thus be worst for the gas-rich systems which are most likely to produce electromagnetic counterparts. Fortunately, higher harmonics of the waveform can make up for this degradation. By including harmonics beyond the quadrupole in our waveform model, we find that the accuracy with which most of the binary's parameters are measured can be substantially improved. In some cases, parameters can be measured as well in partially aligned binaries as they can be when the binary spins are random.

The coalescence of massive black holes is one of the primary sources of gravitational waves (GWs) for LISA. Measurements of the GWs can localize the source on the sky to an ellipse with a major axis of a few tens of arcminutes to a few degrees, depending on source redshift, and a minor axis which is 2--4 times smaller. The distance (and thus an approximate redshift) can be determined to better than a per cent for the closest sources we consider, although weak lensing degrades this performance. It will be of great interest to search this three-dimensional `pixel' for an electromagnetic counterpart to the GW event. The presence of a counterpart allows unique studies which combine electromagnetic and GW information, especially if the counterpart is found prior to final merger of the holes. To understand the feasibility of early counterpart detection, we calculate the evolution of the GW pixel with time. We find that the greatest improvement in pixel size occurs in the final day before merger, when spin precession effects are maximal. The source can be localized to within 10 square degrees as early as a month before merger at $z = 1$; for higher redshifts, this accuracy is only possible in the last few days.

Recent observations support the hypothesis that a large fraction of "short-hard" gamma-ray bursts (SHBs) are associated with compact binary inspiral. Since gravitational-wave (GW) measurements of well-localized inspiraling binaries can measure absolute source distances, simultaneous observation of a binary's GWs and SHB would allow us to independently determine both its luminosity distance and redshift. Such a "standard siren" (the GW analog of a standard candle) would provide an excellent probe of the relatively nearby universe's expansion, complementing other standard candles. In this paper, we examine binary measurement using a Markov Chain Monte Carlo technique to build the probability distributions describing measured parameters. We assume that each SHB observation gives both sky position and the time of coalescence, and we take both binary neutron stars and black hole-neutron star coalescences as plausible SHB progenitors. We examine how well parameters particularly distance) can be measured from GW observations of SHBs by a range of ground-based detector networks. We find that earlier estimates overstate how well distances can be measured, even at fairly large signal-to-noise ratio. The fundamental limitation to determining distance proves to be a degeneracy between distance and source inclination. Overcoming this limitation requires that we either break this degeneracy, or measure enough sources to broadly sample the inclination distribution. (Abridged)

We constrain the parameterized post-Einsteinian framework with binary pulsar observations of orbital period decay due to gravitational wave emission. This framework proposes to enhance the amplitude and phase of gravitational waveform templates through post-Einsteinian parameters to search for generic deviations from General Relativity in gravitational wave data. Such enhancements interpolate between General Relativity and alternative theory predictions, but their magnitude must be such as to satisfy all current experiments and observations. The data that currently constrains the parameterized post-Einsteinian framework the most is the orbital period decay of binary pulsars. We use such observations to place upper limits on the magnitude of post-Einsteinian parameters, which will be critical when gravitational waves are detected and this framework is implemented.

Spectacular breakthroughs in numerical relativity now make it possible to compute spacetime dynamics in almost complete generality, allowing us to model the coalescence and merger of binary black holes with essentially no approximations. The primary limitation of these calculations is now computational. In particular, it is difficult to model systems with large mass ratio and large spins, since one must accurately resolve the multiple lengthscales which play a role in such systems. Perturbation theory can play an important role in extending the reach of computational modeling for binary systems. In this paper, we present first results of a code which allows us to model the gravitational waves generated by the inspiral, merger, and ringdown of a binary system in which one member of the binary is much more massive than the other. This allows us to accurately calibrate binary dynamics in the large mass ratio regime. We focus in this analysis on the recoil imparted to the merged remnant by these waves. We closely examine the "antikick", an anti-phase cancellation of the recoil arising from the plunge and ringdown waves, described in detail by Schnittman et al. We find that, for orbits aligned with the black hole spin, the antikick grows as a function of spin. The total recoil is smallest for prograde coalescence into a rapidly rotating black hole, and largest for retrograde coalescence. Amusingly, this completely reverses the predicted trend for kick versus spin from analyses that only include inspiral information.

Gravitational wave observations will be excellent tools for making precise measurements of processes that occur in very strong-field regions of spacetime. Extreme mass ratio systems, formed by the capture of a stellar mass body compact by a massive black hole, will be targets for planned space-based interferometers such as LISA and DECIGO. These systems will be especially powerful tools for testing the spacetime nature of black hole candidates. In this writeup of the talk I gave at JGRG19, I describe how the properties of black holes are imprinted on their waveforms, and how measurements can be used to study these properties and thereby learn about the astrophysics of black holes and about strong-field gravity.

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between analytical EOB and numerical Teukolsky-based gravitational waveforms can be reduced to less than 10^(-1) rad and less than 2 x 10^(-3), respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of roughly 6-27 rad after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

Our universe contains a great number of extremely compact and massive objects which are generally accepted to be black holes. Precise observations of orbital motion near candidate black holes have the potential to determine if they have the spacetime structure that general relativity demands. As a means of formulating measurements to test the black hole nature of these objects, Collins and Hughes introduced "bumpy black holes": objects that are almost, but not quite, general relativity's black holes. The spacetimes of these objects have multipoles that deviate slightly from the black hole solution, reducing to black holes when the deviation is zero. In this paper, we extend this work in two ways. First, we show how to introduce bumps which are smoother and lead to better behaved orbits than those in the original presentation. Second, we show how to make bumpy Kerr black holes -- objects which reduce to the Kerr solution when the deviation goes to zero. This greatly extends the astrophysical applicability of bumpy black holes. Using Hamilton-Jacobi techniques, we show how a spacetime's bumps are imprinted on orbital frequencies, and thus can be determined by measurements which coherently track a small orbiting body's orbital phase. We find that weak-field orbits of bumpy black holes are modified exactly as expected from a Newtonian analysis of a body with a prescribed multipolar structure, reproducing well-known results from the celestial mechanics literature. The impact of bumps on strong-field orbits is especially strong, suggesting that this framework will allow observations to set robust limits on the extent to which a spacetime's multipoles deviate from the black hole expectation.

Largely motivated by the development of highly sensitive gravitational-wave detectors, our understanding of merging compact binaries and the gravitational waves they generate has improved dramatically in recent years. Breakthroughs in numerical relativity now allow us to model the coalescence of two black holes with no approximations or simplifications. There has also been outstanding progress in our analytical understanding of binaries. We review these developments, examining merging binaries using black hole perturbation theory, post-Newtonian expansions, and direct numerical integration of the field equations. We summarize these approaches and what they have taught us about gravitational waves from compact binaries. We place these results in the context of gravitational-wave generating systems, analyzing the impact gravitational wave emission has on their sources, as well as what we can learn about them from direct gravitational-wave measurements.

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

Increasing the sensitivity of a gravitational-wave (GW) detector improves our ability to measure the characteristics of detected sources. It also increases the number of weak signals that contribute to the data. Because GW detectors have nearly all-sky sensitivity, they can be subject to a confusion limit: Many sources which cannot be distinguished may be measured simultaneously, defining a stochastic noise floor to the sensitivity. For GW detectors operating at present and for their planned upgrades, the projected event rate is sufficiently low that we are far from the confusion-limited regime. However, some detectors currently under discussion may have large enough reach to binary inspiral that they enter the confusion-limited regime. In this paper, we examine the binary inspiral confusion limit for terrestrial detectors. We consider a broad range of inspiral rates in the literature, several planned advanced gravitational-wave detectors, and the highly advanced "Einstein Telescope" design. Though most advanced detectors will not be impacted by this limit, the Einstein Telescope with a very low frequency "seismic wall" may be subject to confusion noise. At a minimum, careful data analysis will be require to separate signals which will appear confused. This result should be borne in mind when designing highly advanced future instruments.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

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.

Following an earlier analysis which examined the effect of the self-gravity of pressure on big-bang nucleosynthesis (BBN), we explore the effect of pressure's self-gravity on the structure of neutron stars. We construct an ad hoc modification of the Tolman-Oppenheimer-Volkoff equation wherein pressure's self-gravity is parameterized by a constant, $\chi$, with $0 \le \chi \le 1$. The full general relativistic contribution to the gravity of pressure is recovered with $\chi = 1$, and is eliminated when $\chi = 0$. This formulation is not proposed as an alternative theory of gravity, but is merely used to quantify the extent to which the self-gravity of pressure contributes to the structure of dense objects. As can be surmised qualitatively, neutron star masses can be quite sensitive to $\chi$, with higher values of neutron-star mass (by $\sim$20--25%) allowed for smaller values of $\chi$. However, for a given equation of state, neither the range of neutron star radii nor the radii at fixed central density depend sensitively on $\chi$. Over the neutron star mass range measured so far, the presence or absence of pressure's self-gravity yields a nearly immeasurable change in radius -- much smaller than the variations in radius due to the uncertainty in the equation of state. In contrast to the result for BBN, we thus find that neutron stars are not likely to be useful testbeds for examining the self-gravity of pressure.

Massive black hole binary coalescences are prime targets for space-based gravitational wave (GW) observatories such as {\it LISA}. GW measurements can localize the position of a coalescing binary on the sky to an ellipse with a major axis of a few tens of arcminutes to a few degrees, depending on source redshift, and a minor axis which is $2 - 4$ times smaller. Neglecting weak gravitational lensing, the GWs would also determine the source's luminosity distance to better than percent accuracy for close sources, degrading to several percent for more distant sources. Weak lensing cannot, in fact, be neglected and is expected to limit the accuracy with which distances can be fixed to errors no less than a few percent. Assuming a well-measured cosmology, the source's redshift could be inferred with similar accuracy. GWs alone can thus pinpoint a binary to a three-dimensional ``pixel'' which can help guide searches for the hosts of these events. We examine the time evolution of this pixel, studying it at merger and at several intervals before merger. One day before merger, the major axis of the error ellipse is typically larger than its final value by a factor of $\sim 1.5-6$. The minor axis is larger by a factor of $\sim 2-9$, and, neglecting lensing, the error in the luminosity distance is larger by a factor of $\sim 1.5-7$. This large change over a short period of time is due to spin-induced precession, which is strongest in the final days before merger. The evolution is slower as we go back further in time. For $z = 1$, we find that GWs will localize a coalescing binary to within $\sim 10\ \mathrm{deg}^2$ as early as a month prior to merger and determine distance (and hence redshift) to several percent.

The coalescence of massive black holes generates gravitational waves (GWs) that will be measurable by space-based detectors such as LISA to large redshifts. The spins of a binary's black holes have an important impact on its waveform. Specifically, geodetic and gravitomagnetic effects cause the spins to precess; this precession then modulates the waveform, adding periodic structure which encodes useful information about the binary's members. Following pioneering work by Vecchio, we examine the impact upon GW measurements of including these precession-induced modulations in the waveform model. We find that the additional periodicity due to spin precession breaks degeneracies among certain parameters, greatly improving the accuracy with which they may be measured. In particular, mass measurements are improved tremendously, by one to several orders of magnitude. Localization of the source on the sky is also improved, though not as much -- low redshift systems can be localized to an ellipse which is roughly ${a few} \times 10$ arcminutes in the long direction and a factor of 2-4 smaller in the short direction. Though not a drastic improvement relative to analyses which neglect spin precession, even modest gains in source localization will greatly facilitate searches for electromagnetic counterparts to GW events. Determination of distance to the source is likewise improved: We find that relative error in measured luminosity distance is commonly $\sim 0.2%-0.7%$ at $z \sim 1$. Finally, with the inclusion of precession, we find that the magnitude of the spins themselves can typically be determined for low redshift systems with an accuracy of about $0.1%-10 %$, depending on the spin value, allowing accurate surveys of mass and spin evolution over cosmic time.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy gamma-rays. The system was observed over several orbital cycles (between September 2006 and February 2007) with the VERITAS array of imaging air-Cherenkov telescopes. A signal of gamma-rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of Gamma=2.40+-0.16_stat+-0.2_sys and a flux above 300 GeV corresponding to 15-20% of the flux from the Crab Nebula.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.

One of the most exciting potential sources of gravitational waves for low-frequency, space-based gravitational wave (GW) detectors such as the proposed Laser Interferometer Space Antenna (LISA) is the inspiral of compact objects into massive black holes in the centers of galaxies. The detection of waves from such "extreme mass ratio inspiral" systems (EMRIs) and extraction of information from those waves require template waveforms. The systems' extreme mass ratio means that their waveforms can be determined accurately using black hole perturbation theory. Such calculations are computationally very expensive. There is a pressing need for families of approximate waveforms that may be generated cheaply and quickly but which still capture the main features of true waveforms. In this paper, we introduce a family of such "kludge" waveforms and describe ways to generate them. We assess performance of the introduced approximations by comparing "kludge" waveforms to accurate waveforms obtained by solving the Teukolsky equation in the adiabatic limit (neglecting GW backreaction). We find that the kludge waveforms do extremely well at approximating the true gravitational waveform, having overlaps with the Teukolsky waveforms of 95% or higher over most of the parameter space for which comparisons can currently be made. Indeed, we find these kludges to be of such high quality (despite their ease of calculation) that it is possible they may play some role in the final search of LISA data for EMRIs.

We revisit the problem of the emission of gravitational waves from a test mass orbiting and thus perturbing a Kerr black hole. The source term of the Teukolsky perturbation equation contains a Dirac delta function which represents a point particle. We present a technique to effectively model the delta function and its derivatives using as few as four points on a numerical grid. The source term is then incorporated into a code that evolves the Teukolsky equation in the time domain as a (2+1) dimensional PDE. The waveforms and energy fluxes are extracted far from the black hole. Our comparisons with earlier work show an order of magnitude gain in performance (speed) and numerical errors less than 1% for a large fraction of parameter space. As a first application of this code, we analyze the effect of finite extraction radius on the energy fluxes. This paper is the first in a series whose goal is to develop adiabatic waveforms describing the inspiral of a small compact body into a massive Kerr black hole.

We have performed a detailed analysis of orbital motion in the vicinity of a nearly extremal Kerr black hole. For very rapidly rotating black holes (spin a=J/M>0.9524M) we have found a class of very strong field eccentric orbits whose angular momentum L_z increases with the orbit's inclination with respect to the equatorial plane, while keeping latus rectum and eccentricity fixed. This behavior is in contrast with Newtonian intuition, and is in fact opposite to the "normal" behavior of black hole orbits. Such behavior was noted previously for circular orbits; since it only applies to orbits very close to the black hole, they were named "nearly horizon-skimming orbits". Our analysis generalizes this result, mapping out the full generic (inclined and eccentric) family of nearly horizon-skimming orbits. The earlier work on circular orbits reported that, under gravitational radiation emission, nearly horizon-skimming orbits tend to evolve to smaller orbit inclination, toward prograde equatorial configuration. Normal orbits, by contrast, always demonstrate slowly growing orbit inclination (orbits evolve toward the retrograde equatorial configuration). Using up-to-date Teukolsky-fluxes, we have concluded that the earlier result was incorrect: all circular orbits, including nearly horizon-skimming ones, exhibit growing orbit inclination. Using kludge fluxes based on a Post-Newtonian expansion corrected with fits to circular and to equatorial Teukolsky-fluxes, we argue that the inclination grows also for eccentric nearly horizon-skimming orbits. We also find that the inclination change is, in any case, very small. As such, we conclude that these orbits are not likely to have a clear and peculiar imprint on the gravitational waveforms expected to be measured by the space-based detector LISA.

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 100 GeV) can impose tight constraints on some GRB emission models. Many GRB afterglow models predict a VHE component similar to that seen in blazars and plerions, in which the GRB spectral energy distribution has a double-peaked shape extending into the VHE regime. VHE emission coincident with delayed X-ray flare emission has also been predicted. GRB follow-up observations have had high priority in the observing program at the Whipple 10m Gamma-ray Telescope and GRBs will continue to be high priority targets as the next generation observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at late times (>~4 hours), from seven GRBs observed with the Whipple Telescope are reported here.

We report on observations of the sky region around the unidentified TeV gamma-ray source TeV J2032+4130 carried out with the Whipple Observatory 10 m atmospheric Cherenkov telescope for a total of 65.5 hrs between 2003 and 2005. The standard two-dimensional analysis developed by the Whipple collaboration for a stand-alone telescope reveals an excess in the field of view at a pre-trials significance level of 6.1 standard deviations. The measured position of this excess is alpha(2000) =20 h 32 m 27 s, delta(2000) = 41 deg 39 min 17 s. The estimated integral flux for this gamma-ray source is about 8% of the Crab-Nebula flux. The data are consistent with a point-like source. Here we present a detailed description of the standard two-dimensional analysis technique used for the analysis of data taken with the Whipple Observatory 10 m telescope and the results for the TeV J2032+4130 campaign. We include a short discussion of the physical mechanisms that may be responsible for the observed gamma-ray emission, based on possible association with known astrophysical objects, in particular Cygnus OB2.

LISA is a planned space-based gravitational-wave (GW) detector that would be sensitive to waves from low-frequency sources, in the band of roughly $(0.03 - 0.1) {\rm mHz} \lesssim f \lesssim 0.1 {\rm Hz}$. This is expected to be an extremely rich chunk of the GW spectrum -- observing these waves will provide a unique view of dynamical processes in astrophysics. Here we give a quick survey of some key LISA sources and what GWs can uniquely teach us about these sources. Particularly noteworthy science which is highlighted here is the potential for LISA to track the moderate to high redshift evolution of black hole masses and spins through the measurement of GWs generated from massive black hole binaries (which in turn form by the merger of galaxies and protogalaxies). Measurement of these binary black hole waves has the potential to determine the masses and spins of the constituent black holes with percent-level accuracy or better, providing a unique high-precision probe of an aspect of early structure growth. This article is based on the ``Astrophysics Tutorial'' talk given by the author at the Sixth International LISA Symposium.

The inspirals of ``small'' ($1 - 100 M_\odot$) compact bodies through highly relativistic orbits of massive (several $\times 10^5 M_\odot -$ several $\times 10^6 M_\odot$) black holes are among the most anticipated sources for the LISA gravitational-wave antenna. The measurement of these waves is expected to map the spacetime of the larger body with high precision, allowing us to test in detail the hypothesis that black hole candidates are described by the Kerr metric of general relativity. In this article, we will briefly describe how these sources can be used to perform such a test. These proposed measurements are often described as ``testing relativity''. This description is at best somewhat glib: Because -- at least to date -- all work related to these measurements assumes general relativity as the theoretical framework in which these tests are performed, the measurements cannot be said to ``test relativity'' in a fundamental way. More accurately, they test the {\it nature of massive compact bodies within general relativity}. A surprising result for such a test could point to deviations from general relativity, and would provide an experimentally motivated direction in which to pursue tests of gravity theories beyond GR.

The first atmospheric Cherenkov telescope of VERITAS (the Very Energetic Radiation Imaging Telescope Array System) has been in operation since February 2005. We present here a technical description of the instrument and a summary of its performance. The calibration methods are described, along with the results of Monte Carlo simulations of the telescope and comparisons between real and simulated data. The analysis of TeV $\gamma$-ray observations of the Crab Nebula, including the reconstructed energy spectrum, is shown to give results consistent with earlier measurements. The telescope is operating as expected and has met or exceeded all design specifications.

In April-May 2003, the blazar 1ES 1959+650 showed an increased level of X-ray activity. This prompted a multiwavelength observation campaign with the Whipple 10 m gamma-ray telescope, the Rossi X-ray Timing Explorer, the Bordeaux Optical Observatory, and the University of Michigan Radio Astrophysical Observatory. We present the multiwavelength data taken from May 2, 2003 to June 7, 2003 and compare the source characteristics with those measured during observations taken during the years 2000 and 2002. The X-ray observations gave a data set with high signal-to-noise light curves and energy spectra; however, the gamma-ray observations did not reveal a major TeV gamma-ray flare. Furthermore, we find that the radio and optical fluxes do not show statistically significant deviations from those measured during the 2002 flaring periods. While the X-ray flux and X-ray photon index appear correlated during subsequent observations, the apparent correlation evolved significantly between the years 2000, 2002, and 2003. We discuss the implications of this finding for the mechanism that causes the flaring activity.

Galaxy clusters might be sources of TeV gamma rays emitted by high-energy protons and electrons accelerated by large scale structure formation shocks, galactic winds, or active galactic nuclei. Furthermore, gamma rays may be produced in dark matter particle annihilation processes at the cluster cores. We report on observations of the galaxy clusters Perseus and Abell 2029 using the 10 m Whipple Cherenkov telescope during the 2003-2004 and 2004-2005 observing seasons. We apply a two-dimensional analysis technique to scrutinize the clusters for TeV emission. In this paper we first determine flux upper limits on TeV gamma-ray emission from point sources within the clusters. Second, we derive upper limits on the extended cluster emission. We subsequently compare the flux upper limits with EGRET upper limits at 100 MeV and theoretical models. Assuming that the gamma-ray surface brightness profile mimics that of the thermal X-ray emission and that the spectrum of cluster cosmic rays extends all the way from thermal energies to multi-TeV energies with a differential spectral index of -2.1, our results imply that the cosmic ray proton energy density is less than 7.9% of the thermal energy density for the Perseus cluster.

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.

Observations of the gravitational radiation from well-localized, inspiraling compact object binaries can measure absolute source distances with high accuracy. When coupled with an independent determination of redshift through an electromagnetic counterpart, these standard sirens can provide an excellent probe of the expansion history of the Universe and the dark energy. Short gamma-ray bursts, if produced by merging neutron star binaries, would be standard sirens with known redshifts detectable by ground-based GW networks such as LIGO-II, Virgo, and AIGO. Depending upon the collimation of these GRBs, a single year of observation of their gravitational waves can measure the Hubble constant to about 2%. When combined with measurement of the absolute distance to the last scattering surface of the cosmic microwave background, this determines the dark energy equation of state parameter w to 9%. Similarly, supermassive binary black hole inspirals will be standard sirens detectable by LISA. Depending upon the precise redshift distribution, 100 sources could measure w at the 4% level.

This article is based on a pair of lectures given at the 2005 SLAC Summer Institute. Our goal is to motivate why most physicists and astrophysicists accept the hypothesis that the most massive, compact objects seen in many astrophysical systems are described by the black hole solutions of general relativity. We describe the nature of the most important black hole solutions, the Schwarzschild and the Kerr solutions. We discuss gravitational collapse and stability in order to motivate why such objects are the most likely outcome of realistic astrophysical collapse processes. Finally, we discuss some of the observations which -- so far at least -- are totally consistent with this viewpoint, and describe planned tests and observations which have the potential to falsify the black hole hypothesis, or sharpen still further the consistency of data with theory.

Einstein's special theory of relativity revolutionized physics by teaching us that space and time are not separate entities, but join as ``spacetime''. His general theory of relativity further taught us that spacetime is not just a stage on which dynamics takes place, but is a participant: The field equation of general relativity connects matter dynamics to the curvature of spacetime. Curvature is responsible for gravity, carrying us beyond the Newtonian conception of gravity that had been in place for the previous two and a half centuries. Much research in gravitation since then has explored and clarified the consequences of this revolution; the notion of dynamical spacetime is now firmly established in the toolkit of modern physics. Indeed, this notion is so well established that we may now contemplate using spacetime as a tool for other science. One aspect of dynamical spacetime -- its radiative character, ``gravitational radiation'' -- will inaugurate entirely new techniques for observing violent astrophysical processes. Over the next one hundred years, much of this subject's excitement will come from learning how to exploit spacetime as a tool for astronomy. This article is intended as a tutorial in the basics of gravitational radiation physics.

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.

The gravitational waves generated in the coalescence of massive binary black holes will be measurable by LISA to enormous distances. Redshifts z~10 or larger (depending somewhat on the mass of the binary) can potentially be probed by such measurements, suggesting that binary coalescences can be made into cosmological tools. We discuss two particularly interesting types of probes. First, by combining gravitational-wave measurements with information about the universe's cosmography, we can study the evolution of black hole masses and merger rates as a function of redshift, providing information about the growth of structures at high redshift and possibly constraining hierarchical merger scenarios. Second, if it is possible to associate an ``electromagnetic'' counterpart with a coalescence, it may be possible to measure both redshift and luminosity distance to an event with less than ~1% error. Such a measurement would constitute an amazingly precise cosmological standard candle. Unfortunately, gravitational lensing uncertainties will reduce the quality of this candle significantly. Though not as amazing as might have been hoped, such a candle would nonetheless very usefully complement other distance-redshift probes, in particular providing a valuable check on systematic effects in such measurements.

Gravitational waves (GWs) from supermassive binary black hole (BBH) inspirals are potentially powerful standard sirens (the GW analog to standard candles) (Schutz 1986, 2002). Because these systems are well-modeled, the space-based GW observatory LISA will be able to measure the luminosity distance (but not the redshift) to some distant massive BBH systems with 1-10% accuracy. This accuracy is largely limited by pointing error: GW sources generally are poorly localized on the sky. Localizing the binary independently (e.g., through association with an electromagnetic counterpart) greatly reduces this positional error. An electromagnetic counterpart may also allow determination of the event's redshift. In this case, BBH coalescence would constitute an extremely precise (better than 1%) standard candle visible to high redshift. In practice, gravitational lensing degrades this precision, though the candle remains precise enough to provide useful information about the distance-redshift relation. Even if very rare, these GW standard sirens would complement, and increase confidence in, other standard candles.

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.

Space-born gravitational-wave interferometers such as {\it LISA} will detect the gravitational wave (GW) signal from the inspiral, plunge and ringdown phases of massive black hole binary mergers at cosmological distances. From the inspiral waves, we will be able to measure the masses of the binaries' members; from the ringdown waves, we will be able to measure the mass of the final merged remnant. A subset of detected events allow the identification of both the inspiral and the ringdown waveforms in a given source, and thus allow a measurement of the total mass-energy lost to GWs over the coalescence, $M_{\rm GW}$. We define ``golden'' binary mergers to be those with measurement errors likely to be small enough for a physically useful determination of $M_{\rm GW}$. A detailed sensitivity study, combined with simple black hole population models, suggests that a few golden binary mergers may be detected during a 3-year {\it LISA} mission lifetime. Any such mass deficit measurement would constitute a robust and valuable observational test of strong-field relativistic gravity. An extension of this concept to include spin measurements may allow a direct empirical test of the black hole area theorem.

Gravitational waves from the coalescence of binary black holes carry linear momentum, causing center of mass recoil. This ``radiation rocket'' has important implications for systems with escape speeds of order the recoil velocity. We describe new recoil calculations using high precision black hole perturbation theory to estimate the magnitude of the recoil for the slow ``inspiral'' coalescence phase; coupled with a cruder calculation for the final ``plunge'', we estimate the total recoil imparted to a merged black hole. We find that velocities of many tens to a few hundred km/sec can be achieved fairly easily. The recoil probably never exceeds about 500 km/sec.

Gravitational waves from the coalescence of binary black holes carry away linear momentum, causing center of mass recoil. This "radiation rocket" effect has important implications for systems with escape speeds of order the recoil velocity. We revisit this problem using black hole perturbation theory, treating the binary as a test mass spiraling into a spinning hole. For extreme mass ratios (q = m1/m2 << 1) we compute the recoil for the slow inspiral epoch of binary coalescence very accurately; these results can be extrapolated to q ~ 0.4 with modest accuracy. Although the recoil from the final plunge contributes significantly to the final recoil, we are only able to make crude estimates of its magnitude. We find that the recoil can easily reach ~ 100-200 km/s, but most likely does not exceed ~ 500 km/s. Though much lower than previous estimates, this recoil is large enough to have important astrophysical consequences. These include the ejection of black holes from globular clusters, dwarf galaxies, and high-redshift dark matter halos.

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular clusters, which may explain the apparent absence of massive black holes in these systems. Ejection from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass precursors.

The collapse of massive stars not only produces observable outbursts across the entire electromagnetic spectrum but, for Galactic (or near-Galactic) supernovae, detectable signals for ground-based neutrino and gravitational wave detectors. Gravitational waves and neutrinos provide the only means to study the actual engine behind the optical outbursts: the collapsed stellar core. While the neutrinos are most sensitive to details of the equation of state, gravitational waves provide a means to study the mass asymmetries in this central core. We present gravitational wave signals from a series of 3-dimensional core-collapse simulations with asymmetries derived from initial perturbations caused by pre-collapse convection, core rotation, and low-mode convection in the explosion engine itself. A Galactic supernovae will allow us to differentiate these different sources of asymmetry. Combining this signal with other observations of the supernova, from neutrinos to gamma-rays to the compact remnant, dramatically increases the predictive power of the gravitational wave signal. We conclude with a discussion of the gravitational wave signal arising from collapsars, the leading engine for long-duration gamma-ray bursts.

Observations have established that extremely compact, massive objects are common in the universe. It is generally accepted that these objects are black holes. As observations improve, it becomes possible to test this hypothesis in ever greater detail. In particular, it is or will be possible to measure the properties of orbits deep in the strong field of a black hole candidate (using x-ray timing or with gravitational-waves) and to test whether they have the characteristics of black hole orbits in general relativity. Such measurements can be used to map the spacetime of a massive compact object, testing whether the object's multipoles satisfy the strict constraints of the black hole hypothesis. Such a test requires that we compare against objects with the ``wrong'' multipole structure. In this paper, we present tools for constructing bumpy black holes: objects that are almost black holes, but that have some multipoles with the wrong value. The spacetimes which we present are good deep into the strong field of the object -- we do not use a large r expansion, except to make contact with weak field intuition. Also, our spacetimes reduce to the black hole spacetimes of general relativity when the ``bumpiness'' is set to zero. We propose bumpy black holes as the foundation for a null experiment: if black hole candidates are the black holes of general relativity, their bumpiness should be zero. By comparing orbits in a bumpy spacetime with those of an astrophysical source, observations should be able to test this hypothesis, stringently testing whether they are the black holes of general relativity. (Abridged)

Following the detection of strong TeV gamma-ray flares from the BL Lac object 1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray and TeV gamma-ray observations from May 18, 2002 to August 14, 2002. Observations with the X-ray telescope RXTE and the Whipple and HEGRA gamma-ray telescopes revealed several strong flares, enabling us to sensitively test the X-ray/gamma-ray flux correlation properties. Although the X-ray and gamma-ray fluxes seemed to be correlated in general, we found an ``orphan'' gamma-ray flare that was not accompanied by an X-ray flare. After describing in detail the radio (UMRAO, VLA), optical (Boltwood, Abastumani), X-ray (RXTE) and gamma-ray (Whipple, HEGRA) light curves and Spectral Energy Distributions (SEDs) we present initial modeling of the SED with a simple Synchrotron Self-Compton (SSC) model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and spectral X-ray and gamma-ray emission characteristics of TeV blazars are very similar, even though the masses estimates of their central black holes differ by up to one order of magnitude.

Massive black holes appear to be present in the nuclei of almost all galaxies, but their genesis and evolution are not well understood. As astrophysical black holes are completely characterized by their masses and spins, the observed joint distribution of these quantities contains important clues to their history. We examine the coevolution of mass and spin in binary merger growth scenarios. We find that holes are typically spun down by mergers. Rapid rotation results only if the binary's larger member already spins quickly and the merger with the smaller hole is consistently near prograde; or, if the binary's mass ratio approaches one. If, as some observations have suggested, observed black holes spin rapidly, then this limits the importance of merger scenarios for the growth of black holes.

The new generation of gravitational wave (GW) detectors have the potential to open a novel window onto the violent dynamics of core collapse. Although it is certain that core collapse events generate gravitational radiation, understanding the characteristics of the radiation -- whether it can be measured with these detectors, and the best way to go about doing so -- is a challenging problem. In this chapter we review the promise of GWs as observational probes, including a discussion of the current state of GW detectors, and discuss the status of work to understand the waves generated by stellar core collapse.

The LIGO (Laser Interferometer Gravitational-Wave Observatory) detectors have just completed their first science run, following many years of planning, research, and development. LIGO is a member of what will be a worldwide network of gravitational-wave observatories, with other members in Europe, Japan, and -- hopefully -- Australia. Plans are rapidly maturing for a low frequency, space-based gravitational-wave observatory: LISA, the Laser Interferometer Space Antenna, to be launched around 2011. The goal of these instruments is to inaugurate the field of {\it gravitational-wave astronomy}: using gravitational-waves as a means of listening to highly relativistic dynamical processes in astrophysics. This review discusses the promise of this field, outlining why gravitational waves are worth pursuing, and what they are uniquely suited to teach us about astrophysical phenomena. We review the current state of the field, both theoretical and experimental, and then highlight some aspects of gravitational-wave science that are particularly exciting (at least to this author).

Gravitational waves generated by the final merger of compact binary systems depend on the structure of the binary's members. If the binary contains neutron stars, measuring such waves can teach us about the properties of matter at extreme densities. Unfortunately, these waves are typically at high frequency where the sensitivity of broad-band detectors is not good. Learning about dense matter from these waves will require networks of broad-band detectors combined with narrow-band detectors that have good sensitivity at high frequencies. This paper presents an algorithm by which a network can be ``tuned'', in accordance with the best available information, in order to most effectively measure merger waves. The algorithm is presented in the context of a toy model that captures the qualitative features of narrow-band detectors and of certain binary neutron star merger wave models. By using what is learned from a sequence of merger measurements, the network can be gradually tuned in order to accurately measure the waves. The number of measurements needed to reach this stage depends upon the waves' signal strength, the number of narrow-band detectors available for the measurement, and the detailed characteristics of the waves that carry the merger information. Future studies will go beyond this toy model, encompassing a more realistic description of both the detectors and the gravitational waves.

We present a new approximate method for constructing gravitational radiation driven inspirals of test-bodies orbiting Kerr black holes. Such orbits can be fully described by a semi-latus rectum $p$, an eccentricity $e$, and an inclination angle $\iota$; or, by an energy $E$, an angular momentum component $L_z$, and a third constant $Q$. Our scheme uses expressions that are exact (within an adiabatic approximation) for the rates of change ($\dot{p}$, $\dot{e}$, $\dot{\iota}$) as linear combinations of the fluxes ($\dot{E}$, $\dot{L_z}$, $\dot{Q}$), but uses quadrupole-order formulae for these fluxes. This scheme thus encodes the exact orbital dynamics, augmenting it with approximate radiation reaction. Comparing inspiral trajectories, we find that this approximation agrees well with numerical results for the special cases of eccentric equatorial and circular inclined orbits, far more accurate than corresponding weak-field formulae for ($\dot{p}$, $\dot{e}$, $\dot{\iota}$). We use this technique to study the inspiral of a test-body in inclined, eccentric Kerr orbits. Our results should be useful tools for constructing approximate waveforms that can be used to study data analysis problems for the future LISA gravitational-wave observatory, in lieu of waveforms from more rigorous techniques that are currently under development.

Gravitational-wave detectors with sensitivities sufficient to measure the radiation from astrophysical sources are rapidly coming into existence. By the end of this decade, there will exist several ground-based instruments in North America, Europe, and Japan, and the joint American-European space-based antenna LISA should be either approaching orbit or in final commissioning in preparation for launch. The goal of these instruments will be to open the field of gravitational-wave astronomy: using gravitational radiation as an observational window on astrophysics and the universe. In this article, we summarize the current status of the various detectors currently being developed, as well as future plans. We also discuss the scientific reach of these instruments, outlining what gravitational-wave astronomy is likely to teach us about the universe.

We derive estimates for the characteristics of gravitational radiation from stellar collapse, using recent models of the core-collapse of Chandrasekhar mass white dwarfs (accretion induced collapse), core-collapse supernovae and collapsars, and the collapse of very massive stars (~> 300 Msun). We study gravitational-wave emission mechanisms using several estimation techniques, including two-dimensional numerical computation of quadrupole wave emission, estimates of bar-mode strength, estimates of r-mode emission, and estimates of waves from black hole ringing. We also review the rate at which the relevant collapses are believed to occur, which has a major impact on their relevance as astrophysical sources. Although the latest supernova progenitor simulations produce cores rotating much slower than those used in the past, we find that bar-mode and r-mode instabilities from core-collapse supernovae remain among the leading candidate sources for LIGO-II. Accretion induced collapse (AIC) of a white dwarf could produce gravitational-wave signals similar to those from core-collapse. In the models that we examine, such collapses are not unstable to bar modes; we note that models recently examined by Liu and Lindblom, which have slightly more angular momentum, are certainly unstable to bar formation. Because AIC events are probably 1,000 times less common than core-collapse supernovae, the typical AIC event will be much further away, and thus the observed waves will be much weaker. In the most optimistic circumstances, we find it may be possible to detect gravitational waves from the collapse of 300 Msun Population III stars.

Binary black hole coalescences emit gravitational waves that will be measurable by the space-based detector it LISA to large redshifts. This suggests that LISA may be able to observe black holes grow and evolve as the universe evolves, mapping the distribution of black hole masses as a function of redshift. An immediate difficulty with this idea is that LISA measures certain redshifted combinations of masses with good accuracy: if a system has some mass parameter m, then LISA measures (1+z)m. This mass-redshift degeneracy makes it difficult to follow the mass evolution. In many cases, LISA will also measure the luminosity distance D of a coalescence accurately. Since cosmological parameters (particularly the mean density, the cosmological constant, and the Hubble constant) are now known with moderate precision, we can obtain z from D and break the degeneracy. This makes it possible to untangle the mass and redshift and to study the mass and merger history of black holes. Mapping the black hole mass distribution could open a window onto an early epoch of structure formation.

We present first results of astrophysically relevant experiments where highly supersonic plasma jets are generated via conically convergent flows. The convergent flows are created by electrodynamic acceleration of plasma in a conical array of fine metallic wires (a modification of the wire array Z-pinch). Stagnation of plasma flow on the axis of symmetry forms a standing conical shock effectively collimating the flow in the axial direction. This scenario is essentially similar to that discussed by Canto\' ~and collaborators as a purely hydrodynamic mechanism for jet formation in astrophysical systems. Experiments using different materials (Al, Fe and W) show that a highly supersonic ($M\sim 20$), well-collimated jet is generated when the radiative cooling rate of the plasma is significant. We discuss scaling issues for the experiments and their potential use for numerical code verification. The experiments also may allow direct exploration of astrophysically relevant issues such as collimation, stability and jet-cloud interactions.

The inspiral of a ``small'' ($\mu \sim 1-100 M_\odot$) compact body into a ``large'' ($M \sim 10^{5-7} M_\odot$) black hole is a key source of gravitational radiation for the space-based gravitational-wave observatory LISA. The waves from such inspirals will probe the extreme strong-field nature of the Kerr metric. In this paper, I investigate the properties of a restricted family of such inspirals (the inspiral of circular, inclined orbits) with an eye toward understanding observable properties of the gravitational waves that they generate. Using results previously presented to calculate the effects of radiation reaction, I assemble the inspiral trajectories (assuming that radiation reacts adiabatically, so that over short timescales the trajectory is approximately geodesic) and calculate the wave generated as the compact body spirals in. I do this analysis for several black hole spins, sampling a range that should be indicative of what spins we will encounter in nature. The spin has a very strong impact on the waveform. In particular, when the hole rotates very rapidly, tidal coupling between the inspiraling body and the event horizon has a very strong influence on the inspiral time scale, which in turn has a big impact on the gravitational wave phasing. The gravitational waves themselves are very usefully described as ``multi-voice chirps'': the wave is a sum of ``voices'', each corresponding to a different harmonic of the fundamental orbital frequencies. Each voice has a rather simple phase evolution. Searching for extreme mass ratio inspirals voice-by-voice may be more effective than searching for the summed waveform all at once.

We present the results of near-infrared observations of extragalactic Cepheids made with the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope. The variables are located in the galaxies IC 1613, IC 4182, M 31, M 81, M 101, NGC 925, NGC 1365, NGC 2090, NGC 3198, NGC 3621, NGC 4496A and NGC 4536. All fields were observed in the F160W bandpass; additional images were obtained in the F110W and F205W filters. Photometry was performed using the DAOPHOT II/ALLSTAR package. Self-consistent distance moduli and color excesses were obtained by fitting Period-Luminosity relations in the H, I and V bands. Our results support the assumption of a standard reddening law adopted by the HST Key Project on the Extragalactic Distance Scale. A companion paper will determine true distance moduli and explore the effects of metallicity on the Cepheid distance scale.

We extend the work developed in previous papers on microlensing with a selection of variable stars. We use the Pixel Method to select variable stars on a set of 2.5 x 10**6 pixel light curves in the LMC Bar presented elsewhere. The previous treatment was done in order to optimise the detection of long timescale variations (larger than a few days) and we further optimise our analysis for the selection of Long Timescale and Long Period Variables (LT&LPV). We choose to perform a selection of variable objects as comprehensive as possible, independent of periodicity and of their position on the colour magnitude diagram. We detail the different thresholds successively applied to the light curves, which allow to produce a catalogue of 632 variable objects. We present a table with the coordinate of each variable, its EROS magnitudes at one epoch and an indicator of blending in both colours, together with a finding chart. A cross-correlation with various catalogues shows that 90% of those variable objects were undetected before, thus enlarging the sample of LT&LPV previously known in this area by a factor of 10. Due to the limitations of both the Pixel Method and the data set, additional data -- namely a longer baseline and near infrared photometry -- are required to further characterise these variable stars, as will be addressed in subsequent papers.

A major focus of much current research in gravitation theory is on understanding how radiation reaction drives the evolution of a binary system, particularly in the extreme mass ratio limit. Such research is of direct relevance to gravitational-wave sources for space-based detectors (such as LISA). We present here a study of the radiative evolution of circular (i.e., constant Boyer-Lindquist coordinate radius), non-equatorial Kerr black hole orbits. Recent theorems have shown that, at least in an adiabatic evolution, such orbits evolve from one circular configuration into another, changing only their radius and inclination angle. This constrains the system's evolution in such a way that the change in its Carter constant can be deduced from knowledge of gravitational wave fluxes propagating to infinity and down the black hole's horizon. Thus, in this particular case, a local radiation reaction force is not needed. In accordance with post-Newtonian weak-field predictions, we find that inclined orbits radiatively evolve to larger inclination angles (although the post-Newtonian prediction overestimates the rate of this evolution in the strong field by a factor $\lesssim 3$). We also find that the gravitational waveforms emitted by these orbits are rather complicated, particularly when the hole is rapidly spinning, as the radiation is influenced by many harmonics of the orbital frequencies.

In this comment, I argue that chaotic effects in binary black hole inspiral will not strongly impact the detection of gravitational waves from such systems.

An unusual set of orbits about extreme Kerr black holes resides at the Boyer-Lindquist radius $r = M$, the coordinate of the hole's event horizon. These ``horizon skimming'' orbits have the property that their angular momentum $L_z$ {\it increases} with inclination angle, opposite to the familiar behavior one encounters at larger radius. In this paper, I show that this behavior is characteristic of a larger family of orbits, the ``nearly horizon skimming'' (NHS) orbits. NHS orbits exist in the very strong field of any black hole with spin $a\agt 0.952412M$. Their unusual behavior is due to the locking of particle motion near the event horizon to the hole's spin, and is therefore a signature of the Kerr metric's extreme strong field. An observational hallmark of NHS orbits is that a small body spiraling into a Kerr black hole due to gravitational-wave emission will be driven into orbits of progressively smaller inclination angle, toward the equator. This is in contrast to the ``normal'' behavior. For circular orbits, the change in inclination is very small, and unlikely to be of observational importance. I argue that the change in inclination may be considerably larger when one considers the evolution of inclined eccentric orbits. If this proves correct, then the gravitational waves produced by evolution through the NHS regime may constitute a very interesting and important probe of the strong-field nature of rotating black holes.

We present here the final results of the Hubble Space Telescope Key Project to measure the Hubble constant. We summarize our method, the results and the uncertainties, tabulate our revised distances, and give the implications of these results for cosmology. The analysis presented here benefits from a number of recent improvements and refinements, including (1) a larger LMC Cepheid sample to define the fiducial period-luminosity (PL) relations, (2) a more recent HST Wide Field and Planetary Camera 2 (WFPC2) photometric calibration, (3) a correction for Cepheid metallicity, and (4) a correction for incompleteness bias in the observed Cepheid PL samples. New, revised distances are given for the 18 spiral galaxies for which Cepheids have been discovered as part of the Key Project, as well as for 13 additional galaxies with published Cepheid data. The new calibration results in a Cepheid distance to NGC 4258 in better agreement with the maser distance to this galaxy. Based on these revised Cepheid distances, we find values (in km/sec/Mpc) of H0 = 71 +/- 2 (random) +/- 6 (systematic) (type Ia supernovae), 71 +/- 2 +/- 7 (Tully-Fisher relation), 70 +/- 5 +/- 6 (surface brightness fluctuations), 72 +/- 9 +/- 7 (type II supernovae), and 82 +/- 6 +/- 9 (fundamental plane). We combine these results for the different methods with 3 different weighting schemes, and find good agreement and consistency with H0 = 72 +/- 8. Finally, we compare these results with other, global methods for measuring the Hubble constant.

In its final year of inspiral, a stellar mass ($1 - 10 M_\odot$) body orbits a massive ($10^5 - 10^7 M_\odot$) compact object about $10^5$ times, spiralling from several Schwarzschild radii to the last stable orbit. These orbits are deep in the massive object's strong field, so the gravitational waves that they produce probe the strong field nature of the object's spacetime. Measuring these waves can, in principle, be used to ``map'' this spacetime, allowing observers to test whether the object is a black hole or something more exotic. Such measurements will require a good theoretical understanding of wave generation during inspiral. In this article, I discuss the major theoretical challenges standing in the way of building such maps from gravitational-wave observations, as well as recent progress in producing extreme mass ratio inspirals and waveforms.

As shown by Teukolsky, the master equation governing the propagation of weak radiation in a black hole spacetime can be separated into four ordinary differential equations, one for each spacetime coordinate. (``Weak'' means the radiation's amplitude is small enough that its own gravitation may be neglected.) Unfortunately, it is difficult to accurately compute solutions to the separated radial equation (the Teukolsky equation), particularly in a numerical implementation. The fundamental reason for this is that the Teukolsky equation's potentials are long ranged. For non-spinning black holes, one can get around this difficulty by applying transformations which relate the Teukolsky solution to solutions of the Regge-Wheeler equation, which has a short-ranged potential. A particularly attractive generalization of this approach to spinning black holes for gravitational radiation (spin weight s = -2) was given by Sasaki and Nakamura. In this paper, I generalize Sasaki and Nakamura's results to encompass radiation fields of arbitrary integer spin weight, and give results directly applicable to scalar (s = 0) and electromagnetic (s = -1) radiation. These results may be of interest for studies of astrophysical radiation processes near black holes, and of programs to compute radiation reaction forces in curved spacetime.

Space-based gravitational-wave interferometers such as LISA will be sensitive to the inspiral of stellar mass compact objects into black holes with masses in the range of roughly 10^5 solar masses to (a few) 10^7 solar masses. During the last year of inspiral, the compact body spends several hundred thousand orbits spiraling from several Schwarzschild radii to the last stable orbit. The gravitational waves emitted from these orbits probe the strong-field region of the black hole spacetime and can make possible high precision tests and measurements of the black hole's properties. Measuring such waves will require a good theoretical understanding of the waves' properties, which in turn requires a good understanding of strong-field radiation reaction and of properties of the black hole's astrophysical environment which could complicate waveform generation. In these proceedings, I review estimates of the rate at which such inspirals occur in the universe, and discuss what is being done and what must be done further in order to calculate the inspiral waveform.

We present a compilation of Cepheid distance moduli and data for four secondary distance indicators that employ stars in the old stellar populations: the planetary nebula luminosity function (PNLF), the globular cluster luminosity function (GCLF), the tip of the red giant branch (TRGB), and the surface brightness fluctuation (SBF) method. The database includes all data published as of July 15, 1999. The main strength of this compilation resides in all data being on a consistent and homogeneous system: all Cepheid distances are derived using the same calibration of the period-luminosity relation, the treatment of errors is consistent for all indicators, measurements which are not considered reliable are excluded. As such, the database is ideal for inter-comparing any of the distance indicators considered, or for deriving a Cepheid calibration to any secondary distance indicator. Specifically, the database includes: 1) Cepheid distances, extinctions and metallicities; 2) apparent magnitudes of the PNLF cutoff; 3) apparent magnitudes and colors of the turnover of the GCLF (both in the V- and B-bands); 4) apparent magnitudes of the TRGB (in the I-band) and V-I colors at and 0.5 magnitudes fainter than the TRGB; 5) apparent surface brightness fluctuation magnitudes I, K', K_short, and using the F814W filter with the HST/WFPC2. In addition, for every galaxy in the database we give reddening estimates from DIRBE/IRAS as well as HI maps, J2000 coordinates, Hubble and T-type morphological classification, apparent total magnitude in B, and systemic velocity. (Abridged)

Since the launch of the Hubble Space Telescope nine years ago, Cepheid distances to 25 galaxies have been determined for the purpose of calibrating secondary distance indicators. A variety of these can now be calibrated, and the accompanying papers by Sakai, Kelson, Ferrarese, and Gibson employ the full set of 25 galaxies to consider the Tully-Fisher relation, the fundamental plane of elliptical galaxies, Type Ia supernovae, and surface brightness fluctuations. When calibrated with Cepheid distances, each of these methods yields a measurement of the Hubble constant and a corresponding measurement uncertainty. We combine these measurements in this paper, together with a model of the velocity field, to yield the best available estimate of the value of H_0 within the range of these secondary distance indicators and its uncertainty. The result is H_0 = 71 +/- 6 km/sec/Mpc. The largest contributor to the uncertainty of this 67% confidence level result is the distance of the Large Magellanic Cloud, which has been assumed to be 50 +/- 3 kpc.

This paper presents the calibration of BVRIH$ Tully-Fisher relations based on Cepheid distances to 21 galaxies within 25 Mpc, and 23 clusters within 10,000 km/s. These relations have been applied to several distant cluster surveys in order to derive a value for the Hubble constant, H0, mainly concentrating on an I-band all-sky survey by Giovanelli and collaborators which consisted of total I magnitudes and 50% linewidth data for ~550 galaxies in 16 clusters. For comparison, we also derive the values of H0 using surveys in B-band and V-band by Bothun and collaborators, and in H-band by Aaronson and collaborators. Careful comparisons with various other databases from literature suggest that the H-band data, whose magnitudes are isophotal magnitudes extrapolated from aperture magnitudes rather than total magnitudes, are subject to systematic uncertainties. Taking a weighted average of the estimates of Hubble constants from four surveys, we obtain H0 = 71 +- 4 (random) +- 7 (systematic) km/s/Mpc. We have also investigated how various systematic uncertainties affect the value of H0 such as the internal extinction correction method used, Tully-Fisher slopes and shapes, a possible metallicity dependence of the Cepheid period-luminosity relation and cluster population incompleteness bias.

Using published photometry and spectroscopy, we construct the fundamental plane and D_n-Sigma relations in Leo I, Virgo and Fornax. The published Cepheid P-L relations to spirals in these clusters fixes the relation between angular size and metric distance for both the fundamental plane and D_n-Sigma relations. Using the locally calibrated fundamental plane, we infer distances to a sample of clusters with a mean redshift of cz \approx 6000 \kms, and derive a value of H_0=78+- 5+- 9 km/s/Mpc (random, systematic) for the local expansion rate. This value includes a correction for depth effects in the Cepheid distances to the nearby clusters, which decreased the deduced value of the expansion rate by 5% +- 5%. If one further adopts the metallicity correction to the Cepheid PL relation, as derived by the Key Project, the value of the Hubble constant would decrease by a further 6%+- 4%. These two sources of systematic error, when combined with a +- 6% error due to the uncertainty in the distance to the Large Magellanic Cloud, a +- 4% error due to uncertainties in the WFPC2 calibration, and several small sources of uncertainty in the fundamental plane analysis, combine to yield a total systematic uncertainty of +- 11%. We find that the values obtained using either the CMB, or a flow-field model, for the reference frame of the distant clusters, agree to within 1%. The Dn-Sigma relation also produces similar results, as expected from the correlated nature of the two scaling relations. A complete discussion of the sources of random and systematic error in this determination of the Hubble constant is also given, in order to facilitate comparison with the other secondary indicators being used by the Key Project.

Using optimal matched filtering, we search 25 hours of data from the LIGO 40-meter prototype laser interferometric gravitational-wave detector for gravitational-wave chirps emitted by coalescing binary systems within our Galaxy. This is the first test of this filtering technique on real interferometric data. An upper limit on the rate R of neutron star binary inspirals in our Galaxy is obtained: with 90% confidence, R< 0.5/hour. Similar experiments with LIGO interferometers will provide constraints on the population of tight binary neutron star systems in the Universe.

The final efforts of the HST Key Project on the Extragalactic Distance Scale are presented. Four distance indicators, the Surface Brightness Fluctuation method, the Fundamental Plane for early-type galaxies, the Tully-Fisher relation and the Type Ia Supernovae, are calibrated using Cepheid distances to galaxies within 25 Mpc. The calibration is then applied to distant samples reaching cz~10000 km/s and (in the case of SNIa) beyond. By combining the constraints imposed on the Hubble constant by the four distance indicators, we obtain H0 = 71+/-6 km/s/Mpc.

A Cepheid-based calibration is derived for four distance indicators that utilize stars in the old stellar populations: the tip of the red giant branch (TRGB), the planetary nebula luminosity function (PNLF), the globular cluster luminosity function (GCLF) and the surface brightness fluctuation method (SBF). The calibration is largely based on the Cepheid distances to 18 spiral galaxies within cz =1500 km/s obtained as part of the HST Key Project on the Extragalactic Distance Scale, but relies also on Cepheid distances from separate HST and ground-based efforts. The newly derived calibration of the SBF method is applied to obtain distances to four Abell clusters in the velocity range between 3800 and 5000 km/s, observed by Lauer et al. (1998) using the HST/WFPC2. Combined with cluster velocities corrected for a cosmological flow model, these distances imply a value of the Hubble constant of H0 = 69 +/- 4 (random) +/- 6 (systematic) km/s/Mpc. This result assumes that the Cepheid PL relation is independent of the metallicity of the variable stars; adopting a metallicity correction as in Kennicutt et al. (1998), would produce a (5 +/- 3)% decrease in H0. Finally, the newly derived calibration allows us to investigate systematics in the Cepheid, PNLF, SBF, GCLF and TRGB distance scales.

Cepheid-based distances to seven Type Ia supernovae (SNe)-host galaxies have been derived using the standard HST Key Project on the Extragalactic Distance Scale pipeline. For the first time, this allows for a transparent comparison of data accumulated as part of three different HST projects, the Key Project, the Sandage et al. Type Ia SNe program, and the Tanvir et al. Leo I Group study. Re-analyzing the Tanvir et al. galaxy and six Sandage et al. galaxies we find a mean (weighted) offset in true distance moduli of 0.12+/-0.07 mag -- i.e., 6% in linear distance -- in the sense of reducing the distance scale, or increasing H0. Adopting the reddening-corrected Hubble relations of Suntzeff et al. (1999), tied to a zero point based upon SNe~1990N, 1981B, 1998bu, 1989B, 1972E and 1960F and the photometric calibration of Hill et al. (1998), leads to a Hubble constant of H0=68+/-2(random)+/-5(systematic) km/s/Mpc. Adopting the Kennicutt et al. (1998) Cepheid period-luminosity-metallicity dependency decreases the inferred H0 by 4%. The H0 result from Type Ia SNe is now in good agreement, to within their respective uncertainties, with that from the Tully-Fisher and surface brightness fluctuation relations.

We report on our numerical implementation of fully relativistic hydrodynamics coupled to Einstein's field equations in three spatial dimensions. We briefly review several steps in our code development, including our recasting of Einstein's equations and several tests which demonstrate its advantages for numerical integrations. We outline our implementation of relativistic hydrodynamics, and present numerical results for the evolution of both stable and unstable Oppenheimer-Volkov equilibrium stars, which represent a very promising first test of our code.

We report on the detection of Cepheids and the first distance measurement to the spiral galaxy NGC 1326-A, a member of the Fornax cluster of galaxies. We have employed data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. Over a 49 day interval, a total of twelve V-band (F555W) and eight I-band (F814W) epochs of observation were obtained. Two photometric reduction packages, ALLFRAME and DoPHOT, have been employed to obtain photometry measures from the three Wide Field CCDs. Variability analysis yields a total of 17 Cepheids in common with both photometry datasets, with periods ranging between 10 and 50 days. Of these 14 Cepheids with high-quality lightcurves are used to fit the V and I period-luminosity relations and derive apparent distance moduli, assuming a Large Magellanic Cloud distance modulus (m-M) (LMC) = 18.50 +- 0.10 mag and color excess E(B-V) = 0.10 mag. Assuming A(V)/E(V-I) = 2.45, the DoPHOT data yield a true distance modulus to NGC 1326-A of (m-M)_0 = 31.36 +- 0.17 (random) +- 0.13 (systematic) mag, corresponding to a distance of 18.7 \pm 1.5 (random) \pm 1.2 (systematic) Mpc. The derived distance to NGC 1326-A is in good agreement with the distance derived previously to NGC 1365, another spiral galaxy member of the Fornax cluster. However the distances to both galaxies are significantly lower than to NGC 1425, a third Cepheid calibrator in the outer parts of the cluster.

The distance to NGC 3319 has been determined from Cepheid variable stars as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen and four epochs of observations, using filters F555W (V) and F814W (I) respectively, were made with the Wide Field Planetary Camera 2. Thirty-three Cepheid variables between periods of 8 and 47 days were discovered. Adopting a Large Magellanic Cloud distance modulus of 18.50 +- 0.10 mag and extinction of E(V-I)=0.13 mag, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.78 +- 0.14 (random) +- 0.10 (systematic) mag and the extinction of E(V-I) = 0.06 mag were determined for NGC 3319. This galaxy is the last galaxy observed for the HST H0 Key Project.

We include matter sources in Einstein's field equations and show that our recently proposed 3+1 evolution scheme can stably evolve strong-field solutions. We insert in our code known matter solutions, namely the Oppenheimer-Volkoff solution for a static star and the Oppenheimer-Snyder solution for homogeneous dust sphere collapse to a black hole, and evolve the gravitational field equations. We find that we can evolve stably static, strong-field stars for arbitrarily long times and can follow dust sphere collapse accurately well past black hole formation. These tests are useful diagnostics for fully self-consistent, stable hydrodynamical simulations in 3+1 general relativity. Moreover, they suggest a successive approximation scheme for determining gravitational waveforms from strong-field sources dominated by longitudinal fields, like binary neutron stars: approximate quasi-equilibrium models can serve as sources for the transverse field equations, which can be evolved without having to re-solve the hydrodynamical equations (``hydro without hydro'').

We report on the discovery of Cepheids in the Virgo spiral galaxy NGC 4535, based on observations made with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. NGC 4535 is one of 18 galaxies observed as a part of The HST Key Project on the Extragalactic Distance Scale which aims to measure the Hubble constant to 10% accuracy. NGC 4535 was observed over 13 epochs using the F555W filter, and over 9 epochs using the F814W filter. The HST F555W and F814W data were transformed to the Johnson V and Kron-Cousins I magnitude systems, respectively. Photometry was performed using two independent programs, DoPHOT and DAOPHOT II/ALLFRAME. Period-luminosity relations in the V and I bands were constructed using 39 high-quality Cepheids present in our set of 50 variable candidates. We obtain a distance modulus of 31.02+/-0.26 mag, corresponding to a distance of 16.0+/-1.9 Mpc. Our distance estimate is based on values of mu = 18.50 +/- 0.10 mag and E(V-I) = 0.13 mag for the distance modulus and reddening of the LMC, respectively.

Using the Hubble Space Telescope (HST) 37 long-period Cepheid variables have been discovered in the Fornax Cluster spiral galaxy NGC 1365. The resulting V and I period-luminosity relations yield a true distance modulus of 31.35 +/- 0.07 mag, which corresponds to a distance of 18.6 +/- 0.6 Mpc. This measurement provides several routes for estimating the Hubble Constant. (1) Assuming this distance for the Fornax Cluster as a whole yields a local Hubble Constant of 70 +/-18_{random} [+/-7]_{systematic} km/s/Mpc. (2) Nine Cepheid-based distances to groups of galaxies out to and including the Fornax and Virgo clusters yield Ho = 73 (+/-16)_r [+/-7]_s km/s/Mpc. (3) Recalibrating the I-band Tully-Fisher relation using NGC 1365 and six nearby spiral galaxies, and applying it to 15 galaxy clusters out to 100 Mpc gives Ho = 76 (+/-3)_r [+/-8]_s km/s/Mpc. (4) Using a broad-based set of differential cluster distance moduli ranging from Fornax to Abell 2147 gives Ho = 72 (+/-)_r [+/-6]_s km/s/Mpc. And finally, (5) Assuming the NGC 1365 distance for the two additional Type Ia supernovae in Fornax and adding them to the SnIa calibration (correcting for light curve shape) gives Ho = 67 (+/-6)_r [+/-7]_s km/s/Mpc out to a distance in excess of 500 Mpc. All five of these Ho determinations agree to within their statistical errors. The resulting estimate of the Hubble Constant combining all these determinations is Ho = 72 (+/-5)_r [+/-12]_s km/s/Mpc.

The distance to NGC 4725 has been derived from Cepheid variables, as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen F555W (V) and four F814W (I) epochs of cosmic-ray-split Wide Field and Planetary Camera 2 observations were obtained. Twenty Cepheids were discovered, with periods ranging from 12 to 49 days. Adopting a Large Magellanic Cloud distance modulus and extinction of 18.50+/-0.10 mag and E(V-I)=0.13 mag, respectively, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.50 +/- 0.16 (random) +/- 0.17 (systematic) mag was determined for NGC 4725. The corresponding of distance of 12.6 +/- 1.0 (random) +/- 1.0 (systematic) Mpc is in excellent agreement with that found with an independent analysis based upon the DoPHOT photometry package. With a foreground reddening of only E(V-I)=0.02, the inferred intrinsic reddening of this field in NGC 4725, E(V-I)=0.19, makes it one of the most highly-reddened, encountered by the HST Key Project, to date.

The Hubble Space Telescope is being used to measure accurate Cepheid distances to nearby galaxies with the ultimate aim of determining the Hubble constant, H_0. For the first time, it has become feasible to use Cepheid variables to derive a distance to a galaxy in the southern hemisphere cluster of Fornax. Based on the discovery of 37 Cepheids in the Fornax galaxy NGC 1365, a distance to this galaxy of 18.6 +/- 0.6 Mpc (statistical error only) is obtained. This distance leads to a value of H_0 = 70 +/- 7 (random) +/- 18 (systematic) km/sec/Mpc in good agreement with estimates of the Hubble constant further afield.

We report the detection of Cepheid variable stars in the barred spiral galaxy NGC 1365, located in the Fornax cluster, using the Hubble Space Telescope Wide Field and Planetary Camera 2. Twelve V (F555W) and four I (F814W) epochs of observation were obtained. The two photometry packages, ALLFRAME and DoPHOT, were separately used to obtain profile-fitting photometry of all the stars in the HST field. The search for Cepheid variable stars resulted in a sample of 52 variables, with periods between 14 and 60 days, in common with both datasets. ALLFRAME photometry and light curves of the Cepheids are presented. A subset of 34 Cepheids were selected on the basis of period, light curve shape, similar ALLFRAME and DoPHOT periods, color, and relative crowding, to fit the Cepheid period-luminosity relations in V and I for both ALLFRAME and DoPHOT. The measured distance modulus to NGC 1365 from the ALLFRAME photometry is 31.31 +/- 0.20 (random) +/- 0.18 (systematic) mag, corresponding to a distance of 18.3 +/- 1.7 (random) +/- 1.6 (systematic) Mpc. The reddening is measured to be E(V-I) = 0.16 +/- 0.08 mag. These values are in excellent agreement with those obtained using the DoPHOT photometry, namely a distance modulus of 31.26 +/- 0.10 mag, and a reddening of 0.15 +/- 0.10 mag (internal errors only).

We report the detection of Cepheids and a new distance to the spiral galaxy NGC 2541, based on data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST). A total of 25 exposures (divided into 13 epochs) are obtained using the F555W filter (transformed to Johnson V), and nine exposures (divided into five epochs) using the F814W filter (transformed to Cousins I). Photometric reduction of the data is performed using two independent packages, DoPHOT and DAOPHOT II/ALLFRAME, which give very good agreement in the measured magnitudes. A total of 34 bona fide Cepheids, with periods ranging from 12 to over 60 days, are identified based on both sets of photometry. By fitting V and I period-luminosity relations, apparent distance moduli are derived assuming a Large Magellanic Cloud distance modulus and mean color excess of 18.50 +/- 0.10 mag and E(B-V) = 0.10 mag respectively. Adopting A(V)/E(V-I)=2.45, we obtain a true distance modulus to NGC 2541 of 30.47 +/- 0.11 (random) +/- 0.12 (systematic) mag (D = 12.4 +/- 0.6 (random) +/- 0.7 (systematic) Mpc), and a total (Galactic plus internal) mean color excess E(B-V) = 0.08 +/- 0.05 (internal error) mag.

The distance to NGC 7331 has been derived from Cepheid variables observed with HST/WFPC2, as part of the Extragalactic Distance Scale Key Project. Multi-epoch exposures in F555W (V) and F814W (I), with photometry derived independently from DoPHOT and DAOPHOT/ALLFRAME programs, were used to detect a total of 13 reliable Cepheids, with periods between 11 and 42 days. The relative distance moduli between NGC 7331 and the LMC, imply an extinction to NGC 7331 of A_V = 0.47+-0.15 mag, and an extinction-corrected distance modulus to NGC 7331 of 30.89+-0.14(random) mag, equivalent to a distance of 15.1 Mpc. There are additional systematic uncertainties in the distance modulus of +-0.12 mag due to the calibration of the Cepheid Period-Luminosity relation, and a systematic offset of +0.05+-0.04 mag if we applied the metallicity correction inferred from the M101 results of Kennicutt et al 1998.

Uncertainty in the metal abundance dependence of the Cepheid variable period- luminosity (PL) relation remains one of the outstanding sources of systematic error in the extragalactic distance scale and Hubble constant. To test for such a metallicity dependence, we have used the WFPC2 camera on the Hubble Space Telescope (HST) to observe Cepheids in two fields in the nearby spiral galaxy M101, which span a range in oxygen abundance of 0.7+-0.15 dex. A differential analysis of the PL relations in V and I in the two fields yields a marginally significant change in the inferred distance modulus on metal abundance, with d(m-M)/d[O/H] = -0.24+-0.16 mag/dex. The trend is in the theoretically predicted sense that metal-rich Cepheids appear brighter and closer than metal-poor stars. External comparisons of Cepheid distances with those derived from three other distance indicators, in particular the tip of the red giant branch method, further constrain the magnitude of any Z-dependence of the PL relation at V and I. The overall effects of any metallicity dependence on the distance scale derived with HST will be of the order of a few percent or less for most applications, though distances to individual galaxies at the extremes of the metal abundance range may be affected at the 10% level.

We estimate the expected signal-to-noise ratios (SNRs) from the three phases (inspiral,merger,ringdown) of coalescing binary black holes (BBHs) for initial and advanced ground-based interferometers (LIGO/VIRGO) and for space-based interferometers (LISA). LIGO/VIRGO can do moderate SNR (a few tens), moderate accuracy studies of BBH coalescences in the mass range of a few to about 2000 solar masses; LISA can do high SNR (of order 10^4) high accuracy studies in the mass range of about 10^5 to 10^8 solar masses. BBHs might well be the first sources detected by LIGO/VIRGO: they are visible to much larger distances (up to 500 Mpc by initial interferometers) than coalescing neutron star binaries (heretofore regarded as the "bread and butter" workhorse source for LIGO/VIRGO, visible to about 30 Mpc by initial interferometers). Low-mass BBHs (up to 50 solar masses for initial LIGO interferometers; 100 for advanced; 10^6 for LISA) are best searched for via their well-understood inspiral waves; higher mass BBHs must be searched for via their poorly understood merger waves and/or their well-understood ringdown waves. A matched filtering search for massive BBHs based on ringdown waves should be capable of finding BBHs in the mass range of about 100 to 700 solar masses out to 200 Mpc (initial LIGO interferometers), and 200 to 3000 solar masses out to about z=1 (advanced interferometers). The required number of templates is of order 6000 or less. Searches based on merger waves could increase the number of detected massive BBHs by a factor of order 10 or more over those found from inspiral and ringdown waves, without detailed knowledge of the waveform shapes, using a "noise monitoring" search algorithm. A full set of merger templates from numerical relativity could further increase the number of detected BBHs by an additional factor of up to 4.

We discuss the extraction of information from detected binary black hole (BBH) coalescence gravitational waves, focusing on the merger phase that occurs after the gradual inspiral and before the ringdown. Our results are: (1) If numerical relativity simulations have not produced template merger waveforms before BBH detections by LIGO/VIRGO, one can band-pass filter the merger waves. For BBHs smaller than about 40 solar masses detected via their inspiral waves, the band pass filtering signal to noise ratio indicates that the merger waves should typically be just barely visible in the noise for initial and advanced LIGO interferometers. (2) We derive an optimized (maximum likelihood) method for extracting a best-fit merger waveform from the noisy detector output; one "perpendicularly projects" this output onto a function space (specified using wavelets) that incorporates our prior knowledge of the waveforms. An extension of the method allows one to extract the BBH's two independent waveforms from outputs of several interferometers. (3) If numerical relativists produce codes for generating merger templates but running the codes is too expensive to allow an extensive survey of the merger parameter space, then a coarse survey of this parameter space, to determine the ranges of the several key parameters and to explore several qualitative issues which we describe, would be useful for data analysis purposes. (4) A complete set of templates could be used to test the nonlinear dynamics of general relativity and to measure some of the binary parameters. We estimate the number of bits of information obtainable from the merger waves (about 10 to 60 for LIGO/VIRGO, up to 200 for LISA), estimate the information loss due to template numerical errors or sparseness in the template grid, and infer approximate requirements on template accuracy and spacing.

We report on the discovery of Cepheids in the field spiral galaxy NGC 3621, based on observations made with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST). NGC 3621 is one of 18 galaxies observed as a part of The HST Key Project on the Extragalactic Distance Scale, which aims to measure the Hubble constant to 10% accuracy. Sixty-nine Cepheids with periods in the range 9--60 days were observed over 12 epochs using the F555W filter, and 4 epochs using the F814W filter. The HST F555W and F814W data were transformed to the Johnson V and Kron-Cousins I magnitude systems, respectively. Photometry was performed using two independent packages, DAOPHOT II/ALLFRAME and DoPHOT. Period-luminosity relations in the V and I bands were constructed using 36 fairly isolated Cepheids present in our set of 69 variables. Extinction-corrected distance moduli relative to the LMC of 10.63 +/- 0.07 mag and 10.56 +/- 0.10 mag were obtained using the ALLFRAME and DoPHOT data, respectively. True distance moduli of 29.13 +/- 0.18 mag and 29.06 +/- 0.18 mag, corresponding to distances of 6.3 Mpc and 6.1 Mpc, were obtained by assuming values of 18.50 +/- 0.10 mag for the distance modulus of the LMC and E(V-I) = 0.13 mag for the reddening of the LMC.

Recent numerical work by Wilson, Mathews, and Marronetti [J. R. Wilson, G. J. Mathews and P. Marronetti, Phys. Rev. D 54, 1317 (1996)] on the coalescence of massive binary neutron stars shows a striking instability as the stars come close together: Each star's central density increases by an amount proportional to 1/(orbital radius). This overwhelms any stabilizing effects of tidal coupling [which are proportional to 1/(orbital radius)^6] and causes the stars to collapse before they merge. Since the claimed increase of density scales with the stars' mass, it should also show up in a perturbation limit where a point particle of mass $\mu$ orbits a neutron star. We prove analytically that this does not happen; the neutron star's central density is unaffected by the companion's presence to linear order in $\mu/R$. We show, further, that the density increase observed by Wilson et. al. could arise as a consequence of not faithfully maintaining boundary conditions.