Showing 1–50 of 183 results for author: Scheel, M A

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.

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Authors: Francois Foucart, Patrick Chi-Kit Cheong, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Binary neutron star mergers play an important role in nuclear astrophysics: their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnet… ▽ More Binary neutron star mergers play an important role in nuclear astrophysics: their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals, and the relative abundances of the produced r-process elements. Existing merger simulations are however limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte-Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant ($(5-10)\,{\rm ms}$). We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates, transport algorithm, as well as recent improvements by Radice et al to the two-moment algorithms can all contribute to changes at the $(10-30)\%$ level in the global properties of the merger remnant and outflows. The limitations of previous moment schemes fixed by Radice et al also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte-Carlo and moment schemes in the context of a low mass binary neutron star system. △ Less

Submitted 22 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.

Comments: 31 pages, 8 figures, comments welcome!

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.

Comments: 20 pages, 8 figures. Submitted to CGQ's focus issue: Gravitational-Wave Memory Effects: From Theory to Observation

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.

Comments: 24 pages, 11 figures, 5 tables

Dynamical ejecta from binary neutron star mergers: Impact of residual eccentricity and equation of state implementation

Authors: Francois Foucart, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Our ability to reliably provide such predictions is however… ▽ More Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a $1.3M_\odot-1.4M_\odot$ binary neutron star system. We find that a residual eccentricity $e\sim 0.01$, as measured $\sim 4-6$ orbits before merger, causes $O(25\%-30\%)$ changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that $O(1\%)$ residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures, to be submitted to PRD

Imprints of Changing Mass and Spin on Black Hole Ringdown

Authors: Hengrui Zhu, Frans Pretorius, Sizheng Ma, Robert Owen, Yitian Chen, Nils Deppe, Lawrence E. Kidder, Maria Okounkova, Harald P. Pfeiffer, Mark A. Scheel, Leo C. Stein

Abstract: We numerically investigate the imprints of gravitational radiation-reaction driven changes to a black hole's mass and spin on the corresponding ringdown waveform. We do so by comparing the dynamics of a perturbed black hole evolved with the full (nonlinear) versus linearized Einstein equations. As expected, we find that the quasinormal mode amplitudes extracted from nonlinear evolution deviate fro… ▽ More We numerically investigate the imprints of gravitational radiation-reaction driven changes to a black hole's mass and spin on the corresponding ringdown waveform. We do so by comparing the dynamics of a perturbed black hole evolved with the full (nonlinear) versus linearized Einstein equations. As expected, we find that the quasinormal mode amplitudes extracted from nonlinear evolution deviate from their linear counterparts at third order in initial perturbation amplitude. For perturbations leading to a change in the black hole mass and spin of $\sim 5\%$, which is reasonable for a remnant formed in an astrophysical merger, we find that nonlinear distortions to the complex amplitudes of some quasinormal modes can be as large as $\sim 50\%$ at the peak of the waveform. Furthermore, the change in the mass and spin results in a drift in the quasinormal mode frequencies, which for large amplitude perturbations causes the nonlinear waveform to rapidly dephase with respect to its linear counterpart. %These two nonlinear effects together create a large distortion in both the amplitude and phase of the ringdown gravitational waveform. Surprisingly, despite these nonlinear effects creating significant deviations in the nonlinear waveform, we show that a linear quasinormal mode model still performs quite well from close to the peak amplitude onwards. △ Less

Submitted 18 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.

Comments: 21 pages, 22 figures

Scalarization of isolated black holes in scalar Gauss-Bonnet theory in the fixing-the-equations approach

Authors: Guillermo Lara, Harald P. Pfeiffer, Nikolas A. Wittek, Nils L. Vu, Kyle C. Nelli, Alexander Carpenter, Geoffrey Lovelace, Mark A. Scheel, William Throwe

Abstract: One of the most promising avenues to perform numerical evolutions in theories beyond General Relativity is the fixing-the-equations approach, a proposal in which new ``driver'' equations are added to the evolution equations in a way that allows for stable numerical evolutions. In this direction, we extend the numerical relativity code SpECTRE to evolve a ``fixed'' version of scalar Gauss-Bonnet th… ▽ More One of the most promising avenues to perform numerical evolutions in theories beyond General Relativity is the fixing-the-equations approach, a proposal in which new ``driver'' equations are added to the evolution equations in a way that allows for stable numerical evolutions. In this direction, we extend the numerical relativity code SpECTRE to evolve a ``fixed'' version of scalar Gauss-Bonnet theory in the decoupling limit, a phenomenologically interesting theory that allows for hairy black hole solutions in vacuum. We focus on isolated black hole systems both with and without linear and angular momentum, and propose a new driver equation to improve the recovery of such stationary solutions. We demonstrate the effectiveness of the latter by numerically evolving black holes that undergo spontaneous scalarization using different driver equations. Finally, we evaluate the accuracy of the obtained solutions by comparing with the original unaltered theory. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 16 pages, 12 figures

Stability of hypermassive neutron stars with realistic rotation and entropy profiles

Authors: Nishad Muhammed, Matthew D. Duez, Pavan Chawhan, Noora Ghadiri, Luisa T. Buchman, Francois Foucart, Patrick Chi-Kit Cheong, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Binary neutron star mergers produce massive, hot, rapidly differentially rotating neutron star remnants; electromagnetic and gravitational wave signals associated with the subsequent evolution depend on the stability of these remnants. Stability of relativistic stars has previously been studied for uniform rotation and for a class of differential rotation with monotonic angular velocity profiles.… ▽ More Binary neutron star mergers produce massive, hot, rapidly differentially rotating neutron star remnants; electromagnetic and gravitational wave signals associated with the subsequent evolution depend on the stability of these remnants. Stability of relativistic stars has previously been studied for uniform rotation and for a class of differential rotation with monotonic angular velocity profiles. Stability of those equilibria to axisymmetric perturbations was found to respect a turning point criterion: along a constant angular momentum sequence, the onset of unstable stars is found at maximum density less than but close to the density of maximum mass. In this paper, we test this turning point criterion for non-monotonic angular velocity profiles and non-isentropic entropy profiles, both chosen to more realistically model post-merger equilibria. Stability is assessed by evolving perturbed equilibria in 2D using the Spectral Einstein Code. We present tests of the code's new capability for axisymmetric metric evolution. We confirm the turning point theorem and determine the region of our rotation law parameter space that provides highest maximum mass for a given angular momentum. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 12 pages, 9 figures

High angular momentum hot differentially rotating equilibrium star evolutions in conformally flat spacetime

Authors: Patrick Chi-Kit Cheong, Nishad Muhammed, Pavan Chawhan, Matthew D. Duez, Francois Foucart, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: The conformal flatness approximation to the Einstein equations has been successfully used in many astrophysical applications such as initial data constructions and dynamical simulations. Although it has been shown that full general relativistic strongly differentially rotating equilibrium models deviate by at most a few percent from their conformally flat counterparts, whether those conformally fl… ▽ More The conformal flatness approximation to the Einstein equations has been successfully used in many astrophysical applications such as initial data constructions and dynamical simulations. Although it has been shown that full general relativistic strongly differentially rotating equilibrium models deviate by at most a few percent from their conformally flat counterparts, whether those conformally flat solutions remain stable has not been fully addressed. To further understand the limitations of the conformal flatness approximation, in this work, we construct spatially-conformally-flat hot hypermassive neutron stars with post-merger-like rotation laws, and perform conformally flat evolutions and analysis over dynamical timescales. We find that enforcing conformally-flat spacetime could change the equilibrium of quasi-toroidal models with high angular momentum for $J \gtrsim 9 ~G M_{\odot}^2 / c$ compared to fully general relativistic cases. In contrast, all the quasi-spherical models considered in this work remain stable even with high angular momentum $J=9~G M_{\odot}^2 / c$. Our investigation suggests that the quasi-spherical models are suitable initial data for long-lived hypermassive neutron star modeling in conformally flat spacetime. △ Less

Submitted 22 July, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions

Authors: Luisa T. Buchman, Matthew D. Duez, Marlo Morales, Mark A. Scheel, Tim M. Kostersitz, Andrew M. Evans, Keefe Mitman

Abstract: Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for subdominant modes. A system of higher order absorbing boundary conditions which greatly reduces this problem was introduced in earlier work [arXiv:gr-qc/0608051]. In this… ▽ More Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for subdominant modes. A system of higher order absorbing boundary conditions which greatly reduces this problem was introduced in earlier work [arXiv:gr-qc/0608051]. In this paper, we devise two new implementations of this boundary condition system in the Spectral Einstein Code (SpEC), and test them in both linear multipolar gravitational wave and inspiralling mass ratio 7:1 binary black hole simulations. One of our implementations in particular is shown to be extremely robust and to produce accuracy superior to the standard freezing-Psi_0 boundary condition usually used by SpEC. △ Less

Submitted 17 July, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: 24 pages, 11 figures (Figs. 9 and 11 showing CCE and extrapolation results are new), accepted by Classical and Quantum Gravity

Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling

Authors: Hengrui Zhu, Justin L. Ripley, Frans Pretorius, Sizheng Ma, Keefe Mitman, Robert Owen, Michael Boyle, Yitian Chen, Nils Deppe, Lawrence E. Kidder, Jordan Moxon, Kyle C. Nelli, Harald P. Pfeiffer, Mark A. Scheel, William Throwe, Nils L. Vu

Abstract: We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and un… ▽ More We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to $5\%$) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counter-rotating and co-rotating quasinormal modes as black hole spin increases. △ Less

Submitted 1 January, 2024; originally announced January 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.

Comments: Data Release and Analysis Scripts: https://github.com/HengruiPrinceton/precession_ringdown

Fully relativistic three-dimensional Cauchy-characteristic matching for physical degrees of freedom

Authors: Sizheng Ma, Jordan Moxon, Mark A. Scheel, Kyle C. Nelli, Nils Deppe, Marceline S. Bonilla, Lawrence E. Kidder, Prayush Kumar, Geoffrey Lovelace, William Throwe, Nils L. Vu

Abstract: A fully relativistic three-dimensional Cauchy-characteristic matching (CCM) algorithm is implemented for physical degrees of freedom in a numerical relativity code SpECTRE. The method is free of approximations and can be applied to any physical system. We test the algorithm with various scenarios involving smooth data, including the propagation of Teukolsky waves within a flat background, the pert… ▽ More A fully relativistic three-dimensional Cauchy-characteristic matching (CCM) algorithm is implemented for physical degrees of freedom in a numerical relativity code SpECTRE. The method is free of approximations and can be applied to any physical system. We test the algorithm with various scenarios involving smooth data, including the propagation of Teukolsky waves within a flat background, the perturbation of a Kerr black hole with a Teukolsky wave, and the injection of a gravitational-wave pulse from the characteristic grid. Our investigations reveal no numerical instabilities in the simulations. In addition, the tests indicate that the CCM algorithm effectively directs characteristic information into the inner Cauchy system, yielding higher precision in waveforms and smaller violations of Bondi-gauge constraints, especially when the outer boundary of the Cauchy evolution is at a smaller radius. △ Less

Submitted 11 June, 2024; v1 submitted 20 August, 2023; originally announced August 2023.

Journal ref: Phys. Rev. D 109, 124027 (2024)

Extending black-hole remnant surrogate models to extreme mass ratios

Authors: Matteo Boschini, Davide Gerosa, Vijay Varma, Cristobal Armaza, Michael Boyle, Marceline S. Bonilla, Andrea Ceja, Yitian Chen, Nils Deppe, Matthew Giesler, Lawrence E. Kidder, Prayush Kumar, Guillermo Lara, Oliver Long, Sizheng Ma, Keefe Mitman, Peter James Nee, Harald P. Pfeiffer, Antoni Ramos-Buades, Mark A. Scheel, Nils L. Vu, Jooheon Yoo

Abstract: Numerical-relativity surrogate models for both black-hole merger waveforms and remnants have emerged as important tools in gravitational-wave astronomy. While producing very accurate predictions, their applicability is limited to the region of the parameter space where numerical-relativity simulations are available and computationally feasible. Notably, this excludes extreme mass ratios. We presen… ▽ More Numerical-relativity surrogate models for both black-hole merger waveforms and remnants have emerged as important tools in gravitational-wave astronomy. While producing very accurate predictions, their applicability is limited to the region of the parameter space where numerical-relativity simulations are available and computationally feasible. Notably, this excludes extreme mass ratios. We present a machine-learning approach to extend the validity of existing and future numerical-relativity surrogate models toward the test-particle limit, targeting in particular the mass and spin of post-merger black-hole remnants. Our model is trained on both numerical-relativity simulations at comparable masses and analytical predictions at extreme mass ratios. We extend the gaussian-process-regression model NRSur7dq4Remnant, validate its performance via cross validation, and test its accuracy against additional numerical-relativity runs. Our fit, which we dub NRSur7dq4EmriRemnant, reaches an accuracy that is comparable to or higher than that of existing remnant models while providing robust predictions for arbitrary mass ratios. △ Less

Submitted 24 October, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 10 pages, 3 figures. Published in PRD. Model publicly available at https://pypi.org/project/surfinBH

Journal ref: Phys.Rev.D 108 (2023) 8, 084015

Gravitational Waves from Binary Neutron Star Mergers with a Spectral Equation of State

Authors: Alexander Knight, Francois Foucart, Matthew D. Duez, Mike Boyle, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are $C^2$ or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the fle… ▽ More In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are $C^2$ or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than $Δ\tildeΛ = 55$ can be captured by the SpEC code for these equations of state. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: 13 pages, 8 figures

Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole-Neutron Star Mergers

Authors: Ore Gottlieb, Danat Issa, Jonatan Jacquemin-Ide, Matthew Liska, Francois Foucart, Alexander Tchekhovskoy, Brian D. Metzger, Eliot Quataert, Rosalba Perna, Daniel Kasen, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly a… ▽ More We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as $L_j\sim t^{-2}$. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For post-merger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization $σ_0>100$ retain significant magnetization ($σ\gg1$) at $r>10^{10}\,{\rm cm}$, emphasizing the importance of magnetic processes in the prompt emission. The jet-wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper. △ Less

Submitted 18 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: For movies of the simulations, see https://oregottlieb.com/bhns.html

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.

Comments: 48 pages, 17 figures. Matches published version, minor changes only

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.

Comments: 14 pages, 11 figures. Accepted for publication in PRD

Journal ref: Phys. Rev. D 108, 064027 (2023)

Numerical simulations of black hole-neutron star mergers in scalar-tensor gravity

Authors: Sizheng Ma, Vijay Varma, Leo C. Stein, Francois Foucart, Matthew D. Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We present a numerical-relativity simulation of a black hole - neutron star merger in scalar-tensor (ST) gravity with binary parameters consistent with the gravitational wave event GW200115. In this exploratory simulation, we consider the Damour-Esposito-Farese extension to Brans-Dicke theory, and maximize the effect of spontaneous scalarization by choosing a soft equation of state and ST theory p… ▽ More We present a numerical-relativity simulation of a black hole - neutron star merger in scalar-tensor (ST) gravity with binary parameters consistent with the gravitational wave event GW200115. In this exploratory simulation, we consider the Damour-Esposito-Farese extension to Brans-Dicke theory, and maximize the effect of spontaneous scalarization by choosing a soft equation of state and ST theory parameters at the edge of known constraints. We extrapolate the gravitational waves, including tensor and scalar (breathing) modes, to future null-infinity. The numerical waveforms undergo ~ 22 wave cycles before the merger, and are in good agreement with predictions from post-Newtonian theory during the inspiral. We find the ST system evolves faster than its general-relativity (GR) counterpart due to dipole radiation, merging a full gravitational-wave cycle before the GR counterpart. This enables easy differentiation between the ST waveforms and GR in the context of parameter estimation. However, we find that dipole radiation's effect may be partially degenerate with the NS tidal deformability during the late inspiral stage, and a full Bayesian analysis is necessary to fully understand the degeneracies between ST and binary parameters in GR. △ Less

Submitted 13 June, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Journal ref: Phys. Rev. D 107, 124051 (2023)

Worldtube excision method for intermediate-mass-ratio inspirals: scalar-field model in 3+1 dimensions

Authors: Nikolas A. Wittek, Mekhi Dhesi, Leor Barack, Harald P. Pfeiffer, Adam Pound, Hannes R. Rüter, Marceline S. Bonilla, Nils Deppe, Lawrence E. Kidder, Prayush Kumar, Mark A. Scheel, William Throwe, Nils L. Vu

Abstract: Binary black hole simulations become increasingly more computationally expensive with smaller mass ratios, partly because of the longer evolution time, and partly because the lengthscale disparity dictates smaller time steps. The program initiated by Dhesi et al. (arXiv:2109.03531) explores a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysic… ▽ More Binary black hole simulations become increasingly more computationally expensive with smaller mass ratios, partly because of the longer evolution time, and partly because the lengthscale disparity dictates smaller time steps. The program initiated by Dhesi et al. (arXiv:2109.03531) explores a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysical range ($10^{-4} \lesssim q \lesssim 10^{-2}$), where purely perturbative methods may not be adequate. A region ("worldtube") much larger than the small black hole is excised from the numerical domain, and replaced with an analytical model approximating a tidally deformed black hole. Here we apply this idea to a toy model of a scalar charge in a fixed circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field. This is a first implementation of the worldtube excision method in full 3+1 dimensions. We demonstrate the accuracy and efficiency of the method, and discuss the steps towards applying it for evolving orbits and, ultimately, in the binary black-hole scenario. Our implementation is publicly accessible in the SpECTRE numerical relativity code. △ Less

Submitted 16 February, 2024; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: 19 pages, 10 figures; remove typo in Eq. (9c)

Journal ref: Phys. Rev. D 108, 024041 (2023)

Laying the foundation of the effective-one-body waveform models SEOBNRv5: improved accuracy and efficiency for spinning non-precessing binary black holes

Authors: Lorenzo Pompili, Alessandra Buonanno, Héctor Estellés, Mohammed Khalil, Maarten van de Meent, Deyan P. Mihaylov, Serguei Ossokine, Michael Pürrer, Antoni Ramos-Buades, Ajit Kumar Mehta, Roberto Cotesta, Sylvain Marsat, Michael Boyle, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel, Hannes R. Rüter, Nils Vu, Reetika Dudi, Sizheng Ma, Keefe Mitman, Denyz Melchor, Sierra Thomas, Jennifer Sanchez

Abstract: We present SEOBNRv5HM, a more accurate and faster inspiral-merger-ringdown gravitational waveform model for quasi-circular, spinning, nonprecessing binary black holes within the effective-one-body (EOB) formalism. Compared to its predecessor, SEOBNRv4HM, the waveform model i) incorporates recent high-order post- Newtonian results in the inspiral, with improved resummations, ii) includes the gravit… ▽ More We present SEOBNRv5HM, a more accurate and faster inspiral-merger-ringdown gravitational waveform model for quasi-circular, spinning, nonprecessing binary black holes within the effective-one-body (EOB) formalism. Compared to its predecessor, SEOBNRv4HM, the waveform model i) incorporates recent high-order post- Newtonian results in the inspiral, with improved resummations, ii) includes the gravitational modes (l, |m|) = (3, 2), (4, 3), in addition to the (2, 2), (3, 3), (2, 1), (4, 4), (5, 5) modes already implemented in SEOBNRv4HM, iii) is calibrated to larger mass-ratios and spins using a catalog of 442 numerical-relativity (NR) simulations and 13 additional waveforms from black-hole perturbation theory, iv) incorporates information from second-order gravitational self-force (2GSF) in the nonspinning modes and radiation-reaction force. Computing the unfaithfulness against NR simulations, we find that for the dominant (2, 2) mode the maximum unfaithfulness in the total mass range $10-300 M_{\odot}$ is below $10^{-3}$ for 90% of the cases (38% for SEOBNRv4HM). When including all modes up to l = 5 we find 98% (49%) of the cases with unfaithfulness below $10^{-2} (10^{-3})$, while these numbers reduce to 88% (5%) when using SEOBNRv4HM. Furthermore, the model shows improved agreement with NR in other dynamical quantities (e.g., the angular momentum flux and binding energy), providing a powerful check of its physical robustness. We implemented the waveform model in a high-performance Python package (pySEOBNR), which leads to evaluation times faster than SEOBNRv4HM by a factor 10 to 50, depending on the configuration, and provides the flexibility to easily include spin-precession and eccentric effects, thus making it the starting point for a new generation of EOBNR waveform models (SEOBNRv5) to be employed for upcoming observing runs of the LIGO-Virgo-KAGRA detectors. △ Less

Submitted 31 March, 2023; originally announced March 2023.

General relativistic simulations of collapsing binary neutron star mergers with Monte-Carlo neutrino transport

Authors: Francois Foucart, Matthew D. Duez, Roland Haas, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel, Elizabeth Spira-Savett

Abstract: Recent gravitational wave observations of neutron star-neutron star and neutron star-black hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. In this manuscript, we present a first set of evolution of massive neutron star binaries using Monte-Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symme… ▽ More Recent gravitational wave observations of neutron star-neutron star and neutron star-black hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. In this manuscript, we present a first set of evolution of massive neutron star binaries using Monte-Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting post-merger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the WhiskyTHC code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the post-merger remnant. One of the main potential use of our simulations is thus as improved initial conditions for longer evolutions of such remnants. △ Less

Submitted 12 April, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 21p, 10 figures, under review by PRD

Eccentric binary black holes: Comparing numerical relativity and small mass-ratio perturbation theory

Authors: Antoni Ramos-Buades, Maarten van de Meent, Harald P. Pfeiffer, Hannes R. Rüter, Mark A. Scheel, Michael Boyle, Lawrence E. Kidder

Abstract: The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a functi… ▽ More The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a function of mass ratio for eccentric non-spinning binary black holes. We produce $52$ NR simulations with mass ratios between $1:10$ and $1:1$ and initial eccentricities up to $0.7$. From these we extract quantities like gravitational wave energy and angular momentum fluxes and periastron advance, and assess their accuracy. To facilitate comparison, we develop tools to map between NR and SMR inspiral evolutions of eccentric binary black holes. We derive post-Newtonian accurate relations between different definitions of eccentricity. Based on these analyses, we introduce a new definition of eccentricity based on the (2,2)-mode of the gravitational radiation, which reduces to the Newtonian definition of eccentricity in the Newtonian limit. From the comparison between NR simulations and SMR results, we quantify the unknown next-to-leading order SMR contributions to the gravitational energy and angular momentum fluxes, and periastron advance. We show that in the comparable mass regime these contributions are subdominant and higher order SMR contributions are negligible. △ Less

Submitted 7 September, 2022; originally announced September 2022.

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.

Comments: 6+2 pages, 4 figures, 1 table. Matches PRL version

Journal ref: Phys. Rev. Lett. 130, 081402 (2023)

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.

Comments: 16+4 pages with 2 appendices, 4 figures, 1 table

Journal ref: Phys. Rev. D 106, 084029 (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.

Comments: 23 pages, 17 figures, 2 tables

Numerical-relativity surrogate modeling with nearly extremal black-hole spins

Authors: Marissa Walker, Vijay Varma, Geoffrey Lovelace, Mark A. Scheel

Abstract: Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost -- especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of Reduced Order Modeling (ROM) has enabled the construction of surrogate models… ▽ More Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost -- especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of Reduced Order Modeling (ROM) has enabled the construction of surrogate models trained on an existing set of NR waveforms. Surrogate models enable the rapid computation of the gravitational waves emitted by BBHs. Typically these models are used for interpolation to compute gravitational waveforms for BBHs with mass ratios and spins within the bounds of the training set. Because simulations with nearly extremal spins are so technically challenging, surrogate models almost always rely on training sets with only moderate spins. In this paper, we explore how well surrogate models can extrapolate to nearly extremal spins when the training set only includes moderate spins. For simplicity, we focus on one-dimensional surrogate models trained on NR simulations of BBHs with equal masses and equal, aligned spins. We assess the performance of the surrogate models at higher spin magnitudes by calculating the mismatches between extrapolated surrogate model waveforms and NR waveforms, by calculating the differences between extrapolated and NR measurements of the remnant black-hole mass, and by testing how the surrogate model improves as the training set extends to higher spins. We find that while extrapolation in this one-dimensional case is viable for current detector sensitivities, surrogate models for next-generation detectors should use training sets that extend to nearly extremal spins. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport effects

Authors: Milad Haddadi, Matthew D. Duez, Francois Foucart, Teresita Ramirez, Rodrigo Fernandez, Alexander L. Knight, Jerred Jesse, Francois Hebert, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first ten seconds after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximat… ▽ More Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first ten seconds after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10,000\,km, so that thermal wind velocities only asymptote beyond this radius and at somewhat higher values than previously reported. △ Less

Submitted 15 March, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 23 pages, 7 figures, version accepted to Classical and Quantum Gravity

Success of the small mass ratio approximation during the final orbits of binary black hole simulations

Authors: Sergi Navarro Albalat, Aaron Zimmerman, Matthew Giesler, Mark A. Scheel

Abstract: Recent studies have shown the surprising effectiveness of the small mass-ratio approximation (SMR) in modeling the relativistic two-body problem even at comparable masses. Up to now this effectiveness has been demonstrated only during inspiral, before the binary transitions into plunge and merger. Here we examine the binding energy of nonspinning binary black hole simulations with mass ratios from… ▽ More Recent studies have shown the surprising effectiveness of the small mass-ratio approximation (SMR) in modeling the relativistic two-body problem even at comparable masses. Up to now this effectiveness has been demonstrated only during inspiral, before the binary transitions into plunge and merger. Here we examine the binding energy of nonspinning binary black hole simulations with mass ratios from 20:1 to equal mass. We show for the first time that the binaries undergo a transition to plunge as predicted by analytic theory, and estimate the size of the transition region, which is $\sim 10$ gravitational wave cycles for equal mass binaries. By including transition, the SMR expansion of the binding energy is accurate until the last cycle of gravitational wave emission. This is true even for comparable mass binaries such as those observed by current gravitational wave detectors, where the transition often makes up much of the observed signal. Our work provides further evidence that the SMR approximation can be directly applied to current gravitational wave observations. △ Less

Submitted 30 May, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 12 pages (5 pages plus appendix and references), 8 figures, 1 table. Reflects published version

Journal ref: Phys. Rev. D 107, 084021 (2023)

Surrogate model for gravitational wave signals from non-spinning, comparable- to large-mass-ratio black hole binaries built on black hole perturbation theory waveforms calibrated to numerical relativity

Authors: Tousif Islam, Scott E. Field, Scott A. Hughes, Gaurav Khanna, Vijay Varma, Matthew Giesler, Mark A. Scheel, Lawrence E. Kidder, Harald P. Pfeiffer

Abstract: We present a reduced-order surrogate model of gravitational waveforms from non-spinning binary black hole systems with comparable to large mass-ratio configurations. This surrogate model, \texttt{BHPTNRSur1dq1e4}, is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) with mass ratios varying from 2.5 to 10,000. \texttt{BHPTNRSur1dq1e4} extends an earlier w… ▽ More We present a reduced-order surrogate model of gravitational waveforms from non-spinning binary black hole systems with comparable to large mass-ratio configurations. This surrogate model, \texttt{BHPTNRSur1dq1e4}, is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) with mass ratios varying from 2.5 to 10,000. \texttt{BHPTNRSur1dq1e4} extends an earlier waveform model, \texttt{EMRISur1dq1e4}, by using an updated transition-to-plunge model, covering longer durations up to 30,500 $m_1$ (where $m_1$ is the mass of the primary black hole), includes several more spherical harmonic modes up to $\ell=10$, and calibrates subdominant modes to numerical relativity (NR) data. In the comparable mass-ratio regime, including mass ratios as low as $2.5$, the gravitational waveforms generated through ppBHPT agree surprisingly well with those from NR after this simple calibration step. We also compare our model to recent SXS and RIT NR simulations at mass ratios ranging from $15$ to $32$, and find the dominant quadrupolar modes agree to better than $\approx 10^{-3}$. We expect our model to be useful to study intermediate-mass-ratio binary systems in current and future gravitational-wave detectors. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 20 pages, 15 figures

Targeted large mass ratio numerical relativity surrogate waveform model for GW190814

Authors: Jooheon Yoo, Vijay Varma, Matthew Giesler, Mark A. Scheel, Carl-Johan Haster, Harald P. Pfeiffer, Lawrence E. Kidder, Michael Boyle

Abstract: Gravitational wave observations of large mass ratio compact binary mergers like GW190814 highlight the need for reliable, high-accuracy waveform templates for such systems. We present NRHybSur2dq15, a new surrogate model trained on hybridized numerical relativity (NR) waveforms with mass ratios $q\leq15$, and aligned spins $|χ_{1z}|\leq0.5$ and $χ_{2z}=0$. We target the parameter space of GW190814… ▽ More Gravitational wave observations of large mass ratio compact binary mergers like GW190814 highlight the need for reliable, high-accuracy waveform templates for such systems. We present NRHybSur2dq15, a new surrogate model trained on hybridized numerical relativity (NR) waveforms with mass ratios $q\leq15$, and aligned spins $|χ_{1z}|\leq0.5$ and $χ_{2z}=0$. We target the parameter space of GW190814-like events as large mass ratio NR simulations are very expensive. The model includes the (2,2), (2,1), (3,3), (4,4), and (5,5) spin-weighted spherical harmonic modes, and spans the entire LIGO bandwidth (with $f_{\mathrm{low}}=20$ Hz) for total masses $M \gtrsim 9.5 \, M_{\odot}$. NRHybSur2dq15 accurately reproduces the hybrid waveforms, with mismatches below $\sim 2 \times 10^{-3}$ for total masses $10 \, M_{\odot} \leq M \leq 300 \, M_{\odot}$. This is at least an order of magnitude improvement over existing semi-analytical models for GW190814-like systems. Finally, we reanalyze GW190814 with the new model and obtain source parameter constraints consistent with previous work. △ Less

Submitted 15 July, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 14 pages, 9 figures; Accepted for publication in PRD

Redshift factor and the small mass-ratio limit in binary black hole simulations

Authors: Sergi Navarro Albalat, Aaron Zimmerman, Matthew Giesler, Mark A. Scheel

Abstract: We present a calculation of the Detweiler redshift factor in binary black hole simulations based on its relation to the surface gravity. The redshift factor has far-reaching applications in analytic approximations, gravitational self-force calculations, and conservative two-body dynamics. By specializing to non-spinning, quasi-circular binaries with mass ratios ranging from $m_A/m_B = 1$ to… ▽ More We present a calculation of the Detweiler redshift factor in binary black hole simulations based on its relation to the surface gravity. The redshift factor has far-reaching applications in analytic approximations, gravitational self-force calculations, and conservative two-body dynamics. By specializing to non-spinning, quasi-circular binaries with mass ratios ranging from $m_A/m_B = 1$ to $m_A/m_B = 9.5$ we are able to recover the leading small-mass-ratio (SMR) prediction with relative differences of order $10^{-5}$ from simulations alone. The next-to-leading order term that we extract agrees with the SMR prediction arising from self-force calculations, with differences of a few percent. These deviations from the first-order conservative prediction are consistent with non-adiabatic effects that can be accommodated in an SMR expansion. This fact is also supported by a comparison to the conservative post-Newtonian prediction of the redshifts. For the individual redshifts, a re-expansion in terms of the symmetric mass ratio $ν$ does not improve the convergence of the series. However we find that when looking at the sum of the redshift factors of both back holes, $z_A + z_B$, which is symmetric under the exchange of the masses, a re-expansion in $ν$ accelerates its convergence. Our work provides further evidence of the surprising effectiveness of SMR approximations in modeling even comparable mass binary black holes. △ Less

Submitted 22 April, 2022; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 22 pages, 23 figures

Gravitational-wave echoes from numerical-relativity waveforms via space-time construction near merging compact objects

Authors: Sizheng Ma, Qingwen Wang, Nils Deppe, François Hébert, Lawrence E. Kidder, Jordan Moxon, William Throwe, Nils L. Vu, Mark A. Scheel, Yanbei Chen

Abstract: We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbat… ▽ More We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbations of the final Kerr black hole, and a nonlinear region. At late times, the boundary between the two regions is an infalling shell. The BHP region contains late-time gravitational-waves emitted toward the future horizon, as well as those emitted toward future null infinity. In this region, by imposing no-ingoing wave conditions at past null infinity, and matching out-going waves at future null infinity with waveforms computed from numerical relativity, we can obtain waves that travel toward the future horizon. In particular, the Newman-Penrose $ψ_0$ associated with the in-going wave on the horizon is related to tidal deformations measured by fiducial observers floating above the horizon. We further determine the boundary of the BHP region on the future horizon by imposing that $ψ_0$ inside the BHP region can be faithfully represented by quasi-normal modes. Using a physically-motivated way to impose boundary conditions near the horizon, and applying the so-called Boltzmann reflectivity, we compute the quasi-normal modes of non-rotating ECOs, as well as gravitational-wave echoes. We also investigate the detectability of these echoes in current and future detectors, and prospects for parameter estimation. △ Less

Submitted 21 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Journal ref: Phys. Rev. D 105, 104007 (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.

Comments: 25 pages, 20 figures, published version. Results are reproducible with the ancillary input files

Journal ref: Phys. Rev. D 105, 084027 (2022)

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.

Comments: 17 + 2 pages, 11 figures, 2 tables

Journal ref: Phys. Rev. D 105, 104015 (2022)

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.

Comments: matches published version. 21 pages, 12 figures. Added KH instability, TOV runs with classical limiters. arXiv admin note: text overlap with arXiv:2109.11645

Universal features of gravitational waves emitted by superkick binary black hole systems

Authors: Sizheng Ma, Matthew Giesler, Vijay Varma, Mark A. Scheel, Yanbei Chen

Abstract: We use numerical relativity to study the merger and ringdown stages of "superkick" binary black hole systems (those with equal mass and anti-parallel spins). We find a universal way to describe the mass and current quadrupole gravitational waves emitted by these systems during the merger and ringdown stage: (i) The time evolutions of these waves are insensitive to the progenitor's parameters (spin… ▽ More We use numerical relativity to study the merger and ringdown stages of "superkick" binary black hole systems (those with equal mass and anti-parallel spins). We find a universal way to describe the mass and current quadrupole gravitational waves emitted by these systems during the merger and ringdown stage: (i) The time evolutions of these waves are insensitive to the progenitor's parameters (spins) after being normalized by their own peak values. (ii) The peak values, which encode all the spin information of the progenitor, can be consistently fitted to formulas inspired by post-Newtonian theory. We find that the universal evolution of the mass quadrupole wave can be accurately modeled by the so-called Backwards One-Body (BOB) model. However, the BOB model, in its present form, leads to a lower waveform match and a significant parameter-estimation bias for the current quadrupole wave. We also decompose the ringdown signal into seven overtones, and study the dependence of mode amplitudes on the progenitor's parameters. Such dependence is found to be insensitive to the overtone index (up to a scaling factor). Finally, we use the Fisher matrix technique to investigate how the ringdown waveform can be at least as important for parameter estimation as the inspiral stage. Assuming the Cosmic Explorer, we find the contribution of ringdown portion dominates as the total mass exceeds ~ 250 solar mass. For massive BBH systems, the accuracy of parameter measurement is improved by incorporating the information of ringdown -- the ringdown sector gives rise to a different parameter correlation from inspiral stage, hence the overall parameter correlation is reduced in the full signal. △ Less

Submitted 10 July, 2021; originally announced July 2021.

Journal ref: Phys. Rev. D 104, 084003 (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.

Comments: 18 pages, 11 figures; Published in Physical Review D

Journal ref: Phys. Rev. D 104, 024051 (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.

Comments: 14 pages, 4 figures, 1 table; published Physical Review D

Journal ref: Phys. Rev. D 103, 124029 (2021)

Implementation of Monte-Carlo transport in the general relativistic SpEC code

Authors: Francois Foucart, Matthew D. Duez, Francois Hebert, Lawrence E. Kidder, Phillip Kovarik, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers, and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of… ▽ More Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers, and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable, but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte-Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the latter is a subdominant source of error at the accuracy reached by current simulations, and for the interactions currently included in our code. We also provide tests of the most important features of this code. △ Less

Submitted 23 July, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: 32p, 1 table and 14 figures; Accepted by ApJ

Journal ref: ApJ 920 82 (2021)

Eccentric binary black hole surrogate models for the gravitational waveform and remnant properties: comparable mass, nonspinning case

Authors: Tousif Islam, Vijay Varma, Jackie Lodman, Scott E. Field, Gaurav Khanna, Mark A. Scheel, Harald P. Pfeiffer, Davide Gerosa, Lawrence E. Kidder

Abstract: We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, \texttt{NRSur2dq1Ecc}, using 47 nonspinning, equal-mass waveforms with eccentricities up to $0.2$ when measured at a reference time… ▽ More We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, \texttt{NRSur2dq1Ecc}, using 47 nonspinning, equal-mass waveforms with eccentricities up to $0.2$ when measured at a reference time of $5500M$ before merger. This is the first waveform model that is directly trained on eccentric numerical relativity simulations and does not require that the binary circularizes before merger. The model includes the $(2,2)$, $(3,2)$, and $(4,4)$ spin-weighted spherical harmonic modes. We also build a final black hole model, \texttt{NRSur2dq1EccRemnant}, which models the mass, and spin of the remnant black hole. We show that our waveform model can accurately predict numerical relativity waveforms with mismatches $\approx 10^{-3}$, while the remnant model can recover the final mass and dimensionless spin with absolute errors smaller than $\approx 5 \times 10^{-4}M$ and $\approx 2 \times10^{-3}$ respectively. We demonstrate that the waveform model can also recover subtle effects like mode-mixing in the ringdown signal without any special ad-hoc modeling steps. Finally, we show that despite being trained only on equal-mass binaries, \texttt{NRSur2dq1Ecc} can be reasonably extended up to mass ratio $q\approx3$ with mismatches $\simeq 10^{-2}$ for eccentricities smaller than $\sim 0.05$ as measured at a reference time of $2000M$ before merger. The methods developed here should prove useful in the building of future eccentric surrogate models over larger regions of the parameter space. △ Less

Submitted 24 March, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 19 pages, 15 figures

Journal ref: Phys. Rev. D 103, 064022 (2021)

Up-down instability of binary black holes in numerical relativity

Authors: Vijay Varma, Matthew Mould, Davide Gerosa, Mark A. Scheel, Lawrence E. Kidder, Harald P. Pfeiffer

Abstract: Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechani… ▽ More Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechanism for the formation of precessing binaries in environments where sources are preferentially formed with (anti) aligned spins. In this paper, we present the first full numerical relativity simulations capturing this instability. These simulations span $\sim 100$ orbits and $\sim 3$-$5$ precession cycles before merger, making them some of the longest numerical relativity simulations to date. Initialized with a small perturbation of $1^{\circ}$-$10^{\circ}$, the instability causes a dramatic growth of the spin misalignments, which can reach $\sim 90^{\circ}$ near merger. We show that this leaves a strong imprint on the subdominant modes of the gravitational wave signal, which can potentially be used to distinguish up-down binaries from other sources. Finally, we show that post-Newtonian and effective-one-body approximants are able to reproduce the unstable dynamics of up-down binaries extracted from numerical relativity. △ Less

Submitted 1 March, 2021; v1 submitted 13 December, 2020; originally announced December 2020.

Comments: Matches PRD version. Animations available at https://davidegerosa.com/spinprecession

Journal ref: Phys. Rev. D 103, 064003 (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.

Comments: 3 pages, 3 figures (journal version)

Report number: LIGO-P2000507

Journal ref: Phys. Rev. Lett. 127, 011103 (2021)

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.

Comments: 13 pages, 9 figures; Corrected a bug affecting the mismatch plot in Figure 7

Journal ref: Phys. Rev. D 103, 024031 (2021)

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.

Comments: 18 pages, 6 figures, 1 table, published in Physical Review D

Journal ref: Phys. Rev. D 103, 024039 (2021)

High-accuracy waveforms for black hole-neutron star systems with spinning black holes

Authors: Francois Foucart, Alexander Chernoglazov, Michael Boyle, Tanja Hinderer, Max Miller, Jordan Moxon, Mark A. Scheel, Nils Deppe, Matthew D. Duez, Francois Hebert, Lawrence E. Kidder, William Throwe, Harald P. Pfeiffer

Abstract: The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-pr… ▽ More The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-precessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to non-spinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy Characteristic Extraction. The new waveforms have $<0.1\,{\rm rad}$ phase errors during inspiral, rising to $\sim (0.2-0.4)\,{\rm rad}$ errors at merger, and $\lesssim 1\%$ error in their amplitude. We compute the faithfulness of recent analytical models to these numerical results, and find that models specifically designed for BHNS binaries perform well ($F>0.99$) for binaries seen face-on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 16p, 14 figs, 2 tables; To be submitted to PRD

Journal ref: Phys. Rev. D 103, 064007 (2021)

Monte-Carlo neutrino transport in neutron star merger simulations

Authors: Francois Foucart, Matthew D. Duez, Francois Hebert, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An importan… ▽ More Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An important limitation of current simulations is the use of approximate methods for neutrino transport that do not converge to a solution of the transport equations as numerical resolution increases, and thus have errors that are impossible to quantify. Here, we report on a first simulation of a binary neutron star merger that uses Monte-Carlo techniques to directly solve the transport equations in low-density regions. In high-density regions, we use approximations inspired by implicit Monte-Carlo to greatly reduce the cost of simulations, while only introducing errors quantifiable through more expensive convergence studies. We simulate an unequal mass neutron star binary merger up to $5\,{\rm ms}$ past merger, and report on the properties of the matter and neutrino outflows. Finally, we compare our results to the output of our best approximate `M1' transport scheme, demonstrating that an M1 scheme that carefully approximates the neutrino energy spectrum only leads to $\sim 10\%$ uncertainty in the composition and velocity of the ejecta, and $\sim20\%$ uncertainty in the $ν_e$ and $\barν_e$ luminosities and energies. The most significant disagreement found between M1 and Monte-Carlo results is a factor of $\sim 2$ difference in the luminosity of heavy-lepton neutrinos. △ Less

Submitted 22 October, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: 9p, 4 figures, 1 table; version accepted by ApJL

Journal ref: ApJL 902 L27 (2020)