Search

Modern epidemiological analyses to understand and combat the spread of disease depend critically on access to, and use of, data. Rapidly evolving data, such as data streams changing during a disease outbreak, are particularly challenging. Data management is further complicated by data being imprecisely identified when used. Public trust in policy decisions resulting from such analyses is easily damaged and is often low, with cynicism arising where claims of "following the science" are made without accompanying evidence. Tracing the provenance of such decisions back through open software to primary data would clarify this evidence, enhancing the transparency of the decision-making process. Here, we demonstrate a Findable, Accessible, Interoperable and Reusable (FAIR) data pipeline developed during the COVID-19 pandemic that allows easy annotation of data as they are consumed by analyses, while tracing the provenance of scientific outputs back through the analytical source code to data sources. Such a tool provides a mechanism for the public, and fellow scientists, to better assess the trust that should be placed in scientific evidence, while allowing scientists to support policy-makers in openly justifying their decisions. We believe that tools such as this should be promoted for use across all areas of policy-facing research.

As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM., Comment: 24 pages, 18 figures. Abstract abridged

Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms ., Comment: 33+18 pages, 13 figures, 4 tables, 2,018 binaries. Catalog metadata in ancillary JSON file. v2: Matches version accepted by CQG. Catalog available at https://www.black-holes.org/waveforms

We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario., Comment: 14 pages, 8 figures, final published version, samples available at https://git.ligo.org/katerina.chatziioannou/gw170729hm_datarelease

Observation and characterisation of gravitational waves from binary black holes requires accurate knowledge of the expected waveforms. The late inspiral and merger phase of the waveform is obtained through direct numerical integration of the full 3-dimensional Einstein equations. The Spectral Einstein Code (SpEC) utilizes a multi-domain pseudo-spectral method tightly adapted to the geometry of the black holes; it is computationally efficient and accurate, but--for high mass-ratios and large spins--sometimes requires manual fine-tuning for the merger-phase of binaries. The Einstein Toolkit (ET) employs finite difference methods and the moving puncture technique; it is less computationally efficient, but highly robust. For some mergers with high mass ratio and large spins, the efficient numerical algorithms used in SpEC have failed, whereas the simpler algorithms used in the ET were successful. Given the urgent need of testing the accuracy of waveform models currently used in LIGO and Virgo inference analyses for high mass ratios and spins, we present here a synergistic approach to numerical-relativity: We combine SpEC and ET waveforms into complete inspiral-merger-ringdown waveforms, taking advantage of the computational efficiency of the pseudo-spectral code during the inspiral, and the robustness of the finite-difference code at the merger. We validate our method against a case where complete waveforms from both codes are available, compute three new hybrid numerical-relativity waveforms, and compare them with analytical waveform models currently used in LIGO and Virgo science. All the waveforms and the hybridization code are publicly available., Comment: 7 pages, 4 figures

For the first time, we construct an inspiral-merger-ringdown waveform model within the effective-one-body formalism for spinning, nonprecessing binary black holes that includes gravitational modes beyond the dominant $(\ell,|m|) = (2,2)$ mode, specifically $(\ell,|m|)=(2,1),(3,3),(4,4),(5,5)$. Our multipolar waveform model incorporates recent (resummed) post-Newtonian results for the inspiral and information from 157 numerical-relativity simulations, and 13 waveforms from black-hole perturbation theory for the (plunge-)merger and ringdown. We quantify the improved accuracy including higher-order modes by computing the faithfulness of the waveform model against the numerical-relativity waveforms used to construct the model. We define the faithfulness as the match maximized over time, phase of arrival, gravitational-wave polarization and sky position of the waveform model, and averaged over binary orientation, gravitational-wave polarization and sky position of the numerical-relativity waveform. When the waveform model contains only the $(2,2)$ mode, we find that the averaged faithfulness to numerical-relativity waveforms containing all modes with $\ell \leq$ 5 ranges from $90\%$ to $99.9\%$ for binaries with total mass $20-200 M_\odot$ (using the Advanced LIGO's design noise curve). By contrast, when the $(2,1),(3,3),(4,4),(5,5)$ modes are also included in the model, the faithfulness improves to $99\%$ for all but four configurations in the numerical-relativity catalog, for which the faithfulness is greater than $98.5\%$. Using our results, we also develop also a (stand-alone) waveform model for the merger-ringdown signal, calibrated to numerical-relativity waveforms, which can be used to measure multiple quasi-normal modes. The multipolar waveform model can be extended to include spin-precession, and will be employed in upcoming observing runs of Advanced LIGO and Virgo., Comment: 28 pages. Version that matches published article

We present a prescription for computing gravitational waveforms for the inspiral, merger and ringdown of non-spinning eccentric binary black hole systems. The inspiral waveform is computed using the post-Newtonian expansion and the merger waveform is computed by interpolating a small number of quasi-circular NR waveforms. The use of circular merger waveforms is possible because eccentric binaries circularize in the last few cycles before the merger, which we demonstrate up to mass ratio $q = m_1/m_2 = 3$. The complete model is calibrated to 23 numerical relativity (NR) simulations starting ~20 cycles before the merger with eccentricities $e_\text{ref} \le 0.08$ and mass ratios $q \le 3$, where $e_\text{ref}$ is the eccentricity ~7 cycles before the merger. The NR waveforms are long enough that they start above 30 Hz (10 Hz) for BBH systems with total mass $M \ge 80 M_\odot$ ($230 M_\odot$). We find that, for the sensitivity of advanced LIGO at the time of its first observing run, the eccentric model has a faithfulness with NR of over 97% for systems with total mass $M \ge 85 M_\odot$ across the parameter space ($e_\text{ref} \le 0.08, q \le 3$). For systems with total mass $M \ge 70 M_\odot$, the faithfulness is over 97% for $e_\text{ref} \lesssim 0.05$ and $q \le 3$. The NR waveforms and the Mathematica code for the model are publicly available.

We present and assess a Bayesian method to interpret gravitational wave signals from binary black holes. Our method directly compares gravitational wave data to numerical relativity simulations. This procedure bypasses approximations used in semi-analytical models for compact binary coalescence. In this work, we use only the full posterior parameter distribution for generic nonprecessing binaries, drawing inferences away from the set of NR simulations used, via interpolation of a single scalar quantity (the marginalized log-likelihood, $\ln {\cal L}$) evaluated by comparing data to nonprecessing binary black hole simulations. We also compare the data to generic simulations, and discuss the effectiveness of this procedure for generic sources. We specifically assess the impact of higher order modes, repeating our interpretation with both $l\le2$ as well as $l\le3$ harmonic modes. Using the $l\le3$ higher modes, we gain more information from the signal and can better constrain the parameters of the gravitational wave signal. We assess and quantify several sources of systematic error that our procedure could introduce, including simulation resolution and duration; most are negligible. We show through examples that our method can recover the parameters for equal mass, zero spin; GW150914-like; and unequal mass, precessing spin sources. Our study of this new parameter estimation method demonstrates we can quantify and understand the systematic and statistical error. This method allows us to use higher order modes from numerical relativity simulations to better constrain the black hole binary parameters., Comment: 30 pages, 22 figures; submitted to PRD

The apparent properties of distant objects encode information about the way the light they emit propagates to an observer, and therefore about the curvature of the underlying spacetime. Measuring the relationship between the redshift $z$ and the luminosity distance $D_{\rm L}$ of a standard candle, for example, yields information on the Universe's matter content. In practice, however, in order to decode this information the observer needs to make an assumption about the functional form of the $D_{\rm L}(z)$ relation; in other words, a cosmological model needs to be assumed. In this work, we use numerical-relativity simulations, equipped with a new ray-tracing module, to numerically obtain this relation for a few black-hole--lattice cosmologies and compare it to the well-known Friedmann-Lema\^itre-Robertson-Walker case, as well as to other relevant cosmologies and to the Empty-Beam Approximation. We find that the latter provides the best estimate of the luminosity distance and formulate a simple argument to account for this agreement. We also find that a Friedmann-Lema\^itre-Robertson-Walker model can reproduce this observable exactly, as long as a time-dependent cosmological constant is included in the fit. Finally, the dependence of these results on the lattice mass-to-spacing ratio $\mu$ is discussed: we discover that, unlike the expansion rate, the $D_{\rm L}(z)$ relation in a black-hole lattice does not tend to that measured in the corresponding continuum spacetime as $\mu \to 0$., Comment: 32 pages, 10 figures, matches published version

We improve the accuracy of the effective-one-body (EOB) waveforms that were employed during the first observing run of Advanced LIGO for binaries of spinning, nonprecessing black holes by calibrating them to a set of 141 numerical-relativity (NR) waveforms. The NR simulations expand the domain of calibration towards larger mass ratios and spins, as compared to the previous EOBNR model. Merger-ringdown waveforms computed in black-hole perturbation theory for Kerr spins close to extremal provide additional inputs to the calibration. For the inspiral-plunge phase, we use a Markov-chain Monte Carlo algorithm to efficiently explore the calibration space. For the merger-ringdown phase, we fit the NR signals with phenomenological formulae. After extrapolation of the calibrated model to arbitrary mass ratios and spins, the (dominant-mode) EOBNR waveforms have faithfulness --- at design Advanced-LIGO sensitivity --- above $99\%$ against all the NR waveforms, including 16 additional waveforms used for validation, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling for these binary configurations. We find that future NR simulations at mass ratios $\gtrsim 4$ and double spin $\gtrsim 0.8$ will be crucial to resolve discrepancies between different ways of extrapolating waveform models. We also find that some of the NR simulations that already exist in such region of parameter space are too short to constrain the low-frequency portion of the models. Finally, we build a reduced-order version of the EOBNR model to speed up waveform generation by orders of magnitude, thus enabling intensive data-analysis applications during the upcoming observation runs of Advanced LIGO., Comment: 27 pages, 15 figures

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers., Comment: 28 pages, 13 figures, 7 tables

We discuss the possibility of a dimensional reduction of the Einstein equations in S3 black-hole lattices. It was reported in previous literature that the evolution of spaces containing curves of local, discrete rotational and reflection Symmetry (LDRRS) can be carried out via a system of ODEs along these curves. However, 3+1 Numerical Relativity computations demonstrate that this is not the case, and we show analytically that this is due to the presence of a tensorial quantity which is not suppressed by the symmetry. We calculate the term analytically, and verify numerically for an 8-black-hole lattice that it fully accounts for the anomalous results, and thus quantify its magnitude in this specific case. The presence of this term prevents the exact evolution of these spaces via previously-reported methods which do not involve a full 3+1 integration of Einstein's equation., Comment: 18 pages, 4 figures

We probe the gravitational interaction of two black holes in the strong-field regime by computing the scattering angle $\chi$ of hyperbolic-like, close binary-black-hole encounters as a function of the impact parameter. The fully general-relativistic result from numerical relativity is compared to two analytic approximations: post-Newtonian theory and the effective-one-body formalism. As the impact parameter decreases, so that black holes pass within a few times their Schwarzschild radii, we find that the post-Newtonian prediction becomes quite inaccurate, while the effective-one-body one keeps showing a good agreement with numerical results. Because we have explored a regime which is very different from the one considered so far with binaries in quasi-circular orbits, our results open a new avenue to improve analytic representations of the general-relativistic two-body Hamiltonian., Comment: 5 pages, 3 figures. Submitted to Physical Review Letters

Starting from a high-level problem description in terms of partial differential equations using abstract tensor notation, the Chemora framework discretizes, optimizes, and generates complete high performance codes for a wide range of compute architectures. Chemora extends the capabilities of Cactus, facilitating the usage of large-scale CPU/GPU systems in an efficient manner for complex applications, without low-level code tuning. Chemora achieves parallelism through MPI and multi-threading, combining OpenMP and CUDA. Optimizations include high-level code transformations, efficient loop traversal strategies, dynamically selected data and instruction cache usage strategies, and JIT compilation of GPU code tailored to the problem characteristics. The discretization is based on higher-order finite differences on multi-block domains. Chemora's capabilities are demonstrated by simulations of black hole collisions. This problem provides an acid test of the framework, as the Einstein equations contain hundreds of variables and thousands of terms., Comment: 18 pages, 4 figures, accepted for publication in Scientific Programming

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters., Comment: 51 pages, 10 figures; published version

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 and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations., Comment: Presented at Amaldi 9

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems., Comment: 62 pages, 20 figures

The peeling theorem of general relativity predicts that the Weyl curvature scalars Psi_n (n=0...4), when constructed from a suitable null tetrad in an asymptotically flat spacetime, fall off asymptotically as r^(n-5) along outgoing radial null geodesics. This leads to the interpretation of Psi_4 as outgoing gravitational radiation at large distances from the source. We have performed numerical simulations in full general relativity of a binary black hole inspiral and merger, and have computed the Weyl scalars in the standard tetrad used in numerical relativity. In contrast with previous results, we observe that all the Weyl scalars fall off according to the predictions of the theorem., Comment: 7 pages, 3 figures, published version

Suitable gauge conditions are fundamental for stable and accurate numerical-relativity simulations of inspiralling compact binaries. A number of well-studied conditions have been developed over the last decade for both the lapse and the shift and these have been successfully used both in vacuum and non-vacuum spacetimes when simulating binaries with comparable masses. At the same time, recent evidence has emerged that the standard "Gamma-driver" shift condition requires a careful and non-trivial tuning of its parameters to ensure long-term stable evolutions of unequal-mass binaries. We present a novel gauge condition in which the damping constant is promoted to be a dynamical variable and the solution of an evolution equation. We show that this choice removes the need for special tuning and provides a shift damping term which is free of instabilities in our simulations and dynamically adapts to the individual positions and masses of the binary black-hole system. Our gauge condition also reduces the variations in the coordinate size of the apparent horizon of the larger black hole and could therefore be useful when simulating binaries with very small mass ratios., Comment: 11 pages, 8 figures

Since the breakthroughs in 2005 which have led to long term stable solutions of the binary black hole problem in numerical relativity, much progress has been made. I present here a short summary of the state of the field, including the capabilities of numerical relativity codes, recent physical results obtained from simulations, and improvements to the methods used to evolve and analyse binary black hole spacetimes., Comment: 14 pages; minor changes and corrections in response to referees

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

Approximate solutions to the Einstein field equations are a valuable tool to investigate gravitational phenomena. An important aspect of any approximation is to investigate and quantify its regime of validity. We present a study that evaluates the effects that approximate puncture initial data, based on "skeleton" solutions to the Einstein constraints as proposed by Faye et al. [PRD 69, 124029 (2004)], have on numerical evolutions. Using data analysis tools, we assess the effectiveness of these constraint-violating initial data and show that the matches of waveforms from skeleton data with the corresponding waveforms from constraint-satisfying initial data are > 0.97 when the total mass of the binary is > 40M(solar). In addition, we demonstrate that the differences between the skeleton and the constraint-satisfying initial data evolutions, and thus waveforms, are due to negative Hamiltonian constraint violations present in the skeleton initial data located in the vicinity of the punctures. During the evolution, the skeleton data develops both Hamiltonian and momentum constraint violations that decay with time, with the binary system relaxing to a constraint-satisfying solution with black holes of smaller mass and thus different dynamics.

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 quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases., Comment: 17 pages, 9 figures. Version accepted by PRD

Generic inspirals and mergers of binary black holes produce beamed emission of gravitational radiation that can lead to a gravitational recoil or kick of the final black hole. The kick velocity depends on the mass ratio and spins of the binary as well as on the dynamics of the binary configuration. Studies have focused so far on the most astrophysically relevant configuration of quasi-circular inspirals, for which kicks as large as 3,300 km/s have been found. We present the first study of gravitational recoil in hyperbolic encounters. Contrary to quasi-circular configurations, in which the beamed radiation tends to average during the inspiral, radiation from hyperbolic encounters is plunge dominated, resulting in an enhancement of preferential beaming. As a consequence, it is possible to achieve kick velocities as large as 10,000 km/s., Comment: 4 pages, 5 figures, 1 table

We present the first comparison between numerical relativity (NR) simulations of an eccentric binary black hole system with corresponding post-Newtonian (PN) results. We evolve an equal-mass, non-spinning configuration with an initial eccentricity e ~ 0.1 for 21 gravitational wave cycles before merger, and find agreement in the gravitational wave phase with an adiabatic eccentric PN model with 2 PN radiation reaction within 0.1 radians for 10 cycles. The NR and PN phase difference grows to 0.7 radians by 5 cycles before merger. We find that these results can be obtained by expanding the eccentric PN expressions in terms of the frequency-related variable x = (omega M)^{2/3} with M the total mass of the binary. When using instead the mean motion n = 2 \pi /P, where P is the orbital period, the comparison leads to significant disagreements with NR., Comment: 15 pages, 7 figures. Replaced extrapolation of fit parameters with values at earliest time - results unchanged. Corrected equation typos in Section IIa. Other minor changes. Corresponds to PRD version

We present a study of the gravitational waveforms from a series of spinning, equal-mass black hole binaries focusing on the harmonic content of the waves and the contribution of the individual harmonics to the signal-to-noise ratio. The gravitational waves were produced from two series of evolutions with black holes of initial spins equal in magnitude and anti-aligned with each other. In one series the magnitude of the spin is varied; while in the second, the initial angle between the black-hole spins and the orbital angular momentum varies. We also conduct a preliminary investigation into using these waveforms as templates for detecting spinning binary black holes. Since these runs are relativity short, containing about two to three orbits, merger and ringdown, we limit our study to systems of total mass greater than 50 solar masses. This choice ensures that our waveforms are present in the ground-based detector band without needing addition gravitational wave cycles. We find that while the mode contribution to the signal-to-noise ratio varies with the initial angle, the total mass of the system caused greater variations in the match., Comment: Invited submission to Class.Quant.Grav. for Amaldi07 proceedings

We present results obtained by scattering a scalar field off the curved background of a coalescing binary black hole system. A massless scalar field is evolved on a set of fixed backgrounds, each provided by a spatial hypersurface generated numerically during a binary black hole merger. We show that the scalar field scattered from the merger region exhibits quasinormal ringing once a common apparent horizon surrounds the two black holes. This occurs earlier than the onset of the perturbative regime as measured by the start of the quasinormal ringing in the gravitational waveforms. We also use the scalar quasinormal frequencies to associate a mass and a spin with each hypersurface, and observe the compatibility of this measure with the horizon mass and spin computed from the dynamical horizon framework., Comment: 10 Pages and 6 figures

We present an investigation into how sensitive the last orbits and merger of binary black hole systems are to the presence of spurious radiation in the initial data. Our numerical experiments consist of a binary black hole system starting the last couple of orbits before merger with additional spurious radiation centered at the origin and fixed initial angular momentum. As the energy in the added spurious radiation increases, the binary is invariably hardened for the cases we tested, i.e. the merger of the two black holes is hastened. The change in merger time becomes significant when the additional energy provided by the spurious radiation increases the Arnowitt-Deser-Misner (ADM) mass of the spacetime by about 1%. While the final masses of the black holes increase due to partial absorption of the radiation, the final spins remain constant to within our numerical accuracy. We conjecture that the spurious radiation is primarily increasing the eccentricity of the orbit and secondarily increasing the mass of the black holes while propagating out to infinity., Comment: 12 pages, 12 figures

We present results from numerical relativity simulations of equal mass, non-spinning binary black hole inspirals and mergers with initial eccentricities e <= 0.8 and coordinate separations D >= 12 M of up to 9 orbits (18 gravitational wave cycles). We extract the mass M_f and spin a_f of the final black hole and find, for eccentricities e < 0.4, that a_f/M_f = 0.69 and M_f/M_adm = 0.96 are independent of the initial eccentricity, suggesting that the binary has circularized by the merger time. For e > 0.5, the black holes plunge rather than orbit, and we obtain a maximum spin parameter a_f/M_f = 0.72 around e = 0.5., Comment: 4 pages, 4 figures, 2 tables. Final version in PRD

We present a study of spinning black hole binaries focusing on the spin dynamics of the individual black holes as well as on the gravitational recoil acquired by the black hole produced by the merger. We consider two series of initial spin orientations away from the binary orbital plane. In one of the series, the spins are anti-aligned; for the second series, one of the spins points away from the binary along the line separating the black holes. We find a remarkable agreement between the spin dynamics predicted at 2nd post-Newtonian order and those from numerical relativity. For each configuration, we compute the kick of the final black hole. We use the kick estimates from the series with anti-aligned spins to fit the parameters in the \KKF{,} and verify that the recoil along the direction of the orbital angular momentum is $\propto \sin\theta$ and on the orbital plane $\propto \cos\theta$, with $\theta$ the angle between the spin directions and the orbital angular momentum. We also find that the black hole spins can be well estimated by evaluating the isolated horizon spin on spheres of constant coordinate radius., Comment: 15 pages, 10 figures, replaced with version accepted for publication in PRD

Tremendous progress has been made towards the solution of the binary-black-hole problem in numerical relativity. The waveforms produced by numerical relativity will play a role in gravitational wave detection as either test-beds for analytic template banks or as template banks themselves. As the parameter space explored by numerical relativity expands, the importance of quantifying the effect that each parameter has on first the detection of gravitational waves and then the parameter estimation of their sources increases. In light of this, we present a study of equal-mass, spinning binary-black-hole evolutions through matched filtering techniques commonly used in data analysis. We study how the match between two numerical waveforms varies with numerical resolution, initial angular momentum of the black holes and the inclination angle between the source and the detector. This study is limited by the fact that the spinning black-hole-binaries are oriented axially and the waveforms only contain approximately two and a half orbits before merger. We find that for detection purposes, spinning black holes require the inclusion of the higher harmonics in addition to the dominant mode, a condition that becomes more important as the black-hole-spins increase. In addition, we conduct a preliminary investigation of how well a template of fixed spin and inclination angle can detect target templates of arbitrary spin and inclination for the axial case considered here.

The inspiral and merger of binary black holes will likely involve black holes with both unequal masses and arbitrary spins. The gravitational radiation emitted by these binaries will carry angular as well as linear momentum. A net flux of emitted linear momentum implies that the black hole produced by the merger will experience a recoil or kick. Previous studies have focused on the recoil velocity from unequal mass, non-spinning binaries. We present results from simulations of equal mass but spinning black hole binaries and show how a significant gravitational recoil can also be obtained in these situations. We consider the case of black holes with opposite spins of magnitude $a$ aligned/anti-aligned with the orbital angular momentum, with $a$ the dimensionless spin parameters of the individual holes. For the initial setups under consideration, we find a recoil velocity of $V = 475 \KMS a$. Supermassive black hole mergers producing kicks of this magnitude could result in the ejection from the cores of dwarf galaxies of the final hole produced by the collision., Comment: 8 pages, 8 figures, replaced with version accepted for publication in ApJ

We show that by adding suitable lower-order terms to the Z4 formulation of the Einstein equations, all constraint violations except constant modes are damped. This makes the Z4 formulation a particularly simple example of a lambda-system as suggested by Brodbeck et al. We also show that the Einstein equations in harmonic coordinates can be obtained from the Z4 formulation by a change of variables that leaves the implied constraint evolution system unchanged. Therefore the same method can be used to damp all constraints in the Einstein equations in harmonic gauge., Comment: Extended derivation of damping mechanism

We extend the notion of numerical stability of finite difference approximations to include hyperbolic systems that are first order in time and second order in space, such as those that appear in Numerical Relativity. By analyzing the symbol of the second order system, we obtain necessary and sufficient conditions for stability in a discrete norm containing one-sided difference operators. We prove stability for certain toy models and the linearized Nagy-Ortiz-Reula formulation of Einstein's equations. We also find that, unlike in the fully first order case, standard discretizations of some well-posed problems lead to unstable schemes and that the Courant limits are not always simply related to the characteristic speeds of the continuum problem. Finally, we propose methods for testing stability for second order in space hyperbolic systems., Comment: 18 pages, 9 figures

We study the absorption of massive spin-half particles by a small Schwarzschild black hole by numerically solving the single-particle Dirac equation in Painleve-Gullstrand coordinates. We calculate the absorption cross section for a range of gravitational couplings Mm/m_P^2 and incident particle energies E. At high couplings, where the Schwarzschild radius R_S is much greater than the wavelength lambda, we find that the cross section approaches the classical result for a point particle. At intermediate couplings we find oscillations around the classical limit whose precise form depends on the particle mass. These oscillations give quantum violations of the equivalence principle. At high energies the cross section converges on the geometric-optics value of 27 \pi R_S^2/4, and at low energies we find agreement with an approximation derived by Unruh. When the hole is much smaller than the particle wavelength we confirm that the minimum possible cross section approaches \pi R_S^2/2., Comment: 11 pages, 3 figures

We present a suite of Mathematica-based computer-algebra packages, termed "Kranc", which comprise a toolbox to convert (tensorial) systems of partial differential evolution equations to parallelized C or Fortran code. Kranc can be used as a "rapid prototyping" system for physicists or mathematicians handling very complicated systems of partial differential equations, but through integration into the Cactus computational toolkit we can also produce efficient parallelized production codes. Our work is motivated by the field of numerical relativity, where Kranc is used as a research tool by the authors. In this paper we describe the design and implementation of both the Mathematica packages and the resulting code, we discuss some example applications, and provide results on the performance of an example numerical code for the Einstein equations., Comment: 24 pages, 1 figure. Corresponds to journal version