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 = 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 ∼ 10 −4 . The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% (90 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 . black holes and of the surrounding spacetime [31,32]. Simulations have also been used for visualizations of curved spacetime [33][34][35][36][37][38][39][40], investigations of spin quantities [41], and the relaxation of spacetime to the Kerr solution following merger [42][43][44]. The motion of the black hole horizons and horizon curvature quantities have been used to explore eccentric dynamics [45][46][47][48], spin precession [49][50][51][52], and the first law of binary black hole mechanics [53][54][55][56][57]. These in turn have been compared to analytic post-Newtonian and self-force approximations (see also [58][59][60]), mapping out the bounds of validity of these approximations.A key application of BBH simulations is the accurate modeling of gravitational waves emitted by these systems during their late inspiral, merger, and final ringdown. Waveforms extracted from BBH simulations are essential for analyzing observed gravitational-wave signals from black hole binaries. Indeed, all BBH observations by LIGO and Virgo were analyzed using waveform families that rely on numerical relativity for their construction, most notably effective-one-body waveform models [61-65] and phenomenological waveform models [66][67][68]. Numerical simulations are also central in validating such waveform models [69][70][71][72][73][74][75][76], and were used to validate GW searches [77][78][79]. Waveforms from numerical relativity are also used directly in parameter estimation [80,81], to construct template banks [82], and to construct waveform families without intermediate analytical models, using methods such as reduced order modeling [83][84][85][86]. Today's simula...

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 taskbased 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. variety of astrophysics codes (e.g., Refs. [6,[9][10][11][12][13][14]) have been designed based on these fundamental building blocks.These strategies work well when the computations are reasonably homogeneous or when one seeks good parallelization to only a few thousand cores. As the number of MPI processes increases, so does the cost of communication which, together with non-uniform workload typical of astrophysics problems, limits the maximum number of useful cores that codes can run on. Efficient core utilization becomes non-trivial, often requiring careful optimization by hand to achieve good scalability [15]. Standard finite-volume and finitedifference methods achieve higher order accuracy with increasingly large (overlapping) stencil sizes, and may require additional effort to achieve scalability on massively parallel machines.As one looks ahead to the arrival of exascale computing, it will become increasingly important to focus on developing algorithms that can take full advantage of these very large machines.Discontinuous Galerkin (DG) methods [16-21], together with a task-based parallelization strategy, have the potential to tackle many of these problems. DG methods offer high-order accuracy in smooth regions (although, for stability, increasing the scheme's order requires decreasing the timestep, which restricts the largest usable order in practice), robustness for shocks and other discontinuities, and grid flexibility including a formulation that allows for comparatively straightforward hp-adaptivity and local timestepping. DG methods can be combined with positivity preserving strategies [22][23][24] or "atmosphere treatments" [25] which seek to maintain non-negative values of the pressure and density in challenging regions such as those containing high-speed astrophysical flow. DG methods are also well suited for parallelization: Their formulation in terms of l...

We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS (Simulating eXtreme Spacetimes) Collaboration's SpEC code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using SpECTRE's version of a Cauchycharacteristic extraction. We find that we can accurately resolve the strain's traditional m ¼ 0 memory modes and some of the m ≠ 0 oscillatory memory modes that have previously only been theorized. We also perform a separate calculation of the memory using equations for the Bondi-Metzner-Sachs charges as well as the energy and angular momentum fluxes at asymptotic infinity. Our new calculation uses only the gravitational wave strain and two of the Weyl scalars at infinity. Also, this computation shows that the memory modes can be understood as a combination of a memory signal throughout the binary's inspiral and merger phases, and a quasinormal mode signal near the ringdown phase. Additionally, we find that the magnetic memory, up to numerical error, is indeed zero as previously conjectured. Last, we find that signalto-noise ratios of memory for LIGO, the Einstein Telescope, and the Laser Interferometer Space Antenna with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or minimal waveform models.

Recently it has been argued that in Einstein gravity anti-de Sitter spacetime is unstable against the formation of black holes for a large class of arbitrarily small perturbations. We examine the effects of including a Gauss-Bonnet term. In five dimensions, spherically symmetric Einstein-Gauss-Bonnet gravity has two key features: Choptuik scaling exhibits a radius gap, and the mass function goes to a finite value as the horizon radius vanishes. These suggest that black holes will not form dynamically if the total mass-energy content of the spacetime is too small, thereby restoring the stability of anti-de Sitter spacetime in this context. We support this claim with numerical simulations and uncover a rich structure in horizon radii and formation times as a function of perturbation amplitude.

Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration's catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Lastly, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package sxs, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws.