Numerical relativity simulations provide a full description of the dynamics of binary systems, including gravitational radiation. The waveforms produced by these simulations have a number of applications in gravitational-wave detection and inference. In this work, we revisit the accuracy of the waveforms produced by the Spectral Einstein Code. Motivated by the wide range of waveform applications, we propose and explore three accuracy metrics between simulation resolutions: (i) the generalized frequency-weighted mismatch, (ii) the relative amplitude difference, and (iii) the phase difference at different times. We confirm that numerical errors accumulate over the binary evolution, but the error is not intrinsically larger during the latest, more dynamical stages. Studying errors across the parameter space, we identify a positive correlation between both the mismatch and the phase difference with precessing spin, but little correlation with aligned spin or eccentricity. Lastly, amplitude and phases differences are symmetric upon exchanging resolutions across the catalog, suggesting that the dominant source of error is random, rather than something systematic that affects all waveforms similarly.

Multimessenger observations of neutron star mergers are unique opportunities to constrain the properties of dense matter and the production site of heavy nuclei. To leverage these observations, we require reliable models of the electromagnetic signals powered by mergers. An important limitation to our ability to develop such models is the use of approximate neutrino physics in simulations. Here, we present simulations using an improved version of our Monte Carlo transport algorithm specifically designed to allow for more advanced on-the-fly calculations of reaction rates that use the simulated energy distribution of neutrinos, including in blocking factors, while still relying on approximations for the angular distribution of neutrinos. We use these new methods to include in simulations inelastic scattering of neutrinos on electrons, and to improve our treatment of neutrino-antineutrino pair annihilation. We find that, without increasing the cost of simulations, we can marginally get to the point when the addition of a single packet represents a change Δf_ν<1 in the angle-integrated distribution function, at the cost of increased shot noise in the coupling to the fluid. With inelastic scattering and a better treatment of pair processes, we find a reduction in the average energy and total luminosity of heavy-lepton neutrinos, and an increase in the amount of mass ejected – here by 50%, although on a relatively low amount of total ejected mass <0.005M_⊙. In a separate set of simulations varying the total mass of the binary away from its prompt collapse threshold, we find rapid variations in the amount of ejected matter and in the geometry and composition of the outflows with the total mass of the system. Finally, we use the simulations with our more advanced transport scheme to study in more detail the energy spectrum of neutrinos across the merger remnant.

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

The spectroscopic study of black hole quasinormal modes in gravitational-wave ringdown observations is hindered by our ignorance of which modes should dominate astrophysical signals for different binary configurations, limiting tests of general relativity and astrophysics. In this work, we present a description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasicircular precessing binary black hole coalescences-a key region of parameter space that has yet to be fully explored within the framework of black hole spectroscopy. We suggest that the remnant perturbation for precessing and nonprecessing systems is approximately the same up to a rotation, which implies that the relative amplitudes of the quasinormal modes in both systems are also related by a rotation. We present evidence for this by analyzing an extensive catalog of numerical relativity simulations. Additional structure in the amplitudes is connected to the system's kick velocity and other asymmetries in the orbital dynamics. We find that the ringdowns of precessing systems need not be dominated by the (l, m, n) = (2, +2,0) quasinormal modes, and that instead the (2, +1, 0) or (2, 0,0) quasinormal modes can dominate. Our results are consistent with a ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521, and may also help in understanding phenomenological inspiral-merger-ringdown waveform model systematics.

Cauchy-characteristic evolution (CCE) is a powerful method for accurately extracting gravitational waves at future null infinity. In this work, we extend the previously implemented CCE system within the numerical relativity code SpECTRE by incorporating a scalar field. This allows the system to capture features of beyond-general-relativity theories. We derive scalar contributions to the equations of motion, Weyl scalar computations, Bianchi identities, and balance laws at future null infinity. Our algorithm, tested across various scenarios, accurately reveals memory effects induced by both scalar and tensor fields and captures Price's power-law tail ( u-l-2) in scalar fields at future null infinity, in contrast to the t-2l-3 tail at future timelike infinity.

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 (NR) code, we present the first simulations of a binary black hole inspiral, merger, and ringdown using discontinuous Galerkin (DG) methods. The efficiency of DG methods allows us to evolve the binary through similar to 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 NR community.

The measurement of orbital eccentricity in gravitational-wave (GW) signals will provide unique insights into the astrophysical origin of binary systems, while ignoring eccentricity in waveform models could introduce significant biases in parameter estimation and tests of general relativity. Upcoming LIGO-Virgo-KAGRA observing runs are expected to detect a subpopulation of eccentric signals, making it vital to develop accurate waveform models for eccentric orbits. Here, employing recent analytical results through the third post-Newtonian order, we develop SEOBNRv5EHM: a new time-domain, effective-one-body, multipolar waveform model for eccentric binary black holes with spins aligned (or antialigned) with the orbital angular momentum. Besides the dominant (2, 2) mode, the model includes the (2, 1), (3, 3), (3, 2), (4, 4), and (4, 3) modes. We validate the model's accuracy by computing its unfaithfulness against 99 (28 public and 71 private) eccentric numerical-relativity (NR) simulations, produced by the Simulating eXtreme Spacetimes Collaboration. Importantly, for NR waveforms with initial GW eccentricities below 0.5, the maximum (2, 2)-mode unfaithfulness across the total mass range 20-200M circle dot is consistently below or close to 1%, with a median value of similar to 0.02%, reflecting an accuracy improvement of approximately an order of magnitude compared to the previous-generation SEOBNRv4EHM and the state-of-the-art TEOBResumS-Dal & iacute; eccentric model. In the quasi-circular-orbit limit, SEOBNRv5EHM is in excellent agreement with the highly accurate SEOBNRv5HM model. The accuracy, robustness, and speed of SEOBNRv5EHM make it suitable for data analysis and astrophysical studies. We demonstrate this by performing a set of recovery studies of synthetic NR-signal injections, and parameter-estimation analyses of the events GW150914 and GW190521, which we find to have no eccentricity signatures.

Using high-accuracy numerical relativity waveforms, we confirm the presence of numerous overtones of the l = 2, m = 2 quasinormal mode early in the ringdown of binary black hole mergers. We do this by demonstrating the stability of the mode amplitudes at different fit times, ruling out the possibility that a linear superposition of modes unphysically fits a highly nonlinear part of the waveform. We also find a number of previously unidentified subdominant second-order quasinormal modes in the (2, 2) mode. Even though these modes are mathematically nonlinear, they nevertheless confirm the validity of perturbation theory as a good approximation for describing much of the ringdown.

We present the results of a search for gravitational-wave transients associated with core-collapse supernova SN 2023ixf, which was observed in the galaxy Messier 101 via optical emission on 2023 May 19, during the LIGO–Virgo–KAGRA 15th Engineering Run. We define a five-day on-source window during which an accompanying gravitational-wave signal may have occurred. No gravitational waves have been identified in data when at least two gravitational-wave observatories were operating, which covered ∼14% of this five-day window. We report the search detection efficiency for various possible gravitational-wave emission models. Considering the distance to M101 (6.7 Mpc), we derive constraints on the gravitational-wave emission mechanism of core-collapse supernovae across a broad frequency spectrum, ranging from 50 Hz to 2 kHz, where we assume the gravitational-wave emission occurred when coincident data are available in the on-source window. Considering an ellipsoid model for a rotating proto-neutron star, our search is sensitive to gravitational-wave energy 1 × 10 −4 M ⊙ c 2 and luminosity 2.6 × 10 −4 M ⊙ c 2 s −1 for a source emitting at 82 Hz. These constraints are around an order of magnitude more stringent than those obtained so far with gravitational-wave data. The constraint on the ellipticity of the proto-neutron star that is formed is as low as 1.08, at frequencies above 1200 Hz, surpassing past results.

We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received with low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum-likelihood Non-imaging Transient Reconstruction and Temporal Search pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10(-3) Hz, we compute the GW-BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers.

Cauchy-characteristic matching (CCM) is a numerical-relativity technique that solves Einstein's equations on an effectively infinite computational domain, thereby eliminating systematic errors associated with artificial boundary conditions. Whether CCM can robustly handle fully nonlinear, dynamical spacetimes, such as binary black hole (BBH) mergers, has remained an open question. In this work, we provide a positive answer by presenting nine successful CCM simulations of BBHs; and demonstrate a key application of this method: computing late-time tails. Our results pave the path for systematic studies of late-time tails in BBH systems, and for producing highly accurate waveforms essential to next-generation gravitational-wave detectors.

When numerically solving Einstein’s equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ‘junk radiation.’ This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produced. Thus, an automated procedure for identifying junk radiation is required. To this end, we present a new algorithm based on the empirical mode decomposition (EMD) from the Hilbert–Huang transform. This approach allows us to isolate and measure the high-frequency oscillations present in the measured irreducible masses of the black holes. The decay of these oscillations allows us to estimate the time from which the junk radiation can be ignored. To make this procedure more precise, we propose three distinct threshold criteria that specify how small the contribution of junk radiation has to be before it can be considered negligible. We apply this algorithm to 3403 BBH simulations from the Simulating eXtreme Spacetime catalog to find appropriate values for the thresholds in the three criteria. We find that this approach yields reliable decay time estimates, i.e. when to consider the simulation physical, for 98.5% of the simulations studied. This demonstrates the efficacy of the EMD as a suitable tool to automatically isolate and characterize junk radiation in the simulation of binary black hole systems.

Highly precise and robust waveform models are required as improvements in detector sensitivity enable us to test general relativity with more precision than ever before. In this work, we introduce a spin-aligned surrogate ringdown model. This ringdown surrogate, NRSur3dq8_RD, is built with numerical waveforms produced using Cauchy-characteristic evolution. In addition, these waveforms are in the superrest frame of the remnant black hole allowing us to do a correct analysis of the ringdown spectrum. The novel prediction of our surrogate model is complex-valued quasinormal mode (QNM) amplitudes, with median relative errors of 10^-2-10^-3 over the parameter space. Like previous remnant surrogates, we also predict the remnant black hole's mass and spin. The QNM mode amplitude errors translate into median errors on ringdown waveforms of 10^-4. The high accuracy and QNM mode content provided by our surrogate will enable high-precision ringdown analyses such as tests of general relativity. Our ringdown model is publicly available through the python package surfinBH.

Numerical relativity simulations of merging black holes provide the most accurate description of the binary dynamics and the emitted gravitational wave signal. However, practical considerations such as imperfect initial data and initial parameters mean that achieving target parameters, such as the orbital eccentricity or the black hole spin directions, at the beginning of the usable part of the simulation is challenging. In this paper, we devise a method to produce simulations with specific target parameters, namely the Keplerian orbital parameters-eccentricity, semimajor axis, mean anomaly-and the black hole spin vectors using SpEC. The method is an extension of the current process for achieving vanishing eccentricity and it is based on a parameter control loop that iteratively numerically evolves the system, fits the orbit with analytical post-Newtonian equations, and calculates updated input parameters. Through SpEC numerical simulations, we demonstrate less than or similar to 10-3 and O(degree) convergence for the orbital eccentricity and the spin directions respectively in <= 7 iterations. These tests extend to binaries with mass ratios q <= 3, eccentricities e <= 0.65, and spin magnitudes j chi j <= 0.75. Our method for controlling the orbital and spin parameters of numerical simulations can be used to produce targeted simulations in sparsely covered regions of the parameter space or study the dynamics of relativistic binaries.

In linear perturbation theory, the ringdown of a gravitational wave (GW) signal is described by a linear combination of quasinormal modes (QNMs). Detecting QNMs from GW signals is a promising way to test GR, central to the developing field of black-hole spectroscopy. More robust black hole spectroscopy tests could also consider the ringdown amplitude-phase consistency. That requires an accurate understanding of the excitation and stability of the QNM expansion coefficients. In this paper, we investigate the robustness of the extracted m = 2 QNM coefficients obtained from a high-accuracy numerical relativity waveform. We explore a framework to assess the robustness of QNM coefficients. Within this framework, we not only consider the traditional criterion related to the constancy of a QNM's expansion coefficients over a window in time, but also emphasize the importance of consistency among fitting models. In addition, we implement an iterative greedy approach within which we fix certain QNM coefficients. We apply this approach to linear fitting, and to nonlinear fitting where the properties of the remnant black hole are treated as unknown variables. We find that the robustness of overtone coefficients is enhanced by our greedy approach, particularly for the (2, 2, 2, +) overtone. Based on our robustness criteria applied to the m = 2 signal modes, we find the (2 4, 2, 0, +) and (2, 2, 1 2, +) modes are robust, while the (3, 2, 1, +) subdominant mode is only marginally robust. After we subtract the contributions of the (2 4, 2, 0, +) and (2 3, 2, 1, +) QNMs from signal mode (4, 2), we also find evidence for the quadratic QNM (2, 1, 0, +) x (2, 1, 0, +).

Simulation of quasicircular compact binaries is a major goal in numerical relativity, as they are expected to constitute most gravitational wave observations. However, given that orbital eccentricity is not well defined in general relativity, providing initial data for such binaries is a challenge for numerical simulations. Most numerical relativity codes obtain initial conditions for low-eccentricity binary simulations by iterating over a sequence of short simulations-measuring eccentricity mid-evolution and correcting the initial data parameters accordingly. Eccentricity measurement depends on a numerically challenging nonlinear fit to an estimator model, and the resulting eccentricity estimate is extremely sensitive to small changes in how the fit is performed. We have developed an improved algorithm that produces more consistent measurements of eccentricity relative to the time window chosen for fitting. The primary innovations are the use of the nonlinear optimization algorithm, variable projection, in place of more conventional routines, an initial fit parameter guess taken from the trajectory frequency spectrum, and additional frequency processing of the trajectory data prior to fitting.

Numerical relativists can now produce gravitational waveforms with memory effects routinely and accurately. The gravitational-wave memory effect contains very low-frequency components, including a persistent offset. The presence of these components violates basic assumptions about time-shift behavior underpinning standard data-analysis techniques in gravitational-wave astronomy. This poses a challenge to the analysis of waveform spectra: How to preserve the low-frequency characteristics when transforming a time-domain waveform to the frequency domain. To tackle this challenge, we revisit the preprocessing procedures applied to the waveforms that contain memory effects. We find inconsistency between the zero-frequency limit of displacement memory and the low- frequency spectrum of the same memory preprocessed using the common scheme in literature. To resolve the inconsistency, we propose a new robust preprocessing scheme that produces the spectra of memory waveforms more faithfully. Using this new scheme, we inspect several characteristics of the spectrum of a memory waveform. In particular, we find a discernible beating pattern formed by the dominant oscillatory mode and the displacement memory. This pattern is absent in the spectrum of a waveform without memory. The difference between the memory and no-memory waveforms is too small to be observed by current-generation detectors in a single binary-black-hole event. Detecting the memory in a single event is likely to occur in the era of next-generation detectors.

Gravitational memory effects and the BMS freedoms exhibited at future null infinity have recently been resolved and utilized in numerical relativity simulations. With this, gravitational wave models and our understanding of the fundamental nature of general relativity have been vastly improved. In this paper, we review the history and intuition behind memory effects and BMS symmetries, how they manifest in gravitational waves, and how controlling the infinite number of BMS freedoms of numerical relativity simulations can crucially improve the waveform models that are used by gravitational wave detectors. We reiterate the fact that, with memory effects and BMS symmetries, not only can these next-generation numerical waveforms be used to observe never-before-seen physics, but they can also be used to test GR and learn new astrophysical information about our universe.

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.

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.