Source author record

James Healy

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

Numerical-relativity validation of effective-one-body waveforms in the intermediate-mass-ratio regime

One of the open problems in developing binary black hole (BBH) waveforms for gravitational wave astronomy is to model the intermediate mass ratio regime and connect it to the extreme mass ratio regime. A natural approach is to employ the effective one body (EOB) approach to the two-body dynamics that, by design, can cover the entire mass ratio range and naturally incorporates the extreme mass ratio limit. Here we use recently obtained numerical relativity (NR) data with mass ratios $m_1/m_2=(7,15,\,32,\,64,\,128)$ to test the accuracy of the state-of-the-art EOB model TEOBResumS in the intermediate mass ratio regime. We generally find an excellent EOB/NR consistency around merger and ringdown for all mass ratios and for all available subdominant multipoles, except for the $\ell=m=5$ one. This mode can be crucially improved using the new large-mass ratio NR data of this paper. The EOB/NR inspirals are also consistent with the estimated NR uncertainties. We also use several NR datasets taken by different public catalogs to probe the universal behavior of the multipolar hierarchy of waveform amplitudes at merger, that smoothly connects the equal-mass BBH to the test-mass result. Interestingly, the universal behavior is strengthened if the nonoscillatory memory contribution is included in the NR waveform. Future NR simulations with improved accuracy will be necessary to further probe, and possibly quantitatively refine, the TEOBResumS transition from late inspiral to plunge in the intermediate mass ratio regime.

preprint2022arXiv

The Fourth RIT binary black hole simulations catalog: Extension to Eccentric Orbits

This fourth release of the RIT public catalog of numerical relativity black-hole-binary waveforms \url{http://ccrg.rit.edu/~RITCatalog} consists of 1881 accurate simulations that include 446 precessing and 611 nonprecessing quasicircular/inspiraling binary systems with mass ratios $q=m_1/m_2$ in the range $1/128\leq q\leq1$ and individual spins up to $s/m^2=0.95$; and 824 in eccentric orbits in the range $0<e\leq1$. The catalog also provides initial parameters of the binary, trajectory information, peak radiation, and final remnant black hole properties. The waveforms are corrected for the center of mass drifting and are extrapolated to future null infinity. As an application of this waveform catalog we reanalyze all of the peak radiation and remnant properties to find new, simple, correlations among them, valid in the presence of eccentricity, for practical astrophysical usage.

preprint2020arXiv

Adapted gauge to a quasilocal measure of the black holes recoil

We explore different gauge choices in the moving puncture formulation in order to improve the accuracy of a linear momentum measure evaluated on the horizon of the remnant black hole produced by the merger of a binary. In particular, motivated by constant values studies, we design a gauge via a variable shift parameter $mη(\vec{r}(t))$ such that it takes a low asymptotic (and at the orbiting punctures) value, while about the standard value of 2 at the final hole horizon. This choice then follows the remnant black hole as it moves due to its net recoil velocity. We find that this choice keeps the accuracy of the binary evolution and, once the asymptotic value of the parameter $mη$ is chosen about or below 1.0, it produces more accurate results for the recoil velocity than the corresponding evaluation of the radiated linear momentum at infinity, for typical numerical resolutions. We also find that the choice of the $\partial_t$-gauge (at our working resolutions) is more accurate in this regard of computing recoil velocities than the $\partial_0$-gauge. Detailed studies of an unequal mass $q=m_1/m_2=1/3$ nonspinning binary are provided and then verified for other mass ratios $(q=1/2,1/5)$ and spinning $(q=1)$ binary black hole mergers.

preprint2020arXiv

Application of the third RIT binary black hole simulations catalog to parameter estimation of gravitational waves signals from the LIGO-Virgo O1/O2 observational runs

Using exclusively the 777 full numerical waveforms of the third Binary Black Holes RIT catalog, we reanalyze the ten black hole merger signals reported in LIGO/Virgo's O1/O2 observation runs. We obtain binary parameters, extrinsic parameters, and the remnant properties of these gravitational waves events which are consistent with, but not identical to previously presented results. We have also analyzed three additional events (GW170121, GW170304, GW170727) reported in Venumadhav et al. 2019, and found closely matching parameters. We finally assess the accuracy of our waveforms with convergence studies applied to O1/O2 events and found them adequate for current estimation of parameters.

preprint2016arXiv

High Energy Collisions of Black Holes Numerically Revisited

We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly reduces the spurious radiation content, allowing for initial speeds much closer to the speed of light, i.e. $v\sim0.99c$. Using these new techniques, We estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{\text{ADM}}=0.13\pm0.01$. This value differs from the second-order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and by previous numerical estimates $(0.14\pm0.03)$.

preprint2016arXiv

Modeling the source of GW150914 with targeted numerical-relativity simulations

In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predicted by general relativity, can be computed only by full numerical relativity, because analytic approximations fail near the time of merger. Moreover, the nearly-equal masses, moderate spins, and small number of orbits of GW150914 are especially straightforward and efficient to simulate with modern numerical-relativity codes. In this paper, we report the modeling of GW150914 with numerical-relativity simulations, using black-hole masses and spins consistent with those inferred from LIGO's measurement. In particular, we employ two independent numerical-relativity codes that use completely different analytical and numerical methods to model the same merging black holes and to compute the emitted gravitational waveform; we find excellent agreement between the waveforms produced by the two independent codes. These results demonstrate the validity, impact, and potential of current and future studies using rapid-response, targeted numerical-relativity simulations for better understanding gravitational-wave observations.

preprint2016arXiv

Spin flips in generic black hole binaries

We study the spin dynamics of individual black holes in a binary system. In particular we focus on the polar precession of spins and the possibility of a complete flip of spins with respect to the orbital plane. We perform a full numerical simulation that displays these characteristics. We evolve equal mass binary spinning black holes for $t=20,000M$ from an initial proper separation of $d=25M$ down to merger after 48.5 orbits. We compute the gravitational radiation from this system and compare it to 3.5 post-Newtonian generated waveforms finding close agreement. We then further use 3.5 post-Newtonian evolutions to show the extension of this spin {flip-flop} phenomenon to unequal mass binaries. We also provide analytic expressions to approximate the maximum {flip-flop} angle and frequency in terms of the binary spins and mass ratio parameters at a given orbital radius. Finally we discuss the effect this spin {flip-flop} would have on accreting matter and other potential observational effects.

preprint2016arXiv

Unstable flip-flopping spinning binary black holes

We give a unified description of the flip-flop effect in spinning binary black holes and the anti-alignment instability in terms of real and imaginary flip-flop frequencies. We find that this instability is only effective for mass ratios $0.5<q<1$. We provide analytic expressions that determine the region of parameter space for which the instability occurs in terms of maps of the mass ratio and spin magnitudes $(q,α_1,α_2)$. This restricts the priors of parameter estimation techniques for the observation of gravitational waves from binary black holes and it is relevant for astrophysical modeling and final recoil computations of such binary systems.

preprint2015arXiv

Flip-flopping binary black holes

We study binary spinning black holes to display the long term individual spin dynamics. We perform a full numerical simulation starting at an initial proper separation of $d\approx25M$ between equal mass holes and evolve them down to merger for nearly 48 orbits, 3 precession cycles, and half of a flip-flop cycle. The simulation lasts for $t=20000M$ and displays a total change in the orientation of the spin of one of the black holes from initially aligned with the orbital angular momentum to a complete anti-alignment after half of a flip-flop cycle. We compare this evolution with an integration of the 3.5 Post-Newtonian equations of motion and spin evolution to show that this process continuously flip-flops the spin during the lifetime of the binary until merger. We also provide lower order analytic expressions for the maximum flip-flop angle and frequency. We discuss the effects this dynamics may have on spin growth in accreting binaries and on the observational consequences for galactic and supermassive binary black holes.

preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

preprint2014arXiv

Critical Collapse of Scalar Fields Beyond Axisymmetry

We investigate non-spherically symmetric, scalar field collapse of a family of initial data consisting of a spherically symmetric profile with a deformation proportional to the real part of the spherical harmonic $Y_{21}(θ,φ)$. Independent of the strength of the anisotropy in the data, we find that supercritical collapse yields a black hole mass scaling $M_h \propto (p-p^*)^γ$ with $γ\approx 0.37$, a value remarkably close to the critical exponent obtained by Choptuik in his pioneering study in spherical symmetry. We also find hints of discrete self-similarity. However, the collapse experiments are not sufficiently close to the critical solution to unequivocally claim that the detected periodicity is from critical collapse echoing.

preprint2014arXiv

Investigating Binary Black Hole Mergers with Principal Component Analysis

Despite recent progress in numerical simulations of the coalescence of binary black hole systems, highly asymmetric spinning systems and the construction of accurate physical templates remain challenging and computationally expensive. We explore the feasibility of a prompt and robust test of whether the signals exhibit evidence for generic features that can educate new simulations. We form catalogs of numerical relativity waveforms with distinct physical effects and compute the relative probability that a gravitational wave signal belongs to each catalog. We introduce an algorithm designed to perform this task for coalescence signals using principal component analysis of waveform catalogs and Bayesian model selection and demonstrate its effectiveness.

preprint2014arXiv

Remnant mass, spin, and recoil from spin aligned black-hole binaries

We perform a set of 36 nonprecessing black-hole binary simulations with spins either aligned or counteraligned with the orbital angular momentum in order to model the final mass, spin, and recoil of the merged black hole as a function of the individual black hole spin magnitudes and the mass ratio of the progenitors. We find that the maximum recoil for these configurations is $V_{max}=526\pm23\,km/s$, which occurs when the progenitor spins are maximal, the mass ratio is $q_{max}=m_1/m_2=0.623\pm0.038$, the smaller black-hole spin is aligned with the orbital angular momentum, and the larger black-hole spin is counteraligned ($α_1=-α_2=1$). This maximum recoil is about $80\,km/s$ larger than previous estimates, but most importantly, because the maximum occurs for smaller mass ratios, the probability for a merging binary to recoil faster than $400\,km/s$ can be as large as $17\%$, while the probability for recoils faster than $250\, km/s$ can be as large as $45\%$. We provide explicit phenomenological formulas for the final mass, spin, and recoil as a function of the individual BH spins and the mass difference between the two black holes. Here we include terms up through fourth-order in the initial spins and mass difference, and find excellent agreement (within a few percent) with independent results available in the literature. The maximum radiated energy is $E_{\rm rad}/m\approx11.3\%$ and final spin $α_{\rm rem}^{\rm max}\approx0.952$ for equal mass, aligned maximally spinning binaries.

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.

preprint2013arXiv

Template Mode Hierarchies for Binary Black Hole Mergers

Matched filtering is a popular data analysis framework used to search for gravitational wave signals emitted by compact object binaries. The templates used in matched filtering searches are constructed predominantly from the quadrupolar mode because this mode is the energetically most dominant channel. However, for highly precessing binaries or binaries with moderately large mass ratios, significant power is also carried by higher-order modes. We investigate how the inclusion of higher modes in the templates increases the prospects for detecting gravitational waves. Specifically, we use numerical relativity waveforms from the late inspiral and coalescence of binary black holes to identify mode hierarchies that cover the sky of binary orientations. We show that the ordering in these hierarchies depends on the characteristics of the binary system and the mode strengths. Our study demonstrates that detecting moderately high precessing or unequal mass binaries requires the inclusion of higher modes in the templates banks.

preprint2012arXiv

Impact of Higher-order Modes on the Detection of Binary Black Hole Coalescences

The inspiral and merger of black-hole binary systems are a promising source of gravitational waves. The most effective method to look for a signal with a well understood waveform, such as the binary black hole signal, is matched filtering against a library of model waveforms. Current model waveforms are comprised solely of the dominant radiation mode, the quadrupole mode, although it is known that there can be significant power in the higher-order modes. The binary black hole waveforms produced by numerical relativity are accurate through late inspiral, merger, and ringdown and include the higher-order modes. The available numerical-relativity waveforms span an increasing portion of the physical parameter space of unequal mass, spin and precession. In this paper, we investigate the degree to which gravitational-wave searches could be improved by the inclusion of higher modes in the model waveforms, for signals with a variety of initial mass ratios and generic spins. Our investigation studies how well the quadrupole-only waveform model matches the signal as a function of the inclination and orientation of the source and how the modes contribute to the distance reach into the Universe of Advanced LIGO for a fixed set of internal source parameters. The mismatch between signals and quadrupole-only waveform can be large, dropping below 0.97 for up to 65% of the source-sky for the non-precessing cases we studied, and over a larger area in one precessing case. There is a corresponding 30% increase in detection volume that could be achieved by adding higher modes to the search; however, this is mitigated by the fact that the mismatch is largest for signals which radiate the least energy and to which the search is therefore least sensitive. Likewise, the mismatch is largest in directions from the source along which the least energy is radiated.

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.

preprint2011arXiv

Late Inspiral and Merger of Binary Black Holes in Scalar-Tensor Theories of Gravity

Gravitational wave observations will probe non-linear gravitational interactions and thus enable strong tests of Einstein's theory of general relativity. We present a numerical relativity study of the late inspiral and merger of binary black holes in scalar-tensor theories of gravity. We consider black hole binaries in an inhomogeneous scalar field, specifically binaries inside a scalar field bubble, in some cases with a potential. We calculate the emission of dipole radiation. We also show how these configurations trigger detectable differences between gravitational waves in scalar-tensor gravity and the corresponding waves in general relativity. We conclude that, barring an external mechanism to induce dynamics in the scalar field, scalar-tensor gravity binary black holes alone are not capable of awaking a dormant scalar field, and are thus observationally indistinguishable from their general relativistic counterparts.

preprint2011arXiv

Mergers of Supermassive Black Holes in Astrophysical Environments

Modeling the late inspiral and merger of supermassive black holes is central to understanding accretion processes and the conditions under which electromagnetic emission accompanies gravitational waves. We use fully general relativistic, hydrodynamics simulations to investigate how electromagnetic signatures correlate with black hole spins, mass ratios, and the gaseous environment in this final phase of binary evolution. In all scenarios, we find some form of characteristic electromagnetic variability whose pattern depends on the spins and binary mass ratios. Binaries in hot accretion flows exhibit a flare followed by a sudden drop in luminosity associated with the plunge and merger, as well as quasi-periodic oscillations correlated with the gravitational waves during the inspiral. Conversely, circumbinary disk systems are characterized by a low luminosity of variable emission, suggesting challenging prospects for their detection.

preprint2010arXiv

Exploring the Use of Numerical Relativity Waveforms in Burst Analysis of Precessing Black Hole Mergers

Recent years have witnessed tremendous progress in numerical relativity and an ever improving performance of ground-based interferometric gravitational wave detectors. In preparation for Advanced LIGO and a new era in gravitational wave astronomy, the numerical relativity and gravitational wave data analysis communities are collaborating to ascertain the most useful role for numerical relativity waveforms in the detection and characterization of binary black hole coalescences. In this paper, we explore the detectability of equal mass, merging black hole binaries with precessing spins and total mass M_T in [80,350]Msol, using numerical relativity waveforms and template-less search algorithms designed for gravitational wave bursts. In particular, we present a systematic study using waveforms produced by the MAYAKRANC code that are added to colored, Gaussian noise and analyzed with the Omega burst search algorithm. Detection efficiency is weighed against the orientation of one of the black-hole's spin axes. We find a strong correlation between the detection efficiency and the radiated energy and angular momentum, and that the inclusion of the l=2, m=+/-1,0 modes, at a minimum, is necessary to account for the full dynamics of precessing systems.

preprint2010arXiv

Final Mass and Spin of Merged Black Holes and the Golden Black Hole

We present results on the mass and spin of the final black hole from mergers of equal mass, spinning black holes. The study extends over a broad range of initial orbital configurations, from direct plunges to quasi-circular inspirals to more energetic orbits (generalizations of Newtonian elliptical orbits). It provides a comprehensive search of those configurations that maximize the final spin of the remnant black hole. We estimate that the final spin can reach a maximum spin $a/M_h \approx 0.99\pm 0.01$ for extremal black hole mergers. In addition, we find that, as one increases the orbital angular momentum from small values, the mergers produce black holes with mass and spin parameters $\lbrace M_h/M, a/M_h \rbrace$ ~spiraling around the values $\lbrace \hat M_h/M, \hat a/M_h \rbrace$ of a {\it golden} black hole. Specifically, $(M_h-\hat M_h)/M \propto e^{\pm B\,ϕ}\cosϕ$ and $(a-\hat a)/M_h \propto e^{\pm C\,ϕ}\sinϕ$, with $ϕ$ a monotonically growing function of the initial orbital angular momentum. We find that the values of the parameters for the \emph{golden} black hole are those of the final black hole obtained from the merger of a binary with the corresponding spinning black holes in a quasi-circular inspiral.

preprint2010arXiv

Intrinsic selection biases of ground-based gravitational wave searches for high-mass BH-BH mergers

The next generation of ground-based gravitational wave detectors may detect a few mergers of comparable-mass M\simeq 100-1000 Msun ("intermediate-mass'', or IMBH) spinning black holes. Black hole spin is known to have a significant impact on the orbit, merger signal, and post-merger ringdown of any binary with non-negligible spin. In particular, the detection volume for spinning binaries depends significantly on the component black hole spins. We provide a fit to the single-detector and isotropic-network detection volume versus (total) mass and arbitrary spin for equal-mass binaries. Our analysis assumes matched filtering to all significant available waveform power (up to l=6 available for fitting, but only l<= 4 significant) estimated by an array of 64 numerical simulations with component spins as large as S_{1,2}/M^2 <= 0.8. We provide a spin-dependent estimate of our uncertainty, up to S_{1,2}/M^2 <= 1. For the initial (advanced) LIGO detector, our fits are reliable for $M\in[100,500]M_\odot$ ($M\in[100,1600]M_\odot$). In the online version of this article, we also provide fits assuming incomplete information, such as the neglect of higher-order harmonics. We briefly discuss how a strong selection bias towards aligned spins influences the interpretation of future gravitational wave detections of IMBH-IMBH mergers.