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

Hiroyuki Nakano

Hiroyuki Nakano 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.COastro-ph.SRhep-thhep-ph

Catalog footprint

What is connected

38works
7topics
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

38 published item(s)

preprint2022arXiv

Adiabatic waveforms from extreme-mass-ratio inspirals: an analytical approach

Scientific analysis for the gravitational-wave detector LISA will require theoretical waveforms from extreme-mass-ratio inspirals (EMRIs) that extensively cover all possible orbital and spin configurations around astrophysical Kerr black holes. However, on-the-fly calculations of these waveforms have not yet overcome the high dimensionality of the parameter space. To confront this challenge, we present a user-ready EMRI waveform model for generic (eccentric and inclined) orbits in Kerr spacetime, using an analytical self-force approach. Our model accurately covers all EMRIs with arbitrary inclination and black hole spin, up to modest eccentricity ($\lesssim 0.3$) and separation ($\gtrsim2$--$10M$ from the last stable orbit). In that regime, our waveforms are accurate at the leading `adiabatic' order, and they approximately capture transient self-force resonances that significantly impact the gravitational-wave phase. The model fills an urgent need for extensive waveforms in ongoing data-analysis studies, and its individual components will continue to be useful in future science-adequate waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraints on Population I/II neutron star-black hole binary formation by gravitational wave and radio observations

Two neutron star (NS)-black hole (BH) binaries, GW200105 and GW200115 found in the LIGO/Virgo O3b run have smaller BH mass of 6--9\,$M_{\odot}$ which is consistent with Population I and II origin. Our population synthesis simulations using $10^6$ Population I and II binaries with appropriate initial parameters show consistent binary mass, event rate, and no detection of radio pulsar (PSR) and BH binaries in our galaxy so far. Especially, we found possible progenitors of GW200105 and GW200115 which were formed at redshift $z=0.15$ and $z=1.6$ with binary mass of $(34M_{\odot},\, 9.2M_{\odot})$ and $(23.7M_{\odot},\, 10.6M_{\odot})$, respectively. The final masses of these binaries are $(6.85M_{\odot},\,2.14M_{\odot})$ and $(6.04M_{\odot},\,1.31M_{\odot})$ which look like $(9.0_{-1.7}^{+1.7}M_{\odot},\, 1.91_{-0.24}^{+0.33}M_{\odot})$ of GW200105 and $(5.9_{-2.5}^{+2.0}M_{\odot},\,1.44_{-0.29}^{+0.85}M_{\odot})$ of GW200115, respectively. We also estimate that 2.68-19.7 PSR-BH binaries in our galaxy will be observed by SKA. The existence of NS-BHs in our galaxy can be confirmed in future SKA era. Using the GW observation of NS-BH mergers and the radio observation of PSR-BHs in future, we can get more severe constraints on the NS-BH formation process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Black hole spectroscopy for KAGRA future prospect in O5

Ringdown gravitational waves of compact binary mergers are an important target to test general relativity. The main components of the ringdown waveform after merger are black hole quasinormal modes. In general relativity, all multipolar quasinormal modes of a black hole should give the same values of black hole parameters. Although the observed binary black hole events so far are not significant enough to perform the test with ringdown gravitational waves, it is expected that the test will be achieved in third generation detectors. The Japanese gravitational wave detector KAGRA, called bKAGRA for the current configuration, has started observation, and discussions for the future upgrade plans have also started. In this study, we consider which KAGRA upgrade plan is the best to detect the subdominant quasinormal modes of black holes in the aim of testing general relativity. We use a numerical relativity waveform as injected signals that contains two multipolar modes and analyze each mode by matched filtering. Our results suggest that the plan FDSQZ, which improves the sensitivity of KAGRA for broad frequency range, is the most suitable configuration for black hole spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accuracy of the Post-Newtonian Approximation for Extreme-Mass Ratio Inspirals from Black-hole Perturbation Approach

We revisit the accuracy of the post-Newtonian (PN) approximation and its region of validity for quasi-circular orbits of a point particle in the Kerr spacetime, by using an analytically known highest post-Newtonian order gravitational energy flux and accurate numerical results in the black hole perturbation approach. It is found that regions of validity become larger for higher PN order results although there are several local maximums in regions of validity for relatively low-PN order results. This might imply that higher PN order calculations are also encouraged for comparable-mass binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational wave quasinormal mode from Population III massive black hole binaries in various models of population synthesis

Focusing on the remnant black holes after merging binary black holes, we show that ringdown gravitational waves of Population III binary black holes mergers can be detected with the rate of $5.9-500~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for various parameters and functions. This rate is estimated for the events with SNR$>8$ for the second generation gravitational wave detectors such as KAGRA. Here, ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Population III star formation rate and the fraction of binaries, respectively. When we consider only the events with SNR$>35$, the event rate becomes $0.046-4.21~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$. This suggest that for remnant black hole's spin $q_f>0.95$ we have the event rate with SNR$>35$ less than $0.037~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$, while it is $3-30~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for the third generation detectors such as Einstein Telescope. If we detect many Population III binary black holes merger, it may be possible to constrain the Population III binary evolution paths not only by the mass distribution but also by the spin distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational waves from dark matter collapse in a star

We investigate the collapse of clusters of weakly interacting massive particles (WIMPs) in the core of a Sun-like star and the possible formation of mini-black holes and the emission of gravitational waves. When the number of WIMPs is small, thermal pressure balances the WIMP cluster's self gravity. If the number of WIMPs is larger than a critical number, thermal pressure cannot balance gravity and the cluster contracts. If WIMPs are collisionless and bosonic, the cluster collapses directly to form a mini-black hole. For fermionic WIMPs, the cluster contracts until it is sustained by Fermi pressure, forming a small compact object. If the fermionic WIMP mass is smaller than $4\times 10^2$ GeV, the radius of the compact object is larger than its Schwarzschild radius and Fermi pressure temporally sustains its self gravity, halting the formation of a black hole. If the fermionic WIMP mass is larger than $4\times 10^2$ GeV, the radius is smaller than its Schwarzschild radius and the compact object becomes a mini-black hole. If the WIMP mass is 1 TeV, the size of the black hole will be approximately 2.5 cm and ultra high frequency gravitational waves will be emitted during black hole formation. The central frequency $f_c$ of ringdown gravitational waves emitted from the black hole will be approximately 2 GHz. To detect the ringdown gravitational waves, the detector's noise must be below $\sqrt{S_h(f_c)}\approx 10^{-30}/\sqrt{\rm Hz}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How close can we approach the event horizon of the Kerr black hole from the detection of the gravitational quasinormal modes?

Using the WKB method, we show that the peak location ($r_{\rm peak}$) of the potential, which determines the quasinormal mode frequency of the Kerr black hole, obeys an accurate empirical relation as a function of the specific angular momentum $a$ and the gravitational mass $M$. If the quasinormal mode with $a/M \sim 1$ is observed by gravitational wave detectors, we can confirm the black-hole space-time around the event horizon, $r_{\rm peak}=r_+ +O(\sqrt{1-q})$ where $r_+$ is the event horizon radius. While if the quasinormal mode is different from that of general relativity, we are forced to seek the true theory of gravity and/or face to the existence of the naked singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transition. Such a transition was previously explored at smaller separations. Here, we perform this transition at a separation of D=20M (large enough that the analytical metric is expected to be accurate), and evolve for six orbits. We find that small constraint violations can have large dynamical effects, but these can be removed by using a constraint-damping system like the conformal covariant formulation of the Z4 system. We find agreement between the subsequent numerical spacetime and the predictions of post-Newtonian theory for the waveform and inspiral rate that is within the post-Newtonian truncation error.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inspiralling, Non-Precessing, Spinning Black Hole Binary Spacetime via Asymptotic Matching

We construct a new global, fully analytic, approximate spacetime which accurately describes the dynamics of non-precessing, spinning black hole binaries during the inspiral phase of the relativistic merger process. This approximate solution of the vacuum Einstein's equations can be obtained by asymptotically matching perturbed Kerr solutions near the two black holes to a post-Newtonian metric valid far from the two black holes. This metric is then matched to a post-Minkowskian metric even farther out in the wave zone. The procedure of asymptotic matching is generalized to be valid on all spatial hypersurfaces, instead of a small group of initial hypersurfaces discussed in previous works. This metric is well suited for long term dynamical simulations of spinning black hole binary spacetimes prior to merger, such as studies of circumbinary gas accretion which requires hundreds of binary orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Possible confirmation of the existence of ergoregion by the Kerr quasinormal mode in gravitational waves from Pop III massive black hole binary

The existence of the ergoregion of the Kerr space-time has not been confirmed observationally yet. We show that the confirmation would be possible by observing the quasinormal mode in gravitational waves. As an example, using the recent population synthesis results of Pop III binary black holes, we find that the peak of the final merger mass ($M_f$) is about $50~\rm M_{\odot}$, while the fraction of the final spin $q_f = a_f/M_f > 0.7$ needed for the confirmation of a part of ergoregion is $\sim 77\%$. To confirm the frequency of the quasinormal mode, ${\rm SNR} > 35$ is needed. The standard model of Pop III population synthesis tells us that the event rate for the confirmation of more than $50\%$ of the ergoregion by the second generation gravitational wave detectors is $\sim 2.3$ ${\rm events\ yr^{-1}\ (SFR_p/(10^{-2.5}\ M_\odot yr^{-1}\ Mpc^{-3}))} \cdot (\rm [f_b/(1+f_b)]/0.33)$ where ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Pop III star formation rate and the fraction of binaries, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pre-DECIGO can get the smoking gun to decide the astrophysical or cosmological origin of GW150914-like binary black holes

Pre-DECIGO consists of three spacecraft arranged in an equilateral triangle with 100km arm lengths orbiting 2000km above the surface of the earth. It is hoped that the launch date will be in the late 2020s. Pre-DECIGO has one clear target: binary black holes (BBHs) like GW150914 and GW151226. Pre-DECIGO can detect $\sim 30M_\odot-30M_\odot$ BBH mergers up to redshift $z\sim 30$. The cumulative event rate is $\sim 1.8\times 10^{5}\,{\rm events~yr^{-1}}$ in the Pop III origin model of BBHs like GW150914, and it saturates at $z\sim 10$, while in the primordial BBH (PBBH) model, the cumulative event rate is $ \sim 3\times 10^{4}\,{\rm events~ yr^{-1}}$ at $z=30$ even if only $0.1\%$ of the dark matter consists of PBHs, and it is still increasing at $z=30$. In the Pop I/II model of BBHs, the cumulative event rate is $(3-10)\times10^{5}\,{\rm events~ yr^{-1}}$ and it saturates at $z \sim 6$. We present the requirements on orbit accuracy, drag free techniques, laser power, frequency stability, and interferometer test mass. For BBHs like GW150914 at 1Gpc, SNR$\sim 90$ is achieved with the definition of Pre-DECIGO in the $0.01-100$Hz band. Pre-DECIGO can measure the mass spectrum and the $z$-dependence of the merger rate to distinguish various models of BBHs like GW150914. Pre-DECIGO can also predict the direction of BBHs at $z=0.1$ with an accuracy of $\sim 0.3\,{\rm deg}^2$ and a merging time accuracy of $\sim 1$s at about a day before the merger so that ground-based GW detectors further developed at that time as well as electromagnetic follow-up observations can prepare for the detection of merger in advance. For intermediate mass BBHs at a large redshift $z > 10$, the QNM frequency after the merger can be within the Pre-DECIGO band so that the ringing tail can also be detectable to confirm the Einstein theory of general relativity with SNR$\sim 35$. [abridged]

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

Spinning, Precessing, Black Hole Binary Spacetime via Asymptotic Matching

We briefly discuss a method to construct a global, analytic, approximate spacetime for precessing, spinning binary black holes. The spacetime construction is broken into three parts: the inner zones are the spacetimes close to each black hole, and are approximated by perturbed Kerr solutions; the near zone is far from the two black holes, and described by the post-Newtonian metric; and finally the wave (far) zone, where retardation effects need to be taken into account, is well modeled by the post-Minkowskian metric. These individual spacetimes are then stitched together using asymptotic matching techniques to obtain a global solution that approximately satisfies the Einstein field equations. Precession effects are introduced by rotating the black hole spin direction according to the precessing equations of motion, in a way that is consistent with the global spacetime construction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The detection of Quasinormal Mode of Pop III Binary Black Hole merger with final $M \sim 60M_\odot$ and $a/M \sim 0.7$ would confirm the strong gravity space-time around $\sim 2M$ which is only $\sim 1.17$ times the event horizon radius

Recent population synthesis simulations of Pop III stars suggest that the event rate of coalescence of $\sim 30M_\odot$--$30M_\odot$ binary black holes can be high enough for the detection by the second generation gravitational wave detectors. The frequencies of chirp signal as well as quasinormal modes are near the best sensitivity of these detectors so that it would be possible to confirm Einstein's general relativity. Using the WKB method, we suggest that for the typical value of spin parameter $a/M\sim 0.7$ from numerical relativity results of the coalescence of binary black holes, the strong gravity of the black hole space-time at around the radius $2M$, which is just $\sim 1.17$ times the event horizon radius, would be confirmed as predicted by general relativity. The expected event rate with the signal-to-noise ratio $> 35$ needed for the determination of the quasinormal mode frequency with the meaningful accuracy is $0.17$--$7.2$~${\rm events~yr^{-1}~(SFR_p/(10^{-2.5}~M_\odot~yr^{-1}~Mpc^{-3}))} \cdot (\rm [f_b/(1+f_b)]/0.33)$ where ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Pop III star formation rate and the fraction of binaries, respectively. As for the possible optical counter part, if the merged black hole of mass $M\sim 60M_\odot$ is in the interstellar matter with $n\sim 100~{\rm cm^{-3}}$ and the proper motion of black hole is $\sim 1~{\rm km~s^{-1}}$, the luminosity is $\sim 10^{40}~{\rm erg~ s^{-1}}$ which can be detected up to $\sim 300~{\rm Mpc}$, for example, by Subaru-HSC and LSST with the limiting magnitude 26.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The detection of quasinormal mode with $a/M \sim 0.95$ would prove a sphere $99\%$ soaking in the ergoregion of the Kerr space-time

Recent numerical relativity simulations of mergers of binary black holes suggest that the maximum final value of $a/M$ is $\sim 0.95$ for the coalescence of two equal mass black holes with aligned spins of the same magnitude $a/M=0.994$ which is close to the upper limit $a/M=0.998$ of accretion spin-up shown by Thorne. Using the WKB method, we suggest that if quasinormal modes with $a/M \sim 0.95$ are detected by the second generation gravitational wave detectors, we could confirm the strong gravity space-time based on Einstein's general relativity up to $1.33M$ which is only $\sim 1.014$ times the event horizon radius and within the ergoregion. One more message about black hole geometry is expected here. If the quasinormal mode is different from that of general relativity, we need to find the true theory of gravity which deviates from general relativity only near the black hole horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A note on gravitational wave extraction from binary simulations

In previous works, we developed a perturbative extrapolation formula to obtain gravitational waves at infinity from binary black hole simulations, based on the Regge-Wheeler-Zerilli formalism. In practice, this formula was basically derived in a background flat spacetime. We derive here an estimation of errors and an improved extrapolation formula by using the Teukolsky perturbative equation for a background Kerr spacetime. The improved formula requires knowledge about a mass and spin of the background spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole-neutron star binary merger: Dependence on black hole spin orientation and equations of state

We systematically performed numerical-relativity simulations for black hole (BH) - neutron star (NS) binary mergers with a variety of the BH spin orientation and equations of state (EOS) of the NS. The initial misalignment angles of the BH spin are chosen in the range of i_tilt,0 = 30--90[deg.]. We employed four models of NS EOS with which the compactness of the NS is in the range of C = M_NS/R_NS = 0.138--0.180, where M_NS and R_NS are the mass and the radius of the NS, respectively. The mass ratio of the BH to the NS, Q = M_BH/M_NS, and the dimensionless spin parameter of the BH, chi, are chosen to be Q = 5 and chi = 0.75, together with M_NS = 1.35 M_sun. We obtain the following results: (i) The inclination angle of i_tilt,0 < 70[deg.] and i_tilt,0 < 50[deg.] are required for the formation of a remnant disk with its mass larger than 0.1 M_sun for the case C = 0.140 and C = 0.160, respectively, while the disk mass is always smaller than 0.1M_sun for C = 0.175. The ejecta with its mass larger than 0.01 M_sun is obtained for i_tilt,0 < 85[deg.] with C = 0.140, for i_tilt,0 < 65[deg.] with C = 0.160, and for i_tilt,0 < 30[deg.] with C = 0.175. (ii) The rotational axis of the dense part of the remnant disk is approximately aligned with the remnant BH spin for i_tilt,0 = 30[deg.]. On the other hand, the disk axis is misaligned initially with ~ 30[deg.] for i_tilt,0 = 60[deg.], and the alignment with the remnant BH spin is achieved at ~ 50--60 ms after the onset of merger. The accretion time scale of the remnant disk is typically ~ 100 ms and depends only weakly on the misalignment angle and the EOS. (iii) The ejecta velocity is typically ~ 0.2--0.3c and depends only weakly on i_tilt,0 and the EOS of the NS, while the morphology of the ejecta depends on its mass. (iv) The gravitational-wave spectra contains the information of the NS compactness in the cutoff frequency for i_tilt,0 < 60[deg.].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Possible golden events for ringdown gravitational waves

There is a forbidden region in the parameter space of quasinormal modes of black holes in general relativity. Using both inspiral and ringdown phases of gravitational waves from binary black holes, we propose two methods to test general relativity. We also evaluate how our methods will work when we apply them to Pop III black-hole binaries with typical masses. Adopting the simple mean of the estimated range of the event rate, we have the expected rate of 500 ${\rm yr^{-1}}$. Then, the rates of events with signal-to-noise ratios greater than 20 and greater than 50 are 32 ${\rm yr^{-1}}$ and 2 ${\rm yr^{-1}}$, respectively. Therefore, there is a good chance to confirm (or refute) the Einstein theory in the strong gravity region by observing the expected quasinormal modes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Approximate black hole binary spacetime via asymptotic matching

We construct a fully analytic, general relativistic, nonspinning black hole binary spacetime that approximately solves the vacuum Einstein equations everywhere in space and time for black holes sufficiently well separated. The metric is constructed by asymptotically matching perturbed Schwarzschild metrics near each black hole to a two-body post-Newtonian metric far from them, and a two-body post-Minkowskian metric farther still. Asymptotic matching is done without linearizing about a particular time slice, and thus it is valid dynamically and for all times, provided the binary is sufficiently well separated. This approximate global metric can be used for long dynamical evolutions of relativistic magnetohydrodynamical, circumbinary disks around inspiraling supermassive black holes to study a variety of phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational Self-Force Correction to the Innermost Stable Circular Equatorial Orbit of a Kerr Black Hole

For a self-gravitating particle of mass μin orbit around a Kerr black hole of mass M >> μ, we compute the O(μ/M) shift in the frequency of the innermost stable circular equatorial orbit (ISCEO) due to the conservative piece of the gravitational self-force acting on the particle. Our treatment is based on a Hamiltonian formulation of the dynamics in terms of geodesic motion in a certain locally-defined effective smooth spacetime. We recover the same result using the so-called first law of binary black-hole mechanics. We give numerical results for the ISCEO frequency shift as a function of the black hole's spin amplitude, and compare with predictions based on the post-Newtonian approximation and the effective one-body model. Our results provide an accurate strong-field benchmark for spin effects in the general relativistic two-body problem.

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

Evolution of the Carter constant for a resonant inspiral into a Kerr black hole: I. The scalar case

We discuss the inspiral of a small body around a Kerr black hole. When the time scale of the radiation reaction is sufficiently longer than its orbital period, the leading order orbital evolution is described only by the knowledge of the averaged evolution of the constants of motion, i.e., the energy, azimuthal angular momentum and the Carter constant. Although there is no conserved current composed of the perturbation field corresponding to the Carter constant, it has been shown that the averaged rate of change of the Carter constant can be given by a simple formula, when there exists a simultaneous turning point of the radial and polar oscillations. However, an inspiralling orbit may cross a "resonance" point, where the frequencies of the radial and polar orbital oscillations are in a rational ratio. At the resonance point, one cannot find a simultaneous turning point in general. Hence, even for the averaged rate of change of the Carter constant, a direct computation of the 'self-force", which is still challenging especially in the case of the Kerr background, seems to be necessary. In this paper, we present a method of computing the averaged rate of change of the Carter constant in a relatively simple manner at the "resonance" point.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotically Matched Spacetime Metric for Non-Precessing, Spinning Black Hole Binaries

We construct a closed-form, fully analytical 4-metric that approximately represents the spacetime evolution of non-precessing, spinning black hole binaries from infinite separations up to a few orbits prior to merger. We employ the technique of asymptotic matching to join a perturbed Kerr metric in the neighborhood of each spinning black hole to a near-zone, post-Newtonian metric farther out. The latter is already naturally matched to a far-zone, post-Minkowskian metric that accounts for full temporal retardation. The result is a 4-metric that is approximately valid everywhere in space and in a small bundle of spatial hypersurfaces. We here restrict our attention to quasi- circular orbits, but the method is valid for any orbital motion or physical scenario, provided an overlapping region of validity or buffer zone exists. A simple extension of such a metric will allow for future studies of the accretion disk and jet dynamics around spinning back hole binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Circumbinary MHD Accretion into Inspiraling Binary Black Holes

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

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

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
preprint2010arXiv

Hybrid black-hole binary initial data

Traditional black-hole binary puncture initial data is conformally flat. This unphysical assumption is coupled with a lack of radiation signature from the binary's past life. As a result, waveforms extracted from evolutions of this data display an abrupt jump. In Kelly et al. [Class.Quant.Grav.27:114005,2010], a new binary black-hole initial data with radiation contents derived in the post-Newtonian (PN) calculation was adapted to puncture evolutions in numerical relativity. This data satisfies the constraint equations to the 2.5PN order, and contains a transverse-traceless "wavy" metric contribution, violating the standard assumption of conformal flatness. Although the evolution contained less spurious radiation, there were undesired features; the unphysical horizon mass loss and the large initial orbital eccentricity. Introducing a hybrid approach to the initial data evaluation, we significantly reduce these undesired features.

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

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
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

Adiabatic Evolution of three 'Constants' of Motion for Greatly Inclined Orbits in Kerr spacetime

General orbits of a particle of small mass $μ$ around a Kerr black hole of mass $M$ are characterized by three parameters: the energy, the angular momentum and the Carter constant. The time-averaged rates of change of the energy and the angular momentum can be obtained by computing the corresponding fluxes of gravitational waves emitted by the particle. By contrast, the time-averaged rate of change of the Carter constant cannot be expressed as a flux of gravitational waves. Recently a method to compute this rate of change was proposed by Mino, and we refined it into a simplified form. In this paper we further extend our previous work to give a new formulation without the aid of expansion in terms of a small inclination angle.

0 citations0 reviewsNo code linkedOpen access