Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.

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 $\chi_1,\chi_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$.

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 .

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.

Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper we construct the first reduced order models of gravitational-wave signals that include the effects of spin-precession, inspiral, merger, and ringdown in compact object binaries, and which are valid for component masses describing binary neutron star, binary black hole and mixed binary systems. This work utilizes the waveform model known as "IMRPhenomPv2". Our ROM enables the use of a fast reduced order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to a half a day, for the longest duration/lowest mass signals. The ROM and ROQ rule are available with the main inference library of the LIGO Scientific Collaboration, LALInference.

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).

One of the main bottlenecks in gravitational wave (GW) astronomy is the high cost of performing parameter estimation and GW searches on the fly. We propose a novel technique based on Reduced Order Quadratures (ROQs), an application and data-specific quadrature rule, to perform fast and accurate likelihood evaluations. These are the dominant cost in Markov chain Monte Carlo (MCMC) algorithms, which are widely employed in parameter estimation studies, and so ROQs offer a new way to accelerate GW parameter estimation. We illustrate our approach using a four dimensional GW burst model embedded in noise. We build an ROQ for this model, and perform four dimensional MCMC searches with both the standard and ROQs quadrature rules, showing that, for this model, the ROQ approach is around 25 times faster than the standard approach with essentially no loss of accuracy. The speed-up from using ROQs is expected to increase for more complex GW signal models and therefore has significant potential to accelerate parameter estimation of GW sources such as compact binary coalescences.

We construct compact and high accuracy Reduced Basis (RB) representations of single and multiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small set of the most relevant waveforms. We find that the exponential convergence of the method allows for a dramatic compression of template banks used for ringdown searches. Compressing a catalog with a minimal match $\MMm=0.99$, we find that the selected RB waveforms are able to represent {\em any} QNM, including those not in the original bank, with extremely high accuracy, typically less than $10^{-13}$. We then extend our studies to two-mode QNMs. Inclusion of a second mode is expected to help with detection, and might make it possible to infer details of the progenitor of the final black hole. We find that the number of RB waveforms needed to represent any two-mode ringdown waveform with the above high accuracy is {\em smaller} than the number of metric-based, one-mode templates with $\MMm=0.99$. For unconstrained two-modes, which would allow for consistency tests of General Relativity, our high accuracy RB has around $10^4$ {\em fewer} waveforms than the number of metric-based templates for $\MMm=0.99$. The number of RB elements grows only linearly with the number of multipole modes versus exponentially with the standard approach, resulting in very compact representations even for many multiple modes. The results of this paper open the possibility of searches of multi-mode ringdown gravitational waves.

We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99% and 99.999% the method generates a factor of about $10-10^5$ fewer templates than standard placement methods. The continuum of gravitational waves can be represented by a finite and comparatively compact basis. The method is robust under variations in the noise of detectors, implying that only a single catalog needs to be generated.