Search

Binary black holes are the most abundant source of gravitational-wave observations. Gravitational-wave observatories in the next decade will require tremendous increases in the accuracy of numerical waveforms modeling binary black holes, compared to today's state of the art. One approach to achieving the required accuracy is using spectral-type methods that scale to many processors. Using the SpECTRE numerical-relativity code, we present the first simulations of a binary black hole inspiral, merger, and ringdown using discontinuous Galerkin methods. The efficiency of discontinuous Galerkin methods allows us to evolve the binary through ~18 orbits at reasonable computational cost. We then use SpECTRE's Cauchy Characteristic Evolution (CCE) code to extract the gravitational waves at future null infinity. The open-source nature of SpECTRE means this is the first time a spectral-type method for simulating binary black hole evolutions is available to the entire numerical-relativity community., Comment: 32 pages, 12 figures, 28 ancillary input files

We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods., Comment: 31 pages, 8 figures, minor wording updates

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

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., Comment: 10 pages, 3 figures. Published in PRD. Model publicly available at https://pypi.org/project/surfinBH

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., Comment: 48 pages, 17 figures. Matches published version, minor changes only

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