Source author record

Scott A. Hughes

Scott A. Hughes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2022arXiv

Measuring quasi-normal mode amplitudes with misaligned binary black hole ringdowns

In recent work, we examined how different modes in the ringdown phase of a binary coalescence are excited as a function of the final plunge geometry. At least in the large mass ratio limit, we found a clean mapping between angles describing the plunge and the amplitude of different quasi-normal modes (QNMs) which constitute the ringdown. In this study, we use that mapping to construct a waveform model expressed as a sum of QNMs where the mode amplitudes and phases are determined by the source plunge parameters. We first generate a large number of calibration waveforms and interpolate between fits of each mode amplitude and phase up to $\ell \leq 8$ and $\ell - |m| \leq 4$. The density of our calibration data allows us to resolve important features such as phase transition discontinuities at large misalignments. Using our ringdown waveform model, we then perform Bayesian parameter estimation with added white Gaussian noise to demonstrate that, in principle, the mode amplitudes can be measured and used to constrain the plunge geometry. We find that inferences are substantially improved by incorporating prior information constraining mode excitation, which motivates work to understand and characterize how the QNM excitation depends on the coalescence geometry. These results are part of a broader effort to map the mode excitation from arbitrary masses and spins, which will be useful for characterizing ringdown waves in upcoming gravitational-wave measurements.

preprint2022arXiv

Precisely computing bound orbits of spinning bodies around black holes I: General framework and results for nearly equatorial orbits

Very large mass ratio binary black hole systems are of interest both as a clean limit of the two-body problem in general relativity, as well as for their importance as sources of low-frequency gravitational waves. At lowest order, the smaller body moves along a geodesic of the larger black hole's spacetime. Post-geodesic effects include the gravitational self force, which incorporates the backreaction of gravitational-wave emission, and the spin-curvature force, which arises from coupling of the small body's spin to the black hole's spacetime curvature. In this paper, we describe a method for precisely computing bound orbits of spinning bodies about black holes. Our analysis builds off of pioneering work by Witzany which demonstrated how to describe the motion of a spinning body to linear order in the small body's spin. Exploiting the fact that in the large mass-ratio limit spinning-body orbits are close to geodesics and using closed-form results due to van de Meent describing precession of the small body's spin along black hole orbits, we develop a frequency-domain formulation of the motion which can be solved very precisely. We examine a range of orbits with this formulation, focusing in this paper on orbits which are eccentric and nearly equatorial (i.e., the orbit's motion is $\mathcal{O}(S)$ out of the equatorial plane), but for which the small body's spin is arbitrarily oriented. We discuss generic orbits with general small-body spin orientation in a companion paper. We characterize the behavior of these orbits and show how the small body's spin shifts the frequencies $Ω_r$ and $Ω_ϕ$ which affect orbital motion. These frequency shifts change accumulated phases which are direct gravitational-wave observables, illustrating the importance of precisely characterizing these quantities for gravitational-wave observations. (Abridged)

preprint2022arXiv

Precisely computing bound orbits of spinning bodies around black holes II: Generic orbits

In this paper, we continue our study of the motion of spinning test bodies orbiting Kerr black holes. Non-spinning test bodies follow geodesics of the spacetime in which they move. A test body's spin couples to the curvature of that spacetime, introducing a "spin-curvature force" which pushes the body's worldline away from a geodesic trajectory. The spin-curvature force is an important example of a post-geodesic effect which must be modeled carefully in order to accurately characterize the motion of bodies orbiting black holes. One motivation for this work is to understand how to include such effects in models of gravitational waves produced from the inspiral of stellar mass bodies into massive black holes. In this paper's predecessor, we describe a technique for computing bound orbits of spinning bodies around black holes with a frequency-domain description which can be solved very precisely. In that paper, we present an overview of our methods, as well as present results for orbits which are eccentric and nearly equatorial (i.e., the orbit's motion is no more than $\mathcal{O}(S)$ out of the equatorial plane). In this paper, we apply this formulation to the fully generic case -- orbits which are inclined and eccentric, with the small body's spin arbitrarily oriented. We compute the trajectories which such orbits follow, and compute how the small body's spin affects important quantities such as the observable orbital frequencies $Ω_r$, $Ω_θ$ and $Ω_ϕ$.

preprint2021arXiv

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI waveforms. Because the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, we also introduce an improved Augmented Analytic Kludge (AAK) model that describes generic Kerr inspirals. Both waveform models can be accelerated using graphics processing unit (GPU) hardware. With the GPU-accelerated waveforms in hand, a variety of studies are performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation. We find relativistic EMRI waveform templates can be generated with fewer harmonic modes ($\sim10-100$) without biasing signal extraction. However, we show for the first time that extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space.

preprint2020arXiv

Prospects for Fundamental Physics with LISA

In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.

preprint2020arXiv

Rapid generation of fully relativistic extreme-mass-ratio-inspiral waveform templates for LISA data analysis

The future space mission LISA will observe a wealth of gravitational-wave sources at millihertz frequencies. Of these, the extreme-mass-ratio inspirals of compact objects into massive black holes are the only sources that combine the challenges of strong-field complexity with that of long-lived signals. Such signals are found and characterized by comparing them against a large number of accurate waveform templates during data analysis, but the rapid generation of such templates is hindered by computing the $\sim10^3$-$10^5$ harmonic modes in a fully relativistic waveform. We use order-reduction and deep-learning techniques to derive a global fit for these modes, and implement it in a complete waveform framework with hardware acceleration. Our high-fidelity waveforms can be generated in under $1\,\mathrm{s}$, and achieve a mismatch of $\lesssim 5\times 10^{-4}$ against reference waveforms that take $\gtrsim 10^4$ times longer. This marks the first time that analysis-length waveforms with full harmonic content can be produced on timescales useful for direct implementation in LISA analysis algorithms.

preprint2020arXiv

Tidal heating as a discriminator for horizons in extreme mass ratio inspirals

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.

preprint2016arXiv

Adiabatic and post-adiabatic approaches to extreme mass ratio inspiral

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.

preprint2016arXiv

Gyroscopes orbiting black holes: A frequency-domain approach to precession and spin-curvature coupling for spinning bodies on generic Kerr 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.

preprint2016arXiv

Inspiral into Gargantua

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.

preprint2016arXiv

Strong-field tidal distortions of rotating black holes: II. Horizon dynamics from eccentric and inclined orbits

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.

preprint2014arXiv

A census of transient orbital resonances encountered during binary inspiral

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, $Ω_θ$ and $Ω_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.

preprint2014arXiv

Gravitational wave astronomy and cosmology

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.

preprint2014arXiv

Research Update on Extreme-Mass-Ratio Inspirals

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.

preprint2014arXiv

Resonantly enhanced and diminished strong-field gravitational-wave fluxes

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 $θ$ and $r$ motions are in an integer ratio: $Ω_θ/Ω_r = β_θ/β_r$, with $β_θ$ and $β_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.

preprint2014arXiv

Small mass plunging into a Kerr black hole: Anatomy of the inspiral-merger-ringdown waveforms

We numerically solve the Teukolsky equation in the time domain to obtain the gravitational-wave emission of a small mass inspiraling and plunging into the equatorial plane of a Kerr black hole. We account for the dissipation of orbital energy using the Teukolsky frequency-domain gravitational-wave fluxes for circular, equatorial orbits, down to the light-ring. We consider Kerr spins $-0.99 \leq q \leq 0.99$, and compute the inspiral-merger-ringdown (2,2), (2,1), (3,3), (3,2), (4,4), and (5,5) modes. We study the large-spin regime, and find a great simplicity in the merger waveforms, thanks to the extremely circular character of the plunging orbits. We also quantitatively examine the mixing of quasinormal modes during the ringdown, which induces complicated amplitude and frequency modulations in the waveforms. Finally, we explain how the study of small mass-ratio black-hole binaries helps extending effective-one-body models for comparable-mass, spinning black-hole binaries to any mass ratio and spin magnitude.

preprint2014arXiv

Strong-field tidal distortions of rotating black holes: Formalism and results for circular, equatorial orbits

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.

preprint2013arXiv

Black hole binary inspiral and trajectory dominance

Gravitational waves emitted during the inspiral, plunge and merger of a black hole binary carry linear momentum. This results in an astrophysically important recoil to the final merged black hole, a ``kick'' that can eject it from the nucleus of a galaxy. In a previous paper we showed that the puzzling partial cancellation of an early kick by a late antikick, and the dependence of the cancellation on black hole spin, can be understood from the phenomenology of the linear momentum waveforms. Here we connect that phenomenology to its underlying cause, the spin-dependence of the inspiral trajectories. This insight suggests that the details of plunge can be understood more broadly with a focus on inspiral trajectories.

preprint2013arXiv

Determining the Hubble constant from gravitational wave observations of merging compact binaries

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.

preprint2013arXiv

Modeling the horizon-absorbed gravitational flux for equatorial-circular orbits in Kerr spacetime

We propose an improved analytical model for the horizon-absorbed gravitational-wave energy flux of a small body in circular orbit in the equatorial plane of a Kerr black hole. Post-Newtonian (PN) theory provides an analytical description of the multipolar components of the absorption flux through Taylor expansions in the orbital frequency. Building on previous work, we construct a mode-by-mode factorization of the absorbed flux whose Taylor expansion agrees with current PN results. This factorized form significantly improves the agreement with numerical results obtained with a frequency-domain Teukolsky code, which evolves through a sequence of circular orbits up to the photon orbit. We perform the comparison between model and numerical data for dimensionless Kerr spins $-0.99 \leq q \leq 0.99$ and for frequencies up to the light ring of the Kerr black hole. Our proposed model enforces the presence of a zero in the flux at an orbital frequency equal to the frequency of the horizon, as predicted by perturbation theory. It also reproduces the expected divergence of the flux close to the light ring. Neither of these features are captured by the Taylor-expanded PN flux. Our proposed absorption flux can also help improve models for the inspiral, merger, ringdown of small mass-ratio binary systems.

preprint2012arXiv

Modeling multipolar gravitational-wave emission from small mass-ratio mergers

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.

preprint2011arXiv

Advanced localization of massive black hole coalescences with LISA

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.

preprint2011arXiv

Measuring parameters of massive black hole binaries with partially aligned spins

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.

preprint2011arXiv

Systematics of black hole binary inspiral kicks and the slowness approximation

During the inspiral and merger of black holes, the interaction of gravitational wave multipoles carries linear momentum away, thereby providing an astrophysically important recoil, or "kick" to the system and to the final black hole remnant. It has been found that linear momentum during the last stage (quasinormal ringing) of the collapse tends to provide an "antikick" that in some cases cancels almost all the kick from the earlier (quasicircular inspiral) emission. We show here that this cancellation is not due to peculiarities of gravitational waves, black holes, or interacting multipoles, but simply to the fact that the rotating flux of momentum changes its intensity slowly. We show furthermore that an understanding of the systematics of the emission allows good estimates of the net kick for numerical simulations started at fairly late times, and is useful for understanding qualitatively what kinds of systems provide large and small net kicks.

preprint2010arXiv

Binary black hole merger gravitational waves and recoil in the large mass ratio limit

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.

preprint2010arXiv

Binary Pulsar Constraints on the Parameterized post-Einsteinian Framework

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.

preprint2010arXiv

Exploring short gamma-ray bursts as gravitational-wave standard sirens

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)

preprint2010arXiv

Extreme Mass-Ratio Inspirals in the Effective-One-Body Approach: Quasi-Circular, Equatorial Orbits around a Spinning Black Hole

We construct effective-one-body waveform models suitable for data analysis with LISA for extreme-mass ratio inspirals in quasi-circular, equatorial orbits about a spinning supermassive black hole. The accuracy of our model is established through comparisons against frequency-domain, Teukolsky-based waveforms in the radiative approximation. The calibration of eight high-order post-Newtonian parameters in the energy flux suffices to obtain a phase and fractional amplitude agreement of better than 1 radian and 1 % respectively over a period between 2 and 6 months depending on the system considered. This agreement translates into matches higher than 97 % over a period between 4 and 9 months, depending on the system. Better agreements can be obtained if a larger number of calibration parameters are included. Higher-order mass ratio terms in the effective-one-body Hamiltonian and radiation-reaction introduce phase corrections of at most 30 radians in a one year evolution. These corrections are usually one order of magnitude larger than those introduced by the spin of the small object in a one year evolution. These results suggest that the effective-one-body approach for extreme mass ratio inspirals is a good compromise between accuracy and computational price for LISA data analysis purposes.

preprint2010arXiv

Modeling Extreme Mass Ratio Inspirals within the Effective-One-Body Approach

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.

preprint2010arXiv

On the falloff of radiated energy in black hole spacetimes

The goal of much research in relativity is to understand gravitational waves generated by a strong-field dynamical spacetime. Quantities of particular interest for many calculations are the Weyl scalar $ψ_4$, which is simply related to the flux of gravitational waves far from the source, and the flux of energy carried to distant observers, $\dot E$. Conservation laws guarantee that, in asympotically flat spacetimes, $ψ_4 \propto 1/r$ and $\dot E \propto 1/r^2$ as $r \to \infty$. Most calculations extract these quantities at some finite extraction radius. An understanding of finite radius corrections to $ψ_4$ and $\dot E$ allows us to more accurately infer their asymptotic values from a computation. In this paper, we show that, if the final state of the system is a black hole, then the leading correction to $ψ_4$ is ${\cal O}(1/r^3)$, and that to the energy flux is ${\cal O}(1/r^4)$ --- not ${\cal O}(1/r^2)$ and ${\cal O}(1/r^3)$ as one might naively guess. Our argument only relies on the behavior of the curvature scalars for black hole spacetimes. Using black hole perturbation theory, we calculate the corrections to the leading falloff, showing that it is quite easy to correct for finite extraction radius effects.

preprint2010arXiv

Probing strong-field gravity and black holes with gravitational waves

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.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

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?

preprint2009arXiv

Confusion background from compact binaries

Double neutron stars are one of the most promizing sources for terrestrial gravitational wave interferometers. For actual interferometers and their planned upgrades, the probability of having a signal present in the data is small, but as the sensitivity improves, the detection rate increases and the waveforms may start to overlap, creating a confusion background, ultimately limiting the capabilities of future detectors. The third generation Einstein Telescope, with an horizon of $z > 1$ and very low frequency ``seismic wall'' may be affected by such 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.

preprint2009arXiv

Gravitational waves from merging compact binaries

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.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

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.

preprint2009arXiv

Spacetime and orbits of bumpy black holes

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.

preprint2008arXiv

"Kludge" gravitational waveforms for a test-body orbiting a Kerr black hole

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.

preprint2008arXiv

Measuring coalescing massive binary black holes with gravitational waves: The impact of spin-induced precession

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.

preprint2005arXiv

Rotating black hole orbit functionals in the frequency domain

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 $θ$ 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 $λ$. In particular, the observer's time $t$ is decomposed with $λ$. The frequency domain description is then built from the $λ$-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 $ψ_4$ used in black hole perturbation theory.

preprint2005arXiv

Using gravitational-wave standard sirens

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.

preprint2004arXiv

Consequences of gravitational radiation recoil

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.

preprint2004arXiv

Gravitational Waves from Stellar Collapse: Correlations to Explosion Asymmetries

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.

preprint2002arXiv

Listening to the Universe with Gravitational-Wave Astronomy

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

preprint2001arXiv

Gravitational Wave Emission From Core-Collapse of Massive Stars

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.

preprint2001arXiv

The evolution of circular, non-equatorial orbits of Kerr black holes due to gravitational-wave emission

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.

preprint2001arXiv

Untangling the merger history of massive black holes with LISA

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.

preprint2000arXiv

Computing radiation from Kerr black holes: Generalization of the Sasaki-Nakamura equation

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.

preprint1997arXiv

Measuring gravitational waves from binary black hole coalescences: II. the waves' information and its extraction, with and without templates

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.

preprint1997arXiv

The central density of a neutron star is unaffected by a binary companion at linear order in $μ/R$

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 $μ$ 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 $μ/R$. We show, further, that the density increase observed by Wilson et. al. could arise as a consequence of not faithfully maintaining boundary conditions.