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

Carlos O. Lousto

Carlos O. Lousto 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.GAastro-ph.COhep-thastro-phastro-ph.SRhep-ph

Catalog footprint

What is connected

43works
8topics
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)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modeling the remnant mass, spin, and recoil from unequal-mass, precessing black-hole binaries: The Intermediate Mass Ratio Regime

We revisit the modeling of the properties of the remnant black hole resulting the merger of a black-hole binary as a function of the parameters of the binary. We provide a set of empirical formulas for the final mass, spin and recoil velocity of the final black hole as a function of the mass ratio and individual spins of the progenitor. In order to determine the fitting coefficients for these formulas, we perform a set of 128 new numerical evolutions of precessing, unequal-mass black-hole binaries, and fit to the resulting remnant mass, spin, and recoil. In order to reduce the complexity of the analysis, we chose configurations that have one of the black holes spinning, with dimensionless spin alpha=0.8, at different angles with respect to the orbital angular momentum, and the other non-spinning. In addition to evolving families of binaries with different spin-inclination angles, we also evolved binaries with mass ratios as small as q=1/6. We use the resulting empirical formulas to predict the probabilities of black hole mergers leading to a given recoil velocity, total radiated gravitational energy, and final black hole spin.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black hole binary remnant mass and spin: A new phenomenological formula

We perform a set of 38 numerical simulations of equal-mass BH binaries in a configuration where the BH spins in the binary are equal in both magnitude and direction, to study precession effects. We vary the initial direction of the total spin S with respect to the orbital angular momentum L, covering the 2 dimensional space of orientation angles with 38 configurations consisting of 36 configurations distributed in the azimuthal angle phi and polar angle theta, and two configurations on the poles. In all cases, we set the initial dimensionless BH spins to 0.8. We observe that during the late-inspiral stage, the total angular momentum of the system J remains within 5 deg of its original direction, with the largest changes in direction occurring when the spins are nearly counter-aligned with the orbital angular momentum. We also observe that the angle between S and L is nearly conserved during the inspiral phase. These two dynamical properties allow us to propose a new phenomenological formula for the final mass and spin of merged BHs in terms of the individual masses and spins of the progenitor binary at far separations. We determine coefficients of this formula (in the equal-mass limit) using a least-squares fit to the results of this new set of 38 runs, an additional set of five new configurations with spins aligned/counteraligned with the orbital angular momentum, and over 100 recent simulations. We find that our formulas reproduce the remnant mass and spin of these simulations to within a relative error of 2.5%. We discuss the region of validity of this dynamical picture for precessing unequal-mass binaries. Finally, we perform a statistical study to see the consequence of this new formula for distributions of spin-magnitudes and remnant masses with applications to BH-spin distributions and gravitational radiation in cosmological scenarios involving several mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear Gravitational Recoil from the Mergers of Precessing Black-Hole Binaries

We present results from an extensive study of 88 precessing, equal-mass black-hole binaries with large spins (83 with intrinsic spins of 0.8 and 5 with intrinsic spins of 0.9)and use these data to model new nonlinear contributions to the gravitational recoil imparted to the merged black hole. We find a new effect, the cross kick, that enhances the recoil for partially aligned binaries beyond the hangup kick effect. This has the consequence of increasing the probabilities of recoils larger than 2000 km/s by nearly a factor two, and, consequently, of black holes getting ejected from galaxies, as well as the observation of large differential redshifts/blueshifts in the cores of recently merged galaxies.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Where angular momentum goes in a precessing black hole binary

We evolve a set of 32 equal-mass black-hole binaries with collinear spins (with intrinsic spin magnitudes $|\vec{S}_{1,2}/m^2_{1,2}|=0.8$) to study the effects of precession in the highly nonlinear plunge and merger regimes. We compare the direction of the instantaneous radiated angular momentum, $\hat{δJ}_{\rm rad}(t)$, to the directions of the total angular momentum, $\hat{J}(t)$, and the orbital angular momentum, $\hat{L}(t)$. We find that $\hat{δJ}_{\rm rad}(t)$ approximately follows $\hat{L}$ throughout the evolution. During the orbital evolution and merger, we observe that the angle between $\vec{L}$ and total spin $\vec{S}$ is approximately conserved to within $1^\circ$, which allows us to propose and test models for the merger remnant's mass and spin. For instance, we verify that the \hangup effect is the dominant effect and largely explains the observed total energy and angular momentum radiated by these precessing systems. We also verify that the total angular momentum, which significantly decreases in magnitude during the inspiral, varies in direction by less than $\sim 5^\circ$. The maximum variation in the direction of $\vec J$ occurs when the spins are nearly antialigned with the orbital angular momentum. Based on our results, we conjecture that transitional precession, which would lead to large variations in the direction of $\vec J$, is not possible for similar-mass binaries and would require a mass ratio $m_1/m_2\lesssim1/4$.

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

Exploring the Outer Limits of Numerical Relativity

We perform several black-hole binary evolutions using fully nonlinear numerical relativity techniques at separations large enough that low-order post-Newtonian expansions are expected to be accurate. As a case study, we evolve an equal-mass nonspinning black-hole binary from a quasicircular orbit at an initial coordinate separation of D=100M for three different resolutions. We find that the orbital period of this binary (in the numerical coordinates) is T=6422M. The orbital motion agrees with post-Newtonian predictions to within 1%. Interestingly, we find that the time derivative of the coordinate separation is dominated by a purely gauge effect leading to an apparent contraction and expansion of the orbit at twice the orbital frequency. Based on these results, we improved our evolution techniques and studied a set of black hole binaries in quasi-circular orbits starting at D=20M, D=50M, and D=100M for ~ 5, 3, and 2 orbits, respectively. We find good agreement between the numerical results and post-Newtonian predictions for the orbital frequency and radial decay rate, radiated energy and angular momentum, and waveform amplitude and phases. The results are relevant for the future computation of long-term waveforms to assist in the detection and analysis of gravitational waves by the next generation of detectors as well as the long-term simulations of black-hole binaries required to accurately model astrophysically realistic circumbinary accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accuracy Issues for Numerical Waveforms

We study the convergence properties of our implementation of the 'moving punctures' approach at very high resolutions for an equal-mass, non-spinning, black-hole binary. We find convergence of the Hamiltonian constraint on the horizons and the L2 norm of the Hamiltonian constraint in the bulk for sixth and eighth-order finite difference implementations. The momentum constraint is more sensitive, and its L2 norm shows clear convergence for a system with consistent sixth-order finite differencing, while the momentum and BSSN constraints on the horizons show convergence for both sixth and eighth-order systems. We analyze the gravitational waveform error from the late-inspiral, merger, and ringdown. We find that using several lower-order techniques for increasing the speed of numerical relativity simulations actually lead to apparently non-convergent errors. Even when using standard high-accuracy techniques, rather than seeing clean convergence, where the waveform phase is a monotonic function of grid resolution, we find that the phase tends to oscillate with resolution, possibly due to stochastic errors induced by grid refinement boundaries. Our results seem to indicate that one can obtain gravitational waveform phases to within 0.05 rad. (and possibly as small as 0.015 rad.), while the amplitude error can be reduced to 0.1%. We then compare with the waveforms obtained using the cZ4 formalism. We find that the cZ4 waveforms have larger truncation errors for a given resolution, but the Richardson extrapolation phase of the cZ4 and BSSN waveforms agree to within 0.01 rad., even during the ringdown.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Recoil From Accretion-Aligned Black-Hole Binaries

We explore the newly discovered "hangup-kick" effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Study of Conformally Flat Initial Data for Highly Spinning Black Holes and their Early Evolutions

We study conformally-flat initial data for an arbitrary number of spinning black holes with exact analytic solutions to the momentum constraints constructed from a linear combination of the classical Bowen-York and conformal Kerr extrinsic curvatures. The solution leading to the largest intrinsic spin, relative to the ADM mass of the spacetime epsilon_S=S/M^2_{ADM}, is a superposition with relative weights of Lambda=0.783 for conformal Kerr and (1-Lambda)=0.217 for Bowen-York. In addition, we measure the spin relative to the initial horizon mass M_{H_0}, and find that the quantity chi=S/M_{H_0}^2 reaches a maximum of χ^{max}=0.9856 for Lambda=0.753. After equilibration, the final black-hole spin should lie in the interval 0.9324<chi_{final}<0.9856. We perform full numerical evolutions to compute the energy radiated and the final horizon mass and spin. We find that the black hole settles to a final spin of chi_{final}^{max}=0.935 when Lambda=0.783. We also study the evolution of the apparent horizon structure of this "maximal" black hole in detail.

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

Hangup Kicks: Still Larger Recoils by Partial Spin/Orbit Alignment of Black-Hole Binaries

We revisit the scenario of the gravitational radiation recoil acquired by the final remnant of a black-hole-binary merger by studying a set of configurations that have components of the spin both aligned with the orbital angular momentum and in the orbital plane. We perform a series of 42 new full numerical simulations for equal-mass and equal-spin-magnitude binaries. We extend previous recoil fitting formulas to include nonlinear terms in the spins and successfully include both the new and known results. The new predicted maximum velocity approaches 5000km/s for spins partially aligned with the orbital angular momentum, which leads to an important increase of the probabilities of large recoils in generic astrophysical mergers. We find non-negligible probabilities for recoils of several thousand km/s from accretion-aligned binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Intermediate-mass-ratio black hole binaries II: Modeling Trajectories and Gravitational Waveforms

We revisit the scenario of small-mass-ratio (q) black-hole binaries; performing new, more accurate, simulations of mass ratios 10:1 and 100:1 for initially nonspinning black holes. We propose fitting functions for the trajectories of the two black holes as a function of time and mass ratio (in the range 1/100 < q < 1/10$) that combine aspects of post-Newtonian trajectories at smaller orbital frequencies and plunging geodesics at larger frequencies. We then use these trajectories to compute waveforms via black hole perturbation theory. Using the advanced LIGO noise curve, we see a match of ~99.5% for the leading (l,m)=(2,2) mode between the numerical relativity and perturbative waveforms. Nonleading modes have similarly high matches. We thus prove the feasibility of efficiently generating a bank of gravitational waveforms in the intermediate-mass-ratio regime using only a sparse set of full numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orbital Evolution of Extreme-Mass-Ratio Black-Hole Binaries with Numerical Relativity

We perform the first fully nonlinear numerical simulations of black-hole binaries with mass ratios 100:1. Our technique for evolving such extreme mass ratios is based on the moving puncture approach with a new gauge condition and an optimal choice of the mesh refinement (plus large computational resources). We achieve a convergent set of results for simulations starting with a small nonspinning black hole just outside the ISCO that then performs over two orbits before plunging into the 100 times more massive black hole. We compute the gravitational energy and momenta radiated as well as the final remnant parameters and compare these quantities with the corresponding perturbative estimates. The results show a close agreement. We briefly discuss the relevance of this simulations for Advanced LIGO, third-generation ground based detectors, and LISA observations, and self-force computations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intermediate Mass Ratio Black Hole Binaries: Numerical Relativity meets Perturbation Theory

We study black-hole binaries in the intermediate-mass-ratio regime 0.01 < q < 0.1 with a new technique that makes use of nonlinear numerical trajectories and efficient perturbative evolutions to compute waveforms at large radii for the leading and nonleading modes. As a proof-of-concept, we compute waveforms for q=1/10. We discuss applications of these techniques for LIGO/VIRGO data analysis and the possibility that our technique can be extended to produce accurate waveform templates from a modest number of fully-nonlinear numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intermediate-mass-ratio black hole binaries: intertwining numerical and perturbative techniques

We describe in detail full numerical and perturbative techniques to compute the gravitational radiation from intermediate-mass-ratio black-hole-binary inspirals and mergers. We perform a series of full numerical simulations of nonspinning black holes with mass ratios q=1/10 and q=1/15 from different initial separations and for different finite-difference resolutions. In order to perform those full numerical runs, we adapt the gauge of the moving punctures approach with a variable damping term for the shift. We also derive an extrapolation (to infinite radius) formula for the waveform extracted at finite radius. For the perturbative evolutions we use the full numerical tracks, transformed into the Schwarzschild gauge, in the source terms of the Regge-Wheller-Zerilli Schwarzschild perturbations formalism. We then extend this perturbative formalism to take into account small intrinsic spins of the large black hole, and validate it by computing the quasinormal mode frequencies, where we find good agreement for spins |a/M|<0.3. Including the final spins improves the overlap functions when comparing full numerical and perturbative waveforms, reaching 99.5% for the leading (l,m)=(2,2) and (3,3) modes, and 98.3% for the nonleading (2,1) mode in the q=1/10 case, which includes 8 orbits before merger. For the q=1/15 case, we obtain overlaps near 99.7% for all three modes. We discuss the modeling of the full inspiral and merger based on a combined matching of post-Newtonian, full numerical, and geodesic trajectories.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling Gravitational Recoil Using Numerical Relativity

We review the developments in modeling gravitational recoil from merging black-hole binaries and introduce a new set of 20 simulations to test our previously proposed empirical formula for the recoil. The configurations are chosen to represent generic binaries with unequal masses and precessing spins. Results of these simulations indicate that the recoil formula is accurate to within a few km/s in the similar mass-ratio regime for the out-of-plane recoil.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling maximum astrophysical gravitational recoil velocities

We measure the recoil velocity as a function of spin for equal-mass, highly-spinning black-hole binaries, with spins in the orbital plane, equal in magnitude and opposite in direction. We confirm that the leading-order effect is linear in the spin and the cosine of angle between the spin direction and the infall direction at merger. We find higher-order corrections that are proportional to the odd powers in both the spin and cosine of this angle. Taking these corrections into account, we predict that the maximum recoil will be 3680+-130 km/s.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Perturbative effects of spinning black holes with applications to recoil velocities

Recently, we proposed an enhancement of the Regge-Wheeler-Zerilli formalism for first-order perturbations about a Schwarzschild background that includes first-order corrections due to the background black-hole spin. Using this formalism, we investigate gravitational wave recoil effects from a spinning black-hole binary system analytically. This allows us to better understand the origin of the large recoils observed in full numerical simulation of spinning black hole binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Remnant Masses, Spins and Recoils from the Merger of Generic Black-Hole Binaries

We obtain empirical formulae for the final remnant black hole mass, spin, and recoil velocity from merging black-hole binaries with arbitrary mass ratios and spins. Our formulae are based on the mass ratio and spin dependence of the post-Newtonian expressions for the instantaneous radiated energy, linear momentum, and angular momentum, as well as the ISCO binding energy and angular momentum. The relative weight between the different terms is fixed by amplitude parameters chosen through a least-squares fit of recently available fully nonlinear numerical simulations. These formulae can be used for statistical studies of N-body simulations of galaxy cores and clusters, and the cosmological growth of supermassive black holes. As an example, we use these formulae to obtain a universal spin magnitude distribution of merged black holes and recoil velocity distributions for dry and hot/cold wet mergers. We also revisit the long term orbital precession and resonances and discuss how they affect spin distributions before the merging regime.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical studies of Spinning Black-Hole Binaries

We study the statistical distributions of the spins of generic black-hole binaries during the inspiral and merger, as well as the distributions of the remnant mass, spin, and recoil velocity. For the inspiral regime, we start with a random uniform distribution of spin directions S1 and S2 and magnitudes S1=S2=0.97 for different mass ratios. Starting from a fiducial initial separation of ri=50m, we perform 3.5PN evolutions down to rf=5m. At this final fiducial separation, we compute the angular distribution of the spins with respect to the final orbital angular momentum, L. We perform 16^4 simulations for six mass ratios between q=1 and q=1/16 and compute the distribution of the angles between L and Delta and L and S, directly related to recoil velocities and total angular momentum. We find a small but statistically significant bias of the distribution towards counter-alignment of both scalar products. To study the merger of black-hole binaries, we turn to full numerical techniques. We introduce empirical formulae to describe the final remnant black hole mass, spin, and recoil velocity for merging black-hole binaries with arbitrary mass ratios and spins. We then evaluate those formulae for randomly chosen directions of the individual spins and magnitudes as well as the binary's mass ratio. We found that the magnitude of the recoil velocity distribution decays as P(v) \exp(-v/2500km/s), <v>=630km/s, and sqrt{<v^2> - <v>^2}= 534km/s, leading to a 23% probability of recoils larger than 1000km/s, and a highly peaked angular distribution along the final orbital axis. The final black-hole spin magnitude show a universal distribution highly peaked at Sf/mf^2=0.73 and a 25 degrees misalignment with respect to the final orbital angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Study of multi black hole and ring singularity apparent horizons

We study critical black hole separations for the formation of a common apparent horizon in systems of $N$ - black holes in a time symmetric configuration. We study in detail the aligned equal mass cases for $N=2,3,4,5$, and relate them to the unequal mass binary black hole case. We then study the apparent horizon of the time symmetric initial geometry of a ring singularity of different radii. The apparent horizon is used as indicative of the location of the event horizon in an effort to predict a critical ring radius that would generate an event horizon of toroidal topology. We found that a good estimate for this ring critical radius is $20/(3π) M$. We briefly discuss the connection of this two cases through a discrete black hole 'necklace' configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Algebraic Classification of Numerical Spacetimes and Black-Hole-Binary Remnants

In this paper we develop a technique for determining the algebraic classification of a numerical spacetime, possibly resulting from a generic black-hole-binary merger, using the Newman-Penrose Weyl scalars. We demonstrate these techniques for a test case involving a close binary with arbitrarily oriented spins and unequal masses. We find that, post merger, the spacetime quickly approaches Petrov type II, and only approaches type D on much longer timescales. These techniques allow us to begin to explore the validity of the "no-hair theorem" for generic merging-black-hole spacetimes.

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
preprint2007arXiv

Maximum gravitational recoil

Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but anti-aligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of ~4000 km/s for maximally-rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Numerical integration of the Teukolsky Equation in the time domain

We present a fourth order convergent (2+1) numerical code to solve the Teukolsky equation in the time domain. Our approach is to rewrite the Teukolsky equation as a system of first order differential equations. In this way we get a system that has the form of an advection equation. This is used in combination with a series expansion of the solution in powers of time. To obtain a fourth order scheme we kept terms up to fourth derivative in time and use the advection-like system of differential equations to substitute the temporal derivatives by spatial derivatives. A local stability study leads to a Courant factor of 1.5 for the nonrotating case. This scheme is used to evolve gravitational perturbations in Schwarzschild and Kerr backgrounds. Our numerical method proved to be fourth order convergent in r* and theta directions. The correct power-law tail, ~1/t^{2\ell+3}, for general initial data, and ~1/t^{2\ell+4}, for time symmetric data, was found in the simulations where the duration in time of the tail phase was long enough. We verified that it is crucial to resolve accurately the angular dependence of the mode at late times in order to obtain these values of the exponents in the power-law decay. In other cases, when the decay was too fast and round-off error was reached before a tail was developed, the quasinormal modes frequencies provided a test to determine the validity of our code.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The imposition of Cauchy data to the Teukolsky equation II: Numerical comparison with the Zerilli-Moncrief approach to black hole perturbations

We revisit the question of the imposition of initial data representing astrophysical gravitational perturbations of black holes. We study their dynamics for the case of nonrotating black holes by numerically evolving the Teukolsky equation in the time domain. In order to express the Teukolsky function Psi explicitly in terms of hypersurface quantities, we relate it to the Moncrief waveform phi_M through a Chandrasekhar transformation in the case of a nonrotating black hole. This relation between Psi and phi_M holds for any constant time hypersurface and allows us to compare the computation of the evolution of Schwarzschild perturbations by the Teukolsky and by the Zerilli and Regge-Wheeler equations. We explicitly perform this comparison for the Misner initial data in the close limit approach. We evolve numerically both, the Teukolsky (with the recent code of Ref. [1]) and the Zerilli equations, finding complete agreement in resulting waveforms within numerical error. The consistency of these results further supports the correctness of the numerical code for evolving the Teukolsky equation as well as the analytic expressions for Psi in terms only of the three-metric and the extrinsic curvature.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The imposition of Cauchy data to the Teukolsky equation III: The rotating case

We solve the problem of expressing the Weyl scalars $ψ$ that describe gravitational perturbations of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express $ψ$ and $\partial _tψ$ as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out of it) of a generic spacelike slice of the spacetime. These results provide the link between initial data and $ψ$ to be evolved by the Teukolsky equation, and can be used to compute the gravitational radiation generated by two orbiting black holes in the close limit approximation. They can also be used to extract waveforms from spacetimes completely generated by numerical methods.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

An exact self-consistent gravitational shock wave in semiclassical gravity

We find a self-consistent pp-gravitational shock wave solution to the semiclassical Einstein equations resulting from the $1/N$ approach to the effective action. We model the renormalized matter stress-energy-momentum tensor by $N$ massless scalar fields in the Minkowski vacuum plus a classical particle. We show that quantum effects generate a milder singularity at the position of the particle than the classical solution, but the singularity does not disappear. At large distances from the particle, the quantum correction decreases slowly, as $1/ρ^2$ ($ρ$ being the distance to the particle in the shock wave plane). We argue that this large distance correction is a necessary consequence of quantum gravity.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Infall of a particle into a black hole as a model for gravitational radiation from the galactic center

I present here the results of the study of the gravitational radiation generated by the infall (from rest at radius $r_0$) of a point particle of mass $m_0$ into a Schwarzschild black hole of mass $M$. We use Laplace's transform methods and find that the spectra of radiation for $\sim5M<r_0<\infty$ presents a series of evenly spaced bumps. The total radiated energy is not monotonically decreasing with $r_0$, but presents a joroba (hunch-back) at around $r_0\approx4.5M$. I finally discuss the detectability of the gravitational radiation coming from the black hole in the center of our galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Exact Gravitational Shock Wave Solution of Higher Order Theories

We find an {\it exact} pp--gravitational wave solution of the fourth order gravity field equations. Outside the (delta--like) source this {\it not} a vacuum solution of General Relativity. It represents the contribution of the massive, $m=(-β)^{-1/2}$, spin--two field associated to the Ricci squared term in the gravitational Lagrangian. The fourth order terms tend to make milder the singularity of the curvature at the point where the particle is located. We generalize this analysis to $D$--dimensions, extended sources, and higher than fourth order theories. We also briefly discuss the scattering of fields by this kind of plane gravitational waves.

0 citations0 reviewsNo code linkedOpen access