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

Bela Szilagyi

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COastro-phastro-ph.SRComputational Engineering, Finance, and Sciencephysics.comp-phphysics.data-anphysics.geo-ph

Catalog footprint

What is connected

39works
9topics
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

39 published item(s)

preprint2020arXiv

Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GipsyX/RTGx, A New Tool Set for Space Geodetic Operations and Research

GipsyX/RTGx is the Jet Propulsion Laboratory's (JPL) next generation software package for positioning, navigation, timing, and Earth science using measurements from three geodetic techniques: Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS); with Very Long Baseline Interferometry (VLBI) under development. The software facilitates combined estimation of geodetic and geophysical parameters using a Kalman filter approach on real or simulated data in both post-processing and in real-time. The estimated parameters include station coordinates and velocities, satellite orbits and clocks, Earth orientation, ionospheric and tropospheric delays. The software is also capable of full realization of a dynamic terrestrial reference through analysis and combination of time series of ground station coordinates. We present some key aspects of its new architecture, and describe some of its major applications, including Real-time orbit determination and ephemeris predictions in the U.S. Air Force Next Generation GPS Operational Control Segment (OCX), as well as in JPL's Global Differential GPS (GDGPS) System, supporting User Range Error (URE) of $<$ 5 cm RMS; precision post-processing GNSS orbit determination, including JPL's contributions to the International GNSS Service (IGS) with URE in the 2 cm RMS range; Precise point positioning (PPP) with ambiguity resolution, both statically and kinematically, for geodetic applications with 2 mm horizontal, and 6.5 mm vertical repeatability for static positioning; Operational orbit and clock determination for Low Earth Orbiting (LEO) satellites, such as NASA's Gravity Recovery and Climate Experiment (GRACE) mission with GRACE relative clock alignment at the 20 ps level.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accuracy of binary black hole waveform models for aligned-spin binaries

Coalescing binary black holes are among the primary science targets for second generation ground-based gravitational wave (GW) detectors. Reliable GW models are central to detection of such systems and subsequent parameter estimation. This paper performs a comprehensive analysis of the accuracy of recent waveform models for binary black holes with aligned spins, utilizing a new set of $84$ high-accuracy numerical relativity simulations. Our analysis covers comparable mass binaries ($1\le m_1/m_2\le 3$), and samples independently both black hole spins up to dimensionless spin-magnitude of $0.9$ for equal-mass binaries and $0.85$ for unequal mass binaries. Furthermore, we focus on the high-mass regime (total mass $\gtrsim 50M_\odot$). The two most recent waveform models considered (PhenomD and SEOBNRv2) both perform very well for signal detection, losing less than 0.5\% of the recoverable signal-to-noise ratio $ρ$, except that SEOBNRv2's efficiency drops mildly for both black hole spins aligned with large magnitude. For parameter estimation, modeling inaccuracies of SEOBNRv2 are found to be smaller than systematic uncertainties for moderately strong GW events up to roughly $ρ\lesssim 15$. PhenomD's modeling errors are found to be smaller than SEOBNRv2's, and are generally irrelevant for $ρ\lesssim 20$. Both models' accuracy deteriorates with increased mass-ratio, and when at least one black hole spin is large and aligned. The SEOBNRv2 model shows a pronounced disagreement with the numerical relativity simulation in the merger phase, for unequal masses and simultaneously both black hole spins very large and aligned. Two older waveform models (PhenomC and SEOBNRv1) are found to be distinctly less accurate than the more recent PhenomD and SEOBNRv2 models. Finally, we quantify the bias expected from all GW models during parameter estimation for recovery of binary's masses and spins.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Black Hole-Neutron Star Mergers with a Hot Nuclear Equation of State: Outflow and Neutrino-Cooled Disk for a Low-Mass, High-Spin Case

Neutrino emission significantly affects the evolution of the accretion tori formed in black hole-neutron star mergers. It removes energy from the disk, alters its composition, and provides a potential power source for a gamma-ray burst. To study these effects, simulations in general relativity with a hot microphysical equation of state and neutrino feedback are needed. We present the first such simulation, using a neutrino leakage scheme for cooling to capture the most essential effects and considering a moderate mass (1.4 M_{\odot} neutron star, 5.6 M_{\odot} black hole), high spin (black hole J/M^2=0.9) system with the K_0=220 MeV Lattimer-Swesty equation of state. We find that about 0.08 M_{\odot} of nuclear matter is ejected from the system, while another 0.3 M_{\odot} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (i) to make the disk much denser and more compact, (ii) to cause the average electron fraction Y_e of the disk to rise to about 0.2 and then gradually decrease again, and (iii) to gradually cool the disk. The disk is initially hot (T~6 MeV) and luminous in neutrinos (L_ν~10^{54} erg s^{-1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effects of neutron-star dynamic tides on gravitational waveforms within the effective-one-body approach

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Initial data for black hole-neutron star binaries, with rotating stars

The coalescence of a neutron star with a black hole is a primary science target of ground-based gravitational wave detectors. Constraining or measuring the neutron star spin directly from gravitational wave observations requires knowledge of the dependence of the emission properties of these systems on the neutron star spin. This paper lays foundations for this task, by developing a numerical method to construct initial data for black hole--neutron star binaries with arbitrary spin on the neutron star. We demonstrate the robustness of the code by constructing initial-data sets in large regions of the parameter space. In addition to varying the neutron star spin-magnitude and spin-direction, we also explore neutron star compactness, mass-ratio, black hole spin, and black hole spin-direction. Specifically, we are able to construct initial data sets with neutron stars spinning near centrifugal break-up, and with black hole spins as large as $S_{\rm BH}/M_{\rm BH}^2=0.99$.

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

On the accuracy and precision of numerical waveforms: Effect of waveform extraction methodology

We present a new set of 95 numerical relativity simulations of non-precessing binary black holes (BBHs). The simulations sample comprehensively both black-hole spins up to spin magnitude of 0.9, and cover mass ratios 1 to 3. The simulations cover on average 24 inspiral orbits, plus merger and ringdown, with low initial orbital eccentricities $e<10^{-4}$. A subset of the simulations extends the coverage of non-spinning BBHs up to mass ratio $q=10$. Gravitational waveforms at asymptotic infinity are computed with two independent techniques, extrapolation, and Cauchy characteristic extraction. An error analysis based on noise-weighted inner products is performed. We find that numerical truncation error, error due to gravitational wave extraction, and errors due to the finite length of the numerical waveforms are of similar magnitude, with gravitational wave extraction errors somewhat dominating at noise-weighted mismatches of $\sim 3\times 10^{-4}$. This set of waveforms will serve to validate and improve aligned-spin waveform models for gravitational wave science.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulations of inspiraling and merging double neutron stars using the Spectral Einstein Code

We present results on the inspiral, merger, and post-merger evolution of a neutron star - neutron star (NSNS) system. Our results are obtained using the hybrid pseudospectral-finite volume Spectral Einstein Code (SpEC). To test our numerical methods, we evolve an equal-mass system for $\approx 22$ orbits before merger. This waveform is the longest waveform obtained from fully general-relativistic simulations for NSNSs to date. Such long (and accurate) numerical waveforms are required to further improve semi-analytical models used in gravitational wave data analysis, for example the effective one body models. We discuss in detail the improvements to SpEC's ability to simulate NSNS mergers, in particular mesh refined grids to better resolve the merger and post-merger phases. We provide a set of consistency checks and compare our results to NSNS merger simulations with the independent BAM code. We find agreement between them, which increases confidence in results obtained with either code. This work paves the way for future studies using long waveforms and more complex microphysical descriptions of neutron star matter in SpEC.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast and accurate prediction of numerical relativity waveforms from binary black hole coalescences using surrogate models

Simulating a binary black hole (BBH) coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. Using reduced order modeling techniques, we construct an accurate surrogate model, which is evaluated in a millisecond to a second, for numerical relativity (NR) waveforms from non-spinning BBH coalescences with mass ratios in $[1, 10]$ and durations corresponding to about $15$ orbits before merger. We assess the model's uncertainty and show that our modeling strategy predicts NR waveforms {\em not} used for the surrogate's training with errors nearly as small as the numerical error of the NR code. Our model includes all spherical-harmonic ${}_{-2}Y_{\ell m}$ waveform modes resolved by the NR code up to $\ell=8.$ We compare our surrogate model to Effective One Body waveforms from $50$-$300 M_\odot$ for advanced LIGO detectors and find that the surrogate is always more faithful (by at least an order of magnitude in most cases).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nearly extremal apparent horizons in simulations of merging black holes

The spin angular momentum $S$ of an isolated Kerr black hole is bounded by the surface area $A$ of its apparent horizon: $8πS \le A$, with equality for extremal black holes. In this paper, we explore the extremality of individual and common apparent horizons for merging, rapidly spinning binary black holes. We consider simulations of merging black holes with equal masses $M$ and initial spin angular momenta aligned with the orbital angular momentum, including new simulations with spin magnitudes up to $S/M^2 = 0.994$. We measure the area and (using approximate Killing vectors) the spin on the individual and common apparent horizons, finding that the inequality $8πS < A$ is satisfied in all cases but is very close to equality on the common apparent horizon at the instant it first appears. We also introduce a gauge-invariant lower bound on the extremality by computing the smallest value that Booth and Fairhurst's extremality parameter can take for any scaling. Using this lower bound, we conclude that the common horizons are at least moderately close to extremal just after they appear. Finally, following Lovelace et al. (2008), we construct quasiequilibrium binary-black-hole initial data with "overspun" marginally trapped surfaces with $8πS > A$ and for which our lower bound on their Booth-Fairhurst extremality exceeds unity. These superextremal surfaces are always surrounded by marginally outer trapped surfaces (i.e., by apparent horizons) with $8πS<A$. The extremality lower bound on the enclosing apparent horizon is always less than unity but can exceed the value for an extremal Kerr black hole. (Abstract abbreviated.)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical relativity reaching into post-Newtonian territory: a compact-object binary simulation spanning 350 gravitational-wave cycles

We present the first numerical-relativity simulation of a compact-object binary whose gravitational waveform is long enough to cover the entire frequency band of advanced gravitational-wave detectors, such as LIGO, Virgo and KAGRA, for mass ratio 7 and total mass as low as $45.5\,M_\odot$. We find that effective-one-body models, either uncalibrated or calibrated against substantially shorter numerical-relativity waveforms at smaller mass ratios, reproduce our new waveform remarkably well, with a negligible loss in detection rate due to modeling error. In contrast, post-Newtonian inspiral waveforms and existing calibrated phenomenological inspiral-merger-ringdown waveforms display greater disagreement with our new simulation. The disagreement varies substantially depending on the specific post-Newtonian approximant used.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Post-merger evolution of a neutron star-black hole binary with neutrino transport

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A sparse representation of gravitational waves from precessing compact binaries

Many relevant applications in gravitational wave physics share a significant common problem: the seven-dimensional parameter space of gravitational waveforms from precessing compact binary inspirals and coalescences is large enough to prohibit covering the space of waveforms with sufficient density. We find that by using the reduced basis method together with a parametrization of waveforms based on their phase and precession, we can construct ultra-compact yet high-accuracy representations of this large space. As a demonstration, we show that less than $100$ judiciously chosen precessing inspiral waveforms are needed for $200$ cycles, mass ratios from $1$ to $10$ and spin magnitudes $\le 0.9$. In fact, using only the first $10$ reduced basis waveforms yields a maximum mismatch of $0.016$ over the whole range of considered parameters. We test whether the parameters selected from the inspiral regime result in an accurate reduced basis when including merger and ringdown; we find that this is indeed the case in the context of a non-precessing effective-one-body model. This evidence suggests that as few as $\sim 100$ numerical simulations of binary black hole coalescences may accurately represent the seven-dimensional parameter space of precession waveforms for the considered ranges.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved methods for simulating nearly extremal binary black holes

Astrophysical black holes could be nearly extremal (that is, rotating nearly as fast as possible); therefore, nearly extremal black holes could be among the binaries that current and future gravitational-wave observatories will detect. Predicting the gravitational waves emitted by merging black holes requires numerical-relativity simulations, but these simulations are especially challenging when one or both holes have mass $m$ and spin $S$ exceeding the Bowen-York limit of $S/m^2=0.93$. We present improved methods that enable us to simulate merging, nearly extremal black holes more robustly and more efficiently. We use these methods to simulate an unequal-mass, precessing binary black hole coalescence, where the larger black hole has $S/m^2=0.99$. We also use these methods to simulate a non-precessing binary black hole coalescence, where both black holes have $S/m^2=0.994$, nearly reaching the Novikov-Thorne upper bound for holes spun up by thin accretion disks. We demonstrate numerical convergence and estimate the numerical errors of the waveforms; we compare numerical waveforms from our simulations with post-Newtonian and effective-one-body waveforms; we compare the evolution of the black-hole masses and spins with analytic predictions; and we explore the effect of increasing spin magnitude on the orbital dynamics (the so-called "orbital hangup" effect).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Key Elements of Robustness in Binary Black Hole Evolutions using Spectral Methods

As a network of advanced-era gravitational wave detectors is nearing its design sensitivity, efficient and accurate waveform modeling becomes more and more relevant. Understanding of the nature of the signal being sought can have an order unity effect on the event rates seen in these instruments. The paper provides a description of key elements of the Spectral Einstein Code ({\tt SpEC}), with details of our spectral adaptive mesh refinement (AMR) algorithm that has been optimized for binary black hole (BBH) evolutions. We expect that the gravitational waveform catalog produced by our code will have a central importance in both the detection and parameter estimation of gravitational waves in these instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron star-black hole mergers with a nuclear equation of state and neutrino cooling: Dependence in the binary parameters

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($χ_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_ν\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solving Einstein's Equation Numerically on Manifolds With Arbitrary Spatial Topologies

This paper develops a method for solving Einstein's equation numerically on multi-cube representations of manifolds with arbitrary spatial topologies. This method is designed to provide a set of flexible, easy to use computational procedures that make it possible to explore the never before studied properties of solutions to Einstein's equation on manifolds with arbitrary toplogical structures. A new covariant, first-order symmetric-hyperbolic representation of Einstein's equation is developed for this purpose, along with the needed boundary conditions at the interfaces between adjoining cubic regions. Numerical tests are presented that demonstrate the long-term numerical stability of this method for evolutions of a complicated, time-dependent solution of Einstein's equation coupled to a complex scalar field on a manifold with spatial topology S^3. The accuracy of these numerical test solutions is evaluated by performing convergence studies and by comparing the full non-linear numerical results to the analytical perturbation solutions, which are also derived here.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A catalog of 174 binary black-hole simulations for gravitational-wave astronomy

This paper presents a publicly available catalog of 174 numerical binary black-hole simulations following up to 35 orbits. The catalog includes 91 precessing binaries, mass ratios up to 8:1, orbital eccentricities from a few percent to $10^{-5}$, black-hole spins up to 98% of the theoretical maximum, and radiated energies up to 11.1% of the initial mass. We establish remarkably good agreement with post-Newtonian precession of orbital and spin directions for two new precessing simulations, and we discuss other applications of this catalog. Formidable challenges remain: e.g., precession complicates the connection of numerical and approximate analytical waveforms, and vast regions of the parameter space remain unexplored.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black hole-neutron star mergers at realistic mass ratios: Equation of state and spin orientation effects

Black hole-neutron star mergers resulting in the disruption of the neutron star and the formation of an accretion disk and/or the ejection of unbound material are prime candidates for the joint detection of gravitational-wave and electromagnetic signals when the next generation of gravitational-wave detectors comes online. However, the disruption of the neutron star and the properties of the post-merger remnant are very sensitive to the parameters of the binary. In this paper, we study the impact of the radius of the neutron star and the alignment of the black hole spin for systems within the range of mass ratio currently deemed most likely for field binaries (M_BH ~ 7 M_NS) and for black hole spins large enough for the neutron star to disrupt (J/M^2=0.9). We find that: (i) In this regime, the merger is particularly sensitive to the radius of the neutron star, with remnant masses varying from 0.3M_NS to 0.1M_NS for changes of only 2 km in the NS radius; (ii) 0.01-0.05M_sun of unbound material can be ejected with kinetic energy >10^51 ergs, a significant increase compared to low mass ratio, low spin binaries. This ejecta could power detectable optical and radio afterglows. (iii) Only a small fraction (<3%) of the Advanced LIGO events in this parameter range have gravitational-wave signals which could offer constraints on the equation of state of the neutron star. (iv) A misaligned black hole spin works against disk formation, with less neutron star material remaining outside of the black hole after merger, and a larger fraction of that material remaining in the tidal tail instead of the forming accretion disk. (v) Large kicks (v>300 km/s) can be given to the final black hole as a result of a precessing BHNS merger, when the disruption of the neutron star occurs just outside or within the innermost stable spherical orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparing Gravitational Waveform Extrapolation to Cauchy-Characteristic Extraction in Binary Black Hole Simulations

We extract gravitational waveforms from numerical simulations of black hole binaries computed using the Spectral Einstein Code. We compare two extraction methods: direct construction of the Newman-Penrose (NP) scalar $Ψ_4$ at a finite distance from the source and Cauchy-characteristic extraction (CCE). The direct NP approach is simpler than CCE, but NP waveforms can be contaminated by near-zone effects---unless the waves are extracted at several distances from the source and extrapolated to infinity. Even then, the resulting waveforms can in principle be contaminated by gauge effects. In contrast, CCE directly provides, by construction, gauge-invariant waveforms at future null infinity. We verify the gauge invariance of CCE by running the same physical simulation using two different gauge conditions. We find that these two gauge conditions produce the same CCE waveforms but show differences in extrapolated-$Ψ_4$ waveforms. We examine data from several different binary configurations and measure the dominant sources of error in the extrapolated-$Ψ_4$ and CCE waveforms. In some cases, we find that NP waveforms extrapolated to infinity agree with the corresponding CCE waveforms to within the estimated error bars. However, we find that in other cases extrapolated and CCE waveforms disagree, most notably for $m=0$ "memory" modes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective-one-body model for black-hole binaries with generic mass ratios and spins

Gravitational waves emitted by black-hole binary systems have the highest signal-to-noise ratio in LIGO and Virgo detectors when black-hole spins are aligned with the orbital angular momentum and extremal. For such systems, we extend the effective-one-body inspiral-merger-ringdown waveforms to generic mass ratios and spins calibrating them to 38 numerical-relativity nonprecessing waveforms produced by the SXS Collaboration. The numerical-relativity simulations span mass ratios from 1 to 8, spin magnitudes up to 98% of extremality, and last for 40 to 60 gravitational-wave cycles. When the total mass of the binary is between 20Msun and 200Msun, the effective-one-body nonprecessing (dominant mode) waveforms have overlaps above 99% (using the advanced-LIGO design noise spectral density) with all of the 38 nonprecessing numerical waveforms, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling. Moreover, without further calibration, we show that the precessing effective-one-body (dominant mode) waveforms have overlaps above 97% with two very long, strongly precessing numerical-relativity waveforms, when maximizing only on the initial phase and time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First direct comparison of non-disrupting neutron star-black hole and binary black hole merger simulations

We present the first direct comparison of numerical simulations of neutron star-black hole and black hole-black hole mergers in full general relativity. We focus on a configuration with non spinning objects and within the most likely range of mass ratio for neutron star-black hole systems (q=6). In this region of the parameter space, the neutron star is not tidally disrupted prior to merger, and we show that the two types of mergers appear remarkably similar. The effect of the presence of a neutron star on the gravitational wave signal is not only undetectable by the next generation of gravitational wave detectors, but also too small to be measured in the numerical simulations: even the plunge, merger and ringdown signals appear in perfect agreement for both types of binaries. The characteristics of the post-merger remnants are equally similar, with the masses of the final black holes agreeing within dM< 5 10^{-4}M_BH and their spins within da< 10^{-3}M_BH. The rate of periastron advance in the mixed binary agrees with previously published binary black hole results, and we use the inspiral waveforms to place constraints on the accuracy of our numerical simulations independent of algorithmic choices made for each type of binary. Overall, our results indicate that non-disrupting neutron star-black hole mergers are exceptionally well modeled by black hole-black hole mergers, and that given the absence of mass ejection, accretion disk formation, or differences in the gravitational wave signals, only electromagnetic precursors could prove the presence of a neutron star in low-spin systems of total mass ~10Msun, at least until the advent of gravitational wave detectors with a sensitivity comparable to that of the proposed Einstein Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inspiral-merger-ringdown waveforms of spinning, precessing black-hole binaries in the effective-one-body formalism

We describe a general procedure to generate spinning, precessing waveforms that include inspiral, merger and ringdown stages in the effective-one-body (EOB) approach. The procedure uses a precessing frame in which precession-induced amplitude and phase modulations are minimized, and an inertial frame, aligned with the spin of the final black hole, in which we carry out the matching of the inspiral-plunge to merger-ringdown waveforms. As a first application, we build spinning, precessing EOB waveforms for the gravitational modes l=2 such that in the nonprecessing limit those waveforms agree with the EOB waveforms recently calibrated to numerical-relativity waveforms. Without recalibrating the EOB model, we then compare EOB and post-Newtonian precessing waveforms to two numerical-relativity waveforms produced by the Caltech-Cornell-CITA collaboration. The numerical waveforms are strongly precessing and have 35 and 65 gravitational-wave cycles. We find a remarkable agreement between EOB and numerical-relativity precessing waveforms and spins' evolutions. The phase difference is ~ 0.2 rad at merger, while the mismatches, computed using the advanced-LIGO noise spectral density, are below 2% when maximizing only on the time and phase at coalescence and on the polarization angle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Massive disk formation in the tidal disruption of a neutron star by a nearly extremal black hole

Black hole-neutron star (BHNS) binaries are important sources of gravitational waves for second-generation interferometers, and BHNS mergers are also a proposed engine for short, hard gamma-ray bursts. The behavior of both the spacetime (and thus the emitted gravitational waves) and the neutron star matter in a BHNS merger depend strongly and nonlinearly on the black hole's spin. While there is a significant possibility that astrophysical black holes could have spins that are nearly extremal (i.e. near the theoretical maximum), to date fully relativistic simulations of BHNS binaries have included black-hole spins only up to $S/M^2$=0.9, which corresponds to the black hole having approximately half as much rotational energy as possible, given the black hole's mass. In this paper, we present a new simulation of a BHNS binary with a mass ratio $q=3$ and black-hole spin $S/M^2$=0.97, the highest simulated to date. We find that the black hole's large spin leads to the most massive accretion disk and the largest tidal tail outflow of any fully relativistic BHNS simulations to date, even exceeding the results implied by extrapolating results from simulations with lower black-hole spin. The disk appears to be remarkably stable. We also find that the high black-hole spin persists until shortly before the time of merger; afterwards, both merger and accretion spin down the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Periastron advance in spinning black hole binaries: comparing effective-one-body and Numerical Relativity

We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasi-circular orbits and having spins collinear with the orbital angular momentum. We compare the predictions with the periastron advance recently computed in accurate numerical-relativity simulations and find remarkable agreement for a wide range of spins and mass ratios. These results do not use any numerical-relativity calibration of the effective-one-body model, and stem from two key ingredients in the effective-one-body Hamiltonian: (i) the mapping of the two-body dynamics of spinning particles onto the dynamics of an effective spinning particle in a (deformed) Kerr spacetime, fully symmetrized with respect to the two-body masses and spins, and (ii) the resummation, in the test-particle limit, of all post-Newtonian (PN) corrections linear in the spin of the particle. In fact, even when only the leading spin PN corrections are included in the effective-one-body spinning Hamiltonian but all the test-particle corrections linear in the spin of the particle are resummed we find very good agreement with the numerical results (within the numerical error for equal-mass binaries and discrepancies of at most 1% for larger mass ratios). Furthermore, we specialize to the extreme mass-ratio limit and derive, using the equations of motion in the gravitational skeleton approach, analytical expressions for the periastron advance, the meridional Lense-Thirring precession and spin precession frequency in the case of a spinning particle on a nearly circular equatorial orbit in Kerr spacetime, including also terms quadratic in the spin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stability of nonspinning effective-one-body model in approximating two-body dynamics and gravitational-wave emission

The detection of gravitational waves and the extraction of physical information from them requires the prediction of accurate waveforms to be used in template banks. For that purpose, the accuracy of effective-one-body (EOB) waveforms has been improved over the last years by calibrating them to numerical-relativity (NR) waveforms. So far, the calibration has employed a handful of NR waveforms with a total length of ~30 cycles, the length being limited by the computational cost of NR simulations. Here we address the outstanding problem of the stability of the EOB calibration with respect to the length of NR waveforms. Performing calibration studies against NR waveforms of nonspinning black-hole binaries with mass ratios 1, 1.5, 5, and 8, and with a total length of ~60 cycles, we find that EOB waveforms calibrated against either 30 or 60 cycles will be indistinguishable by the advanced detectors LIGO and Virgo when the signal-to-noise ratio (SNR) is below 110. When extrapolating to a very large number of cycles, using very conservative assumptions, we can conclude that state-of-the-art nonspinning EOB waveforms of any length are sufficiently accurate for parameter estimation with advanced detectors when the SNR is below 20, the mass ratio is below 5 and total mass is above 20 Msun. The results are not conclusive for the entire parameter space because of current NR errors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Suitability of hybrid gravitational waveforms for unequal-mass binaries

This article studies sufficient accuracy criteria of hybrid post-Newtonian (PN) and numerical relativity (NR) waveforms for parameter estimation of strong binary black-hole sources in second- generation ground-based gravitational-wave detectors. We investigate equal-mass non-spinning binaries with a new 33-orbit NR waveform, as well as unequal-mass binaries with mass ratios 2, 3, 4 and 6. For equal masses, the 33-orbit NR waveform allows us to recover previous results and to extend the analysis toward matching at lower frequencies. For unequal masses, the errors between different PN approximants increase with mass ratio. Thus, at 3.5PN, hybrids for higher-mass-ratio systems would require NR waveforms with many more gravitational-wave (GW) cycles to guarantee no adverse impact on parameter estimation. Furthermore, we investigate the potential improvement in hybrid waveforms that can be expected from 4th order post-Newtonian waveforms, and find that knowledge of this 4th post-Newtonian order would significantly improve the accuracy of hybrid waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Template Banks for Binary black hole searches with Numerical Relativity waveforms

Gravitational waves (GW) from coalescing stellar-mass black hole binaries (BBH) are expected to be detected by the Advanced Laser Interferometer Gravitational-wave Observatory and Advanced Virgo. Detection searches operate by matched-filtering the detector data using a bank of waveform templates. Traditionally, template banks for BBH are constructed from intermediary analytical waveform models which are calibrated against numerical relativity simulations and which can be aluated for any choice of BBH parameters. This paper explores an alternative to the traditional approach, namely the construction of template banks directly from numerical BBH simulations. Using non-spinning BBH systems as an example, we demonstrate which regions of the mass-parameter plane can be covered with existing numerical BBH waveforms. We estimate the required number and required length of BBH simulations to cover the entire non-spinning BBH parameter plane up to mass-ratio 10, thus illustrating that our approach can be used to guide parameter placement of future numerical simulations. We derive error bounds which are independent of analytical waveform models; therefore, our formalism can be used to independently test the accuracy of such waveform models. The resulting template banks are suitable for advanced LIGO searches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Are different approaches to constructing initial data for binary black hole simulations of the same astrophysical situation equivalent?

Initial data for numerical evolutions of binary-black holes have been dominated by "conformally flat" (CF) data (i.e., initial data where the conformal background metric is chosen to be flat) because they are easy to construct. However, CF initial data cannot simulate nearly extremal spins, while more complicated "conformally curved" initial data (i.e., initial data in which the background metric is \emph{not} explicitly chosen to be flat), such as initial data where the spatial metric is chosen to be proportional to a weighted superposition of two Kerr-Schild (SKS) black holes can. Here we establish the consistency between the astrophysical results of these two initial data schemes for nonspinning binary systems. We evolve the inspiral, merger, and ringdown of two equal-mass, nonspinning black holes using SKS initial data and compare with an analogous simulation using CF initial data. We find that the resultant gravitational-waveform phases agree to within $δϕ\lesssim 10^{-2}$ radians and the amplitudes agree to within $δA/A \lesssim 5 \times 10^{-3}$, which are within the numerical errors of the simulations. Furthermore, we find that the final mass and spin of the remnant black hole agree to one part in $10^{5}

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black hole-neutron star mergers for 10 solar mass black holes

General relativistic simulations of black hole-neutron star mergers have currently been limited to low-mass black holes (less than 7 solar mass), even though population synthesis models indicate that a majority of mergers might involve more massive black holes (10 solar mass or more). We present the first general relativistic simulations of black hole-neutron star mergers with 10 solar mass black holes. For massive black holes, the tidal forces acting on the neutron star are usually too weak to disrupt the star before it reaches the innermost stable circular orbit of the black hole. Varying the spin of the black hole in the range a/M = 0.5-0.9, we find that mergers result in the disruption of the star and the formation of a massive accretion disk only for large spins a/M>0.7-0.9. From these results, we obtain updated constraints on the ability of BHNS mergers to be the progenitors of short gamma-ray bursts as a function of the mass and spin of the black hole. We also discuss the dependence of the gravitational wave signal on the black hole parameters, and provide waveforms and spectra from simulations beginning 7-8 orbits before merger.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of non-equal mass black hole binaries with spectral methods

This paper presents techniques and results for simulations of unequal mass, non-spinning black hole binaries with pseudo-spectral methods. Specifically, we develop an efficient root-finding procedure to ensure the black hole initial data have the desired masses and spins, we extend the dual coordinate frame method and eccentricity removal to asymmetric binaries. Furthermore, we describe techniques to simulate mergers of unequal mass black holes. The second part of the paper presents numerical simulations of non-spinning black hole binaries with mass ratios 2, 3, 4 and 6, covering between 15 and 22 orbits, merger and ringdown. We discuss the accuracy of these simulations, the evolution of the (initially zero) black hole spins, and the remnant black hole properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving Partial Differential Equations Numerically on Manifolds with Arbitrary Spatial Topologies

A multi-cube method is developed for solving systems of elliptic and hyperbolic partial differential equations numerically on manifolds with arbitrary spatial topologies. It is shown that any three-dimensional manifold can be represented as a set of non-overlapping cubic regions, plus a set of maps to identify the faces of adjoining regions. The differential structure on these manifolds is fixed by specifying a smooth reference metric tensor. Matching conditions that ensure the appropriate levels of continuity and differentiability across region boundaries are developed for arbitrary tensor fields. Standard numerical methods are then used to solve the equations with the appropriate boundary conditions, which are determined from these inter-region matching conditions. Numerical examples are presented which use pseudo-spectral methods to solve simple elliptic equations on multi-cube representations of manifolds with the topologies T^3, S^2 x S^1 and S^3. Examples are also presented of numerical solutions of simple hyperbolic equations on multi-cube manifolds with the topologies R x T^3, R x S^2 x S^1 and R x S^3.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accurate gravitational waveforms for binary-black-hole mergers with nearly extremal spins

Motivated by the possibility of observing gravitational waves from merging black holes whose spins are nearly extremal (i.e., 1 in dimensionless units), we present numerical waveforms from simulations of merging black holes with the highest spins simulated to date: (1) a 25.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.97 aligned with the orbital angular momentum; and (2) a previously reported 12.5-orbit inspiral, merger, and ringdown of two holes with equal masses and spins of magnitude 0.95 anti-aligned with the orbital angular momentum. First, we consider the horizon mass and spin evolution of the new aligned-spin simulation. During the inspiral, the horizon area and spin evolve in remarkably close agreement with Alvi's analytic predictions, and the remnant hole's final spin agrees reasonably well with several analytic predictions. We also find that the total energy emitted by a real astrophysical system with these parameters---almost all of which is radiated during the time included in this simulation---would be 10.952% of the initial mass at infinite separation. Second, we consider the gravitational waveforms for both simulations. After estimating their uncertainties, we compare the waveforms to several post-Newtonian approximants, finding significant disagreement well before merger, although the phase of the TaylorT4 approximant happens to agree remarkably well with the numerical prediction in the aligned-spin case. We find that the post-Newtonian waveforms have sufficient uncertainty that hybridized waveforms will require far longer numerical simulations (in the absence of improved post-Newtonian waveforms) for accurate parameter estimation of low-mass binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulating merging binary black holes with nearly extremal spins

Astrophysically realistic black holes may have spins that are nearly extremal (i.e., close to 1 in dimensionless units). Numerical simulations of binary black holes are important tools both for calibrating analytical templates for gravitational-wave detection and for exploring the nonlinear dynamics of curved spacetime. However, all previous simulations of binary-black-hole inspiral, merger, and ringdown have been limited by an apparently insurmountable barrier: the merging holes' spins could not exceed 0.93, which is still a long way from the maximum possible value in terms of the physical effects of the spin. In this paper, we surpass this limit for the first time, opening the way to explore numerically the behavior of merging, nearly extremal black holes. Specifically, using an improved initial-data method suitable for binary black holes with nearly extremal spins, we simulate the inspiral (through 12.5 orbits), merger and ringdown of two equal-mass black holes with equal spins of magnitude 0.95 antialigned with the orbital angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Spin Diagrams for Equal-Mass Black-Hole Binaries with Aligned Spins

Binary black-hole systems with spins aligned with the orbital angular momentum are of special interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbitrary. By least-square fitting the results of these simulations we have constructed two "spin diagrams" which provide straightforward information about the recoil velocity |v_kick| and the final black-hole spin a_fin in terms of the dimensionless spins a_1 and a_2 of the two initial black holes. Overall they suggest a maximum recoil velocity of |v_kick|=441.94 km/s, and minimum and maximum final spins a_fin=0.3471 and a_fin=0.9591, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Some mathematical problems in numerical relativity

The main goal of numerical relativity is the long time simulation of highly nonlinear spacetimes that cannot be treated by perturbation theory. This involves analytic, computational and physical issues. At present, the major impasses to achieving global simulations of physical usefulness are of an analytic/computational nature. We present here some examples of how analytic insight can lend useful guidance for the improvement of numerical approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Well-Posed Initial-Boundary Evolution in General Relativity

Maximally dissipative boundary conditions are applied to the initial-boundary value problem for Einstein's equations in harmonic coordinates to show that it is well-posed for homogeneous boundary data and for boundary data that is small in a linearized sense. The method is implemented as a nonlinear evolution code which satisfies convergence tests in the nonlinear regime and is robustly stable in the weak field regime. A linearized version has been stably matched to a characteristic code to compute the gravitational waveform radiated to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Boundary conditions in linearized harmonic gravity

We investigate the initial-boundary value problem for linearized gravitational theory in harmonic coordinates. Rigorous techniques for hyperbolic systems are applied to establish well-posedness for various reductions of the system into a set of six wave equations. The results are used to formulate computational algorithms for Cauchy evolution in a 3-dimensional bounded domain. Numerical codes based upon these algorithms are shown to satisfy tests of robust stability for random constraint violating initial data and random boundary data; and shown to give excellent performance for the evolution of typical physical data. The results are obtained for plane boundaries as well as piecewise cubic spherical boundaries cut out of a Cartesian grid.

0 citations0 reviewsNo code linkedOpen access