Search

The Boundary-to-Bound (B2B) correspondence, which connects orbital and radiative observables between bound and unbound orbits, has recently been introduced and demonstrated in the perturbative regime. We produce a large number of numerical relativity simulations of bound and unbound encounters between two nonspinning equal mass black holes in order to test this correspondence in the non-perturbative regime. We focus on testing the radiated energy and angular momentum, as well as orbital parameters such as the period and periastron advance. We find that, across a wide range of eccentricities, the B2B relationships do not hold in the non-perturbative regime, thereby placing a clear limit on the applicability of these relationships. We also approximate the separatrix between bound and unbound relativistic encounters as a function of their initial energies and angular momenta., Comment: published in PRD

We present an analytical and numerical study of a system composed of a stellar binary pair and a massless, locally isothermal viscous accretion disk that is coplanar to the binary orbital plane. Analytically, we study the effect of the binary's gravitational potential over short timescales through the study of stability for epicyclic orbits, and over long timescales by revisiting the concept of resonant torques. Numerically, we perform two-dimensional Newtonian numerical simulations of the disk-binary system over a range of binary mass ratios. We find that the results of our simulations are consistent with previous numerical studies. We additionally show, by comparison of the analytical and numerical results, that the circumbinary gap is maintained on the orbital timescale through the driving of epicyclic instabilities, and does not depend on resonant torquing, contrary to standard lore. While our results are applicable to any disk-binary system, we highlight the importance of this result in the search for electromagnetic and gravitational-wave signatures from supermassive black-hole binaries.

Black-holes are known to span at least 9 orders of magnitude in mass: from the stellar-mass objects observed by the Laser Interferometer Gravitational-Wave Observatory Scientific Collaboration and Virgo Collaboration, to supermassive black-holes like the one observed by the Event Horizon Telescope at the heart of M87. Regardless of the mass scale, all of these objects are expected to form binaries and eventually emit observable gravitational radiation, with more massive objects emitting at ever lower gravitational-wave frequencies. We present the tool, gwent, for modelling the sensitivities of current and future generations of gravitational wave detectors across the entire gravitational-wave spectrum of coalescing black-hole binaries. We provide methods to generate sensitivity curves for pulsar timing arrays (PTAs) using a novel realistic PTA sensitivity curve generator, space-based interferometers using adaptive models that can represent a wide range of proposed detector designs, and ground-based interferometers using realistic noise models that can reproduce current, second, and third generation designs, as well as novel variations of the essential design parameters. To model the signal from black-hole binaries at any mass scale, we use phenomenological waveforms capable of modelling the inspiral, merger, and ringdown for sources with varying mass ratios and spins. Using this adaptable framework, we produce signal-to-noise ratios for the combination of any modelled parameter, associated with either the detector or the source. By allowing variation across each detector and source parameter, we can pinpoint the most important factors to determining the optimal performance for particular instrument designs. The adaptability of our detector and signal models can easily be extended to new detector designs and other models of gravitational wave signals., Comment: 42 pages, 16 figures, Our code, gwent, is available at https://github.com/ark0015/gwent and on PyPI https://pypi.org/project/gwent/ Accepted by CQG 15 December 2020

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

When two relativistically boosted, nonspinning black holes pass by one another on a scattering trajectory, we might expect the tidal interaction to spin up each black hole. We present the first exploration of this effect, appearing at fourth post-Newtonian order, with full numerical relativity calculations. The basic setup for the calculations involves two free parameters: the initial boost of each black hole and the initial angle between the velocity vectors and a line connecting the centers of the black holes, with zero angle corresponding to a head-on trajectory. To minimize gauge effects, we measure final spins only if the black holes reach a final separation of at least $20M$. Fixing the initial boost, we find that as the initial angle decreases toward the scattering/nonscattering limit, the spin-up grows nonlinearly. In addition, as initial boosts are increased from $0.42c$ to $0.78c$, the largest observed final dimensionless spin on each black hole increases nonlinearly from $0.02$ to $0.20$. Based on these results, we conclude that much higher spin-ups may be possible with larger boosts, although achieving this will require improved numerical techniques., Comment: 8 pages, 5 figures, matches published version

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA., Comment: Astro2020 APC White Paper Primary Thematic Science Area: Multi-Messenger Astronomy and Astrophysics Secondary Areas: Cosmology and Fundamental Physics, Galaxy Evolution, Formation and Evolution of Compact Objects

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps., Comment: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics). v2: fixed a reference

While still in its infancy, the budding field of gravitational-wave astronomy has so far exceeded most expectations, and the achievements that have already been made bode well for the decade to come. While the discoveries made possible by LIGO have captured the imagination of experts and nonexperts alike, it is important when looking ahead to consider those discoveries in the context of the field as a whole. Just as radio, optical, and x-ray radiation probe different physical phenomena occurring on a range of length and energy scales, the future of gravitational-wave astrophysics depends on our ability to open up the entire spectrum. We will describe the scientific prospects for the field of gravitational-wave astronomy as a whole as we enter the coming decade, and we will place the specific contributions from a future space-based gravitational-wave observatory within this context., Comment: submission to the 2020-2030 Astronomy and Astrophysics Decadal Survey (Astro2020)

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center., Comment: Submitted to the Astro2020 Decadal Survey call for science white papers

The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$., Comment: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics)

We present a highly accurate, fully analytical model for the late inspiral, merger, and ringdown of black-hole binaries with arbitrary mass ratios and spin vectors, including the contributions of harmonics beyond the fundamental mode. This model assumes only that nonlinear effects remain small throughout the entire coalescence, and is developed based on a physical understanding of the dynamics of late stage binary evolution, in particular on the tendency of the dynamical binary spacetime to behave like a linear perturbation of the static merger-remnant spacetime, even at times before the merger has occurred. We demonstrate that our model agrees with the most accurate numerical relativity results to within their own uncertainties throughout the merger-ringdown phase, and it does so for example cases spanning the full range of binary parameter space that is currently testable with numerical relativity. Furthermore, our model maintains accuracy back to the innermost stable circular orbit of the merger-remnant spacetime over much of the relevant parameter space, greatly decreasing the need to introduce phenomenological degrees of freedom to describe the late inspiral., Comment: Replaced with version accepted to PRL. Includes Supplementary Material with additional content, in particular an illustration of the physical merger model, that was requested by the referees

When a gravitational wave is detected by Advanced LIGO/Virgo, sophisticated parameter estimation (PE) pipelines spring into action. These pipelines leverage approximants to generate large numbers of theoretical gravitational waveform predictions to characterize the detected signal. One of the most accurate and physically comprehensive classes of approximants in wide use is the "Spinning Effective One Body--Numerical Relativity" (SEOBNR) family. Waveform generation with these approximants can be computationally expensive, which has limited their usefulness in multiple data analysis contexts. In prior work we improved the performance of the aligned-spin approximant SEOBNR version 2 (v2) by nearly 300x. In this work we focus on optimizing the full eight-dimensional, precessing approximant SEOBNR version 3 (v3). While several v2 optimizations were implemented during its development, v3 is far too slow for use in state-of-the-art source characterization efforts for long-inspiral detections. Completion of a PE run after such a detection could take centuries to complete using v3. Here we develop and implement a host of optimizations for v3, calling the optimized approximant v3_Opt. Our optimized approximant is about 340x faster than v3, and generates waveforms that are numerically indistinguishable., Comment: 16 pages, 1 figure, 9 tables; updated to clarify science background and fix typos

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals., Comment: Published in the Astrophysical Journal Supplement Series

We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to make possible large-scale simulations of magnetized plasmas in strong, dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate that GiRaFFE passes a large suite of both flat and curved-spacetime code tests passed by a number of other state-of-the-art GRFFE codes, and is thus ready for production-scale simulations of GRFFE phenomena of key interest to relativistic astrophysics., Comment: 23 pages, 4 figures. Consistent with published version

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA 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 infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz., Comment: Submitted to ESA on January 13th in response to the call for missions for the L3 slot in the Cosmic Vision Programme

The Spinning Effective One Body-Numerical Relativity (SEOBNR) series of gravitational wave approximants are among the best available for Advanced LIGO data analysis. Unfortunately, SEOBNR codes as they currently exist within LALSuite are generally too slow to be directly useful for standard Markov-Chain Monte Carlo-based parameter estimation (PE). Reduced-Order Models (ROMs) of SEOBNR have been developed for this purpose, but there is no known way to make ROMs of the full eight-dimensional intrinsic parameter space more efficient for PE than the SEOBNR codes directly. So as a proof of principle, we have sped up the original LALSuite SEOBNRv2 approximant code, which models waveforms from aligned-spin systems, by nearly 300x. Our optimized code shortens the timescale for conducting PE with this approximant to months, assuming a purely serial analysis, so that even modest parallelization combined with our optimized code will make running the full PE pipeline with SEOBNR codes directly a realistic possibility. A number of our SEOBNRv2 optimizations have already been applied to SEOBNRv3, a new approximant capable of modeling sources with all eight (precessing) intrinsic degrees of freedom. We anticipate that once all of our optimizations have been applied to SEOBNRv3, a similar speed-up may be achieved., Comment: 17 pages, 5 figures. Matches published version

We present a detailed analysis of the expected signal-to-noise ratios of supermassive black hole binaries on eccentric orbits observed by pulsar timing arrays. We derive several analytical relations that extend the results of Peters and Mathews [Phys. Rev. D 131, 435 (1963)] to quantify the impact of eccentricity in the detection of single resolvable binaries in the pulsar timing array band. We present ready-to-use expressions to compute the increase/loss in signal-to-noise ratio of eccentric single resolvable sources whose dominant harmonic is located in the low/high frequency sensitivity regime of pulsar timing arrays. Building upon the work of Phinney (arXiv:astro-ph/0108028) and Enoki and Nagashima [Prog. Theor. Phys. 117, 241 (2007)], we present an analytical framework that enables the construction of rapid spectra for a stochastic gravitational-wave background generated by a cosmological population of eccentric sources. We confirm previous findings which indicate that, relative to a population of quasicircular binaries, the strain of a stochastic, isotropic gravitational-wave background generated by a cosmological population of eccentric binaries will be suppressed in the frequency band of pulsar timing arrays. We quantify this effect in terms of signal-to-noise ratios in a pulsar timing array., Comment: 18 pages, 9 figures, 3 Appendices. Submitted to PRD. v2: typos corrected, references added. Accepted to PRD

The detection of a stochastic gravitational-wave signal from the superposition of many inspiraling supermassive black holes with pulsar timing arrays (PTAs) is likely to occur within the next decade. With this detection will come the opportunity to learn about the processes that drive black-hole-binary systems toward merger through their effects on the gravitational-wave spectrum. We use Bayesian methods to investigate the extent to which effects other than gravitational-wave emission can be distinguished using PTA observations. We show that, even in the absence of a detection, it is possible to place interesting constraints on these dynamical effects for conservative predictions of the population of tightly bound supermassive black-hole binaries. For instance, if we assume a relatively weak signal consistent with a low number of bound binaries and a low black-hole-mass to galaxy-mass correlation, we still find that a non-detection by a simulated array, with a sensitivity that should be reached in practice within a few years, disfavors gravitational-wave-dominated evolution with an odds ratio of $\sim$30:1. Such a finding would suggest either that all existing astrophysical models for the population of tightly bound binaries are overly optimistic, or else that some dynamical effect other than gravitational-wave emission is actually dominating binary evolution even at the relatively high frequencies/small orbital separations probed by PTAs., Comment: 14 pages, 8 figures

Recently, a family of exact force-free electrodynamic (FFE) solutions was given by Brennan, Gralla and Jacobson, which generalizes earlier solutions by Michel, Menon and Dermer, and other authors. These solutions have been proposed as useful models for describing the outer magnetosphere of conducting stars. As with any exact analytical solution that aspires to describe actual physical systems, it is vitally important that the solution possess the necessary stability. In this paper, we show via fully nonlinear numerical simulations that the aforementioned FFE solutions, despite being highly special in their properties, are nonetheless stable under small perturbations. Through this study, we also introduce a three-dimensional pseudospectral relativistic FFE code that achieves exponential convergence for smooth test cases, as well as two additional well-posed FFE evolution systems in the appendix that have desirable mathematical properties. Furthermore, we provide an explicit analysis that demonstrates how propagation along degenerate principal null directions of the spacetime curvature tensor simplifies scattering, thereby providing an intuitive understanding of why these exact solutions are tractable, i.e. why they are not backscattered by spacetime curvature., Comment: 33 pages, 21 figures; V2 updated to match published version

Compact binaries that emit gravitational waves in the sensitivity band of ground-based detectors can have non-negligible eccentricities just prior to merger, depending on the formation scenario. We develop a purely analytic, frequency-domain model for gravitational waves emitted by compact binaries on orbits with small eccentricity, which reduces to the quasi-circular post-Newtonian approximant TaylorF2 at zero eccentricity and to the post-circular approximation of Yunes et al. (2009) at small eccentricity. Our model uses a spectral approximation to the (post-Newtonian) Kepler problem to model the orbital phase as a function of frequency, accounting for eccentricity effects up to ${\cal{O}}(e^8)$ at each post-Newtonian order. Our approach accurately reproduces an alternative time-domain eccentric waveform model for eccentricities $e\in [0, 0.4]$ and binaries with total mass less than 12 solar masses. As an application, we evaluate the signal amplitude that eccentric binaries produce in different networks of existing and forthcoming gravitational waves detectors. Assuming a population of eccentric systems containing black holes and neutron stars that is uniformly distributed in co-moving volume, we estimate that second generation detectors like Advanced LIGO could detect approximately 0.1-10 events per year out to redshift $z\sim 0.2$, while an array of Einstein Telescope detectors could detect hundreds of events per year to redshift $z \sim 2.3$., Comment: 12 pages, 6 figures, 1 appendix. Submitted to Phys. Rev. D. v2: affiliations updated, one reference corrected. Accepted to Phys. Rev. D

In dense stellar regions, highly eccentric binaries of black holes and neutron stars can form through various n-body interactions. Such a binary could emit a significant fraction of its binding energy in a sequence of largely isolated gravitational wave bursts prior to merger. Given expected black hole and neutron star masses, many such systems will emit these repeated bursts at frequencies within the sensitive band of contemporary ground-based gravitational wave detectors. Unfortunately, existing gravitational wave searches are ill-suited to detect these signals. In this work, we adapt a "power stacking" method to the detection of gravitational wave signals from highly eccentric binaries. We implement this method as an extension of the Q-transform, a projection onto a multiresolution basis of windowed complex exponentials that has previously been used to analyze data from the network of LIGO/Virgo detectors. Our method searches for excess power over an ensemble of time-frequency tiles. We characterize the performance of our method using Monte Carlo experiments with signals injected in simulated detector noise. Our results indicate that the power stacking method achieves substantially better sensitivity to eccentric binary signals than existing localized burst searches., Comment: 17 pages, 20 figures

We develop a waveform model to describe the inspiral, merger and ringdown of binary systems with comparable and intermediate mass-ratios. This model incorporates first-order conservative self-force corrections to the energy and angular momentum, which are valid in the strong-field regime [1]. We model the radiative part of the self-force by deriving second-order radiative corrections to the energy flux. These corrections are obtained by minimizing the phase discrepancy between our self-force model and the effective one body model [2, 3] for a variety of mass-ratios. We show that our model performs substantially better than post-Newtonian approximants currently used to model neutron star-black hole mergers from early inspiral to the innermost stable circular orbit. In order to match the late inspiral evolution onto the plunge regime, we extend the 'transition phase' developed by Ori and Thorne [4] by including finite mass-ratio corrections and modelling the orbital phase evolution using an implicit rotating source [5]. We explicitly show that the implicit rotating source approach provides a natural transition from late-time radiation to ringdown that is equivalent to ringdown waveform modelling based on a sum of quasinormal modes., Comment: 22 pages, 12 figures. Submitted to Phys. Rev. D. v2: Accepted to Phys. Rev. D. References updated. No corrections

We respond to a recent criticism from Zaslavskii (arXiv:1301.3429v2 [gr-qc]) of our article, which appeared in Phys. Rev. Lett. 110 011102 (2013). We clarify the points of apparent disagreement, none of which relate to the central result of our original paper, which is the demonstration that the diverging center-of-masss energy of collisions in Kerr spacetime, known as the "BSW effect", can never be observed in nature. To emphasize this result, we present it in a more concise way which completely avoids any of the questions brought up by Zaslavskii., Comment: 2 pages, 1 figure

Gravitational waves at suitable frequencies can resonantly interact with a binary system, inducing changes to its orbit. A stochastic gravitational-wave background causes the orbital elements of the binary to execute a classic random walk, with the variance of orbital elements growing with time. The lack of such a random walk in binaries that have been monitored with high precision over long time-scales can thus be used to place an upper bound on the gravitational-wave background. Using periastron time data from the Hulse-Taylor binary pulsar spanning ~30 years, we obtain a bound of h_c < 7.9*10^(-14) at ~10^(-4) Hz, where h_c is the strain amplitude per logarithmic frequency interval. Our constraint complements those from pulsar timing arrays, which probe much lower frequencies, and ground-based gravitational-wave observations, which probe much higher frequencies. Interesting sources in our frequency band, which overlaps the lower sensitive frequencies of proposed space-based observatories, include white-dwarf/supermassive black-hole binaries in the early/late stages of inspiral, and TeV scale preheating or phase transitions. The bound improves as (time span)^(-2) and (sampling rate)^(-1/2). The Hulse-Taylor constraint can be improved to ~3.8*10^(-15) with a suitable observational campaign over the next decade. Our approach can also be applied to other binaries, including (with suitable care) the Earth-Moon system, to obtain constraints at different frequencies. The observation of additional binary pulsars with the SKA could reach a sensitivity of h_c ~ 3*10^(-17)., Comment: 8 pages, 1 figure. version 2, clarified several points thanks to private communications with Marc Kamionkowski

It has been suggested that maximally spinning black holes can serve as particle accelerators, reaching arbitrarily high center-of-mass energies. Despite several objections regarding the practical achievability of such high energies, and demonstrations past and present that such large energies could never reach a distant observer, interest in this problem has remained substantial. We show that, unfortunately, a maximally spinning black hole can never serve as a probe of high energy collisions, even in principle and despite the correctness of the original diverging energy calculation. Black holes can indeed facilitate dark matter annihilation, but the most energetic photons can carry little more than the rest energy of the dark matter particles to a distant observer, and those photons are actually generated relatively far from the black hole where relativistic effects are negligible. Therefore, any strong gravitational potential could probe dark matter equally well, and an appeal to black holes for facilitating such collisions is unnecessary., Comment: 5 pages, 3 figures, accepted to Phys. Rev. Lett

We assess the detectability of the gravitational wave signals from highly eccentric compact binaries. We use a simple model for the inspiral, merger, and ringdown of these systems. The model is based on mapping the binary to an effective single black hole system described by a Kerr metric, thereby including certain relativistic effects such as zoom-whirl-type behavior. The resultant geodesics source quadrupolar radiation and, in turn, are evolved under its dissipative effects. At the light ring, we attach a merger model that was previously developed for quasicircular mergers but also performs well for eccentric mergers with little modification. We apply this model to determine the detectability of these sources for initial, Enhanced, and Advanced LIGO across the parameter space of nonspinning close capture compact binaries. We conclude that, should these systems exist in nature, the vast majority will be missed by conventional burst searches or by quasicircular waveform templates in the advanced detector era. Other methods, such as eccentric templates or, more practically, a stacked excess power search, must be developed to avoid losing these sources. These systems would also have been missed frequently in the initial LIGO data analysis. Thus, previous null coincidence results with detected gamma-ray bursts cannot exclude the possibility of coincident gravitational wave signals from eccentric binaries., Comment: 18 pages, 18 figures; revised to match accepted version, PRD in press

We present a model for merger-driven evolution of the mass function for massive galaxies and their central supermassive black holes at late times. We discuss the current observational evidence in favor of merger-driven massive galaxy evolution during this epoch, and demonstrate that the observed evolution of the mass function can be reproduced by evolving an initial mass function under the assumption of negligible star formation. We calculate the stochastic gravitational wave signal from the resulting black-hole binary mergers in the low redshift universe (z <= 1) implied by this model, and find that this population has a signal-to-noise ratio as much as ~5x larger than previous estimates for pulsar timing arrays, with an expectation value for the characteristic strain h_c (f=1 yr^{-1}) = 4.1 x 10^{-15} that may already be in tension with observational constraints, and a {2-sigma, 3-sigma} lower limit within this model of h_c (f=1 yr^{-1}) = {1.1 x 10^{-15}, 6.8 x 10^{-16}}. The strength of this signal is sufficient to make it detectable with high probability under conservative assumptions within the next several years, if the principle assumption of merger-driven galaxy evolution since z = 1 holds true. For cases where a galaxy merger fails to lead to a black hole merger, we estimate the probability for a given number of satellite unmerged black holes to remain within a massive host galaxy, and interpret the result in light of ULX observations. In particular, we find that the brightest cluster galaxies should have 1-2 such sources with luminosities above 10^{39} erg/s, which is consistent with the statistics of observed ULXs., Comment: 11 pages, 5 figures, submitted to ApJ, v2 includes the referee's requested changes

Recent observations of massive galaxies indicate that they double in mass and quintuple in size between redshift z = 1 and the present, despite undergoing very little star formation, suggesting that galaxy mergers drive the evolution. Since these galaxies will contain supermassive black holes, this suggests a larger black hole merger rate, and therefore a larger gravitational-wave signal, than previously expected. We calculate the merger-driven evolution of the mass function, and find that merger rates are 10 to 30 times higher and gravitational waves are 3 to 5 times stronger than previously estimated, so that the gravitational-wave signal may already be detectable with existing data from pulsar timing arrays. We also provide an explanation for the disagreement with past estimates that were based on dark matter halo simulations., Comment: 6 pages, 1 figure

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime., Comment: 20 pages, 8 figures, proceedings of the 9th Amaldi Conference on Gravitational Waves. Final journal version. For a longer exposition of the eLISA science case, see http://arxiv.org/abs/1201.3621

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits., Comment: Submitted to GW Notes; abstract abridged due to arXiv requirements

We present a mission concept, the Geostationary Antenna for Disturbance-Free Laser Interferometry (GADFLI), for a space-based gravitational-wave interferometer consisting of three satellites in geostationary orbit around the Earth. Compared to the nominal design of the Laser Interferometer Space Antenna (LISA), this concept has the advantages of significantly decreased requirements on the telescope size and laser power, decreased launch mass, substantially improved shot noise resulting from the shorter 73000 km armlengths, simplified and less expensive communications, and an overall lower cost which we (roughly) estimate at $1.2B. GADFLI preserves much of the science of LISA, particularly the observation of massive black-hole binary coalescences, although the SNR is diminished for all masses in the potential designs we consider., Comment: 6 pages, 3 figures. This document was submitted to NASA in response to a Request For Information, Solicitation #NNH11ZDA019L, Concepts for the NASA Gravitational-Wave Mission

We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission-an implicit rotating source. We further apply the late-time merger-ringdown model for the rotational frequency introduced in [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)], along with an improved amplitude model appropriate for the dominant (2, \pm2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above \sim150M\odot. This model may be directly applicable to gravitational-wave detection of intermediate-mass black-hole mergers., Comment: 19 pages, 17 figures; new version matches published article

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA's sensitivity band, with total rest mass of 2\times10^6 M\odot, a redshift of z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately \sim3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well- determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well., Comment: 9 pages, 5 figures, submitted to Phys. Rev. D

The coalescence of black hole-neutron star binaries is expected to be a principal source of gravitational waves for the next generation of detectors, Advanced LIGO and Advanced Virgo. Ideally, these and other gravitational wave sources would have a distinct electromagnetic counterpart, as significantly more information could be gained through two separate channels. In addition, since these detectors will probe distances with non-negligible redshift, a coincident observation of an electromagnetic counterpart to a gravitational wave signal would facilitate a novel measurement of dark energy [1]. For black hole masses not much larger than the neutron star mass, the tidal disruption and subsequent accretion of the neutron star by the black hole provides one avenue for generating an electromagnetic counterpart [2]. However, in this work, we demonstrate that, for all black hole-neutron star binaries observable by Advanced LIGO/Virgo, the interaction of the black hole with the magnetic field of the neutron star will drive a Poynting flux. This Poynting flux generates synchrotron/curvature radiation as the electron-positron plasma in the neutron star magnetosphere is accel- erated, and thermal radiation as the plasma is focused onto the neutron star magnetic poles, creating a "hot spot" on the neutron star surface. This novel effect will gener- ate copious luminosity, comparable to supernovae and active galactic nuclei, so that black hole-neutron star coalescences detectable with gravitational waves by Advanced LIGO/Virgo could also potentially be detectable electromagnetically., Comment: 17 pages, 2 figures, submitted to Nature

The advent of long-term stability in numerical relativity has yielded a windfall of answers to long-standing questions regarding the dynamics of space-time, matter, and electromagnetic fields in the strong-field regime of black-hole binary mergers. In this review, we will briefly summarize the methodology currently applied to these problems, emphasizing the most recent advancements. We will discuss recent results of astrophysical relevance, and present some novel interpretation. Though we primarily present a review, we also present a simple analytical model for the time-dependent Poynting flux from two orbiting black holes immersed in a magnetic field, which compares favorably with recent numerical results. Finally, we will discuss recent advancements in our theoretical understanding of merger dynamics and gravitational waveforms that have resulted from interpreting the ever-growing body of numerical relativity results., Comment: Invited review for NRDA/CAPRA 2010, CQG special issue - 7 pages, 1 figure, 1 table

Given the absence of observations of black hole binaries, it is critical that the full range of accessible parameter space be explored in anticipation of future observation with gravitational wave detectors. To this end, we compile the Hamiltonian equations of motion describing the conservative dynamics of the most general black hole binaries, as computed by Will and collaborators, and incorporate an effective treatment of dissipation through gravitational radiation. We evolve these equations for systems with orbital eccentricity and precessing spins. We find that, while spin-spin coupling corrections can destroy constant radius orbits in principle, the effect is so small that orbits will reliably tend to quasi-spherical orbits as angular momentum and energy are lost to gravitational radiation. Still, binaries that are initially highly eccentric may retain eccentricity as they pass into the detectable bandwidth of ground-based gravitational wave detectors. We also show that a useful set of natural frequencies for an orbit demonstrating both spin precession and periastron precession is comprised of (1) the frequency of angular motion in the orbital plane, (2) the frequency of the plane precession, and (3) the frequency of radial oscillations. These three natural harmonics shape the observed waveform., Comment: 29 pages, 8 figures, published in CQG

Advances in the field of numerical relativity now make it possible to calculate the final, most powerful merger phase of binary black-hole coalescence for generic binaries. The state of the art has advanced well beyond the equal-mass case into the unequal-mass and spinning regions of parameter space. We present a study of the nonspinning portion of parameter space, primarily using an analytic waveform model tuned to available numerical data, with an emphasis on observational implications. We investigate the impact of varied mass ratio on merger signal-to-noise ratios (SNRs) for several detectors, and compare our results with expectations from the test-mass limit. We note a striking similarity of the waveform phasing of the merger waveform across the available mass ratios. Motivated by this, we calculate the match between our 1:1 (equal mass) and 4:1 mass-ratio waveforms during the merger as a function of location on the source sky, using a new formalism for the match that accounts for higher harmonics. This is an indicator of the amount of degeneracy in mass ratio for mergers of moderate-mass-ratio systems., Comment: 13 pages, 11 figures, submitted to Phys. Rev. D

Some braneworld models may have observable consequences that, if detected, would validate a requisite element of string theory. In the infinite Randall-Sundrum model (RS2), the AdS radius of curvature, L, of the extra dimension supports a single bound state of the massless graviton on the brane, thereby reproducing Newtonian gravity in the weak-field limit. However, using the AdS/CFT correspondence, it has been suggested that one possible consequence of RS2 is an enormous increase in Hawking radiation emitted by black holes. We utilize this possibility to derive two novel methods for constraining L via gravitational wave measurements. We show that the EMRI event rate detected by LISA can constrain L at the ~1 micron level for optimal cases, while the observation of a single galactic black hole binary with LISA results in an optimal constraint of L <= 5 microns., Comment: 4 pages, replaced with version published in Phys. Rev. Lett.

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise., Comment: 12 pages, 2 figures, proceedings of the 8th Edoardo Amaldi Conference on Gravitational Waves, New York, June 21-26, 2009

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin)., Comment: 16 pages, 9 figures, replaced with version published in Phys. Rev. D

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal., Comment: 5 pages, 4 figures

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise., Comment: Proceedings of Numerical Relativity and Data Analysis NRDA 2008 (Syracuse NY, August 2008) - 14 pages, 1 figure

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses., Comment: 56 pages, 25 figures; various clarifications; accepted to CQG

We present a detailed descriptive analysis of the gravitational radiation from black-hole binary mergers of nonspinning black holes, based on numerical simulations of systems varying from equal-mass to a 6:1 mass ratio. Our primary goal is to present relatively complete information about the waveforms, including all the leading multipolar components, to interested researchers. In our analysis, we pursue the simplest physical description of the dominant features in the radiation, providing an interpretation of the waveforms in terms of an {\em implicit rotating source}. This interpretation applies uniformly to the full wave train, from inspiral through ringdown. We emphasize strong relationships among the $\ell=m$ modes that persist through the full wave train. Exploring the structure of the waveforms in more detail, we conduct detailed analytic fitting of the late-time frequency evolution, identifying a key quantitative feature shared by the $\ell=m$ modes among all mass ratios. We identify relationships, with a simple interpretation in terms of the implicit rotating source, among the evolution of frequency and amplitude, which hold for the late-time radiation. These detailed relationships provide sufficient information about the late-time radiation to yield a predictive model for the late-time waveforms, an alternative to the common practice of modeling by a sum of quasinormal mode overtones. We demonstrate an application of this in a new effective-one-body-based analytic waveform model., Comment: 25 pages, 21 figures, 5 tables; published version, with added references

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters., Comment: 5 pages, 1 figure, 3 tables. Version published in Astrophys. J. Lett

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics., Comment: Added reference on large spin-kick configurations. 14 pages, 5 figures, 1 table. To be published in the Proceedings of the 2007 Spring Symposium of the Space Telescope Science Institute (Baltimore, MD), held from Monday April 23, 2007 to Thursday April 26, 2007

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes., Comment: 28 pages, 20 figures, submitted to PRD; v2: minor revisions from referee report

We present an accurate approximation of the full gravitational radiation waveforms generated in the merger of non-eccentric systems of two non-spinning black holes. Utilizing information from recent numerical relativity simulations and the natural flexibility of the effective-one-body (EOB) model, we extend the latter so that it can successfully match the numerical relativity waveforms during the last stages of inspiral, merger and ringdown. By ``successfully'' here, we mean with phase differences < 8% of a gravitational-wave cycle accumulated by the end of the ringdown phase, maximizing only over time of arrival and initial phase. We obtain this result by simply adding a 4-post-Newtonian order correction in the EOB radial potential and determining the (constant) coefficient by imposing high-matching performances with numerical waveforms of mass ratios m1/m2 = 1, 3/2, 2 and 4, m1 and m2 being the individual black-hole masses. The final black-hole mass and spin predicted by the numerical simulations are used to determine the ringdown frequency and decay time of three quasi-normal-mode damped sinusoids that are attached to the EOB inspiral-(plunge) waveform at the EOB light-ring. The EOB waveforms might be tested and further improved in the future by comparison with extremely long and accurate inspiral numerical-relativity waveforms. They may already be employed for coherent searches and parameter estimation of gravitational waves emitted by non-spinning coalescing binary black holes with ground-based laser-interferometer detectors., Comment: 15 pages, 9 figures

We compare waveforms obtained by numerically evolving nonspinning binary black holes to post-Newtonian (PN) template families currently used in the search for gravitational waves by ground-based detectors. We find that the time-domain 3.5PN template family, which includes the inspiral phase, has fitting factors (FFs) >= 0.96 for binary systems with total mass M = 10 ~ 20 Msun. The time-domain 3.5PN effective-one-body template family, which includes the inspiral, merger and ring-down phases, gives satisfactory signal-matching performance with FFs >= 0.96 for binary systems with total mass M = 10 ~ 120 Msun. If we introduce a cutoff frequency properly adjusted to the final black-hole ring-down frequency, we find that the frequency-domain stationary-phase-approximated template family at 3.5PN order has FFs >= 0.96 for binary systems with total mass M = 10 ~ 20 Msun. However, to obtain high matching performances for larger binary masses, we need to either extend this family to unphysical regions of the parameter space or introduce a 4PN order coefficient in the frequency-domain GW phase. Finally, we find that the phenomenological Buonanno-Chen-Vallisneri family has FFs >= 0.97 with total mass M=10 ~ 120Msun. The main analyses use the noise spectral-density of LIGO, but several tests are extended to VIRGO and advanced LIGO noise-spectral densities., Comment: 19 pages, 17 figures, added references, corrected typos, changed figure 16

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%., Comment: 5 pages, 3 figures. Replaced with published version