-
Automated determination of the end time of junk radiation in binary black hole simulations
Authors:
Isabella G. Pretto,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
When numerically solving Einstein's equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ''junk radiation.'' This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produ…
▽ More
When numerically solving Einstein's equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ''junk radiation.'' This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produced. Thus, an automated procedure for identifying junk radiation is required. To this end, we present a new algorithm based on the empirical mode decomposition (EMD) from the Hilbert-Huang transform. This approach allows us to isolate and measure the high-frequency oscillations present in the measured irreducible masses of the black holes. The decay of these oscillations allows us to estimate the time from which the junk radiation can be ignored. To make this procedure more precise, we propose three distinct threshold criteria that specify how small the contribution of junk radiation has to be before it can be considered negligible. We apply this algorithm to 3403 BBH simulations from the SXS catalog to find appropriate values for the thresholds in the three criteria. We find that this approach yields reliable decay time estimates, i.e., when to consider the simulation physical, for over 98.6% of the simulations studied. This demonstrates the efficacy of the EMD as a suitable tool to automatically isolate and characterize junk radiation in the simulation of binary black hole systems.
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method
Authors:
Nils Deppe,
Francois Foucart,
Marceline S. Bonilla,
Michael Boyle,
Nicholas J. Corso,
Matthew D. Duez,
Matthew Giesler,
François Hébert,
Lawrence E. Kidder,
Yoonsoo Kim,
Prayush Kumar,
Isaac Legred,
Geoffrey Lovelace,
Elias R. Most,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybr…
▽ More
We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
A Review of Gravitational Memory and BMS Frame Fixing in Numerical Relativity
Authors:
Keefe Mitman,
Michael Boyle,
Leo C. Stein,
Nils Deppe,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
Gravitational memory effects and the BMS freedoms exhibited at future null infinity have recently been resolved and utilized in numerical relativity simulations. With this, gravitational wave models and our understanding of the fundamental nature of general relativity have been vastly improved. In this paper, we review the history and intuition behind memory effects and BMS symmetries, how they ma…
▽ More
Gravitational memory effects and the BMS freedoms exhibited at future null infinity have recently been resolved and utilized in numerical relativity simulations. With this, gravitational wave models and our understanding of the fundamental nature of general relativity have been vastly improved. In this paper, we review the history and intuition behind memory effects and BMS symmetries, how they manifest in gravitational waves, and how controlling the infinite number of BMS freedoms of numerical relativity simulations can crucially improve the waveform models that are used by gravitational wave detectors. We reiterate the fact that, with memory effects and BMS symmetries, not only can these next-generation numerical waveforms be used to observe never-before-seen physics, but they can also be used to test GR and learn new astrophysical information about our universe.
△ Less
Submitted 14 May, 2024;
originally announced May 2024.
-
Improved frequency spectra of gravitational waves with memory in a binary-black-hole simulation
Authors:
Yitian Chen,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
Keefe Mitman,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
William Throwe,
Nils L. Vu,
Saul A. Teukolsky
Abstract:
Numerical relativists can now produce gravitational waveforms with memory effects routinely and accurately. The gravitational-wave memory effect contains very low-frequency components, including a persistent offset. The presence of these components violates basic assumptions about time-shift behavior underpinning standard data-analysis techniques in gravitational-wave astronomy. This poses a chall…
▽ More
Numerical relativists can now produce gravitational waveforms with memory effects routinely and accurately. The gravitational-wave memory effect contains very low-frequency components, including a persistent offset. The presence of these components violates basic assumptions about time-shift behavior underpinning standard data-analysis techniques in gravitational-wave astronomy. This poses a challenge to the analysis of waveform spectra: How to preserve the low-frequency characteristics when transforming a time-domain waveform to the frequency domain. To tackle this challenge, we revisit the preprocessing procedures applied to the waveforms that contain memory effects. We find inconsistency between the zero-frequency limit of displacement memory and the low- frequency spectrum of the same memory preprocessed using the common scheme in literature. To resolve the inconsistency, we propose a new robust preprocessing scheme that produces the spectra of memory waveforms more faithfully. Using this new scheme, we inspect several characteristics of the spectrum of a memory waveform. In particular, we find a discernible beating pattern formed by the dominant oscillatory mode and the displacement memory. This pattern is absent in the spectrum of a waveform without memory. The difference between the memory and no-memory waveforms is too small to be observed by current-generation detectors in a single binary-black-hole event. Detecting the memory in a single event is likely to occur in the era of next-generation detectors.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method
Authors:
Yoonsoo Kim,
Elias R. Most,
William Throwe,
Saul A. Teukolsky,
Nils Deppe
Abstract:
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical approach for simulating such force-free plasmas, which allows for high accuracy in smooth regions as well as capturing dissipation in current sheets. Using a high…
▽ More
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical approach for simulating such force-free plasmas, which allows for high accuracy in smooth regions as well as capturing dissipation in current sheets. Using a high-order accurate discontinuous Galerkin method augmented with a conservative finite-difference method, we demonstrate efficient global simulations of black hole and neutron star magnetospheres. In addition to a series of challenging test problems, we show that our approach can-depending on the physical properties of the system and the numerical implementation-be up to 10x more efficient than conventional simulations, with a speedup of 2-3x for most problems we consider in practice.
△ Less
Submitted 12 June, 2024; v1 submitted 1 April, 2024;
originally announced April 2024.
-
Optimizing post-Newtonian parameters and fixing the BMS frame for numerical-relativity waveform hybridizations
Authors:
Dongze Sun,
Michael Boyle,
Keefe Mitman,
Mark A. Scheel,
Leo C. Stein,
Saul A. Teukolsky,
Vijay Varma
Abstract:
Numerical relativity (NR) simulations of binary black holes provide precise waveforms, but are typically too computationally expensive to produce waveforms with enough orbits to cover the whole frequency band of gravitational-wave observatories. Accordingly, it is important to be able to hybridize NR waveforms with analytic, post-Newtonian (PN) waveforms, which are accurate during the early inspir…
▽ More
Numerical relativity (NR) simulations of binary black holes provide precise waveforms, but are typically too computationally expensive to produce waveforms with enough orbits to cover the whole frequency band of gravitational-wave observatories. Accordingly, it is important to be able to hybridize NR waveforms with analytic, post-Newtonian (PN) waveforms, which are accurate during the early inspiral phase. We show that to build such hybrids, it is crucial to both fix the Bondi-Metzner-Sachs (BMS) frame of the NR waveforms to match that of PN theory, and optimize over the PN parameters. We test such a hybridization procedure including all spin-weighted spherical harmonic modes with $|m|\leq \ell$ for $\ell\leq 8$, using 29 NR waveforms with mass ratios $q\leq 10$ and spin magnitudes $|χ_1|, |χ_2|\leq 0.8$. We find that for spin-aligned systems, the PN and NR waveforms agree very well. The difference is limited by the small nonzero orbital eccentricity of the NR waveforms, or equivalently by the lack of eccentric terms in the PN waveforms. To maintain full accuracy of the simulations, the matching window for spin-aligned systems should be at least 5 orbits long and end at least 15 orbits before merger. For precessing systems, the errors are larger than for spin-aligned cases. The errors are likely limited by the absence of precession-related spin-spin PN terms. Using $10^5\,M$ long NR waveforms, we find that there is no optimal choice of the matching window within this time span, because the hybridization result for precessing cases is always better if using earlier or longer matching windows. We provide the mean orbital frequency of the smallest acceptable matching window as a function of the target error between the PN and NR waveforms and the black hole spins.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Black Hole Spectroscopy for Precessing Binary Black Hole Coalescences
Authors:
Hengrui Zhu,
Harrison Siegel,
Keefe Mitman,
Maximiliano Isi,
Will M. Farr,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
Sizheng Ma,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Vijay Varma,
Nils L. Vu
Abstract:
The spectroscopic study of black hole quasinormal modes in gravitational-wave ringdown observations is hindered by our ignorance of which modes should dominate astrophysical signals for different binary configurations, limiting tests of general relativity and astrophysics. In this work, we present a description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-ci…
▽ More
The spectroscopic study of black hole quasinormal modes in gravitational-wave ringdown observations is hindered by our ignorance of which modes should dominate astrophysical signals for different binary configurations, limiting tests of general relativity and astrophysics. In this work, we present a description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences -- a key region of parameter space that has yet to be fully explored within the framework of black hole spectroscopy. We suggest that the remnant perturbation for precessing and non-precessing systems is approximately the same up to a rotation, which implies that the relative amplitudes of the quasinormal modes in both systems are also related by a rotation. We present evidence for this by analyzing an extensive catalog of numerical relativity simulations. Additional structure in the amplitudes is connected to the system's kick velocity and other asymmetries in the orbital dynamics. We find that the ringdowns of precessing systems need not be dominated by the ${(\ell,m)=(2,\pm 2)}$ quasinormal modes, and that instead the $(2,\pm 1)$~or~$(2,0)$ quasinormal modes can dominate. Our results are consistent with a ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521, and may also help in understanding phenomenological inspiral-merger-ringdown waveform model systematics.
△ Less
Submitted 18 July, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
A positivity-preserving adaptive-order finite-difference scheme for GRMHD
Authors:
Nils Deppe,
Lawrence E. Kidder,
Saul A. Teukolsky,
Marceline S. Bonilla,
François Hébert,
Yoonsoo Kim,
Mark A. Scheel,
William Throwe,
Nils L. Vu
Abstract:
We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is achieved by monitoring the relative power in the highest mode of the reconstructed polynomial and reducing the order when the polynomial series no longer converges. Ou…
▽ More
We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is achieved by monitoring the relative power in the highest mode of the reconstructed polynomial and reducing the order when the polynomial series no longer converges. Our approach is similar to the multidimensional optimal order detection (MOOD) strategy, but differs in several ways. The approach is a priori and so does not require retaking a time step. It can also readily be combined with positivity-preserving flux limiters that have gained significant traction in computational astrophysics and numerical relativity. This combination ultimately guarantees a physical solution both during reconstruction and time stepping. We demonstrate the capabilities of the method using a standard suite of very challenging 1d, 2d, and 3d general relativistic magnetohydrodynamics test problems.
△ Less
Submitted 18 January, 2024; v1 submitted 7 June, 2023;
originally announced June 2023.
-
Numerical relativity surrogate model with memory effects and post-Newtonian hybridization
Authors:
Jooheon Yoo,
Keefe Mitman,
Vijay Varma,
Michael Boyle,
Scott E. Field,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Leo C. Stein,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
Numerical relativity simulations provide the most precise templates for the gravitational waves produced by binary black hole mergers. However, many of these simulations use an incomplete waveform extraction technique -- extrapolation -- that fails to capture important physics, such as gravitational memory effects. Cauchy-characteristic evolution (CCE), by contrast, is a much more physically accur…
▽ More
Numerical relativity simulations provide the most precise templates for the gravitational waves produced by binary black hole mergers. However, many of these simulations use an incomplete waveform extraction technique -- extrapolation -- that fails to capture important physics, such as gravitational memory effects. Cauchy-characteristic evolution (CCE), by contrast, is a much more physically accurate extraction procedure that fully evolves Einstein's equations to future null infinity and accurately captures the expected physics. In this work, we present a new surrogate model, NRHybSur3dq8$\_$CCE, built from CCE waveforms that have been mapped to the post-Newtonian (PN) BMS frame and then hybridized with PN and effective one-body (EOB) waveforms. This model is trained on 102 waveforms with mass ratios $q\leq8$ and aligned spins $χ_{1z}, \, χ_{2z} \in \left[-0.8, 0.8\right]$. The model spans the entire LIGO-Virgo-KAGRA (LVK) frequency band (with $f_{\text{low}}=20\text{Hz}$) for total masses $M\gtrsim2.25M_{\odot}$ and includes the $\ell\leq4$ and $(\ell,m)=(5,5)$ spin-weight $-2$ spherical harmonic modes, but not the $(3,1)$, $(4,2)$ or $(4,1)$ modes. We find that NRHybSur3dq8$\_$CCE can accurately reproduce the training waveforms with mismatches $\lesssim2\times10^{-4}$ for total masses $2.25M_{\odot}\leq M\leq300M_{\odot}$ and can, for a modest degree of extrapolation, capably model outside of its training region. Most importantly, unlike previous waveform models, the new surrogate model successfully captures memory effects.
△ Less
Submitted 14 September, 2023; v1 submitted 5 June, 2023;
originally announced June 2023.
-
Nonlinearities in Black Hole Ringdowns
Authors:
Keefe Mitman,
Macarena Lagos,
Leo C. Stein,
Sizheng Ma,
Lam Hui,
Yanbei Chen,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary…
▽ More
The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic $(4,4)$ mode's amplitude exhibits quadratic scaling with the fundamental $(2,2)$ mode -- its parent mode. The nonlinear mode's amplitude is comparable to or even larger than that of the linear $(4,4)$ mode. Therefore, correctly modeling the ringdown of higher harmonics -- improving mode mismatches by up to 2 orders of magnitude -- requires the inclusion of nonlinear effects.
△ Less
Submitted 22 February, 2023; v1 submitted 15 August, 2022;
originally announced August 2022.
-
Fixing the BMS frame of numerical relativity waveforms with BMS charges
Authors:
Keefe Mitman,
Leo C. Stein,
Michael Boyle,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results th…
▽ More
The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results that are in different BMS frames. Consequently, to build a model for the waveforms produced during the merging of compact objects, which ideally would be a hybridization of PN, NR, and black hole perturbation theory, one needs a fast and robust method for fixing the BMS freedoms. In this work, we present the first means of fixing the entire BMS freedom of NR waveforms to match the frame of either PN waveforms or black hole perturbation theory. We achieve this by finding the BMS transformations that change certain charges in a prescribed way -- e.g., finding the center-of-mass transformation that maps the center-of-mass charge to a mean of zero. We find that this new method is 20 times faster, and more correct when mapping to the superrest frame, than previous methods that relied on optimization algorithms. Furthermore, in the course of developing this charge-based frame fixing method, we compute the PN expression for the Moreschi supermomentum to 3PN order without spins and 2PN order with spins. This Moreschi supermomentum is effectively equivalent to the energy flux or the null memory contribution at future null infinity $\mathscr{I}^{+}$. From this PN calculation, we also compute oscillatory ($m\not=0$ modes) and spin-dependent memory terms that have not been identified previously or have been missing from strain expressions in the post-Newtonian literature.
△ Less
Submitted 21 October, 2022; v1 submitted 8 August, 2022;
originally announced August 2022.
-
Multipole moments on the common horizon in a binary-black-hole simulation
Authors:
Yitian Chen,
Prayush Kumar,
Neev Khera,
Nils Deppe,
Arnab Dhani,
Michael Boyle,
Matthew Giesler,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We construct the covariantly defined multipole moments on the common horizon of an equal-mass, non-spinning, quasicircular binary-black-hole system. We see a strong correlation between these multipole moments and the gravitational waveform. We find that the multipole moments are well described by the fundamental quasinormal modes at sufficiently late times. For each multipole moment, at least two…
▽ More
We construct the covariantly defined multipole moments on the common horizon of an equal-mass, non-spinning, quasicircular binary-black-hole system. We see a strong correlation between these multipole moments and the gravitational waveform. We find that the multipole moments are well described by the fundamental quasinormal modes at sufficiently late times. For each multipole moment, at least two fundamental modes of different $\ell$ are detectable in the best model. These models provide faithful estimates of the true mass and spin of the remnant black hole. We also show that by including overtones, the $\ell=m=2$ mass multipole moment admits an excellent quasinormal-mode description at all times after the merger. This demonstrates the perhaps surprising power of perturbation theory near the merger.
△ Less
Submitted 4 August, 2022;
originally announced August 2022.
-
Modeling Compact Binary Merger Waveforms Beyond General Relativity
Authors:
Gabriel S. Bonilla,
Prayush Kumar,
Saul A. Teukolsky
Abstract:
The parameterized post-Einsteinian framework modifies inspiral waveform models to incorporate effects beyond General Relativity. We extend the existing model into the merger-ringdown regime. The modification introduced here adds a single degree of freedom that corresponds to a change in the binary coalescence time. Other merger properties remain as predicted by GR. We discuss parameter estimation…
▽ More
The parameterized post-Einsteinian framework modifies inspiral waveform models to incorporate effects beyond General Relativity. We extend the existing model into the merger-ringdown regime. The modification introduced here adds a single degree of freedom that corresponds to a change in the binary coalescence time. Other merger properties remain as predicted by GR. We discuss parameter estimation with this model, and how it can be used to extract information from beyond-GR waveforms.
△ Less
Submitted 26 March, 2022;
originally announced March 2022.
-
A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
Authors:
Nils L. Vu,
Harald P. Pfeiffer,
Gabriel S. Bonilla,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Geoffrey Lovelace,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nikolas A. Wittek,
Tom Włodarczyk
Abstract:
Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple days to solve these problems at high resolution and when matter is involved, because they are either hard to parallelize or require a large amount of computationa…
▽ More
Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple days to solve these problems at high resolution and when matter is involved, because they are either hard to parallelize or require a large amount of computational resources. Here we present a new solver for linear and nonlinear elliptic problems that is designed to scale with resolution and to parallelize on computing clusters. To achieve this we employ a discontinuous Galerkin discretization, an iterative multigrid-Schwarz preconditioned Newton-Krylov algorithm, and a task-based parallelism paradigm. To accelerate convergence of the elliptic solver we have developed novel subdomain-preconditioning techniques. We find that our multigrid-Schwarz preconditioned elliptic solves achieve iteration counts that are independent of resolution, and our task-based parallel programs scale over 200 million degrees of freedom to at least a few thousand cores. Our new code solves a classic initial data problem for binary black holes faster than the spectral code SpEC when distributed to only eight cores, and in a fraction of the time on more cores. It is publicly accessible in the next-generation SpECTRE numerical relativity code. Our results pave the way for highly parallel elliptic solves in numerical relativity and beyond.
△ Less
Submitted 18 April, 2022; v1 submitted 12 November, 2021;
originally announced November 2021.
-
High Precision Ringdown Modeling: Multimode Fits and BMS Frames
Authors:
Lorena Magaña Zertuche,
Keefe Mitman,
Neev Khera,
Leo C. Stein,
Michael Boyle,
Nils Deppe,
François Hébert,
Dante A. B. Iozzo,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils Vu
Abstract:
Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically im…
▽ More
Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), an extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim10^5$ when using multimode fitting as opposed to only fitting the $(2,\pm2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms$-$which exhibit the memory effect$-$to a BMS frame known as the super rest frame, there is an improvement of $\sim10^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim4$ better than using the publicly-available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by NR and comparing our results to previous QNM studies.
△ Less
Submitted 11 May, 2022; v1 submitted 29 October, 2021;
originally announced October 2021.
-
The SpECTRE Cauchy-characteristic evolution system for rapid, precise waveform extraction
Authors:
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
Nils Deppe,
Nils Fischer,
Francois Hébert,
Lawrence E. Kidder,
William Throwe
Abstract:
We give full details regarding the new Cauchy-characteristic evolution (CCE) system in SpECTRE. The implementation is built to provide streamlined flexibility for either extracting waveforms during the process of a SpECTRE binary compact object simulation, or as a standalone module for extracting waveforms from worldtube data provided by another code base. Using our recently presented improved ana…
▽ More
We give full details regarding the new Cauchy-characteristic evolution (CCE) system in SpECTRE. The implementation is built to provide streamlined flexibility for either extracting waveforms during the process of a SpECTRE binary compact object simulation, or as a standalone module for extracting waveforms from worldtube data provided by another code base. Using our recently presented improved analytic formulation, the CCE system is free of pure-gauge logarithms that would spoil the spectral convergence of the scheme. It gracefully extracts all five Weyl scalars, in addition to the news and the strain. The SpECTRE CCE system makes significant improvements on previous implementations in modularity, ease of use, and speed of computation.
△ Less
Submitted 16 October, 2021;
originally announced October 2021.
-
Simulating magnetized neutron stars with discontinuous Galerkin methods
Authors:
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
William Throwe,
Isha Anantpurkar,
Cristóbal Armaza,
Gabriel S. Bonilla,
Michael Boyle,
Himanshu Chaudhary,
Matthew D. Duez,
Nils L. Vu,
Francois Foucart,
Matthew Giesler,
Jason S. Guo,
Yoonsoo Kim,
Prayush Kumar,
Isaac Legred,
Dongjun Li,
Geoffrey Lovelace,
Sizheng Ma,
Alexandra Macedo,
Denyz Melchor,
Marlo Morales,
Jordan Moxon,
Kyle C. Nelli
, et al. (11 additional authors not shown)
Abstract:
Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativisti…
▽ More
Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$ΛΠ^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method.
△ Less
Submitted 28 June, 2022; v1 submitted 24 September, 2021;
originally announced September 2021.
-
A high-order shock capturing discontinuous Galerkin-finite-difference hybrid method for GRMHD
Authors:
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited discontinuous Galerkin candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken…
▽ More
We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited discontinuous Galerkin candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust finite-difference methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the discontinuous Galerkin and finite-difference solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1d, 2d, and 3d relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.
△ Less
Submitted 16 January, 2024; v1 submitted 23 September, 2021;
originally announced September 2021.
-
Efficient simulations of high-spin black holes with a new gauge
Authors:
Yitian Chen,
Nils Deppe,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
We present a new choice of initial data for binary black hole simulations that significantly improves the efficiency of high-spin simulations. We use spherical Kerr-Schild coordinates, where the horizon of a rotating black hole is spherical, for each black hole. The superposed spherical Kerr-Schild initial data reduce the runtime by a factor of 2 compared to standard superposed Kerr-Schild for an…
▽ More
We present a new choice of initial data for binary black hole simulations that significantly improves the efficiency of high-spin simulations. We use spherical Kerr-Schild coordinates, where the horizon of a rotating black hole is spherical, for each black hole. The superposed spherical Kerr-Schild initial data reduce the runtime by a factor of 2 compared to standard superposed Kerr-Schild for an intermediate resolution spin-0.99 binary-black-hole simulation. We also explore the possibility of delaying the transition from the initial data gauge to the evolution gauge, which produces an additional speed-up of 1.3.
△ Less
Submitted 14 October, 2021; v1 submitted 4 August, 2021;
originally announced August 2021.
-
Fixing the BMS Frame of Numerical Relativity Waveforms
Authors:
Keefe Mitman,
Neev Khera,
Dante A. B. Iozzo,
Leo C. Stein,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe
Abstract:
Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools…
▽ More
Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE's Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during postprocessing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincaré charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincaré transformations.
△ Less
Submitted 31 January, 2022; v1 submitted 5 May, 2021;
originally announced May 2021.
-
Comparing Remnant Properties from Horizon Data and Asymptotic Data in Numerical Relativity
Authors:
Dante A. B. Iozzo,
Neev Khera,
Leo C. Stein,
Keefe Mitman,
Michael Boyle,
Nils Deppe,
Francois Hebert,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe
Abstract:
We present a new study of remnant black hole properties from 13 binary black hole systems, numerically evolved using the Spectral Einstein Code. The mass, spin, and recoil velocity of each remnant were determined quasi-locally from apparent horizon data and asymptotically from Bondi data $(h, ψ_4, ψ_3, ψ_2, ψ_1)$ computed at future null infinity using SpECTRE's Cauchy characteristic evolution. We…
▽ More
We present a new study of remnant black hole properties from 13 binary black hole systems, numerically evolved using the Spectral Einstein Code. The mass, spin, and recoil velocity of each remnant were determined quasi-locally from apparent horizon data and asymptotically from Bondi data $(h, ψ_4, ψ_3, ψ_2, ψ_1)$ computed at future null infinity using SpECTRE's Cauchy characteristic evolution. We compare these independent measurements of the remnant properties in the bulk and on the boundary of the spacetime, giving insight into how well asymptotic data are able to reproduce local properties of the remnant black hole in numerical relativity. We also discuss the theoretical framework for connecting horizon quantities to asymptotic quantities and how it relates to our results. This study recommends a simple improvement to the recoil velocities reported in the Simulating eXtreme Spacetimes waveform catalog, provides an improvement to future surrogate remnant models, and offers new analysis techniques for evaluating the physical accuracy of numerical simulations.
△ Less
Submitted 15 June, 2021; v1 submitted 14 April, 2021;
originally announced April 2021.
-
Testing the black-hole area law with GW150914
Authors:
Maximiliano Isi,
Will M. Farr,
Matthew Giesler,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We present observational confirmation of Hawking's black-hole area theorem based on data from GW150914, finding agreement with the prediction with 97% (95%) probability when we model the ringdown including (excluding) overtones of the quadrupolar mode. We obtain this result from a new time-domain analysis of the pre- and postmerger data. We also confirm that the inspiral and ringdown portions of t…
▽ More
We present observational confirmation of Hawking's black-hole area theorem based on data from GW150914, finding agreement with the prediction with 97% (95%) probability when we model the ringdown including (excluding) overtones of the quadrupolar mode. We obtain this result from a new time-domain analysis of the pre- and postmerger data. We also confirm that the inspiral and ringdown portions of the signal are consistent with the same remnant mass and spin, in agreement with general relativity.
△ Less
Submitted 9 June, 2021; v1 submitted 8 December, 2020;
originally announced December 2020.
-
Adding Gravitational Memory to Waveform Catalogs using BMS Balance Laws
Authors:
Keefe Mitman,
Dante A. B. Iozzo,
Neev Khera,
Michael Boyle,
Tommaso De Lorenzo,
Nils Deppe,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe
Abstract:
Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration's catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spaceti…
▽ More
Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration's catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Lastly, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package $\texttt{sxs}$, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws.
△ Less
Submitted 18 January, 2021; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Extending Gravitational Wave Extraction Using Weyl Characteristic Fields
Authors:
Dante A. B. Iozzo,
Michael Boyle,
Nils Deppe,
Jordan Moxon,
Mark A. Scheel,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Saul A. Teukolsky
Abstract:
We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in postprocessing. The…
▽ More
We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in postprocessing. These methods have been used to produce $Ψ_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes (SXS) waveform catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy characteristic extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors.
△ Less
Submitted 21 January, 2021; v1 submitted 28 October, 2020;
originally announced October 2020.
-
Computation of Displacement and Spin Gravitational Memory in Numerical Relativity
Authors:
Keefe Mitman,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
William Throwe
Abstract:
We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS Collaboration's $\texttt{SpEC}$ code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using $\texttt{SpECTRE}$'s version of a Cauchy-characteris…
▽ More
We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS Collaboration's $\texttt{SpEC}$ code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using $\texttt{SpECTRE}$'s version of a Cauchy-characteristic extraction. We find that we can accurately resolve the strain's traditional $m=0$ memory modes and some of the $m\not=0$ oscillatory memory modes that have previously only been theorized. We also perform a separate calculation of the memory using equations for the Bondi-Metzner-Sachs charges as well as the energy and angular momentum fluxes at asymptotic infinity. Our new calculation uses only the gravitational wave strain and two of the Weyl scalars at infinity. Also, this computation shows that the memory modes can be understood as a combination of a memory signal throughout the binary's inspiral and merger phases, and a quasinormal mode signal near the ringdown phase. Additionally, we find that the magnetic memory, up to numerical error, is indeed zero as previously conjectured. Lastly, we find that signal-to-noise ratios of memory for LIGO, the Einstein Telescope (ET), and the Laser Interferometer Space Antenna (LISA) with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or Minimal Waveform models.
△ Less
Submitted 26 January, 2021; v1 submitted 22 July, 2020;
originally announced July 2020.
-
Improved Cauchy-characteristic evolution system for high-precision numerical relativity waveforms
Authors:
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at $\mathcal{I}^+$ from numerical relativity simulations. Cauchy-characteristic evolution combines an interior solution of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to $\mathcal{I}^+$. The foundati…
▽ More
We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at $\mathcal{I}^+$ from numerical relativity simulations. Cauchy-characteristic evolution combines an interior solution of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to $\mathcal{I}^+$. The foundation of our improved algorithm is a comprehensive method of handling the gauge transformations between the arbitrarily specified coordinates of the interior Cauchy evolution and the unique (up to BMS transformations) Bondi-Sachs coordinate system of the exterior characteristic evolution. We present a reformulated set of characteristic evolution equations better adapted to numerical implementation. In addition, we develop a method to ensure that the angular coordinates used in the volume during the characteristic evolution are asymptotically inertial. This provides a direct route to an expanded set of waveform outputs and is guaranteed to avoid pure-gauge logarithmic dependence that has caused trouble for previous spectral implementations of the characteristic evolution equations. We construct a set of Weyl scalars compatible with the Bondi-like coordinate systems used in characteristic evolution, and determine simple, easily implemented forms for the asymptotic Weyl scalars in our suggested set of coordinates.
△ Less
Submitted 20 August, 2021; v1 submitted 2 July, 2020;
originally announced July 2020.
-
Numerical relativity simulation of GW150914 beyond general relativity
Authors:
Maria Okounkova,
Leo C. Stein,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the le…
▽ More
We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $\gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis.
△ Less
Submitted 7 May, 2020; v1 submitted 6 November, 2019;
originally announced November 2019.
-
Numerical binary black hole collisions in dynamical Chern-Simons gravity
Authors:
Maria Okounkova,
Leo C. Stein,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational w…
▽ More
We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational wave detection. Head-on collisions are interesting in their own right, however, as they cleanly probe the quasi-normal mode spectrum of the final black hole. We thus compute the leading-order dynamical Chern-Simons modifications to the complex frequencies of the post-merger gravitational radiation. We consider equal-mass systems, with equal spins oriented along the axis of collision, resulting in remnant black holes with spin. We find modifications to the complex frequencies of the quasi-normal mode spectrum that behave as a power law with the spin of the remnant, and that are not degenerate with the frequencies associated with a Kerr black hole of any mass and spin. We discuss these results in the context of testing general relativity with gravitational wave observations.
△ Less
Submitted 15 November, 2019; v1 submitted 20 June, 2019;
originally announced June 2019.
-
Testing the no-hair theorem with GW150914
Authors:
Maximiliano Isi,
Matthew Giesler,
Will M. Farr,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode ($\ell=m=2$), with $3.6σ$ confidence. A ringdown model inc…
▽ More
We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode ($\ell=m=2$), with $3.6σ$ confidence. A ringdown model including overtones allows us to measure the final mass and spin magnitude of the remnant exclusively from postinspiral data, obtaining an estimate in agreement with the values inferred from the full signal. The mass and spin values we measure from the ringdown agree with those obtained using solely the fundamental mode at a later time, but have smaller uncertainties. Agreement between the postinspiral measurements of mass and spin and those using the full waveform supports the hypothesis that the GW150914 merger produced a Kerr black hole, as predicted by general relativity, and provides a test of the no-hair theorem at the ${\sim}10\%$ level. An independent measurement of the frequency of the first overtone yields agreement with the no-hair hypothesis at the ${\sim 20}\%$ level. As the detector sensitivity improves and the detected population of black hole mergers grows, we can expect that using overtones will provide even stronger tests.
△ Less
Submitted 8 August, 2019; v1 submitted 2 May, 2019;
originally announced May 2019.
-
The SXS Collaboration catalog of binary black hole simulations
Authors:
Michael Boyle,
Daniel Hemberger,
Dante A. B. Iozzo,
Geoffrey Lovelace,
Serguei Ossokine,
Harald P. Pfeiffer,
Mark A. Scheel,
Leo C. Stein,
Charles J. Woodford,
Aaron B. Zimmerman,
Nousha Afshari,
Kevin Barkett,
Jonathan Blackman,
Katerina Chatziioannou,
Tony Chu,
Nicholas Demos,
Nils Deppe,
Scott E. Field,
Nils L. Fischer,
Evan Foley,
Heather Fong,
Alyssa Garcia,
Matthew Giesler,
Francois Hebert,
Ian Hinder
, et al. (15 additional authors not shown)
Abstract:
Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of…
▽ More
Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms .
△ Less
Submitted 11 September, 2019; v1 submitted 9 April, 2019;
originally announced April 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.
-
Evolving Metric Perturbations in dynamical Chern-Simons Gravity
Authors:
Maria Okounkova,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We present a well-posed constraint-preserving scheme for evolving first-order metric perturbations on an arbitrary background with arbitrary source. We use this scheme to evolve the leading-order metric perturbation in order-reduced dynamical Chern-Simons gravity (dCS) on a Kerr background. In particular we test the stability of stationary dCS data on a Kerr background with stationary first-order…
▽ More
We present a well-posed constraint-preserving scheme for evolving first-order metric perturbations on an arbitrary background with arbitrary source. We use this scheme to evolve the leading-order metric perturbation in order-reduced dynamical Chern-Simons gravity (dCS) on a Kerr background. In particular we test the stability of stationary dCS data on a Kerr background with stationary first-order dCS scalar field source. We find that the leading-order metric perturbation numerically exhibits linear growth, but that the level of this growth converges to zero with numerical resolution. This analysis shows that spinning black holes in dCS gravity are numerically stable to leading-order perturbations in the metric.
△ Less
Submitted 31 January, 2019; v1 submitted 26 November, 2018;
originally announced November 2018.
-
Numerical black hole initial data and shadows in dynamical Chern-Simons gravity
Authors:
Maria Okounkova,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We present a scheme for generating first-order metric perturbation initial data for an arbitrary background and source. We then apply this scheme to derive metric perturbations in order-reduced dynamical Chern-Simons gravity (dCS). In particular, we solve for metric perturbations on a black hole background that are sourced by a first-order dCS scalar field. This gives us the leading-order metric p…
▽ More
We present a scheme for generating first-order metric perturbation initial data for an arbitrary background and source. We then apply this scheme to derive metric perturbations in order-reduced dynamical Chern-Simons gravity (dCS). In particular, we solve for metric perturbations on a black hole background that are sourced by a first-order dCS scalar field. This gives us the leading-order metric perturbation to the spacetime in dCS gravity. We then use these solutions to compute black hole shadows in the linearly perturbed spacetime by evolving null geodesics. We present a novel scheme to decompose the shape of the shadow into multipoles parametrized by the spin of the background black hole and the perturbation parameter $\varepsilon^2$. We find that we can differentiate the presence of a pure Kerr spacetime from a spacetime with a dCS perturbation using the shadow, allowing in part for a null-hypothesis test of general relativity. We then consider these results in the context of the Event Horizon Telescope.
△ Less
Submitted 10 January, 2019; v1 submitted 11 October, 2018;
originally announced October 2018.
-
Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates
Authors:
Prayush Kumar,
Jonathan Blackman,
Scott E. Field,
Mark Scheel,
Chad R. Galley,
Michael Boyle,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Bela Szilagyi,
Saul A. Teukolsky
Abstract:
Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy bl…
▽ More
Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants.
We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.
△ Less
Submitted 19 April, 2020; v1 submitted 24 August, 2018;
originally announced August 2018.
-
General-relativistic neutron star evolutions with the discontinuous Galerkin method
Authors:
François Hébert,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
Simulations of relativistic hydrodynamics often need both high accuracy and robust shock-handling properties. The discontinuous Galerkin method combines these features --- a high order of convergence in regions where the solution is smooth and shock-capturing properties for regions where it is not --- with geometric flexibility and is therefore well suited to solve the partial differential equatio…
▽ More
Simulations of relativistic hydrodynamics often need both high accuracy and robust shock-handling properties. The discontinuous Galerkin method combines these features --- a high order of convergence in regions where the solution is smooth and shock-capturing properties for regions where it is not --- with geometric flexibility and is therefore well suited to solve the partial differential equations describing astrophysical scenarios. We present here evolutions of a general-relativistic neutron star with the discontinuous Galerkin method. In these simulations, we simultaneously evolve the spacetime geometry and the matter on the same computational grid, which we conform to the spherical geometry of the problem. To verify the correctness of our implementation, we perform standard convergence and shock tests. We then show results for evolving, in three dimensions, a Kerr black hole; a neutron star in the Cowling approximation (holding the spacetime metric fixed); and, finally, a neutron star where the spacetime and matter are both dynamical. The evolutions show long-term stability, good accuracy, and an improved rate of convergence versus a comparable-resolution finite-volume method.
△ Less
Submitted 28 August, 2018; v1 submitted 5 April, 2018;
originally announced April 2018.
-
Critical behavior in 3-d gravitational collapse of massless scalar fields
Authors:
Nils Deppe,
Lawrence E. Kidder,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
We present results from the first study of critical behavior in 3-d gravitational collapse. The source of the gravitational field is a massless scalar field. This is a well-studied case for spherically symmetric gravitational collapse, allowing us to understand the reliability and accuracy of the simulations. We study both supercritical and subcritical evolutions to see if one provides more accura…
▽ More
We present results from the first study of critical behavior in 3-d gravitational collapse. The source of the gravitational field is a massless scalar field. This is a well-studied case for spherically symmetric gravitational collapse, allowing us to understand the reliability and accuracy of the simulations. We study both supercritical and subcritical evolutions to see if one provides more accurate results than the other. We find that even for non-spherical initial data with 35 percent of the power in the $\ell=2$ spherical harmonic, the critical solution is the same as in spherical symmetry.
△ Less
Submitted 18 December, 2018; v1 submitted 23 February, 2018;
originally announced February 2018.
-
On choosing the start time of binary black hole ringdown
Authors:
Swetha Bhagwat,
Maria Okounkova,
Stefan W. Ballmer,
Duncan A. Brown,
Matthew Giesler,
Mark A. Scheel,
Saul A. Teukolsky
Abstract:
The final stage of a binary black hole merger is ringdown, in which the system is described by a Kerr black hole with quasinormal mode perturbations. It is far from straightforward to identify the time at which the ringdown begins. Yet determining this time is important for precision tests of the general theory of relativity that compare an observed signal with quasinormal mode descriptions of the…
▽ More
The final stage of a binary black hole merger is ringdown, in which the system is described by a Kerr black hole with quasinormal mode perturbations. It is far from straightforward to identify the time at which the ringdown begins. Yet determining this time is important for precision tests of the general theory of relativity that compare an observed signal with quasinormal mode descriptions of the ringdown, such as tests of the no-hair theorem. We present an algorithmic method to analyze the choice of ringdown start time in the observed waveform. This method is based on determining how close the strong field is to a Kerr black hole (Kerrness). Using numerical relativity simulations, we characterize the Kerrness of the strong-field region close to the black hole using a set of local, gauge-invariant geometric and algebraic conditions that measure local isometry to Kerr. We produce a map that associates each time in the gravitational waveform with a value of each of these Kerrness measures; this map is produced by following outgoing null characteristics from the strong and near-field regions to the wave zone. We perform this analysis on a numerical relativity simulation with parameters consistent with GW150914- the first gravitational wave detection. We find that the choice of ringdown start time of $3\,\mathrm{ms}$ after merger used in the GW150914 study to test general relativity corresponds to a high dimensionless perturbation amplitude of $ \sim 7.5 \times 10^{-3}$ in the strong-field region. This suggests that in higher signal-to-noise detections, one would need to start analyzing the signal at a later time for studies that depend on the validity of black hole perturbation theory.
△ Less
Submitted 6 November, 2017; v1 submitted 2 November, 2017;
originally announced November 2017.
-
SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
Authors:
Lawrence E. Kidder,
Scott E. Field,
Francois Foucart,
Erik Schnetter,
Saul A. Teukolsky,
Andy Bohn,
Nils Deppe,
Peter Diener,
François Hébert,
Jonas Lippuner,
Jonah Miller,
Christian D. Ott,
Mark A. Scheel,
Trevor Vincent
Abstract:
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resoluti…
▽ More
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads.
△ Less
Submitted 21 July, 2017; v1 submitted 31 August, 2016;
originally announced September 2016.
-
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.
-
A Parallel Adaptive Event Horizon Finder for Numerical Relativity
Authors:
Andy Bohn,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
With Advanced LIGO detecting the gravitational waves emitted from a pair of merging black holes in late 2015, we have a new perspective into the strong field regime of binary black hole systems. Event horizons are the defining features of such black hole spacetimes. We introduce a new code for locating event horizons in numerical simulations based on a Delaunay triangulation on a topological spher…
▽ More
With Advanced LIGO detecting the gravitational waves emitted from a pair of merging black holes in late 2015, we have a new perspective into the strong field regime of binary black hole systems. Event horizons are the defining features of such black hole spacetimes. We introduce a new code for locating event horizons in numerical simulations based on a Delaunay triangulation on a topological sphere. The code can automatically refine arbitrary regions of the event horizon surface to find and explore features such as the hole in a toroidal event horizon, as discussed in our companion paper. We also investigate various ways of integrating the geodesic equation and find evolution equations that can be integrated efficiently with high accuracy.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
Toroidal Horizons in Binary Black Hole Mergers
Authors:
Andy Bohn,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
We find the first binary black hole event horizon with a toroidal topology. It had been predicted that generically the event horizons of merging black holes should briefly have a toroidal topology, but such a phase has never been seen prior to this work. In all previous binary black hole simulations, in the coordinate slicing used to evolve the black holes, the topology of the event horizon transi…
▽ More
We find the first binary black hole event horizon with a toroidal topology. It had been predicted that generically the event horizons of merging black holes should briefly have a toroidal topology, but such a phase has never been seen prior to this work. In all previous binary black hole simulations, in the coordinate slicing used to evolve the black holes, the topology of the event horizon transitions directly from two spheres during the inspiral to a single sphere as the black holes merge. We present a coordinate transformation to a foliation of spacelike hypersurfaces that "cut a hole" through the event horizon surface, resulting in a toroidal event horizon. A torus could potentially provide a mechanism for violating topological censorship. However, these toroidal event horizons satisfy topological censorship by construction, because we can always trivially apply the inverse coordinate transformation to remove the topological feature.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
Simulations of inspiraling and merging double neutron stars using the Spectral Einstein Code
Authors:
Roland Haas,
Christian D. Ott,
Bela Szilagyi,
Jeffrey D. Kaplan,
Jonas Lippuner,
Mark A. Scheel,
Kevin Barkett,
Curran D. Muhlberger,
Tim Dietrich,
Matthew D. Duez,
Francois Foucart,
Harald P. Pfeiffer,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
We present results on the inspiral, merger, and post-merger evolution of a neutron star - neutron star (NSNS) system. Our results are obtained using the hybrid pseudospectral-finite volume Spectral Einstein Code (SpEC). To test our numerical methods, we evolve an equal-mass system for $\approx 22$ orbits before merger. This waveform is the longest waveform obtained from fully general-relativistic…
▽ More
We present results on the inspiral, merger, and post-merger evolution of a neutron star - neutron star (NSNS) system. Our results are obtained using the hybrid pseudospectral-finite volume Spectral Einstein Code (SpEC). To test our numerical methods, we evolve an equal-mass system for $\approx 22$ orbits before merger. This waveform is the longest waveform obtained from fully general-relativistic simulations for NSNSs to date. Such long (and accurate) numerical waveforms are required to further improve semi-analytical models used in gravitational wave data analysis, for example the effective one body models. We discuss in detail the improvements to SpEC's ability to simulate NSNS mergers, in particular mesh refined grids to better resolve the merger and post-merger phases. We provide a set of consistency checks and compare our results to NSNS merger simulations with the independent BAM code. We find agreement between them, which increases confidence in results obtained with either code. This work paves the way for future studies using long waveforms and more complex microphysical descriptions of neutron star matter in SpEC.
△ Less
Submitted 25 June, 2016; v1 submitted 4 April, 2016;
originally announced April 2016.
-
Formulation of discontinuous Galerkin methods for relativistic astrophysics
Authors:
Saul A. Teukolsky
Abstract:
The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to arbitrary time-dependent curved spacetimes. We show how to formulate the algorithm in such spacetimes for applications in relativistic astrophysics. We also sho…
▽ More
The DG algorithm is a powerful method for solving pdes, especially for evolution equations in conservation form. Since the algorithm involves integration over volume elements, it is not immediately obvious that it will generalize easily to arbitrary time-dependent curved spacetimes. We show how to formulate the algorithm in such spacetimes for applications in relativistic astrophysics. We also show how to formulate the algorithm for equations in non-conservative form, such as Einstein's field equations themselves. We find two computationally distinct formulations in both cases, one of which has seldom been used before for flat space in curvilinear coordinates but which may be more efficient. We also give a new derivation of the ALE algorithm (Arbitrary Lagrangian-Eulerian) using 4-vector methods that is much simpler than the usual derivation and explains why the method preserves the conservation form of the equations. The various formulations are explored with some simple numerical experiments that also explore the effect of the metric identities on the results.
△ Less
Submitted 5 October, 2015;
originally announced October 2015.
-
The Kerr Metric
Authors:
Saul A. Teukolsky
Abstract:
This review describes the events leading up to the discovery of the Kerr metric in 1963 and the enormous impact the discovery has had in the subsequent 50 years. The review discusses the Penrose process, the four laws of black hole mechanics, uniqueness of the solution, and the no-hair theorems. It also includes Kerr perturbation theory and its application to black hole stability and quasi-normal…
▽ More
This review describes the events leading up to the discovery of the Kerr metric in 1963 and the enormous impact the discovery has had in the subsequent 50 years. The review discusses the Penrose process, the four laws of black hole mechanics, uniqueness of the solution, and the no-hair theorems. It also includes Kerr perturbation theory and its application to black hole stability and quasi-normal modes. The Kerr metric's importance in the astrophysics of quasars and accreting stellar-mass black hole systems is detailed. A theme of the review is the "miraculous" nature of the solution, both in describing in a simple analytic formula the most general rotating black hole, and in having unexpected mathematical properties that make many calculations tractable. Also included is a pedagogical derivation of the solution suitable for a first course in general relativity.
△ Less
Submitted 15 January, 2015; v1 submitted 8 October, 2014;
originally announced October 2014.
-
Initial data for high-compactness black hole-neutron star binaries
Authors:
Katherine Henriksson,
François Foucart,
Lawrence E. Kidder,
Saul A. Teukolsky
Abstract:
For highly compact neutron stars, constructing numerical initial data for black hole-neutron star binary evolutions is very difficult. We describe improvements to an earlier method that enable it to handle these more challenging cases. We examine the case of a 6:1 mass ratio system in inspiral close to merger, where the star is governed by a polytropic $Γ=2$, an SLy, or an LS220 equation of state.…
▽ More
For highly compact neutron stars, constructing numerical initial data for black hole-neutron star binary evolutions is very difficult. We describe improvements to an earlier method that enable it to handle these more challenging cases. We examine the case of a 6:1 mass ratio system in inspiral close to merger, where the star is governed by a polytropic $Γ=2$, an SLy, or an LS220 equation of state. In particular, we are able to obtain a solution with a realistic LS220 equation of state for a star with compactness 0.26 and mass 1.98 $M_\odot$, which is representative of the highest reliably determined neutron star masses. For the SLy equation of state, we can obtain solutions with a comparable compactness of 0.25, while for a family of polytropic equations of state, we obtain solutions with compactness up to 0.21, the largest compactness that is stable in this family. These compactness values are significantly higher than any previously published results. We find that improvements in adapting the computational domain to the neutron star surface and in accounting for the center of mass drift of the system are the key ingredients allowing us to obtain these solutions.
△ Less
Submitted 25 September, 2014;
originally announced September 2014.
-
Magnetic effects on the low-T/|W| instability in differentially rotating neutron stars
Authors:
Curran D. Muhlberger,
Fatemeh Hossein Nouri,
Matthew D. Duez,
Francois Foucart,
Lawrence E. Kidder,
Christian D. Ott,
Mark A. Scheel,
Béla Szilágyi,
Saul A. Teukolsky
Abstract:
Dynamical instabilities in protoneutron stars may produce gravitational waves whose observation could shed light on the physics of core-collapse supernovae. When born with sufficient differential rotation, these stars are susceptible to a shear instability (the "low-T/|W| instability"), but such rotation can also amplify magnetic fields to strengths where they have a considerable impact on the dyn…
▽ More
Dynamical instabilities in protoneutron stars may produce gravitational waves whose observation could shed light on the physics of core-collapse supernovae. When born with sufficient differential rotation, these stars are susceptible to a shear instability (the "low-T/|W| instability"), but such rotation can also amplify magnetic fields to strengths where they have a considerable impact on the dynamics of the stellar matter. Using a new magnetohydrodynamics module for the Spectral Einstein Code, we have simulated a differentially-rotating neutron star in full 3D to study the effects of magnetic fields on this instability. Though strong toroidal fields were predicted to suppress the low-T/|W| instability, we find that they do so only in a small range of field strengths. Below 4e13 G, poloidal seed fields do not wind up fast enough to have an effect before the instability saturates, while above 5e14 G, magnetic instabilities can actually amplify a global quadrupole mode (this threshold may be even lower in reality, as small-scale magnetic instabilities remain difficult to resolve numerically). Thus, the prospects for observing gravitational waves from such systems are not in fact diminished over most of the magnetic parameter space.
Additionally, we report that the detailed development of the low-T/|W| instability, including its growth rate, depends strongly on the particular numerical methods used. The high-order methods we employ suggest that growth might be considerably slower than found in some previous simulations.
△ Less
Submitted 9 May, 2014;
originally announced May 2014.
-
Periastron advance in spinning black hole binaries: comparing effective-one-body and Numerical Relativity
Authors:
Tanja Hinderer,
Alessandra Buonanno,
Abdul H. Mroué,
Daniel A. Hemberger,
Geoffrey Lovelace,
Harald P. Pfeiffer,
Lawrence E. Kidder,
Mark A. Scheel,
Bela Szilagyi,
Nicholas W. Taylor,
Saul A. Teukolsky
Abstract:
We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasi-circular orbits and having spins collinear with the orbital angular momentum. We compare the predictions with the periastron advance recently computed in accurate numerical-relativity simulations and find remarkable agreement for a wide range of spins and mass ratios. These results do no…
▽ More
We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasi-circular orbits and having spins collinear with the orbital angular momentum. We compare the predictions with the periastron advance recently computed in accurate numerical-relativity simulations and find remarkable agreement for a wide range of spins and mass ratios. These results do not use any numerical-relativity calibration of the effective-one-body model, and stem from two key ingredients in the effective-one-body Hamiltonian: (i) the mapping of the two-body dynamics of spinning particles onto the dynamics of an effective spinning particle in a (deformed) Kerr spacetime, fully symmetrized with respect to the two-body masses and spins, and (ii) the resummation, in the test-particle limit, of all post-Newtonian (PN) corrections linear in the spin of the particle. In fact, even when only the leading spin PN corrections are included in the effective-one-body spinning Hamiltonian but all the test-particle corrections linear in the spin of the particle are resummed we find very good agreement with the numerical results (within the numerical error for equal-mass binaries and discrepancies of at most 1% for larger mass ratios). Furthermore, we specialize to the extreme mass-ratio limit and derive, using the equations of motion in the gravitational skeleton approach, analytical expressions for the periastron advance, the meridional Lense-Thirring precession and spin precession frequency in the case of a spinning particle on a nearly circular equatorial orbit in Kerr spacetime, including also terms quadratic in the spin.
△ Less
Submitted 16 November, 2013; v1 submitted 2 September, 2013;
originally announced September 2013.
-
Periastron Advance in Spinning Black Hole Binaries: Gravitational Self-Force from Numerical Relativity
Authors:
Alexandre Le Tiec,
Alessandra Buonanno,
Abdul H. Mroué,
Harald P. Pfeiffer,
Daniel A. Hemberger,
Geoffrey Lovelace,
Lawrence E. Kidder,
Mark A. Scheel,
Bela Szilágyi,
Nicholas W. Taylor,
Saul A. Teukolsky
Abstract:
We study the general relativistic periastron advance in spinning black hole binaries on quasi-circular orbits, with spins aligned or anti-aligned with the orbital angular momentum, using numerical-relativity simulations, the post-Newtonian approximation, and black hole perturbation theory. By imposing a symmetry by exchange of the bodies' labels, we devise an improved version of the perturbative r…
▽ More
We study the general relativistic periastron advance in spinning black hole binaries on quasi-circular orbits, with spins aligned or anti-aligned with the orbital angular momentum, using numerical-relativity simulations, the post-Newtonian approximation, and black hole perturbation theory. By imposing a symmetry by exchange of the bodies' labels, we devise an improved version of the perturbative result, and use it as the leading term of a new type of expansion in powers of the symmetric mass ratio. This allows us to measure, for the first time, the gravitational self-force effect on the periastron advance of a non-spinning particle orbiting a Kerr black hole of mass M and spin S = -0.5 M^2, down to separations of order 9M. Comparing the predictions of our improved perturbative expansion with the exact results from numerical simulations of equal-mass and equal-spin binaries, we find a remarkable agreement over a wide range of spins and orbital separations.
△ Less
Submitted 29 November, 2013; v1 submitted 2 September, 2013;
originally announced September 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.
-
Final spin and radiated energy in numerical simulations of binary black holes with equal masses and equal, aligned or anti-aligned spins
Authors:
Daniel A. Hemberger,
Geoffrey Lovelace,
Thomas J. Loredo,
Lawrence E. Kidder,
Mark A. Scheel,
Béla Szilágyi,
Nicholas W. Taylor,
Saul A. Teukolsky
Abstract:
The behavior of merging black holes (including the emitted gravitational waves and the properties of the remnant) can currently be computed only by numerical simulations. This paper introduces ten numerical relativity simulations of binary black holes with equal masses and equal spins aligned or anti-aligned with the orbital angular momentum. The initial spin magnitudes have $|χ_i| \lesssim 0.95$…
▽ More
The behavior of merging black holes (including the emitted gravitational waves and the properties of the remnant) can currently be computed only by numerical simulations. This paper introduces ten numerical relativity simulations of binary black holes with equal masses and equal spins aligned or anti-aligned with the orbital angular momentum. The initial spin magnitudes have $|χ_i| \lesssim 0.95$ and are more concentrated in the aligned direction because of the greater astrophysical interest of this case. We combine this data with five previously reported simulations of the same configuration, but with different spin magnitudes, including the highest spin simulated to date, $χ_i \approx 0.97$. This data set is sufficiently accurate to enable us to offer improved analytic fitting formulae for the final spin and for the energy radiated by gravitational waves as a function of initial spin. The improved fitting formulae can help to improve our understanding of the properties of binary black hole merger remnants and can be used to enhance future approximate waveforms for gravitational wave searches, such as Effective-One-Body waveforms.
△ Less
Submitted 23 September, 2013; v1 submitted 26 May, 2013;
originally announced May 2013.