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

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

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 2 May, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

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: To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulation… ▽ More To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulations, we confirm that the relative fundamental quasinormal mode amplitudes of precessing systems are related to those of non-precessing systems by a simple rotation, and that additional structure in the spectrum 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. These results build upon previous works on waveform modeling, and are consistent with a recent ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521. △ Less

Submitted 13 December, 2023; originally announced December 2023.

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

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

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)

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

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

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)

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

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)

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)

A comparison of momentum transport models for numerical relativity

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

Abstract: The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales have mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system… ▽ More The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales have mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on post-merger simulations. We find that correct handling of turbulent heating can be important for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows. △ Less

Submitted 14 December, 2020; v1 submitted 11 August, 2020; originally announced August 2020.

Comments: 15 pagers, 11 figures

Journal ref: Phys. Rev. D 102, 104050 (2020)

Numerical relativity simulation of GW150914 beyond general relativity

Authors: Maria Okounkova, Leo C. Stein, Jordan Moxon, Mark A. Scheel, Saul A. Teukolsky

Abstract: We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the le… ▽ More We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $\gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis. △ Less

Submitted 7 May, 2020; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: 7 pages + appendices, 8 figures, Updated to match Phys. D. Rev article

Journal ref: Phys. Rev. D 101, 104016 (2020)

Smooth equations of state for high-accuracy simulations of neutron star binaries

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

Abstract: High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational cost, however, remains a significant challenge. One issue is that high-order convergence of the solution requires the use of smooth evolution variables, while many… ▽ More High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational cost, however, remains a significant challenge. One issue is that high-order convergence of the solution requires the use of smooth evolution variables, while many of the equations of state used to model the neutron star matter have discontinuities, typically in the first derivative of the pressure. Spectral formulations of the equation of state have been proposed as a potential solution to this problem. Here, we report on the numerical implementation of spectral equations of state in the Spectral Einstein Code. We show that, in our code, spectral equations of state allow for high-accuracy simulations at a lower computational cost than commonly used `piecewise polytrope' equations state. We also demonstrate that not all spectral equations of state are equally useful: different choices for the low-density part of the equation of state can significantly impact the cost and accuracy of simulations. As a result, simulations of neutron star mergers present us with a trade-off between the cost of simulations and the physical realism of the chosen equation of state. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Comments: Submitted to PRD, 13p, 7 figures, 3 tables

Journal ref: Phys. Rev. D 100, 104048 (2019)

Numerical binary black hole collisions in dynamical Chern-Simons gravity

Authors: Maria Okounkova, Leo C. Stein, Mark A. Scheel, Saul A. Teukolsky

Abstract: We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational w… ▽ More We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational wave detection. Head-on collisions are interesting in their own right, however, as they cleanly probe the quasi-normal mode spectrum of the final black hole. We thus compute the leading-order dynamical Chern-Simons modifications to the complex frequencies of the post-merger gravitational radiation. We consider equal-mass systems, with equal spins oriented along the axis of collision, resulting in remnant black holes with spin. We find modifications to the complex frequencies of the quasi-normal mode spectrum that behave as a power law with the spin of the remnant, and that are not degenerate with the frequencies associated with a Kerr black hole of any mass and spin. We discuss these results in the context of testing general relativity with gravitational wave observations. △ Less

Submitted 15 November, 2019; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: 13 pages (+ appendices), 15 figures

Journal ref: Phys. Rev. D 100, 104026 (2019)

Surrogate models for precessing binary black hole simulations with unequal masses

Authors: Vijay Varma, Scott E. Field, Mark A. Scheel, Jonathan Blackman, Davide Gerosa, Leo C. Stein, Lawrence E. Kidder, Harald P. Pfeiffer

Abstract: Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second mod… ▽ More Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second model, \RemnantModel, predicts the properties of the remnant black hole. These models extend previous 7-dimensional, non-eccentric precessing models to higher mass ratios, and have been trained against 1528 simulations with mass ratios $q\leq4$ and spin magnitudes $χ_1,χ_2 \leq 0.8$, with generic spin directions. The waveform model, NRSur7dq4, which begins about 20 orbits before merger, includes all $\ell \leq 4$ spin-weighted spherical harmonic modes, as well as the precession frame dynamics and spin evolution of the black holes. The final black hole model, \RemnantModel, models the mass, spin, and recoil kick velocity of the remnant black hole. In their training parameter range, both models are shown to be more accurate than existing models by at least an order of magnitude, with errors comparable to the estimated errors in the numerical relativity simulations. We also show that the surrogate models work well even when extrapolated outside their training parameter space range, up to mass ratios $q=6$. △ Less

Submitted 10 October, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: Matches published version. Models publicly available at https://zenodo.org/record/3455886#.XZ9s1-dKjBI and https://pypi.org/project/surfinBH

Journal ref: Phys. Rev. Research 1, 033015 (2019)

Testing the no-hair theorem with GW150914

Authors: Maximiliano Isi, Matthew Giesler, Will M. Farr, Mark A. Scheel, Saul A. Teukolsky

Abstract: We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode ($\ell=m=2$), with $3.6σ$ confidence. A ringdown model inc… ▽ More We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode ($\ell=m=2$), with $3.6σ$ confidence. A ringdown model including overtones allows us to measure the final mass and spin magnitude of the remnant exclusively from postinspiral data, obtaining an estimate in agreement with the values inferred from the full signal. The mass and spin values we measure from the ringdown agree with those obtained using solely the fundamental mode at a later time, but have smaller uncertainties. Agreement between the postinspiral measurements of mass and spin and those using the full waveform supports the hypothesis that the GW150914 merger produced a Kerr black hole, as predicted by general relativity, and provides a test of the no-hair theorem at the ${\sim}10\%$ level. An independent measurement of the frequency of the first overtone yields agreement with the no-hair hypothesis at the ${\sim 20}\%$ level. As the detector sensitivity improves and the detected population of black hole mergers grows, we can expect that using overtones will provide even stronger tests. △ Less

Submitted 8 August, 2019; v1 submitted 2 May, 2019; originally announced May 2019.

Comments: v2: journal version

Report number: LIGO-P1900135

Journal ref: Phys. Rev. Lett. 123, 111102 (2019)

The SXS Collaboration catalog of binary black hole simulations

Authors: Michael Boyle, Daniel Hemberger, Dante A. B. Iozzo, Geoffrey Lovelace, Serguei Ossokine, Harald P. Pfeiffer, Mark A. Scheel, Leo C. Stein, Charles J. Woodford, Aaron B. Zimmerman, Nousha Afshari, Kevin Barkett, Jonathan Blackman, Katerina Chatziioannou, Tony Chu, Nicholas Demos, Nils Deppe, Scott E. Field, Nils L. Fischer, Evan Foley, Heather Fong, Alyssa Garcia, Matthew Giesler, Francois Hebert, Ian Hinder , et al. (15 additional authors not shown)

Abstract: Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of… ▽ More Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms . △ Less

Submitted 11 September, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: 33+18 pages, 13 figures, 4 tables, 2,018 binaries. Catalog metadata in ancillary JSON file. v2: Matches version accepted by CQG. Catalog available at https://www.black-holes.org/waveforms

Journal ref: Class. Quantum Grav. 36, 195006 (2019)

Numerical simulations of neutron star-black hole binaries in the near-equal-mass regime

Authors: F. Foucart, M. D. Duez, L. E. Kidder, S. Nissanke, H. P. Pfeiffer, M. A. Scheel

Abstract: Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power… ▽ More Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $χ_{\rm BH}=0.84$ for our equal mass case. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. D 99, 103025 (2019)

On the properties of the massive binary black hole merger GW170729

Authors: Katerina Chatziioannou, Roberto Cotesta, Sudarshan Ghonge, Jacob Lange, Ken K. -Y. Ng, Juan Calderon Bustillo, James Clark, Carl-Johan Haster, Sebastian Khan, Michael Puerrer, Vivien Raymond, Salvatore Vitale, Nousha Afshari, Stanislav Babak, Kevin Barkett, Jonathan Blackman, Alejandro Bohe, Michael Boyle, Alessandra Buonanno, Manuela Campanelli, Gregorio Carullo, Tony Chu, Eric Flynn, Heather Fong, Alyssa Garcia , et al. (31 additional authors not shown)

Abstract: We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusio… ▽ More We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario. △ Less

Submitted 7 November, 2019; v1 submitted 15 March, 2019; originally announced March 2019.

Comments: 14 pages, 8 figures, final published version, samples available at https://git.ligo.org/katerina.chatziioannou/gw170729hm_datarelease

Journal ref: Phys. Rev. D 100, 104015 (2019)

Gravitational waveforms from SpEC simulations : neutron star-neutron star and low-mass black hole-neutron star binaries

Authors: Francois Foucart, Matthew D. Duez, Tanja Hinderer, Jesus Caro, Andrew R. Williamson, Michael Boyle, Alessandra Buonanno, Roland Haas, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more tha… ▽ More Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more than twenty cycles for low-mass black hole-neutron star binaries, including binaries with non-spinning objects, and binaries with rapidly spinning neutron stars that maximize the impact on the gravitational wave signal of the near-resonant growth of the fundamental excitation mode of the neutron star (f-mode). We also provide for the first time with SpEC a high-accuracy neutron star-neutron star waveform. These waveforms are made publicly available as part of the SxS catalogue. We compare our results to analytical waveform models currently implemented in data analysis pipelines. For most simulations, the models lie outside of the predicted numerical errors in the last few orbits before merger, but do not show systematic deviations from the numerical results: comparing different models appears to provide reasonable estimates of the modeling errors. The sole exception is the equal-mass simulation using a rapidly counter-rotating neutron star to maximize the impact of the excitation of the f-mode, for which all models perform poorly. This is however expected, as even the single model that takes f-mode excitation into account ignores the significant impact of the neutron star spin on the f-mode excitation frequency. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 17p, 10 figures

Journal ref: Phys. Rev. D 99, 044008 (2019)

Numerical black hole initial data and shadows in dynamical Chern-Simons gravity

Authors: Maria Okounkova, Mark A. Scheel, Saul A. Teukolsky

Abstract: We present a scheme for generating first-order metric perturbation initial data for an arbitrary background and source. We then apply this scheme to derive metric perturbations in order-reduced dynamical Chern-Simons gravity (dCS). In particular, we solve for metric perturbations on a black hole background that are sourced by a first-order dCS scalar field. This gives us the leading-order metric p… ▽ More We present a scheme for generating first-order metric perturbation initial data for an arbitrary background and source. We then apply this scheme to derive metric perturbations in order-reduced dynamical Chern-Simons gravity (dCS). In particular, we solve for metric perturbations on a black hole background that are sourced by a first-order dCS scalar field. This gives us the leading-order metric perturbation to the spacetime in dCS gravity. We then use these solutions to compute black hole shadows in the linearly perturbed spacetime by evolving null geodesics. We present a novel scheme to decompose the shape of the shadow into multipoles parametrized by the spin of the background black hole and the perturbation parameter $\varepsilon^2$. We find that we can differentiate the presence of a pure Kerr spacetime from a spacetime with a dCS perturbation using the shadow, allowing in part for a null-hypothesis test of general relativity. We then consider these results in the context of the Event Horizon Telescope. △ Less

Submitted 10 January, 2019; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: 31 pages, 8 figures Updated to match accepted version in Classical and Quantum Gravity

Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

Authors: Kabir Chakravarti, Anuradha Gupta, Sukanta Bose, Matthew D. Duez, Jesus Caro, Wyatt Brege, Francois Foucart, Shaon Ghosh, Koutarou Kyutoku, Benjamin D. Lackey, Masaru Shibata, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the… ▽ More We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Journal ref: Phys. Rev. D 99, 024049 (2019)

Distinguishing the nature of comparable-mass neutron star binary systems with multimessenger observations: GW170817 case study

Authors: Tanja Hinderer, Samaya Nissanke, Francois Foucart, Kenta Hotokezaka, Trevor Vincent, Mansi Kasliwal, Patricia Schmidt, Andrew R. Williamson, David Nichols, Matthew Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM obs… ▽ More The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations. △ Less

Submitted 12 March, 2020; v1 submitted 11 August, 2018; originally announced August 2018.

Comments: 8 pages

Journal ref: Phys. Rev. D 100, 063021 (2019)

Evaluating radiation transport errors in merger simulations using a Monte-Carlo algorithm

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

Abstract: Neutrino-matter interactions play an important role in the post-merger evolution of neutron star-neutron star and black hole-neutron star mergers. Most notably, they determine the properties of the bright optical/infrared transients observable after a merger. Unfortunately, Boltzmann's equations of radiation transport remain too costly to be evolved directly in merger simulations. Simulations rely… ▽ More Neutrino-matter interactions play an important role in the post-merger evolution of neutron star-neutron star and black hole-neutron star mergers. Most notably, they determine the properties of the bright optical/infrared transients observable after a merger. Unfortunately, Boltzmann's equations of radiation transport remain too costly to be evolved directly in merger simulations. Simulations rely instead on approximate transport algorithms with unquantified modeling errors. In this paper, we use for the first time a time-dependent general relativistic Monte-Carlo (MC) algorithm to solve Boltzmann's equations and estimate important properties of the neutrino distribution function ~10ms after a neutron star merger. We do not fully couple the MC algorithm to the fluid evolution, but use a short evolution of the merger remnant to critically assess errors in our approximate gray two-moment transport scheme. We demonstrate that the analytical closure used by the moment scheme is highly inaccurate in the polar regions, but performs well elsewhere. While the average energy of polar neutrinos is reasonably well captured by the two-moment scheme, estimates for the neutrino energy become less accurate at lower latitudes. The two-moment formalism also overestimates the density of neutrinos in the polar regions by ~50%, and underestimates the neutrino pair-annihilation rate at the poles by factors of 2-3. Although the latter is significantly more accurate than one might have expected before this study, our results indicate that predictions for the properties of polar outflows and for the creation of a baryon-free region at the poles are likely to be affected by errors in the two-moment scheme, thus limiting our ability to reliably model kilonovae and gamma-ray bursts. △ Less

Submitted 6 June, 2018; originally announced June 2018.

Comments: 15 pages, 10 figures

Journal ref: Phys. Rev. D 98, 063007 (2018)

Black hole-neutron star mergers using a survey of finite-temperature equations of state

Authors: Wyatt Brege, Matthew D. Duez, Francois Foucart, M. Brett Deaton, Jesus Caro, Daniel A. Hemberger, Lawrence E. Kidder, Evan O'Connor, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However,… ▽ More Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However, most of the simulations on which these formulae are based use equations of state without finite temperature and composition-dependent nuclear physics. In this paper, we simulate black hole-neutron star mergers varying both the neutron star mass and the equation of state, using three finite-temperature nuclear models of varying stiffness. Our simulations largely vindicate formulae for ejecta properties but do not find the expected dependence of disk mass on neutron star compaction. We track the early evolution of the accretion disk, largely driven by shocking and fallback inflow, and do find notable equation of state effects on the structure of this early-time, neutrino-bright disk. △ Less

Submitted 25 April, 2018; originally announced April 2018.

Comments: 10 pages, 11 figures, submitted to Phys. Rev. D

Journal ref: Phys. Rev. D 98, 063009 (2018)

Measuring the properties of nearly extremal black holes with gravitational waves

Authors: Katerina Chatziioannou, Geoffrey Lovelace, Michael Boyle, Matthew Giesler, Daniel A. Hemberger, Reza Katebi, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel, Béla Szilágyi

Abstract: Characterizing the properties of black holes is one of the most important science objectives for gravitational-wave observations. Astrophysical evidence suggests that black holes that are nearly extremal (i.e. spins near the theoretical upper limit) might exist and thus might be among the merging black holes observed with gravitational waves. In this paper, we explore how well current gravitationa… ▽ More Characterizing the properties of black holes is one of the most important science objectives for gravitational-wave observations. Astrophysical evidence suggests that black holes that are nearly extremal (i.e. spins near the theoretical upper limit) might exist and thus might be among the merging black holes observed with gravitational waves. In this paper, we explore how well current gravitational wave parameter estimation methods can measure the spins of rapidly spinning black holes in binaries. We simulate gravitational-wave signals using numerical-relativity waveforms for nearly-extremal, merging black holes. For simplicity, we confine our attention to binaries with spins parallel or antiparallel with the orbital angular momentum. We find that recovering the holes' nearly extremal spins is challenging. When the spins are nearly extremal and parallel to each other, the resulting parameter estimates do recover spins that are large, though the recovered spin magnitudes are still significantly smaller than the true spin magnitudes. When the spins are nearly extremal and antiparallel to each other, the resulting parameter estimates recover the small effective spin but incorrectly estimate the individual spins as nearly zero. We study the effect of spin priors and argue that a commonly used prior (uniform in spin magnitude and direction) hinders unbiased recovery of large black-hole spins. △ Less

Submitted 22 August, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: 8 pages, 7 figures, final published version

Journal ref: Phys. Rev. D 98, 044028 (2018)

On choosing the start time of binary black hole ringdown

Authors: Swetha Bhagwat, Maria Okounkova, Stefan W. Ballmer, Duncan A. Brown, Matthew Giesler, Mark A. Scheel, Saul A. Teukolsky

Abstract: The final stage of a binary black hole merger is ringdown, in which the system is described by a Kerr black hole with quasinormal mode perturbations. It is far from straightforward to identify the time at which the ringdown begins. Yet determining this time is important for precision tests of the general theory of relativity that compare an observed signal with quasinormal mode descriptions of the… ▽ More The final stage of a binary black hole merger is ringdown, in which the system is described by a Kerr black hole with quasinormal mode perturbations. It is far from straightforward to identify the time at which the ringdown begins. Yet determining this time is important for precision tests of the general theory of relativity that compare an observed signal with quasinormal mode descriptions of the ringdown, such as tests of the no-hair theorem. We present an algorithmic method to analyze the choice of ringdown start time in the observed waveform. This method is based on determining how close the strong field is to a Kerr black hole (Kerrness). Using numerical relativity simulations, we characterize the Kerrness of the strong-field region close to the black hole using a set of local, gauge-invariant geometric and algebraic conditions that measure local isometry to Kerr. We produce a map that associates each time in the gravitational waveform with a value of each of these Kerrness measures; this map is produced by following outgoing null characteristics from the strong and near-field regions to the wave zone. We perform this analysis on a numerical relativity simulation with parameters consistent with GW150914- the first gravitational wave detection. We find that the choice of ringdown start time of $3\,\mathrm{ms}$ after merger used in the GW150914 study to test general relativity corresponds to a high dimensionless perturbation amplitude of $ \sim 7.5 \times 10^{-3}$ in the strong-field region. This suggests that in higher signal-to-noise detections, one would need to start analyzing the signal at a later time for studies that depend on the validity of black hole perturbation theory. △ Less

Submitted 6 November, 2017; v1 submitted 2 November, 2017; originally announced November 2017.

Comments: 23+4 pages, 22 figures

Journal ref: Phys. Rev. D 97, 104065 (2018)

Evolution of the Magnetized, Neutrino-Cooled Accretion Disk in the Aftermath of a Black Hole Neutron Star Binary Merger

Authors: Fatemeh Hossein Nouri, Matthew D. Duez, Francois Foucart, M. Brett Deaton, Roland Haas, Milad Haddadi, Lawrence E. Kidder, Christian D. Ott, Harald P. Pfeiffer, Mark A. Scheel, Bela Szilagyi

Abstract: Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-n… ▽ More Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-neutron star binary system using the Spectral Einstein Code (SpEC) from an initial post-merger state provided by previous numerical relativity simulations. We use a finite-temperature nuclear equation of state and incorporate neutrino effects in a leakage approximation. To achieve the needed accuracy, we introduce improvements to SpEC's implementation of general-relativistic magnetohydrodynamics (MHD), including the use of cubed-sphere multipatch grids and an improved method for dealing with supersonic accretion flows where primitive variable recovery is difficult. We find that a seed magnetic field triggers a sustained source of heating, but its thermal effects are largely cancelled by the accretion and spreading of the torus from MHD-related angular momentum transport. The neutrino luminosity peaks at the start of the simulation, and then drops significantly over the first 20\,ms but in roughly the same way for magnetized and nonmagnetized disks. The heating rate and disk's luminosity decrease much more slowly thereafter. These features of the evolution are insensitive to grid structure and resolution, formulation of the MHD equations, and seed field strength, although turbulent effects are not fully converged △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: 17 pages, 18 figures

Journal ref: Phys. Rev. D 97, 083014 (2018)

Numerical binary black hole mergers in dynamical Chern-Simons: I. Scalar field

Authors: Maria Okounkova, Leo C. Stein, Mark A. Scheel, Daniel A. Hemberger

Abstract: Testing general relativity in the non-linear, dynamical, strong-field regime of gravity is one of the major goals of gravitational wave astrophysics. Performing precision tests of general relativity (GR) requires numerical inspiral, merger, and ringdown waveforms for binary black hole (BBH) systems in theories beyond GR. Currently, GR and scalar-tensor gravity are the only theories amenable to num… ▽ More Testing general relativity in the non-linear, dynamical, strong-field regime of gravity is one of the major goals of gravitational wave astrophysics. Performing precision tests of general relativity (GR) requires numerical inspiral, merger, and ringdown waveforms for binary black hole (BBH) systems in theories beyond GR. Currently, GR and scalar-tensor gravity are the only theories amenable to numerical simulations. In this article, we present a well-posed perturbation scheme for numerically integrating beyond-GR theories that have a continuous limit to GR. We demonstrate this scheme by simulating BBH mergers in dynamical Chern-Simons gravity (dCS), to linear order in the perturbation parameter. We present mode waveforms and energy fluxes of the dCS pseudoscalar field from our numerical simulations. We find good agreement with analytic predictions at early times, including the absence of pseudoscalar dipole radiation. We discover new phenomenology only accessible through numerics: a burst of dipole radiation during merger. We also quantify the self-consistency of the perturbation scheme. Finally, we estimate bounds that GR-consistent LIGO detections could place on the new dCS length scale, approximately $\ell \lesssim \mathcal{O}(10)~\mathrm{km}$. △ Less

Submitted 17 August, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 14+4 pages, 8 figures, 1 table; visualization available at http://www.youtube.com/watch?v=WQH-1b_XUM4 . Matches published version

Journal ref: Phys. Rev. D 96, 044020 (2017)

Dynamical ejecta from precessing neutron star-black hole mergers with a hot, nuclear-theory based equation of state

Authors: Francois Foucart, Dhruv Desai, Wyatt Brege, Matthew D. Duez, Daniel Kasen, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Neutron star-black hole binaries are among the strongest sources of gravitational waves detectable by current observatories. They can also power bright electromagnetic signals (gamma-ray bursts, kilonovae), and may be a significant source of production of r-process nuclei. A misalignment of the black hole spin with respect to the orbital angular momentum leads to precession of that spin and of the… ▽ More Neutron star-black hole binaries are among the strongest sources of gravitational waves detectable by current observatories. They can also power bright electromagnetic signals (gamma-ray bursts, kilonovae), and may be a significant source of production of r-process nuclei. A misalignment of the black hole spin with respect to the orbital angular momentum leads to precession of that spin and of the orbital plane, and has a significant effect on the properties of the post-merger remnant and of the material ejected by the merger. We present a first set of simulations of precessing neutron star-black hole mergers using a hot, composition dependent, nuclear-theory based equation of state (DD2). We show that the mass of the remnant and of the dynamical ejecta are broadly consistent with the result of simulations using simpler equations of state, while differences arise when considering the dynamics of the merger and the velocity of the ejecta. We show that the latter can easily be understood from assumptions about the composition of low-density, cold material in the different equations of state, and propose an updated estimate for the ejecta velocity which takes those effects into account. We also present an updated mesh-refinement algorithm which allows us to improve the numerical resolution used to evolve neutron star-black hole mergers. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 23 pages, 8 figures

Journal ref: Classical and Quantum Gravity, Volume 34, Number 4 (2017)

Complete waveform model for compact binaries on eccentric orbits

Authors: E. A. Huerta, Prayush Kumar, Bhanu Agarwal, Daniel George, Hsi-Yu Schive, Harald P. Pfeiffer, Roland Haas, Wei Ren, Tony Chu, Michael Boyle, Daniel A. Hemberger, Lawrence E. Kidder, Mark A. Scheel, Bela Szilagyi

Abstract: We present a time domain waveform model that describes the inspiral-merger-ringdown (IMR) of compact binary systems whose components are non-spinning, and which evolve on orbits with low to moderate eccentricity. The inspiral evolution is described using third order post-Newtonian equations both for the equations of motion of the binary, and its far-zone radiation field. This latter component also… ▽ More We present a time domain waveform model that describes the inspiral-merger-ringdown (IMR) of compact binary systems whose components are non-spinning, and which evolve on orbits with low to moderate eccentricity. The inspiral evolution is described using third order post-Newtonian equations both for the equations of motion of the binary, and its far-zone radiation field. This latter component also includes instantaneous, tails and tails-of-tails contributions, and a contribution due to non-linear memory. This framework reduces to the post-Newtonian approximant TaylorT4 at third post-Newtonian order in the zero eccentricity limit. To improve phase accuracy, we incorporate higher-order post-Newtonian corrections for the energy flux of quasi-circular binaries and gravitational self-force corrections to the binding energy of compact binaries. This enhanced inspiral evolution prescription is combined with an analytical prescription for the merger-ringdown evolution using a catalog of numerical relativity simulations. This IMR waveform model reproduces effective-one-body waveforms for systems with mass-ratios between 1 to 15 in the zero eccentricity limit. Using a set of eccentric numerical relativity simulations, not used during calibration, we show that our eccentric model accurately reproduces the features of eccentric compact binary coalescence throughout the merger. Using this model we show that the gravitational wave transients GW150914 and GW151226 can be effectively recovered with template banks of quasi-circular, spin-aligned waveforms if the eccentricity $e_0$ of these systems when they enter the aLIGO band at a gravitational wave frequency of 14 Hz satisfies $e_0^{\rm GW150914}\leq0.15$ and $e_0^{\rm GW151226}\leq0.1$. △ Less

Submitted 1 February, 2017; v1 submitted 19 September, 2016; originally announced September 2016.

Comments: 31 pages, 19 figures, 3 appendices. Submitted to Phys Rev D. v2: direct comparison to eccentric numerical relativity simulations, not used for calibration of the model, included. References added. Accepted to Phys Rev D

ACM Class: J.2

Journal ref: Phys. Rev. D 95, 024038 (2017)

SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics

Authors: Lawrence E. Kidder, Scott E. Field, Francois Foucart, Erik Schnetter, Saul A. Teukolsky, Andy Bohn, Nils Deppe, Peter Diener, François Hébert, Jonas Lippuner, Jonah Miller, Christian D. Ott, Mark A. Scheel, Trevor Vincent

Abstract: We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resoluti… ▽ More We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads. △ Less

Submitted 21 July, 2017; v1 submitted 31 August, 2016; originally announced September 2016.

Comments: 41 pages, 13 figures, and 7 tables. Ancillary data contains simulation input files

Journal ref: Journal of Computational Physics, Volume 335, 2017, Pages 84-114

Impact of an improved neutrino energy estimate on outflows in neutron star merger simulations

Authors: Francois Foucart, Evan O'Connor, Luke Roberts, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of… ▽ More Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and timescale of the electromagnetic signal depend on the orientation of the binary with respect to an observer's line-of-sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme. △ Less

Submitted 2 December, 2016; v1 submitted 25 July, 2016; originally announced July 2016.

Comments: 19p, 17 figures, Accepted by Phys.Rev.D

Journal ref: Phys. Rev. D 94, 123016 (2016)

Modeling the source of GW150914 with targeted numerical-relativity simulations

Authors: Geoffrey Lovelace, Carlos O. Lousto, James Healy, Mark A. Scheel, Alyssa Garcia, Richard O'Shaughnessy, Michael Boyle, Manuela Campanelli, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Bela Szilagyi, Saul A. Teukolsky, Yosef Zlochower

Abstract: In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predic… ▽ More In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predicted by general relativity, can be computed only by full numerical relativity, because analytic approximations fail near the time of merger. Moreover, the nearly-equal masses, moderate spins, and small number of orbits of GW150914 are especially straightforward and efficient to simulate with modern numerical-relativity codes. In this paper, we report the modeling of GW150914 with numerical-relativity simulations, using black-hole masses and spins consistent with those inferred from LIGO's measurement. In particular, we employ two independent numerical-relativity codes that use completely different analytical and numerical methods to model the same merging black holes and to compute the emitted gravitational waveform; we find excellent agreement between the waveforms produced by the two independent codes. These results demonstrate the validity, impact, and potential of current and future studies using rapid-response, targeted numerical-relativity simulations for better understanding gravitational-wave observations. △ Less

Submitted 18 July, 2016; originally announced July 2016.

Comments: 11 pages, 3 figures, submitted to Classical and Quantum Gravity

Journal ref: Class. Quantum Grav. 33, 244002 (2016)

Properties of the Binary Black Hole Merger GW150914

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca , et al. (980 additional authors not shown)

Abstract: On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in ge… ▽ More On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime. △ Less

Submitted 20 June, 2016; v1 submitted 11 February, 2016; originally announced February 2016.

Comments: 19 pages, 6 figures

Report number: LIGO-P1500218

Journal ref: Phys. Rev. Lett. 116, 241102 (2016)

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

Authors: Tanja Hinderer, Andrea Taracchini, Francois Foucart, Alessandra Buonanno, Jan Steinhoff, Matthew Duez, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel, Bela Szilagyi, Kenta Hotokezaka, Koutarou Kyutoku, Masaru Shibata, Cory W. Carpenter

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

Submitted 9 May, 2016; v1 submitted 1 February, 2016; originally announced February 2016.

Journal ref: Phys. Rev. Lett. 116, 181101 (2016)

Low mass binary neutron star mergers : gravitational waves and neutrino emission

Authors: Francois Foucart, Roland Haas, Matthew D. Duez, Evan O'Connor, Christian D. Ott, Luke Roberts, Lawrence E. Kidder, Jonas Lippuner, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger… ▽ More Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes. △ Less

Submitted 22 January, 2016; v1 submitted 21 October, 2015; originally announced October 2015.

Comments: Accepted by PRD; 23 pages; 24 figures; 4 tables

Journal ref: Phys. Rev. D 93, 044019 (2016)

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

Authors: Jonathan Blackman, Scott E. Field, Chad R. Galley, Bela Szilagyi, Mark A. Scheel, Manuel Tiglio, Daniel A. Hemberger

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

Submitted 7 October, 2015; v1 submitted 26 February, 2015; originally announced February 2015.

Comments: Updated to published version, which includes a section comparing the surrogate and effective-one-body models. The surrogate is publicly available for download at http://www.black-holes.org/surrogates/ . 6 pages, 6 figures

Journal ref: Phys. Rev. Lett. 115, 121102 (2015)