-
AsterX: a new open-source GPU-accelerated GRMHD code for dynamical spacetimes
Authors:
Jay V. Kalinani,
Liwei Ji,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Lucas Timotheo Sanches,
Bing-Jyun Tsao,
Steven R Brandt,
Manuela Campanelli,
Riccardo Ciolfi,
Bruno Giacomazzo,
Roland Haas,
Erik Schnetter,
Yosef Zlochower
Abstract:
We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel…
▽ More
We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel applications. AsterX employs the Valencia formulation for GRMHD, coupled with the 'Z4c' formalism for spacetime evolution, while incorporating high resolution shock capturing schemes to accurately handle the hydrodynamics. AsterX has undergone rigorous testing in both static and dynamic spacetime, demonstrating remarkable accuracy and agreement with other codes in literature. Benchmarking the code through scaling tests on OLCF's Frontier supercomputer, we demonstrate a weak scaling efficiency of about 67%-77% on 4096 nodes compared to an 8-node performance.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes
Authors:
Liwei Ji,
Vassilios Mewes,
Yosef Zlochower,
Lorenzo Ennoggi,
Federico G. Lopez Armengol,
Manuela Campanelli,
Federico Cipolletta,
Zachariah B. Etienne
Abstract:
Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of flui…
▽ More
Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition timestep limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double FFT filter and implementing it within the fully MPI-parallelized SphericalNR framework in the Einstein Toolkit. We demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the timestep significantly compared to unfiltered simulations.
△ Less
Submitted 9 November, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Massively parallel simulations of binary black holes with Dendro-GR
Authors:
Milinda Fernando,
David Neilsen,
Yosef Zlochower,
Eric W. Hirschmann,
Hari Sundar
Abstract:
We present results from the new Dendro-GR code. These include simulations of binary black hole mergers for mass ratios up to q=16. Dendro-GR uses Wavelet Adaptive Multi-Resolution (WAMR) to generate an unstructured grid adapted to the spacetime geometry together with an octree based data structure. We demonstrate good scaling, improved convergence properties and efficient use of computational reso…
▽ More
We present results from the new Dendro-GR code. These include simulations of binary black hole mergers for mass ratios up to q=16. Dendro-GR uses Wavelet Adaptive Multi-Resolution (WAMR) to generate an unstructured grid adapted to the spacetime geometry together with an octree based data structure. We demonstrate good scaling, improved convergence properties and efficient use of computational resources. We validate the code with comparisons to LazEv.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations
Authors:
Federico G. Lopez Armengol,
Zachariah B. Etienne,
Scott C. Noble,
Bernard J. Kelly,
Leonardo R. Werneck,
Brendan Drachler,
Manuela Campanelli,
Federico Cipolletta,
Yosef Zlochower,
Ariadna Murguia-Berthier,
Lorenzo Ennoggi,
Mark Avara,
Riccardo Ciolfi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Julian H. Krolik,
Vassilios Mewes,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril,
Jeremy Schnittman
Abstract:
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg…
▽ More
We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach.
△ Less
Submitted 31 October, 2022; v1 submitted 17 December, 2021;
originally announced December 2021.
-
The population of merging compact binaries inferred using gravitational waves through GWTC-3
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
A. Allocca,
P. A. Altin,
A. Amato
, et al. (1612 additional authors not shown)
Abstract:
We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8…
▽ More
We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8 $\rm{Gpc^{-3}\, yr^{-1}}$ and 140 $\rm{Gpc^{-3} yr^{-1}}$ , assuming a constant rate density versus comoving volume and taking the union of 90% credible intervals for methods used in this work. Accounting for the BBH merger rate to evolve with redshift, we find the BBH merger rate to be between 17.9 $\rm{Gpc^{-3}\, yr^{-1}}$ and 44 $\rm{Gpc^{-3}\, yr^{-1}}$ at a fiducial redshift (z=0.2). We obtain a broad neutron star mass distribution extending from $1.2^{+0.1}_{-0.2} M_\odot$ to $2.0^{+0.3}_{-0.3} M_\odot$. We can confidently identify a rapid decrease in merger rate versus component mass between neutron star-like masses and black-hole-like masses, but there is no evidence that the merger rate increases again before 10 $M_\odot$. We also find the BBH mass distribution has localized over- and under-densities relative to a power law distribution. While we continue to find the mass distribution of a binary's more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above $\sim 60 M_\odot$. The rate of BBH mergers is observed to increase with redshift at a rate proportional to $(1+z)^κ$ with $κ= 2.9^{+1.7}_{-1.8}$ for $z\lesssim 1$. Observed black hole spins are small, with half of spin magnitudes below $χ_i \simeq 0.25$. We observe evidence of negative aligned spins in the population, and an increase in spin magnitude for systems with more unequal mass ratio.
△ Less
Submitted 23 February, 2022; v1 submitted 5 November, 2021;
originally announced November 2021.
-
GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
T. Akutsu,
S. Albanesi,
A. Allocca,
P. A. Altin
, et al. (1637 additional authors not shown)
Abstract:
The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There ar…
▽ More
The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There are 35 compact binary coalescence candidates identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} > 0.5$. Of these, 18 were previously reported as low-latency public alerts, and 17 are reported here for the first time. Based upon estimates for the component masses, our O3b candidates with $p_\mathrm{astro} > 0.5$ are consistent with gravitational-wave signals from binary black holes or neutron star-black hole binaries, and we identify none from binary neutron stars. However, from the gravitational-wave data alone, we are not able to measure matter effects that distinguish whether the binary components are neutron stars or black holes. The range of inferred component masses is similar to that found with previous catalogs, but the O3b candidates include the first confident observations of neutron star-black hole binaries. Including the 35 candidates from O3b in addition to those from GWTC-2.1, GWTC-3 contains 90 candidates found by our analysis with $p_\mathrm{astro} > 0.5$ across the first three observing runs. These observations of compact binary coalescences present an unprecedented view of the properties of black holes and neutron stars.
△ Less
Submitted 23 October, 2023; v1 submitted 5 November, 2021;
originally announced November 2021.
-
GWTC-2.1: Deep Extended Catalog of Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
S. Albanesi,
A. Allocca,
P. A. Altin,
A. Amato,
C. Anand,
S. Anand
, et al. (1407 additional authors not shown)
Abstract:
The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibra…
▽ More
The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the astrophysical probability for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have an astrophysical probability greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. If the 8 additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects $\geq 3M_\odot$) is increased compared to GWTC-2, with total masses from $\sim 14 M_\odot$ for GW190924_021846 to $\sim 182 M_\odot$ for GW190426_190642. The primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than $0.65$ and $0.44$ at $90\%$ probability for GW190403_051519 and GW190917_114630 respectively), and find that 2 of the 8 new events have effective inspiral spins $χ_\mathrm{eff} > 0$ (at $90\%$ credibility), while no binary is consistent with $χ_\mathrm{eff} < 0$ at the same significance.
△ Less
Submitted 10 May, 2022; v1 submitted 2 August, 2021;
originally announced August 2021.
-
HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage
Authors:
Ariadna Murguia-Berthier,
Scott C. Noble,
Luke F. Roberts,
Enrico Ramirez-Ruiz,
Leonardo R. Werneck,
Michael Kolacki,
Zachariah B. Etienne,
Mark Avara,
Manuela Campanelli,
Riccardo Ciolfi,
Federico Cipolletta,
Brendan Drachler,
Lorenzo Ennoggi,
Joshua Faber,
Grace Fiacco,
Bruno Giacomazzo,
Tanmayee Gupte,
Trung Ha,
Bernard J. Kelly,
Julian H. Krolik,
Federico G. Lopez Armengol,
Ben Margalit,
Tim Moon,
Richard O'Shaughnessy,
Jesús M. Rueda-Becerril
, et al. (3 additional authors not shown)
Abstract:
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge…
▽ More
The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.
△ Less
Submitted 31 January, 2022; v1 submitted 9 June, 2021;
originally announced June 2021.
-
GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
Authors:
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva,
S. B. Anderson
, et al. (1327 additional authors not shown)
Abstract:
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave event…
▽ More
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1.
△ Less
Submitted 8 March, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Properties and astrophysical implications of the 150 Msun binary black hole merger GW190521
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
A. Aich,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand
, et al. (1233 additional authors not shown)
Abstract:
The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theo…
▽ More
The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range $65 - 120\,M_{\odot}$. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger $(142^{+28}_{-16}\,M_{\odot})$ classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular binary black hole coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be $0.13^{+0.30}_{-0.11}\,{\rm Gpc}^{-3}\,\rm{yr}^{-1}$. We discuss the astrophysical implications of GW190521 for stellar collapse, and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescence, or via hierarchical merger of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
Authors:
Vassilios Mewes,
Yosef Zlochower,
Manuela Campanelli,
Thomas W. Baumgarte,
Zachariah B. Etienne,
Federico G. Lopez Armengol,
Federico Cipolletta
Abstract:
We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-S…
▽ More
We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.
△ Less
Submitted 5 May, 2020; v1 submitted 14 February, 2020;
originally announced February 2020.
-
PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order
Authors:
Dennis B. Bowen,
Mark Avara,
Vassilios Mewes,
Yosef Zlochower,
Scott C. Noble,
Manuela Campanelli,
Hotaka Shiokawa,
Roseanne M. Cheng,
Julian H. Krolik
Abstract:
We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprise…
▽ More
We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprised of an arbitrary number of moving, local meshes, or patches, which are free to employ their own resolution, coordinate system/topology, physics equations, reference frame, and in our new approach, numerical method. Each local patch exchanges boundary data with a single global patch on which all other patches reside through a client-router-server parallelization model. In generalizing Patchwork to be compatible with arbitrary order time integration, PatchworkMHD and PatchworkWave have significantly improved the interpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in the original Patchwork code. Furthermore, we extend Patchwork to be multimethod by allowing multiple state vectors to be updated simultaneously, with each state vector providing its own interpatch interpolation and transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partial differential equations. We demonstrate our changes through the implementation of a scalar wave toy-model that is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.
△ Less
Submitted 31 January, 2020;
originally announced February 2020.
-
Hybrid waveforms for generic precessing binaries for gravitational-wave data analysis
Authors:
Jam Sadiq,
Yosef Zlochower,
Richard O'Shaughnessy,
Jacob Lange
Abstract:
We construct hybrid binary black holes merger waveforms using analytical model waveforms for the early inspiral phase and numerical relativity waveforms for late inspiral to merger and post merger phases. To hybridize analytical and numerical waveforms, we first perform a 3-dimensional rotation to align the instantaneous orbital planes associated with the two waveforms at some fiducial frequency,…
▽ More
We construct hybrid binary black holes merger waveforms using analytical model waveforms for the early inspiral phase and numerical relativity waveforms for late inspiral to merger and post merger phases. To hybridize analytical and numerical waveforms, we first perform a 3-dimensional rotation to align the instantaneous orbital planes associated with the two waveforms at some fiducial frequency, we then find appropriate phase and time translations that maximize the overlap of the two waveforms in a hybridization interval. We discuss the accuracy and limitations for hybrids constructed by this procedure in the context of LIGO-Virgo-KAGRA observations. Our goal is to hybridize waveforms for more generic precessing binaries and construct longer waveforms that are sufficiently accurate for the parameter estimation techniques for upcoming LIGO observations.
△ Less
Submitted 7 July, 2020; v1 submitted 20 January, 2020;
originally announced January 2020.
-
Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
M. A. Aloy,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1169 additional authors not shown)
Abstract:
GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of stat…
▽ More
GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.
△ Less
Submitted 6 March, 2020; v1 submitted 2 August, 2019;
originally announced August 2019.
-
On the properties of the massive binary black hole merger GW170729
Authors:
Katerina Chatziioannou,
Roberto Cotesta,
Sudarshan Ghonge,
Jacob Lange,
Ken K. -Y. Ng,
Juan Calderon Bustillo,
James Clark,
Carl-Johan Haster,
Sebastian Khan,
Michael Puerrer,
Vivien Raymond,
Salvatore Vitale,
Nousha Afshari,
Stanislav Babak,
Kevin Barkett,
Jonathan Blackman,
Alejandro Bohe,
Michael Boyle,
Alessandra Buonanno,
Manuela Campanelli,
Gregorio Carullo,
Tony Chu,
Eric Flynn,
Heather Fong,
Alyssa Garcia
, et al. (31 additional authors not shown)
Abstract:
We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusio…
▽ More
We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario.
△ Less
Submitted 7 November, 2019; v1 submitted 15 March, 2019;
originally announced March 2019.
-
The second RIT binary black hole simulations catalog and its application to gravitational waves parameter estimation
Authors:
James Healy,
Carlos O. Lousto,
Jacob Lange,
Richard O'Shaughnessy,
Yosef Zlochower,
Manuela Campanelli
Abstract:
The RIT numerical relativity group is releasing the second public catalog of black-hole-binary waveforms \url{http://ccrg.rit.edu/~RITCatalog}. This release consists of 320 accurate simulations that include 46 precessing and 274 nonprecessing binary systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$ and individual spins up to $s/m^2=0.95$. The new catalog contains search and orderi…
▽ More
The RIT numerical relativity group is releasing the second public catalog of black-hole-binary waveforms \url{http://ccrg.rit.edu/~RITCatalog}. This release consists of 320 accurate simulations that include 46 precessing and 274 nonprecessing binary systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$ and individual spins up to $s/m^2=0.95$. The new catalog contains search and ordering tools for the waveforms based on initial parameters of the binary, trajectory information, peak radiation, and final remnant black hole properties. The final black hole remnant properties provided here can be used to model the merger of black-hole binaries from its initial configurations. The waveforms are extrapolated to infinite observer location and can be used to independently interpret gravitational wave signals from laser interferometric detectors. As an application of this waveform catalog we reanalyze the signal of GW150914 implementing parameter estimation techniques that make use of only numerical waveforms without any reference to information from phenomenological models.
△ Less
Submitted 27 March, 2019; v1 submitted 8 January, 2019;
originally announced January 2019.
-
GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
M. A. Aloy,
P. A. Altin,
A. Amato,
A. Ananyeva
, et al. (1126 additional authors not shown)
Abstract:
We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1$\mathrm{M}_\odot$ during the first and second observing runs of the Advanced gravitational-wave detector network. During the first observing run (O1), from September $12^\mathrm{th}$, 2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary black hole mer…
▽ More
We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1$\mathrm{M}_\odot$ during the first and second observing runs of the Advanced gravitational-wave detector network. During the first observing run (O1), from September $12^\mathrm{th}$, 2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run (O2), which ran from November $30^\mathrm{th}$, 2016 to August $25^\mathrm{th}$, 2017, saw the first detection of gravitational waves from a binary neutron star inspiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818 and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between $18.6_{-0.7}^{+3.2}\mathrm{M}_\odot$, and $84.4_{-11.1}^{+15.8} \mathrm{M}_\odot$, and range in distance between $320_{-110}^{+120}$ Mpc and $2840_{-1360}^{+1400}$ Mpc. No neutron star - black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of $110\, -\, 3840$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary neutron stars and $9.7\, -\, 101$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary black holes assuming fixed population distributions, and determine a neutron star - black hole merger rate 90% upper limit of $610$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$.
△ Less
Submitted 4 October, 2019; v1 submitted 30 November, 2018;
originally announced November 2018.
-
Numerical Relativity of Compact Binaries in the 21st Century
Authors:
Matthew D. Duez,
Yosef Zlochower
Abstract:
We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star--neutron s…
▽ More
We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star--neutron star mergers. Modeling such events requires numerical integration of the field equations of general relativity in three spatial dimensions, coupled, in the case of neutron-star containing binaries, with increasingly sophisticated treatment of fluids, electromagnetic fields, and neutrino radiation. However, it was not until 2005 that accurate long-term evolutions of binaries containing black holes were even possible. Since then, there has been an explosion of new results and insights into the physics of strongly-gravitating system. Particular emphasis has been placed on understanding the gravitational wave and electromagnetic signatures from these extreme events. And with the recent dramatic discoveries of gravitational waves from merging black holes by the Laser Interferometric Gravitational Wave Observatory and Virgo, and the subsequent discovery of both electromagnetic and gravitational wave signals from a merging neutron star binary, numerical relativity became an indispensable tool for the new field of multimessenger astronomy.
△ Less
Submitted 13 September, 2018; v1 submitted 17 August, 2018;
originally announced August 2018.
-
GW170817: Measurements of Neutron Star Radii and Equation of State
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
B. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
G. Allen,
A. Allocca,
M. A. Aloy
, et al. (1127 additional authors not shown)
Abstract:
On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of th…
▽ More
On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(ρ)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 \,M_\odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(ρ)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}\times 10^{34} \,\mathrm{dyn}/\mathrm{cm}^{2}$ at the 90% level.
△ Less
Submitted 15 October, 2018; v1 submitted 29 May, 2018;
originally announced May 2018.
-
Properties of the binary neutron star merger GW170817
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
B. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
G. Allen,
A. Allocca,
M. A. Aloy
, et al. (1126 additional authors not shown)
Abstract:
On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spin…
▽ More
On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spins, and tidal parameters, using the known source location, improved modeling, and recalibrated Virgo data. We extend the range of gravitational-wave frequencies considered down to 23 Hz, compared to 30 Hz in the initial analysis. We also compare results inferred using several signal models, which are more accurate and incorporate additional physical effects as compared to the initial analysis. We improve the localization of the gravitational-wave source to a 90% credible region of $16~\mathrm{deg}^2$. We find tighter constraints on the masses, spins, and tidal parameters, and continue to find no evidence for nonzero component spins. The component masses are inferred to lie between 1.00 and 1.89 $M_\odot$ when allowing for large component spins, and to lie between 1.16 and 1.60 $M_\odot$ (with a total mass $2.73^{+0.04}_{-0.01} \, M_\odot$) when the spins are restricted to be within the range observed in Galactic binary neutron stars. Under minimal assumptions about the nature of the compact objects, our constraints for the tidal deformability parameter $\tilde Λ$ are $(0,630)$ when we allow for large component spins, and $300^{+420}_{-230}$ (using a 90% highest posterior density interval) when restricting the magnitude of the component spins, ruling out several equation-of-state models at the 90% credible level. Finally, with LIGO and GEO600 data, we use a Bayesian analysis to place upper limits on the amplitude and spectral energy density of a possible post-merger signal. (Abridged)
△ Less
Submitted 22 January, 2019; v1 submitted 29 May, 2018;
originally announced May 2018.
-
Numerical relativity in spherical coordinates with the Einstein Toolkit
Authors:
Vassilios Mewes,
Yosef Zlochower,
Manuela Campanelli,
Ian Ruchlin,
Zachariah B. Etienne,
Thomas W. Baumgarte
Abstract:
Numerical relativity codes that do not make assumptions on spatial symmetries most commonly adopt Cartesian coordinates. While these coordinates have many attractive features, spherical coordinates are much better suited to take advantage of approximate symmetries in a number of astrophysical objects, including single stars, black holes and accretion disks. While the appearance of coordinate singu…
▽ More
Numerical relativity codes that do not make assumptions on spatial symmetries most commonly adopt Cartesian coordinates. While these coordinates have many attractive features, spherical coordinates are much better suited to take advantage of approximate symmetries in a number of astrophysical objects, including single stars, black holes and accretion disks. While the appearance of coordinate singularities often spoils numerical relativity simulations in spherical coordinates, especially in the absence of any symmetry assumptions, it has recently been demonstrated that these problems can be avoided if the coordinate singularities are handled analytically. This is possible with the help of a reference-metric version of the Baumgarte-Shapiro-Shibata-Nakamura formulation together with a proper rescaling of tensorial quantities. In this paper we report on an implementation of this formalism in the Einstein Toolkit. We adapt the Einstein Toolkit infrastructure, originally designed for Cartesian coordinates, to handle spherical coordinates, by providing appropriate boundary conditions at both inner and outer boundaries. We perform numerical simulations for a disturbed Kerr black hole, extract the gravitational wave signal, and demonstrate that the noise in these signals is orders of magnitude smaller when computed on spherical grids rather than Cartesian grids. With the public release of our new Einstein Toolkit thorns, our methods for numerical relativity in spherical coordinates will become available to the entire numerical relativity community.
△ Less
Submitted 2 May, 2018; v1 submitted 26 February, 2018;
originally announced February 2018.
-
Comparing an analytical spacetime metric for a merging binary to a fully nonlinear numerical evolution using curvature scalars
Authors:
Jam Sadiq,
Yosef Zlochower,
Hiroyuki Nakano
Abstract:
We introduce a new geometrically invariant prescription for comparing two different spacetimes based on geodesic deviation. We use this method to compare a family of recently introduced analytical spacetime representing inspiraling black-hole binaries to fully nonlinear numerical solutions to the Einstein equations. Our method can be used to improve analytical spacetime models by providing a local…
▽ More
We introduce a new geometrically invariant prescription for comparing two different spacetimes based on geodesic deviation. We use this method to compare a family of recently introduced analytical spacetime representing inspiraling black-hole binaries to fully nonlinear numerical solutions to the Einstein equations. Our method can be used to improve analytical spacetime models by providing a local measure of the effects that violations of the Einstein equations will have on timelike geodesics, and indirectly, gas dynamics. We also discuss the advantages and limitations of this method.
△ Less
Submitted 28 March, 2018; v1 submitted 8 February, 2018;
originally announced February 2018.
-
Targeted numerical simulations of binary black holes for GW170104
Authors:
James Healy,
Jacob Lange,
Richard O'Shaughnessy,
Carlos Lousto,
Manuela Campanelli,
Andrew Williamson,
Yosef Zlochower,
Juan Calderon Bustillo,
James Clark,
Christopher Evans,
D. Ferguson,
Sudarshan Ghonge,
Karan Jani,
Bhavesh Khamesra,
Pablo Laguna,
Deirdre M. Shoemaker,
Alyssa Garcia,
Michael Boyle,
Daniel Hemberger,
Lawrence Kidder,
Prayush Kumar,
Geoffrey Lovelace,
Harald Pfeiffer,
Mark Scheel,
Saul Teukolsky
Abstract:
In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at multiple resolutions and with two independent techniques to solve Einstein's equations. For the nonprecessing and precessing simulations, we demonstrate the tw…
▽ More
In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at multiple resolutions and with two independent techniques to solve Einstein's equations. For the nonprecessing and precessing simulations, we demonstrate the two techniques agree mode by mode, at a precision substantially in excess of statistical uncertainties in current LIGO's observations. Conversely, we demonstrate our full numerical solutions contain information which is not accurately captured with the approximate phenomenological models commonly used to infer compact binary parameters. To quantify the impact of these differences on parameter inference for GW170104 specifically, we compare the predictions of our simulations and these approximate models to LIGO's observations of GW170104.
△ Less
Submitted 19 December, 2017; v1 submitted 15 December, 2017;
originally announced December 2017.
-
Evolutions of unequal mass, highly spinning black hole binaries
Authors:
James Healy,
Carlos O. Lousto,
Ian Ruchlin,
Yosef Zlochower
Abstract:
We evolve a binary black hole system bearing a mass ratio of $q=m_1/m_2=2/3$ and individual spins of $S^z_1/m_1^2=0.95$ and $S^z_2/m_2^2=-0.95$ in a configuration where the large black hole has its spin antialigned with the orbital angular momentum, $L^z$, and the small black hole has its spin aligned with $L^z$. This configuration was chosen to measure the maximum recoil of the remnant black hole…
▽ More
We evolve a binary black hole system bearing a mass ratio of $q=m_1/m_2=2/3$ and individual spins of $S^z_1/m_1^2=0.95$ and $S^z_2/m_2^2=-0.95$ in a configuration where the large black hole has its spin antialigned with the orbital angular momentum, $L^z$, and the small black hole has its spin aligned with $L^z$. This configuration was chosen to measure the maximum recoil of the remnant black hole for nonprecessing binaries. We find that the remnant black hole recoils at 500km/s, the largest recorded value from numerical simulations for aligned spin configurations. The remnant mass, spin, and gravitational waveform peak luminosity and frequency also provide a valuable point in parameter space for source modeling.
△ Less
Submitted 22 November, 2017;
originally announced November 2017.
-
Evolutions of Nearly Maximally Spinning Black Hole Binaries Using the Moving Puncture Approach
Authors:
Yosef Zlochower,
James Healy,
Carlos O. Lousto,
Ian Ruchlin
Abstract:
We demonstrate that numerical relativity codes based on the moving punctures formalism are capable of evolving nearly maximally spinning black hole binaries. We compare a new evolution of an equal-mass, aligned-spin binary with dimensionless spin chi=0.99 using puncture-based data with recent simulations of the SXS Collaboration. We find that the overlap of our new waveform with the published resu…
▽ More
We demonstrate that numerical relativity codes based on the moving punctures formalism are capable of evolving nearly maximally spinning black hole binaries. We compare a new evolution of an equal-mass, aligned-spin binary with dimensionless spin chi=0.99 using puncture-based data with recent simulations of the SXS Collaboration. We find that the overlap of our new waveform with the published results of the SXS Collaboration is larger than 0.999. To generate our new waveform, we use the recently introduced HiSpID puncture data, the CCZ4 evolution system, and a modified lapse condition that helps keep the horizon radii reasonably large.
△ Less
Submitted 8 August, 2017; v1 submitted 6 June, 2017;
originally announced June 2017.
-
A Parameter Estimation Method that Directly Compares Gravitational Wave Observations to Numerical Relativity
Authors:
Jacob Lange,
Richard O'Shaughnessy,
Michael Boyle,
Juan Calderón Bustillo,
Manuela Campanelli,
Tony Chu,
James A. Clark,
Nicholas Demos,
Heather Fong,
James Healy,
Daniel Hemberger,
Ian Hinder,
Karan Jani,
Bhavesh Khamesra,
Lawrence E. Kidder,
Prayush Kumar,
Pablo Laguna,
Carlos O. Lousto,
Geoffrey Lovelace,
Serguei Ossokine,
Harald Pfeiffer,
Mark A. Scheel,
Deirdre Shoemaker,
Bela Szilagyi,
Saul Teukolsky
, et al. (1 additional authors not shown)
Abstract:
We present and assess a Bayesian method to interpret gravitational wave signals from binary black holes. Our method directly compares gravitational wave data to numerical relativity simulations. This procedure bypasses approximations used in semi-analytical models for compact binary coalescence. In this work, we use only the full posterior parameter distribution for generic nonprecessing binaries,…
▽ More
We present and assess a Bayesian method to interpret gravitational wave signals from binary black holes. Our method directly compares gravitational wave data to numerical relativity simulations. This procedure bypasses approximations used in semi-analytical models for compact binary coalescence. In this work, we use only the full posterior parameter distribution for generic nonprecessing binaries, drawing inferences away from the set of NR simulations used, via interpolation of a single scalar quantity (the marginalized log-likelihood, $\ln {\cal L}$) evaluated by comparing data to nonprecessing binary black hole simulations. We also compare the data to generic simulations, and discuss the effectiveness of this procedure for generic sources. We specifically assess the impact of higher order modes, repeating our interpretation with both $l\le2$ as well as $l\le3$ harmonic modes. Using the $l\le3$ higher modes, we gain more information from the signal and can better constrain the parameters of the gravitational wave signal. We assess and quantify several sources of systematic error that our procedure could introduce, including simulation resolution and duration; most are negligible. We show through examples that our method can recover the parameters for equal mass, zero spin; GW150914-like; and unequal mass, precessing spin sources. Our study of this new parameter estimation method demonstrates we can quantify and understand the systematic and statistical error. This method allows us to use higher order modes from numerical relativity simulations to better constrain the black hole binary parameters.
△ Less
Submitted 27 May, 2017;
originally announced May 2017.
-
The nonspinning binary black hole merger scenario revisited
Authors:
James Healy,
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We present the results of 14 simulations of nonspinning black hole binaries with mass ratios $q=m_1/m_2$ in the range $1/100\leq q\leq1$. For each of these simulations we perform three runs at increasing resolution to assess the finite difference errors and to extrapolate the results to infinite resolution. For $q\geq 1/6$, we follow the evolution of the binary typically for the last ten orbits pr…
▽ More
We present the results of 14 simulations of nonspinning black hole binaries with mass ratios $q=m_1/m_2$ in the range $1/100\leq q\leq1$. For each of these simulations we perform three runs at increasing resolution to assess the finite difference errors and to extrapolate the results to infinite resolution. For $q\geq 1/6$, we follow the evolution of the binary typically for the last ten orbits prior to merger. By fitting the results of these simulations, we accurately model the peak luminosity, peak waveform frequency and amplitude, and the recoil of the remnant hole for unequal mass nonspinning binaries. We verify the accuracy of these new models and compare them to previously existing empirical formulas. These new fits provide a basis for a hierarchical approach to produce more accurate remnant formulas in the generic precessing case. They also provide input to gravitational waveform modeling.
△ Less
Submitted 19 May, 2017;
originally announced May 2017.
-
Modeling the Black hole Merger of QSO 3C 186
Authors:
Carlos O. Lousto,
Yosef Zlochower,
Manuela Campanelli
Abstract:
Recent detailed observations of the radio-loud quasar 3C 186 indicate the possibility that a supermassive recoiling black hole is moving away from the host galaxy at a speed of nearly 2100km/s. If this is the case, we can model the mass ratio and spins of the progenitor binary black hole using the results of numerical relativity simulations. We find that the black holes in the progenitor must have…
▽ More
Recent detailed observations of the radio-loud quasar 3C 186 indicate the possibility that a supermassive recoiling black hole is moving away from the host galaxy at a speed of nearly 2100km/s. If this is the case, we can model the mass ratio and spins of the progenitor binary black hole using the results of numerical relativity simulations. We find that the black holes in the progenitor must have comparable masses with a mass ratio $q=m_1/m_2>1/4$ and the spin of the primary black hole must be $α_2=S_2/m_2^2>0.4$. The final remnant of the merger is bounded by $α_f>0.45$ and at least $4\%$ of the total mass of the binary system is radiated into gravitational waves. We consider four different pre-merger scenarios that further narrow those values. Assuming, for instance, a cold accretion driven merger model, we find that the binary had comparable masses with $q=0.70^{+0.29}_{-0.21}$ and the normalized spins of the larger and smaller black holes were $α_2=0.94^{+0.06}_{-0.22}$ and $α_1=0.95^{+0.05}_{-0.09}$. We can also estimate the final recoiling black hole spin $α_f=0.93^{+0.02}_{-0.03}$ and that the system radiated $9.6^{+0.8}_{-1.4}\%$ of its total mass, making the merger of those black holes the most energetic event ever observed.
△ Less
Submitted 3 April, 2017;
originally announced April 2017.
-
The RIT binary black hole simulations catalog
Authors:
James Healy,
Carlos O. Lousto,
Yosef Zlochower,
Manuela Campanelli
Abstract:
The RIT numerical relativity group is releasing a public catalog of black-hole-binary waveforms. The initial release of the catalog consists of 126 recent simulations that include precessing and non precessing systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$. The catalog contains information about the initial data of the simulation, the waveforms extrapolated to infinity, as well…
▽ More
The RIT numerical relativity group is releasing a public catalog of black-hole-binary waveforms. The initial release of the catalog consists of 126 recent simulations that include precessing and non precessing systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$. The catalog contains information about the initial data of the simulation, the waveforms extrapolated to infinity, as well as information about the peak luminosity and final remnant black hole properties. These waveforms can be used to independently interpret gravitational wave signals from laser interferometric detectors and
△ Less
Submitted 9 March, 2017;
originally announced March 2017.
-
Post-Newtonian Quasicircular Initial Orbits for Numerical Relativity
Authors:
James Healy,
Carlos O. Lousto,
Hiroyuki Nakano,
Yosef Zlochower
Abstract:
We use post-Newtonian (PN) approximations to determine the initial orbital and spin parameters of black hole binaries that lead to low-eccentricity inspirals when evolved with numerical relativity techniques. In particular, we seek initial configurations that lead to very small eccentricities at small separations, as is expected for astrophysical systems. We consider three cases: (i) quasicircular…
▽ More
We use post-Newtonian (PN) approximations to determine the initial orbital and spin parameters of black hole binaries that lead to low-eccentricity inspirals when evolved with numerical relativity techniques. In particular, we seek initial configurations that lead to very small eccentricities at small separations, as is expected for astrophysical systems. We consider three cases: (i) quasicircular orbits with no radial velocity, (ii) quasicircular orbits with an initial radial velocity determined by radiation reaction, and (iii) parameters obtained form evolution of the PN equations of motion from much larger separations. We study seven cases of spinning, nonprecessing, unequal mass binaries. We then use several definitions of the eccentricity, based on orbital separations and waveform phase and amplitude, and find that using the complete 3PN Hamiltonian for quasicircular orbits to obtain the tangential orbital momentum, and using the highest-known-order radiation reaction expressions to obtain the radial momentum, leads to the lowest eccentricity. The accuracy of this method even exceeds that of inspiral data based on 3PN and 4PN evolutions.
△ Less
Submitted 9 November, 2020; v1 submitted 2 February, 2017;
originally announced February 2017.
-
Modeling the source of GW150914 with targeted numerical-relativity simulations
Authors:
Geoffrey Lovelace,
Carlos O. Lousto,
James Healy,
Mark A. Scheel,
Alyssa Garcia,
Richard O'Shaughnessy,
Michael Boyle,
Manuela Campanelli,
Daniel A. Hemberger,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Bela Szilagyi,
Saul A. Teukolsky,
Yosef Zlochower
Abstract:
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 predic…
▽ More
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.
△ Less
Submitted 18 July, 2016;
originally announced July 2016.
-
Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (957 additional authors not shown)
Abstract:
We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing…
▽ More
We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin χ_eff \in [-0.3, 0.2], where χ_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and χ_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1].
△ Less
Submitted 3 June, 2016;
originally announced June 2016.
-
Tests of general relativity with GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca,
P. A. Altin,
S. B. Anderson
, et al. (963 additional authors not shown)
Abstract:
The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hol…
▽ More
The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hole merger in general relativity. We find that the final remnant's mass and spin, as determined from the low-frequency (inspiral) and high-frequency (post-inspiral) phases of the signal, are mutually consistent with the binary black-hole solution in general relativity. Furthermore, the data following the peak of GW150914 are consistent with the least-damped quasi-normal mode inferred from the mass and spin of the remnant black hole. By using waveform models that allow for parameterized general-relativity violations during the inspiral and merger phases, we perform quantitative tests on the gravitational-wave phase in the dynamical regime and we determine the first empirical bounds on several high-order post-Newtonian coefficients. We constrain the graviton Compton wavelength, assuming that gravitons are dispersed in vacuum in the same way as particles with mass, obtaining a $90\%$-confidence lower bound of $10^{13}$ km. In conclusion, within our statistical uncertainties, we find no evidence for violations of general relativity in the genuinely strong-field regime of gravity.
△ Less
Submitted 23 October, 2018; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Properties of the Binary Black Hole Merger GW150914
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
The LIGO Scientific Collaboration,
the Virgo Collaboration,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
B. Allen,
A. Allocca
, et al. (980 additional authors not shown)
Abstract:
On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in ge…
▽ More
On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.
△ Less
Submitted 20 June, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
High Energy Collisions of Black Holes Numerically Revisited
Authors:
James Healy,
Ian Ruchlin,
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly…
▽ More
We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly reduces the spurious radiation content, allowing for initial speeds much closer to the speed of light, i.e. $v\sim0.99c$. Using these new techniques, We estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{\text{ADM}}=0.13\pm0.01$. This value differs from the second-order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and by previous numerical estimates $(0.14\pm0.03)$.
△ Less
Submitted 24 August, 2016; v1 submitted 19 June, 2015;
originally announced June 2015.
-
Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques
Authors:
Yosef Zlochower,
Hiroyuki Nakano,
Bruno C. Mundim,
Manuela Campanelli,
Scott Noble,
Miguel Zilhao
Abstract:
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 transitio…
▽ More
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.
△ Less
Submitted 23 June, 2016; v1 submitted 1 April, 2015;
originally announced April 2015.
-
Modeling the remnant mass, spin, and recoil from unequal-mass, precessing black-hole binaries: The Intermediate Mass Ratio Regime
Authors:
Yosef Zlochower,
Carlos O. Lousto
Abstract:
We revisit the modeling of the properties of the remnant black hole resulting the merger of a black-hole binary as a function of the parameters of the binary. We provide a set of empirical formulas for the final mass, spin and recoil velocity of the final black hole as a function of the mass ratio and individual spins of the progenitor. In order to determine the fitting coefficients for these form…
▽ More
We revisit the modeling of the properties of the remnant black hole resulting the merger of a black-hole binary as a function of the parameters of the binary. We provide a set of empirical formulas for the final mass, spin and recoil velocity of the final black hole as a function of the mass ratio and individual spins of the progenitor. In order to determine the fitting coefficients for these formulas, we perform a set of 128 new numerical evolutions of precessing, unequal-mass black-hole binaries, and fit to the resulting remnant mass, spin, and recoil. In order to reduce the complexity of the analysis, we chose configurations that have one of the black holes spinning, with dimensionless spin alpha=0.8, at different angles with respect to the orbital angular momentum, and the other non-spinning. In addition to evolving families of binaries with different spin-inclination angles, we also evolved binaries with mass ratios as small as q=1/6. We use the resulting empirical formulas to predict the probabilities of black hole mergers leading to a given recoil velocity, total radiated gravitational energy, and final black hole spin.
△ Less
Submitted 20 July, 2016; v1 submitted 25 March, 2015;
originally announced March 2015.
-
Perturbative extraction of gravitational waveforms generated with Numerical Relativity
Authors:
H. Nakano,
J. Healy,
C. O. Lousto,
Y. Zlochower
Abstract:
We derive an analytical expression for extracting the gravitational waveforms at null infinity using the Weyl scalar $ψ_4$ measured at a finite radius. Our expression is based on a series solution in orders of 1/r to the equations for gravitational perturbations about a spinning black hole. We compute this expression to order $1/r^2$ and include the spin parameter $a$ of the Kerr background. We te…
▽ More
We derive an analytical expression for extracting the gravitational waveforms at null infinity using the Weyl scalar $ψ_4$ measured at a finite radius. Our expression is based on a series solution in orders of 1/r to the equations for gravitational perturbations about a spinning black hole. We compute this expression to order $1/r^2$ and include the spin parameter $a$ of the Kerr background. We test the accuracy of this extraction procedure by measuring the waveform for a merging black-hole binary at ten different extraction radii (in the range r/M=75-190) and for three different resolutions in the convergence regime. We find that the extraction formula provides a set of values for the radiated energy and momenta that at finite extraction radii converges towards the expected values with increasing resolution, which is not the case for the `raw' waveform at finite radius. We also examine the phase and amplitude errors in the waveform as a function of observer location and again observe the benefits of using our extraction formula. The leading corrections to the phase are ${\cal O}(1/r)$ and to the amplitude are ${\cal O}(1/r^2)$. This method provides a simple and practical way of estimating the waveform at infinity, and may be especially useful for scenarios such as well separated binaries, where the radiation zone is far from the sources, that would otherwise require extended simulation grids in order to extrapolate the `raw' waveform to infinity. Thus this method saves important computational resources and provides an estimate of errors.
△ Less
Submitted 28 April, 2015; v1 submitted 2 March, 2015;
originally announced March 2015.
-
Puncture Initial Data for Black-Hole Binaries with High Spins and High Boosts
Authors:
Ian Ruchlin,
James Healy,
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We solve the Hamiltonian and momentum constraints of general relativity for two black holes with nearly extremal spins and relativistic boosts in the puncture formalism. We use a non-conformally-flat ansatz with an attenuated superposition of two Lorentz-boosted, conformally Kerr or conformally Schwarzschild 3-metrics and their corresponding extrinsic curvatures. We compare evolutions of these dat…
▽ More
We solve the Hamiltonian and momentum constraints of general relativity for two black holes with nearly extremal spins and relativistic boosts in the puncture formalism. We use a non-conformally-flat ansatz with an attenuated superposition of two Lorentz-boosted, conformally Kerr or conformally Schwarzschild 3-metrics and their corresponding extrinsic curvatures. We compare evolutions of these data with the standard Bowen-York conformally flat ansatz (technically limited to intrinsic spins $χ=S/M^2_{\text{ADM}}=0.928$ and boosts $P/M_{\text{ADM}}=0.897$), finding, typically, an order of magnitude smaller burst of spurious radiation and agreement with inspiral and merger. As a first case study, we evolve two equal-mass black holes from rest with an initial separation of $d=12M$ and spins $χ_i=S_i/m_i^2=0.99$, compute the waveforms produced by the collision, the energy and angular momentum radiated, and the recoil of the final remnant black hole. We find that the black-hole trajectories curve at close separations, leading to the radiation of angular momentum. We also study orbiting nonspinning and moderate-spin black-hole binaries and compare these with standard Bowen-York data. We find a substantial reduction in the nonphysical initial burst of radiation which leads to cleaner waveforms. Finally, we study the case of orbiting binary black-hole systems with spin magnitude $χ_i=0.95$ in an aligned configuration and compare waveform and final remnant results with those of the SXS Collaboration, finding excellent agreement. This represents the first moving punctures evolution of orbiting and spinning black holes exceeding the Bowen-York limit. Finally, we study different choices of the initial lapse and lapse evolution equation in the moving punctures approach to improve the accuracy and efficiency of the simulations.
△ Less
Submitted 11 January, 2017; v1 submitted 30 October, 2014;
originally announced October 2014.
-
Resolving the relative influence of strong field spacetime dynamics and MHD on circumbinary disk physics
Authors:
Miguel Zilhão,
Scott C. Noble,
Manuela Campanelli,
Yosef Zlochower
Abstract:
In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hy…
▽ More
In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hydrodynamic evolutions between the two spacetimes down to separations of approximately $40GM/c^2$ where $M$ is the total mass of the binary. For smaller separations, the differences due to PN-order become comparable to differences caused by using initial data further from equilibrium. For magnetized gas, MHD stresses, which drives the accretion dynamics, tends to mask all higher order PN effects even at separations of $20GM/c^2$, leading to essentially the same observed electromagnetic luminosity. This implies that our calculations of the EM signal may be robust down to small binary separations. Our investigation is the first to demonstrate how the level of PN accuracy affects a circumbinary disk's evolution and informs us of the range in separation within which to trust the PN approximation for this kind of study. We also address the influence the initial conditions and binary separation have on simulation predictions.
△ Less
Submitted 31 January, 2015; v1 submitted 16 September, 2014;
originally announced September 2014.
-
Remnant mass, spin, and recoil from spin aligned black-hole binaries
Authors:
James Healy,
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We perform a set of 36 nonprecessing black-hole binary simulations with spins either aligned or counteraligned with the orbital angular momentum in order to model the final mass, spin, and recoil of the merged black hole as a function of the individual black hole spin magnitudes and the mass ratio of the progenitors. We find that the maximum recoil for these configurations is…
▽ More
We perform a set of 36 nonprecessing black-hole binary simulations with spins either aligned or counteraligned with the orbital angular momentum in order to model the final mass, spin, and recoil of the merged black hole as a function of the individual black hole spin magnitudes and the mass ratio of the progenitors. We find that the maximum recoil for these configurations is $V_{max}=526\pm23\,km/s$, which occurs when the progenitor spins are maximal, the mass ratio is $q_{max}=m_1/m_2=0.623\pm0.038$, the smaller black-hole spin is aligned with the orbital angular momentum, and the larger black-hole spin is counteraligned ($α_1=-α_2=1$). This maximum recoil is about $80\,km/s$ larger than previous estimates, but most importantly, because the maximum occurs for smaller mass ratios, the probability for a merging binary to recoil faster than $400\,km/s$ can be as large as $17\%$, while the probability for recoils faster than $250\, km/s$ can be as large as $45\%$. We provide explicit phenomenological formulas for the final mass, spin, and recoil as a function of the individual BH spins and the mass difference between the two black holes. Here we include terms up through fourth-order in the initial spins and mass difference, and find excellent agreement (within a few percent) with independent results available in the literature. The maximum radiated energy is $E_{\rm rad}/m\approx11.3\%$ and final spin $α_{\rm rem}^{\rm max}\approx0.952$ for equal mass, aligned maximally spinning binaries.
△ Less
Submitted 27 June, 2014;
originally announced June 2014.
-
The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the NINJA-2 Collaboration,
:,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
A. Ain,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca
, et al. (867 additional authors not shown)
Abstract:
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of…
▽ More
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]
△ Less
Submitted 5 January, 2014;
originally announced January 2014.
-
Approximate black hole binary spacetime via asymptotic matching
Authors:
Bruno C. Mundim,
Hiroyuki Nakano,
Nicolás Yunes,
Manuela Campanelli,
Scott C. Noble,
Yosef Zlochower
Abstract:
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-…
▽ More
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.
△ Less
Submitted 10 May, 2014; v1 submitted 23 December, 2013;
originally announced December 2013.
-
Black hole binary remnant mass and spin: A new phenomenological formula
Authors:
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We perform a set of 38 numerical simulations of equal-mass BH binaries in a configuration where the BH spins in the binary are equal in both magnitude and direction, to study precession effects. We vary the initial direction of the total spin S with respect to the orbital angular momentum L, covering the 2 dimensional space of orientation angles with 38 configurations consisting of 36 configuratio…
▽ More
We perform a set of 38 numerical simulations of equal-mass BH binaries in a configuration where the BH spins in the binary are equal in both magnitude and direction, to study precession effects. We vary the initial direction of the total spin S with respect to the orbital angular momentum L, covering the 2 dimensional space of orientation angles with 38 configurations consisting of 36 configurations distributed in the azimuthal angle phi and polar angle theta, and two configurations on the poles. In all cases, we set the initial dimensionless BH spins to 0.8. We observe that during the late-inspiral stage, the total angular momentum of the system J remains within 5 deg of its original direction, with the largest changes in direction occurring when the spins are nearly counter-aligned with the orbital angular momentum. We also observe that the angle between S and L is nearly conserved during the inspiral phase. These two dynamical properties allow us to propose a new phenomenological formula for the final mass and spin of merged BHs in terms of the individual masses and spins of the progenitor binary at far separations. We determine coefficients of this formula (in the equal-mass limit) using a least-squares fit to the results of this new set of 38 runs, an additional set of five new configurations with spins aligned/counteraligned with the orbital angular momentum, and over 100 recent simulations. We find that our formulas reproduce the remnant mass and spin of these simulations to within a relative error of 2.5%. We discuss the region of validity of this dynamical picture for precessing unequal-mass binaries. Finally, we perform a statistical study to see the consequence of this new formula for distributions of spin-magnitudes and remnant masses with applications to BH-spin distributions and gravitational radiation in cosmological scenarios involving several mergers.
△ Less
Submitted 21 May, 2014; v1 submitted 19 December, 2013;
originally announced December 2013.
-
Where angular momentum goes in a precessing black hole binary
Authors:
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We evolve a set of 32 equal-mass black-hole binaries with collinear spins (with intrinsic spin magnitudes $|\vec{S}_{1,2}/m^2_{1,2}|=0.8$) to study the effects of precession in the highly nonlinear plunge and merger regimes. We compare the direction of the instantaneous radiated angular momentum, $\hat{δJ}_{\rm rad}(t)$, to the directions of the total angular momentum, $\hat{J}(t)$, and the orbita…
▽ More
We evolve a set of 32 equal-mass black-hole binaries with collinear spins (with intrinsic spin magnitudes $|\vec{S}_{1,2}/m^2_{1,2}|=0.8$) to study the effects of precession in the highly nonlinear plunge and merger regimes. We compare the direction of the instantaneous radiated angular momentum, $\hat{δJ}_{\rm rad}(t)$, to the directions of the total angular momentum, $\hat{J}(t)$, and the orbital angular momentum, $\hat{L}(t)$. We find that $\hat{δJ}_{\rm rad}(t)$ approximately follows $\hat{L}$ throughout the evolution. During the orbital evolution and merger, we observe that the angle between $\vec{L}$ and total spin $\vec{S}$ is approximately conserved to within $1^\circ$, which allows us to propose and test models for the merger remnant's mass and spin. For instance, we verify that the \hangup effect is the dominant effect and largely explains the observed total energy and angular momentum radiated by these precessing systems. We also verify that the total angular momentum, which significantly decreases in magnitude during the inspiral, varies in direction by less than $\sim 5^\circ$. The maximum variation in the direction of $\vec J$ occurs when the spins are nearly antialigned with the orbital angular momentum. Based on our results, we conjecture that transitional precession, which would lead to large variations in the direction of $\vec J$, is not possible for similar-mass binaries and would require a mass ratio $m_1/m_2\lesssim1/4$.
△ Less
Submitted 12 January, 2014; v1 submitted 23 July, 2013;
originally announced July 2013.
-
Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration
Authors:
Ian Hinder,
Alessandra Buonanno,
Michael Boyle,
Zachariah B. Etienne,
James Healy,
Nathan K. Johnson-McDaniel,
Alessandro Nagar,
Hiroyuki Nakano,
Yi Pan,
Harald P. Pfeiffer,
Michael Pürrer,
Christian Reisswig,
Mark A. Scheel,
Erik Schnetter,
Ulrich Sperhake,
Bela Szilágyi,
Wolfgang Tichy,
Barry Wardell,
Anıl Zenginoglu,
Daniela Alic,
Sebastiano Bernuzzi,
Tanja Bode,
Bernd Brügmann,
Luisa T. Buchman,
Manuela Campanelli
, et al. (31 additional authors not shown)
Abstract:
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 i…
▽ More
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.
△ Less
Submitted 11 December, 2013; v1 submitted 19 July, 2013;
originally announced July 2013.
-
Exploring the Outer Limits of Numerical Relativity
Authors:
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We perform several black-hole binary evolutions using fully nonlinear numerical relativity techniques at separations large enough that low-order post-Newtonian expansions are expected to be accurate. As a case study, we evolve an equal-mass nonspinning black-hole binary from a quasicircular orbit at an initial coordinate separation of D=100M for three different resolutions. We find that the orbita…
▽ More
We perform several black-hole binary evolutions using fully nonlinear numerical relativity techniques at separations large enough that low-order post-Newtonian expansions are expected to be accurate. As a case study, we evolve an equal-mass nonspinning black-hole binary from a quasicircular orbit at an initial coordinate separation of D=100M for three different resolutions. We find that the orbital period of this binary (in the numerical coordinates) is T=6422M. The orbital motion agrees with post-Newtonian predictions to within 1%. Interestingly, we find that the time derivative of the coordinate separation is dominated by a purely gauge effect leading to an apparent contraction and expansion of the orbit at twice the orbital frequency. Based on these results, we improved our evolution techniques and studied a set of black hole binaries in quasi-circular orbits starting at D=20M, D=50M, and D=100M for ~ 5, 3, and 2 orbits, respectively. We find good agreement between the numerical results and post-Newtonian predictions for the orbital frequency and radial decay rate, radiated energy and angular momentum, and waveform amplitude and phases. The results are relevant for the future computation of long-term waveforms to assist in the detection and analysis of gravitational waves by the next generation of detectors as well as the long-term simulations of black-hole binaries required to accurately model astrophysically realistic circumbinary accretion disks.
△ Less
Submitted 20 May, 2013; v1 submitted 14 April, 2013;
originally announced April 2013.
-
Nonlinear Gravitational Recoil from the Mergers of Precessing Black-Hole Binaries
Authors:
Carlos O. Lousto,
Yosef Zlochower
Abstract:
We present results from an extensive study of 88 precessing, equal-mass black-hole binaries with large spins (83 with intrinsic spins of 0.8 and 5 with intrinsic spins of 0.9)and use these data to model new nonlinear contributions to the gravitational recoil imparted to the merged black hole. We find a new effect, the cross kick, that enhances the recoil for partially aligned binaries beyond the h…
▽ More
We present results from an extensive study of 88 precessing, equal-mass black-hole binaries with large spins (83 with intrinsic spins of 0.8 and 5 with intrinsic spins of 0.9)and use these data to model new nonlinear contributions to the gravitational recoil imparted to the merged black hole. We find a new effect, the cross kick, that enhances the recoil for partially aligned binaries beyond the hangup kick effect. This has the consequence of increasing the probabilities of recoils larger than 2000 km/s by nearly a factor two, and, consequently, of black holes getting ejected from galaxies, as well as the observation of large differential redshifts/blueshifts in the cores of recently merged galaxies.
△ Less
Submitted 12 January, 2014; v1 submitted 29 November, 2012;
originally announced November 2012.
-
Accuracy Issues for Numerical Waveforms
Authors:
Yosef Zlochower,
Marcelo Ponce,
Carlos O. Lousto
Abstract:
We study the convergence properties of our implementation of the 'moving punctures' approach at very high resolutions for an equal-mass, non-spinning, black-hole binary. We find convergence of the Hamiltonian constraint on the horizons and the L2 norm of the Hamiltonian constraint in the bulk for sixth and eighth-order finite difference implementations.
The momentum constraint is more sensitive,…
▽ More
We study the convergence properties of our implementation of the 'moving punctures' approach at very high resolutions for an equal-mass, non-spinning, black-hole binary. We find convergence of the Hamiltonian constraint on the horizons and the L2 norm of the Hamiltonian constraint in the bulk for sixth and eighth-order finite difference implementations.
The momentum constraint is more sensitive, and its L2 norm shows clear convergence for a system with consistent sixth-order finite differencing, while the momentum and BSSN constraints on the horizons show convergence for both sixth and eighth-order systems. We analyze the gravitational waveform error from the late-inspiral, merger, and ringdown. We find that using several lower-order techniques for increasing the speed of numerical relativity simulations actually lead to apparently non-convergent errors. Even when using standard high-accuracy techniques, rather than seeing clean convergence, where the waveform phase is a monotonic function of grid resolution, we find that the phase tends to oscillate with resolution, possibly due to stochastic errors induced by grid refinement boundaries. Our results seem to indicate that one can obtain gravitational waveform phases to within 0.05 rad. (and possibly as small as 0.015 rad.), while the amplitude error can be reduced to 0.1%. We then compare with the waveforms obtained using the cZ4 formalism. We find that the cZ4 waveforms have larger truncation errors for a given resolution, but the Richardson extrapolation phase of the cZ4 and BSSN waveforms agree to within 0.01 rad., even during the ringdown.
△ Less
Submitted 19 November, 2012; v1 submitted 27 August, 2012;
originally announced August 2012.
-
Circumbinary MHD Accretion into Inspiraling Binary Black Holes
Authors:
Scott C. Noble,
Bruno C. Mundim,
Hiroyuki Nakano,
Julian H. Krolik,
Manuela Campanelli,
Yosef Zlochower,
Nicolás Yunes
Abstract:
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 evolut…
▽ More
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.
△ Less
Submitted 4 April, 2012;
originally announced April 2012.
