BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sascha Husa

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COastro-ph.GA

Catalog footprint

What is connected

43works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

43 published item(s)

preprint2025arXiv

Eccentric or circular? A reanalysis of binary black hole gravitational wave events for orbital eccentricity signatures

We present a reanalysis of 17 gravitational-wave events detected with Advanced LIGO and Advanced Virgo in their first three observing runs, using the new IMRPhenomTEHM model -- a phenomenological time-domain multipolar waveform model for aligned-spin black-hole binaries in elliptical orbits with two eccentric parameters: eccentricity and mean anomaly. We also analyze all events with the underlying quasi-circular model IMRPhenomTHM to study the impact of including eccentricity and compare the eccentric and quasi-circular binary hypotheses. The high computational efficiency of IMRPhenomTEHM enables us to explore the impact of two different eccentricity priors -- uniform and log-uniform -- as well as different sampler and data settings. We find evidence for eccentricity in two publicly available LVK events, GW200129 and GW200208_22, with Bayes factors favoring the eccentric hypothesis over the quasi-circular aligned-spin scenario: $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[1.30^{+0.15}_{-0.15}, 5.14^{+0.15}_{-0.15}\right]$ and $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[0.49^{+0.08}_{-0.08}, 1.14^{+0.08}_{-0.08}\right]$, respectively. Additionally, the two high-mass events GW190701 and GW190929 exhibit potential eccentric features. For all four events, we conduct further analyses to study the impact of different sampler settings. We also investigate waveform systematics by exploring the support for spin precession using IMRPhenomTPHM and NRSur7dq4, offering new insights into the formation channels of detected binaries. Our results highlight the importance of considering eccentric waveform models in future observing runs, alongside precessing models, as they can help mitigate potential biases in parameter estimation studies. This will be particularly relevant with the expected increase in the diversity of the binary black hole population with new detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

First eccentric inspiral-merger-ringdown analysis of neutron star-black hole mergers

The gravitational wave event GW200105 was the first confident neutron star-black hole (NSBH) merger identified by the LIGO-Virgo-KAGRA collaboration. A recent analysis by Morras et al. with an eccentric precessing waveform model that describes the inspiral phase of the $l=2$ and $m=\{0,\pm 2\}$ modes has identified this event as the first NSBH merger with strong evidence of orbital eccentricity. In this paper we perform the first analysis of this event with an aligned-spin eccentric waveform model that describes the full inspiral, merger, and ringdown, includes subdominant harmonics, and is partially calibrated to numerical relativity simulations. This analysis confirms the results and finds evidence in favor of eccentricity even with a log-uniform prior in eccentricity. We also analyze the NSBH events GW200115 and GW230529, completing the analysis of all NSBHs with IMRPhenomTEHM, and find that these signal are consistent with vanishing eccentricity. Finally, we briefly discuss computational challenges when performing the analysis with time-domain eccentric waveform models.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Hybrid waveforms for precessing quasi-circular binary systems

The demand for long and accurate gravitational waveforms is increasing as we prepare for the next generation of detectors and seek to improve current waveform models. However, numerical relativity waveforms, while highly accurate, are often too short for these applications due to their high computational cost. Hybrid waveforms, which stitch together gravitational wave signals from different modeling approaches, provide a way to generate complete inspiral-merger-ringdown signals. While hybridization is well-established for aligned-spin systems, precession introduces additional complexities due to gauge ambiguities, frame dependence, or spin dynamics. Here we study the challenges associated with alignment of precessing waveforms and present a systematic approach for constructing hybrid waveforms of precessing quasi-circular systems. Our approach relies on minimal assumptions about the merger waveforms and employs the quadrupole-aligned frame to mitigate mode-mixing. Our method is designed to be robust and broadly applicable, imposing minimal constraints on the input waveforms. This framework expands the applicability of hybridization techniques, facilitating flexible hybrid construction for parameter estimation, model calibration, and gravitational-wave data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Time-domain phenomenological multipolar waveforms for aligned-spin binary black holes in elliptical orbits

We introduce IMRPhenomTEHM, a new phenomenological time-domain model for eccentric aligned-spin binary black holes. Building upon the accurate quasi-circular IMRPhenomTHM model, IMRPhenomTEHM integrates the eccentric post-Newtonian (PN) dynamics and introduces eccentric corrections into the waveform multipoles up to 3PN, including spin effects. The model incorporates the dominant (2, $\pm$2) spherical harmonic mode, as well as the subdominant modes (2, $\pm$1), (3, $\pm$3), (4, $\pm$4), and (5, $\pm$5), assuming the binary has circularized by the time of merger. This approach ensures a smooth transition to the non-eccentric limit, providing an accurate quasi-circular limit against the IMRPhenomTHM model. When comparing against 28 public eccentric numerical relativity simulations from the Simulating eXtreme Spacetimes catalog, IMRPhenomTEHM achieves lower than 2% unfaithfulness, confirming its accurate description without calibration to numerical relativity eccentric datasets. IMRPhenomTEHM provides a reliable description of the evolution of eccentric black hole binaries with aligned spins and eccentricities lower than $e=0.4$ at a frequency of 10 Hz, making it suitable for upcoming gravitational-wave observing runs. We validate the model's accuracy through parameter estimation studies, recovering injected parameters within 90% credible intervals for three numerical relativity eccentric simulations and reanalyzing GW150914 and GW190521, obtaining results consistent with the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A detailed analysis of GW190521 with phenomenological waveform models

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

New twists in compact binary waveform modelling: a fast time domain model for precession

We present IMRPhenomTPHM, a phenomenological model for the gravitational wave signals emitted by the coalescence of quasi-circular precessing binary black holes systems. The model is based on the "twisting up" approximation, which maps non-precessing signals to precessing ones in terms of a time dependent rotation described by three Euler angles, and which has been utilized in several frequency domain waveform models that have become standard tools in gravitational wave data analysis. Our model is however constructed in the time domain, which allows several improvements over the frequency domain models: we do not use the stationary phase approximation, we employ a simple approximation for the precessing Euler angles for the ringdown signal, and we implement a new method for computing the Euler angles through the evolution of the spin dynamics of the system, which is more accurate and also computationally efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A first survey of spinning eccentric black hole mergers: Numerical relativity simulations, hybrid waveforms, and parameter estimation

We analyze a new numerical relativity dataset of spinning but nonprecessing binary black holes on eccentric orbits, with eccentricities from approximately $0.1$ to $0.5$, with dimensionless spins up to $0.75$ included at mass ratios $q=m_1/m_2 = (1, 2)$, and further nonspinning binaries at mass ratios $q = (1.5, 3, 4)$. A comparison of the final mass and spin of these simulations with noneccentric data extends previous results in the literature on circularization of eccentric binaries to the spinning case. For the $(l,m)=(2,2)$ spin-weighted spherical harmonic mode we construct eccentric hybrid waveforms that connect the numerical relativity data to a post-Newtonian description for the inspiral, and we discuss the limitations in the current knowledge about post-Newtonian theory which complicate the generation of eccentric hybrid waveforms. We also perform a Bayesian parameter estimation study, quantifying the parameter biases introduced when using three different quasicircular waveform models to estimate the parameters of highly eccentric binary systems. We find that the used aligned-spin quasicircular model including higher order modes produces lower bias in certain parameters than the nonprecessing quasicircular model without higher order modes and the quasicircular precessing model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating the evaluation of inspiral-merger-ringdown waveforms with adapted grids

This paper presents an algorithm to accelerate the evaluation of inspiral-merger-ringdown waveform models for gravitational wave data analysis. While the idea can also be applied in the time domain, here we focus on the frequency domain, which is most typically used to reduced computational cost in gravitational wave data analysis. Our work extends the idea of multibanding, which has been developed to accelerate frequency domain waveforms, to include the merger and ringdown and spherical harmonics beyond the dominant quadrupole spherical harmonic. The original method is based on a heuristic algorithm based on the inspiral to de-refine the equi-spaced frequency grid used for data analysis where a coarser grid is sufficient for accurate evaluation of a waveform model. Here we use a different criterion, based on the local interpolation error, which is more flexible and can easily be adapted to general waveforms, if their phenomenology is understood. We discuss our implementation in the LIGO Algorithm Library for the PhenomXHM frequency domain model, and report the acceleration in different parts of the parameter space of compact binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of eccentricity on the gravitational wave searches for binary black holes: High mass case

The possible formation of stellar-mass binary black holes through dynamical interactions in dense stellar environments predicts the existence of binaries with non-negligible eccentricity in the frequency band of ground-based gravitational wave detectors; the detection of binary black hole mergers with measurable orbital eccentricity would validate the existence of this formation channel. Waveform templates currently used in the matched-filter gravitational-wave searches of LIGO-Virgo data neglect effects of eccentricity which is expected to reduce their efficiency to detect eccentric binary black holes. Meanwhile, the sensitivity of coherent unmodeled gravitational-wave searches (with minimal assumptions about the signal model) have been shown to be largely unaffected by the presence of even sizable orbital eccentricity. In this paper, we compare the performance of two state-of-the-art search algorithms recently used by LIGO and Virgo to search for binary black holes in the second Observing Run (O2), quantifying their search sensitivity by injecting numerical-relativity simulations of inspiral-merger-ringdown eccentric waveforms into O2 LIGO data. Our results show that the matched-filter search PyCBC performs better than the unmodeled search cWB for the high chirp mass ($>20 M_{\odot}$) and low eccentricity region ($e_{30 Hz} < 0.3$) of parameter space. For moderate eccentricities and low chirp mass, on the other hand, the unmodeled search is more sensitive than the modeled search.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Setting the cornerstone for the IMRPhenomX family of models for gravitational waves from compact binaries: The dominant harmonic for non-precessing quasi-circular black holes

In this paper we present IMRPhenomXAS, a thorough overhaul of the IMRPhenomD [1,2] waveform model, which describes the dominant $l=2, \:| m | = 2$ spherical harmonic mode of non-precessing coalescing black holes in terms of piecewise closed form expressions in the frequency domain. Improvements include in particular the accurate treatment of unequal spin effects, and the inclusion of extreme mass ratio waveforms. IMRPhenomD has previously been extended to approximately include spin precession [3] and subdominant spherical harmonics [4], and with its extensions it has become a standard tool in gravitational wave parameter estimation. Improved extensions of IMRPhenomXAS are discussed in companion papers [5,6].

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simple procedures to reduce eccentricity of binary black hole simulations

We present simple procedures to construct quasi-circular initial data for numerical evolutions of binary black hole spacetimes. Our method consists of using Post-Newtonian theory in three ways: first to provide an initial guess for the initial momenta at 3.5PN order that implies low residual eccentricity, second to measure the resulting eccentricity, and third to calculate corrections to the momenta or initial separation which further reduce the eccentricity. Regarding the initial guess, we compare numerical evolutions in post-Newtonian theory to the post-circular and post-post-circular analytical approximations to quasi-circular data. We discuss a robust fitting procedure to measure eccentricity from numerical simulations using the orbital frequency, and derive from the quasi-Keplerian parametrization at 1PN oder the correction factors for the tangential and radial momentum components required to achieve reduce the measured eccentricity to zero. We first test our procedure integrating PN equations of motion at $3.5$PN where low eccentric initial data is easily obtained, and then apply our method to sets of binary black hole numerical relativity simulations with different mass ratios (q=1,2,...,8), spin configurations and separations. Our set of simulations contains non-spinning, spin-aligned and precessing simulations. We observe that the iterative procedure produces low eccentric simulations with eccentricities of the order (10^{-4}) with only one iteration. The simplicity of the procedure allows to obtain low eccentric NR simulations easily and saving computational resources. Moreover, the analytical PN formulas derived in this paper will be useful to generate eccentric hybrid waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards the routine use of subdominant harmonics in gravitational-wave inference: re-analysis of GW190412 with generation X waveform models

We re-analyse the gravitational-wave event GW190412 with state-of-the-art phenomenological waveform models. This event, which has been associated with a black hole merger, is interesting due to the significant contribution from subdominant harmonics. We use both frequency-domain and time-domain waveform models. The PhenomX waveform models constitute the fourth generation of frequency-domain phenomenological waveforms for black hole binary coalescence; they have more recently been complemented by the time-domain PhenomT models, which open up new strategies to model precession and eccentricity, and to perform tests of general relativity with the phenomenological waveforms approach. Both PhenomX and PhenomT have been constructed with similar techniques and accuracy goals, and due to their computational efficiency this "generation X" model family allows the routine use of subdominant spherical harmonics in Bayesian inference. We show the good agreement between these and other state-of-the-art waveform models for GW190412, and discuss the improvements over the previous generation of phenomenological waveform models. We also discuss practical aspects of Bayesian inference such as run convergence, variations of sampling parameters, and computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Validity of common modelling approximations for precessing binary black holes with higher-order modes

The current paradigm for constructing waveforms from precessing compact binaries is to first construct a waveform in a non-inertial, co-precessing binary source frame followed by a time-dependent rotation to map back to the physical, inertial frame. A key insight in the construction of these models is that the co-precessing waveform can be effectively mapped to some equivalent aligned spin waveform. Secondly, the time-dependent rotation implicitly introduces $m$-mode mixing, necessitating an accurate description of higher-order modes in the co-precessing frame. We assess the efficacy of this modelling strategy in the strong field regime using Numerical Relativity simulations. We find that this framework allows for the highly accurate construction of $(2,\pm 2)$ modes in our data set, while for higher order modes, especially the $(2,|1|), (3,|2|)$ and $(4,|3|)$ modes, we find rather large mismatches. We also investigate a variant of the approximate map between co-precessing and aligned spin waveforms, where we only identify the slowly varying part of the time domain co-precessing waveforms with the aligned-spin one, but find no significant improvement. Our results indicate that the simple paradigm to construct precessing waveforms does not provide an accurate description of higher order modes in the strong-field regime, and demonstrate the necessity for modelling mode asymmetries and mode-mixing to significantly improve the description of precessing higher order modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Free hyperboloidal evolution in spherical symmetry

We address the hyperboloidal initial value problem in the context of Numerical Relativity, motivated by its evolution on hyperboloidal slices: smooth spacelike slices that reach future null infinity, the "location" in spacetime where radiation is to be extracted. Our approach uses the BSSN and Z4 formulations and a time-independent conformal factor. The resulting system of PDEs includes formally diverging terms at null infinity. Here we discuss a regularized numerical scheme in spherical symmetry. A critical ingredient are the gauge conditions, which control the treatment of future null infinity. Stable numerical evolutions have been performed with regular and black hole initial data on a hyperboloidal slice. A sufficiently large scalar field perturbation will create a black hole, whose final stationary state is different from the trumpet initial data derived here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparison of subdominant gravitational wave harmonics between post-Newtonian and numerical relativity calculations and construction of multi-mode hybrids

Gravitational waveforms which describe the inspiral, merger and ringdown of coalescing binaries are usually constructed by synthesising information from perturbative descriptions, in particular post-Newtonian theory and black-hole perturbation theory, with numerical solutions of the full Einstein equations. In this paper we discuss the "glueing" of numerical and post-Newtonian waveforms to produce hybrid waveforms which include subdominant spherical harmonics ("higher order modes"), and focus in particular on the process of consistently aligning the waveforms, which requires a comparison of both descriptions and a discussion of their imprecisions. We restrict to the non-precessing case, and illustrate the process using numerical waveforms of up to mass ratio $q=18$ produced with the BAM code, and publicly available waveforms from the SXS catalogue. The results also suggest new ways of analysing finite radius errors in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Frequency-domain gravitational waves from non-precessing black-hole binaries. I. New numerical waveforms and anatomy of the signal

In this paper we discuss the anatomy of frequency-domain gravitational-wave signals from non-precessing black-hole coalescences with the goal of constructing accurate phenomenological waveform models. We first present new numerical-relativity simulations for mass ratios up to 18 including spins. From a comparison of different post-Newtonian approximants with numerical-relativity data we select the uncalibrated SEOBNRv2 model as the most appropriate for the purpose of constructing hybrid post-Newtonian/numerical-relativity waveforms, and we discuss how we prepare time-domain and frequency-domain hybrid data sets. We then use our data together with results in the literature to calibrate simple explicit expressions for the final spin and radiated energy. Equipped with our prediction for the final state we then develop a simple and accurate merger-ringdown-model based on modified Lorentzians in the gravitational wave amplitude and phase, and we discuss a simple method to represent the low frequency signal augmenting the TaylorF2 post-Newtonian approximant with terms corresponding to higher orders in the post-Newtonian expansion. We finally discuss different options for modelling the small intermediate frequency regime between inspiral and merger-ringdown. A complete phenomenological model based on the present work is presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Frequency-domain gravitational waves from non-precessing black-hole binaries. II. A phenomenological model for the advanced detector era

We present a new frequency-domain phenomenological model of the gravitational-wave signal from the inspiral, merger and ringdown of non-precessing (aligned-spin) black-hole binaries. The model is calibrated to 19 hybrid effective-one-body--numerical-relativity waveforms up to mass ratios of 1:18 and black-hole spins of $|a/m| \sim 0.85$ ($0.98$ for equal-mass systems). The inspiral part of the model consists of an extension of frequency-domain post-Newtonian expressions, using higher-order terms fit to the hybrids. The merger-ringdown is based on a phenomenological ansatz that has been significantly improved over previous models. The model exhibits mismatches of typically less than 1\% against all 19 calibration hybrids, and an additional 29 verification hybrids, which provide strong evidence that, over the calibration region, the model is sufficiently accurate for all relevant gravitational-wave astronomy applications with the Advanced LIGO and Virgo detectors. Beyond the calibration region the model produces physically reasonable results, although we recommend caution in assuming that \emph{any} merger-ringdown waveform model is accurate outside its calibration region. As an example, we note that an alternative non-precessing model, SEOBNRv2 (calibrated up to spins of only 0.5 for unequal-mass systems), exhibits mismatch errors of up to 10\% for high spins outside its calibration region. We conclude that waveform models would benefit most from a larger number of numerical-relativity simulations of high-aligned-spin unequal-mass binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impact of gravitational radiation higher order modes on single aligned-spin gravitational wave searches for binary black holes

Current template-based gravitational wave searches for compact binary coalescences (CBC) use waveform models that neglect the higher order modes content of the gravitational radiation emitted, considering only the quadrupolar $(\ell,|m|)=(2,2)$ modes. We study the effect of such a neglection for the case of aligned-spin CBC searches for equal-spin (and non-spinning) binary black holes in the context of two versions of Advanced LIGO: the upcoming 2015 version, known as early Advanced LIGO (eaLIGO) and its Zero-Detuned High Energy Power version, that we will refer to as Advanced LIGO (AdvLIGO). In addition, we study the case of a non-spinning search for initial LIGO (iLIGO). We do this via computing the effectualness of the aligned-spin SEOBNRv1 ROM waveform family, which only considers quadrupolar modes, towards hybrid post-Newtonian/Numerical Relativity waveforms which contain higher order modes. We find that for all LIGO versions, losses of more than $10\%$ of events occur for mass ratio $q\geq6$ and $M \geq 100M_\odot$ due to the neglection of higher modes. Moreover, for iLIGO and eaLIGO, losses notably increase up to $(39,23)\%$ respectively for the highest mass $(220M_\odot)$ and mass ratio ($q=8$) studied. For the case of early AdvLIGO, losses of $10\%$ occur for $M>50M_\odot$ and $q\geq6$. Neglection of higher modes leads to observation-averaged systematic parameter biases towards lower spin, total mass and chirp mass. For completeness, we perform a preliminar, non-exhaustive comparison of systematic biases to statistical errors. We find that, for a given SNR, systematic biases dominate over statistical errors at much lower total mass for eaLIGO than for AdvLIGO.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spherical symmetry as a test case for unconstrained hyperboloidal evolution

We consider the hyperboloidal initial value problem in numerical relativity, motivated by the goal to evolve radiating compact objects such as black hole binaries with a numerical grid that includes null infinity. Unconstrained evolution schemes promise optimal efficiency, but are difficult to regularize at null infinity, where the compactified Einstein equations are formally singular. In this work we treat the spherically symmetry case, which already poses nontrivial problems and constitutes an important first step. We have carried out stable numerical evolutions with the generalized BSSN and Z4 equations coupled to a scalar field. The crucial ingredients have been to find an appropriate evolution equation for the lapse function and to adapt constraint damping terms to handle null infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A simple model of complete precessing black-hole-binary gravitational waveforms

The construction of a model of the gravitational-wave (GW) signal from generic configurations of spinning-black-hole binaries, through inspiral, merger and ringdown, is one of the most pressing theoretical problems in the build-up to the era of GW astronomy. We present the first such model in the frequency domain, "PhenomP", which captures the basic phenomenology of the seven-dimensional parameter space of binary configurations with only three key physical parameters. Two of these (the binary's mass ratio and an effective total spin parallel to the orbital angular momentum, which determines the inspiral rate) define an underlying non-precessing-binary model. The non-precessing-binary waveforms are then "twisted up" with approximate expressions for the precessional motion, which require only one additional physical parameter, an effective precession spin, $χ_p$. All other parameters (total mass, sky location, orientation and polarisation, and initial phase) can be specified trivially. The model is constructed in the frequency domain, which will be essential for efficient GW searches and source measurements. We have tested the model's fidelity for GW applications by comparison against hybrid post-Newtonian-numerical-relativity waveforms at a variety of configurations --although we did not use these numerical simulations in the construction of the model. Our model can be used to develop GW searches, to study the implications for astrophysical measurements, and as a simple conceptual framework to form the basis of generic-binary waveform modelling in the advanced-detector era.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational-wave observations of binary black holes: Effect of non-quadrupole modes

We study the effect of non-quadrupolar modes in the detection and parameter estimation of gravitational waves (GWs) from non-spinning black-hole binaries. We evaluate the loss of signal-to-noise ratio and the systematic errors in the estimated parameters when one uses a quadrupole-mode template family to detect GW signals with all the relevant modes, for target signals with total masses $20 M_\odot \leq M \leq 250 M_\odot$ and mass ratios $1 \leq q \leq 18$. Target signals are constructed by matching numerical-relativity simulations describing the late inspiral, merger and ringdown of the binary with post-Newtonian/effective-one-body waveforms describing the early inspiral. We find that waveform templates modeling only the quadrupolar modes of the GW signal are sufficient (loss of detection rate $< 10\%$) for the detection of GWs with mass ratios $q\leq4$ using advanced GW observatories. Neglecting the effect of non-quadrupole modes will introduce systematic errors in the estimated parameters. The systematic errors are larger than the expected $1\,σ$ statistical errors for binaries with large, unequal masses ($q\gtrsim4, M \gtrsim 150 M_\odot$), for sky-averaged signal-to-noise ratios larger than $8$. We provide a summary of the regions in the parameter space where neglecting non-quadrupole modes will cause unacceptable loss of detection rates and unacceptably large systematic biases in the estimated parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unconstrained hyperboloidal evolution of black holes in spherical symmetry with GBSSN and Z4c

We consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry already poses nontrivial problems and constitutes an important first step to regularize the resulting singular PDEs. We evolve the Einstein equations in their generalized BSSN and Z4 formulations coupled to a massless self-gravitating scalar field. Stable numerical evolutions are achieved for black hole initial data, and critically rely on the construction of appropriate gauge conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing the validity of the single-spin approximation in inspiral-merger-ringdown waveforms

Gravitational-wave signals from black-hole binaries with non-precessing spins are described by four parameters -- each black hole's mass and spin. It has been shown that the dominant spin effects can be modeled by a \emph{single} spin parameter, leading to the development of several \emph{three-parameter} waveform models. Previous studies indicate that these models should be adequate for gravitational-wave detection. In this paper we focus on the systematic biases that would result from using them to estimate binary parameters, and consider a one-parameter family of configurations at one choice of mass ratio and effective single spin. We find that for low-mass binaries within that family of configurations, where the observable waveform is dominated by the inspiral, the systematic bias in all physical parameters is smaller than the parameter uncertainty due to degeneracies between the mass ratio and the spins, at least up to signal-to-noise ratios (SNRs) of 50. For higher-mass binaries, where the merger and ringdown make a greater contribution to the observed signal, the bias in the mass ratio is comparable to its uncertainty at SNRs of only $\sim$30, and the bias in the measurement of the total spin is \emph{larger} than the uncertainty defined by the 90% confidence region even at an SNR of only 10. Although this bias may be mitigated in future models by a better choice of single-effective-spin parameter, these results suggest that it may be possible to accurately measure \emph{both} black-hole spins in intermediate-mass binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Towards models of gravitational waveforms from generic binaries: A simple approximate mapping between precessing and non-precessing inspiral signals

One of the greatest theoretical challenges in the build-up to the era of second-generation gravitational-wave detectors is the modeling of generic binary waveforms. We introduce an approximation that has the potential to significantly simplify this problem. We show that generic precessing-binary inspiral waveforms (covering a seven-dimensional space of intrinsic parameters) can be mapped to a two-dimensional space of non-precessing binaries, characterized by the mass ratio and a single effective total spin. The mapping consists of a time-dependent rotation of the waveforms into the quadrupole-aligned frame, and is extremely accurate (matches $> 0.99$ with parameter biases in the total spin of $Δχ\leq 0.04$), even in the case of transitional precession. In addition, we demonstrate a simple method to construct hybrid post-Newtonian--numerical-relativity precessing-binary waveforms in the quadrupole-aligned frame, and provide evidence that our approximate mapping can be used all the way to the merger. Finally, based on these results, we outline a general proposal for the construction of generic waveform models, which will be the focus of future work.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An efficient iterative method to reduce eccentricity in numerical-relativity simulations of compact binary inspiral

We present a new iterative method to reduce eccentricity in black-hole-binary simulations. Given a good first estimate of low-eccentricity starting momenta, we evolve puncture initial data for ~4 orbits and construct improved initial parameters by comparing the inspiral with post-Newtonian calculations. Our method is the first to be applied directly to the gravitational-wave (GW) signal, rather than the orbital motion. The GW signal is in general less contaminated by gauge effects, which, in moving-puncture simulations, limit orbital-motion-based measurements of the eccentricity to an uncertainty of $Δe \sim 0.002$, making it difficult to reduce the eccentricity below this value. Our new method can reach eccentricities below $10^{-3}$ in one or two iteration steps; we find that this is well below the requirements for GW astronomy in the advanced detector era. Our method can be readily adapted to any compact-binary simulation with GW emission, including black-hole-binary simulations that use alternative approaches, and neutron-star-binary simulations. We also comment on the differences in eccentricity estimates based on the strain $h$, and the Newman-Penrose scalar $Ψ_4$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black-hole hair loss: learning about binary progenitors from ringdown signals

Perturbed Kerr black holes emit gravitational radiation, which (for the practical purposes of gravitational-wave astronomy) consists of a superposition of damped sinusoids termed quasi-normal modes. The frequencies and time-constants of the modes depend only on the mass and spin of the black hole - a consequence of the no-hair theorem. It has been proposed that a measurement of two or more quasi-normal modes could be used to confirm that the source is a black hole and to test if general relativity continues to hold in ultra-strong gravitational fields. In this paper we propose a practical approach to testing general relativity with quasi-normal modes. We will also argue that the relative amplitudes of the various quasi-normal modes encode important information about the origin of the perturbation that caused them. This helps in inferring the nature of the perturbation from an observation of the emitted quasi-normal modes. In particular, we will show that the relative amplitudes of the different quasi-normal modes emitted in the process of the merger of a pair of nonspinning black holes can be used to measure the component masses of the progenitor binary.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reliability of complete gravitational waveform models for compact binary coalescences

With recent advances in post-Newtonian (PN) theory and numerical relativity (NR) it has become possible to construct inspiral-merger-ringdown waveforms by combining both descriptions into one hybrid signal. While addressing the reliability of such waveforms, previous studies have identified the PN contribution as the dominant source of error, which can be reduced by incorporating longer NR simulations. Here we overcome the two outstanding issues that make it difficult to determine the minimum NR simulation length necessary to produce suitably accurate hybrids: (1) the criteria for a GW search is the mismatch between the true waveform and a set of model waveforms, optimized over all waveforms in the model, but for discrete hybrids this optimization was not yet possible. (2) these calculations typically require that numerical waveforms already exist, while we develop an algorithm to estimate hybrid mismatches errors without numerical data. Our procedure relies on combining supposedly equivalent PN models at highest available order with common data in the NR regime, and their difference serves as a measure of the uncertainty assumed in each waveform. Contrary to some earlier studies, we estimate that ~10 NR orbits before merger should allow for the construction of waveform families that are accurate enough for detection in a broad range of parameters, only excluding highly spinning, unequal-mass systems. Nonspinning systems, even with high mass-ratio (q>=20) are well modeled for astrophysically reasonable component masses. The parameter bias is only of the order of 1% for total mass and symmetric mass-ratio and less than 0.1 for the dimensionless spin magnitude. We take the view that similar NR waveform lengths will remain the state of the art in the advanced detector era, and begin to assess the limits of the science that can be done with them.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational-wave detectability of equal-mass black-hole binaries with aligned spins

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than "three times as loud" as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kranc: a Mathematica application to generate numerical codes for tensorial evolution equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Length requirements for numerical-relativity waveforms

One way to produce complete inspiral-merger-ringdown gravitational waveforms from black-hole-binary systems is to connect post-Newtonian (PN) and numerical-relativity (NR) results to create "hybrid" waveforms. Hybrid waveforms are central to the construction of some phenomenological models for GW search templates, and for tests of GW search pipelines. The dominant error source in hybrid waveforms arises from the PN contribution, and can be reduced by increasing the number of NR GW cycles that are included in the hybrid. Hybrid waveforms are considered sufficiently accurate for GW detection if their mismatch error is below 3% (i.e., a fitting factor about 0.97). We address the question of the length requirements of NR waveforms such that the final hybrid waveforms meet this requirement, considering nonspinning binaries with q = M_2/M_1 \in [1,4] and equal-mass binaries with χ= S_i/M_i^2 \in [-0.5,0.5]. We conclude that for the cases we study simulations must contain between three (in the equal-mass nonspinning case) and ten (the χ= 0.5 case) orbits before merger, but there is also evidence that these are the regions of parameter space for which the least number of cycles will be needed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations of black-hole binaries with unequal masses or non-precessing spins: accuracy, physical properties, and comparison with post-Newtonian results

We present gravitational waveforms for the last orbits and merger of black-hole-binary (BBH) systems along two branches of the BBH parameter space: equal-mass binaries with equal non-precessing spins, and nonspinning unequal-mass binaries. The waveforms are calculated from numerical solutions of Einstein's equations for black-hole binaries that complete between six and ten orbits before merger. Along the equal-mass spinning branch, the spin parameter of each BH is $χ_i = S_i/M_i^2 \in [-0.85,0.85]$, and along the unequal-mass branch the mass ratio is $q =M_2/M_1 \in [1,4]$. We discuss the construction of low-eccentricity puncture initial data for these cases, the properties of the final merged BH, and compare the last 8-10 GW cycles up to $Mω= 0.1$ with the phase and amplitude predicted by standard post-Newtonian (PN) approximants. As in previous studies, we find that the phase from the 3.5PN TaylorT4 approximant is most accurate for nonspinning binaries. For equal-mass spinning binaries the 3.5PN TaylorT1 approximant (including spin terms up to only 2.5PN order) gives the most robust performance, but it is possible to treat TaylorT4 in such a way that it gives the best accuracy for spins $χ_i > -0.75$. When high-order amplitude corrections are included, the PN amplitude of the $(\ell=2,m=\pm2)$ modes is larger than the NR amplitude by between 2-4%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tracking the precession of compact binaries from their gravitational-wave signal

We present a simple method to track the precession of a black-hole-binary system, using only information from the gravitational-wave (GW) signal. Our method consists of locating the frame from which the magnitude of the $(\ell=2,|m|=2)$ modes is maximized, which we denote the "quadrupole-aligned" frame. We demonstrate the efficacy of this method when applied to waveforms from numerical simulations. In the test case of an equal-mass nonspinning binary, our method locates the direction of the orbital angular momentum to within $(Δθ, Δϕ) = (0.05^{\circ},0.2^{\circ})$. We then apply the method to a $q = M_2/M_1 = 3$ binary that exhibits significant precession. In general a spinning binary's orbital angular momentum $\mathbf{L}$ is \emph{not} orthogonal to the orbital plane. Evidence that our method locates the direction of $\mathbf{L}$ rather than the normal of the orbital plane is provided by comparison with post-Newtonian (PN) results. Also, we observe that it accurately reproduces similar higher-mode amplitudes to a comparable non-spinning (and therefore non-precessing) binary, and that the frequency of the $(\ell=2,|m|=2)$ modes is consistent with the "total frequency" of the binary's motion. The simple form of the quadrupole-aligned waveform will be useful in attempts to analytically model the inspiral-merger-ringdown (IMR) signal of precessing binaries, and in standardizing the representation of waveforms for studies of accuracy and consistency of source modelling efforts, both numerical and analytical.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Samurai Project: verifying the consistency of black-hole-binary waveforms for gravitational-wave detection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Where do moving punctures go?

Currently the most popular method to evolve black-hole binaries is the ``moving puncture'' method. It has recently been shown that when puncture initial data for a Schwarzschild black hole are evolved using this method, the numerical slices quickly lose contact with the second asymptotically flat end, and end instead on a cylinder of finite Schwarzschild coordinate radius. These slices are stationary, meaning that their geometry does not evolve further. We will describe these results in the context of maximal slices, and present time-independent puncture-like data for the Schwarzschild spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

News from Critical Collapse: Bondi Mass, Tails and Quasinormal Modes

We discuss critical gravitational collapse on the threshold of apparent horizon formation as a model both for the discussion of global aspects of critical collapse and for numerical studies in a compactified context. For our matter model we choose a self-gravitating massless scalar field in spherical symmetry, which has been studied extensively in the critical collapse literature. Our evolution system is based on Bondi coordinates, the mass function is used as an evolution variable to ensure regularity at null infinity. We compute radiation quantities like the Bondi mass and news function and find that they reflect the DSS behavior. Surprisingly, the period of radiation at null infinity is related to the formal result for the leading quasi-normal mode of a black hole with rapidly decreasing mass. Furthermore, our investigations shed some light on global versus local issues in critical collapse, and the validity and usefulness of the concept of null infinity when predicting detector signals.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

SU(2) Cosmological Solitons

We present a class of numerical solutions to the SU(2) nonlinear $σ$-model coupled to the Einstein equations with cosmological constant $Λ\geq 0$ in spherical symmetry. These solutions are characterized by the presence of a regular static region which includes a center of symmetry. They are parameterized by a dimensionless ``coupling constant'' $β$, the sign of the cosmological constant, and an integer ``excitation number'' $n$. The phenomenology we find is compared to the corresponding solutions found for the Einstein-Yang-Mills (EYM) equations with positive $Λ$ (EYM$Λ$). If we choose $Λ$ positive and fix $n$, we find a family of static spacetimes with a Killing horizon for $0 \leq β< β_{max}$. As a limiting solution for $β= β_{max}$ we find a {\em globally} static spacetime with $Λ=0$, the lowest excitation being the Einstein static universe. To interpret the physical significance of the Killing horizon in the cosmological context, we apply the concept of a trapping horizon as formulated by Hayward. For small values of $β$ an asymptotically de Sitter dynamic region contains the static region within a Killing horizon of cosmological type. For strong coupling the static region contains an ``eternal cosmological black hole''.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Type II Critical Collapse of a Self-Gravitating Nonlinear $σ$-Model

We report on the existence and phenomenology of type II critical collapse within the one-parameter family of SU(2) $σ$-models coupled to gravity. Numerical investigations in spherical symmetry show discretely self-similar (DSS) behavior at the threshold of black hole formation for values of the dimensionless coupling constant $\ccbeta$ ranging from 0.2 to 100; at 0.18 we see small deviations from DSS. While the echoing period $Δ$ of the critical solution rises sharply towards the lower limit of this range, the characteristic mass scaling has a critical exponent $γ$ which is almost independent of $\ccbeta$, asymptoting to $0.1185 \pm 0.0005$ at large $\ccbeta$. We also find critical scaling of the scalar curvature for near-critical initial data. Our numerical results are based on an outgoing-null-cone formulation of the Einstein-matter equations, specialized to spherical symmetry. Our numerically computed initial-data critical parameters $p^*$ show 2nd order convergence with the grid resolution, and after compensating for this variation in $p^*$, our individual evolutions are uniformly 2nd order convergent even very close to criticality.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The Asymmetric Merger of Black Holes

We study event horizons of non-axisymmetric black holes and show how features found in axisymmetric studies of colliding black holes and of toroidal black holes are non-generic and how new features emerge. Most of the details of black hole formation and black hole merger are known only in the axisymmetric case, in which numerical evolution has successfully produced dynamical space-times. The work that is presented here uses a new approach to construct the geometry of the event horizon, not by locating it in a given spacetime, but by direct construction. In the axisymmetric case, our method produces the familiar pair-of-pants structure found in previous numerical simulations of black hole mergers, as well as event horizons that go through a toroidal epoch as discovered in the collapse of rotating matter. The main purpose of this paper is to show how new - substantially different - features emerge in the non-axisymmetric case. In particular, we show how black holes generically go through a toroidal phase before they become spherical, and how this fits together with the merger of black holes.

0 citations0 reviewsNo code linkedOpen access