Search

Brownian coating thermal noise in detector test masses is limiting the sensitivity of current gravitational-wave detectors on Earth. Therefore, accurate numerical models can inform the ongoing effort to minimize Brownian coating thermal noise in current and future gravitational-wave detectors. Such numerical models typically require significant computational resources and time, and often involve closed-source commercial codes. In contrast, open-source codes give complete visibility and control of the simulated physics, enable direct assessment of the numerical accuracy, and support the reproducibility of results. In this article, we use the open-source SpECTRE numerical relativity code and adopt a novel discontinuous Galerkin numerical method to model Brownian coating thermal noise. We demonstrate that SpECTRE achieves significantly higher accuracy than a previous approach at a fraction of the computational cost. Furthermore, we numerically model Brownian coating thermal noise in multiple sub-wavelength crystalline coating layers for the first time. Our new numerical method has the potential to enable fast exploration of realistic mirror configurations, and hence to guide the search for optimal mirror geometries, beam shapes and coating materials for gravitational-wave detectors., Comment: 15 pages, 5 figures. Results are reproducible with the ancillary input files

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., 6+2 pages, 4 figures, 1 table. Matches PRL version

The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results that are in different BMS frames. Consequently, to build a model for the waveforms produced during the merging of compact objects, which ideally would be a hybridization of PN, NR, and black hole perturbation theory, one needs a fast and robust method for fixing the BMS freedoms. In this work, we present the first means of fixing the entire BMS freedom of NR waveforms to match the frame of either PN waveforms or black hole perturbation theory. We achieve this by finding the BMS transformations that change certain charges in a prescribed way -- e.g., finding the center-of-mass transformation that maps the center-of-mass charge to a mean of zero. We find that this new method is 20 times faster, and more correct when mapping to the superrest frame, than previous methods that relied on optimization algorithms. Furthermore, in the course of developing this charge-based frame fixing method, we compute the PN expression for the Moreschi supermomentum to 3PN order without spins and 2PN order with spins. This Moreschi supermomentum is effectively equivalent to the energy flux or the null memory contribution at future null infinity $\mathscr{I}^{+}$. From this PN calculation, we also compute oscillatory ($m\not=0$ modes) and spin-dependent memory terms that have not been identified previously or have been missing from strain expressions in the post-Newtonian literature., 16+4 pages with 2 appendices, 4 figures, 1 table

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/$\Lambda\Pi^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., Comment: matches published version. 21 pages, 12 figures. Added KH instability, TOV runs with classical limiters. arXiv admin note: text overlap with arXiv:2109.11645

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., Comment: 17 + 2 pages, 11 figures, 2 tables

We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbations of the final Kerr black hole, and a nonlinear region. At late times, the boundary between the two regions is an infalling shell. The BHP region contains late-time gravitational-waves emitted toward the future horizon, as well as those emitted toward future null infinity. In this region, by imposing no-ingoing wave conditions at past null infinity, and matching out-going waves at future null infinity with waveforms computed from numerical relativity, we can obtain waves that travel toward the future horizon. In particular, the Newman-Penrose $\psi_0$ associated with the in-going wave on the horizon is related to tidal deformations measured by fiducial observers floating above the horizon. We further determine the boundary of the BHP region on the future horizon by imposing that $\psi_0$ inside the BHP region can be faithfully represented by quasi-normal modes. Using a physically-motivated way to impose boundary conditions near the horizon, and applying the so-called Boltzmann reflectivity, we compute the quasi-normal modes of non-rotating ECOs, as well as gravitational-wave echoes. We also investigate the detectability of these echoes in current and future detectors, and prospects for parameter estimation.

Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple days to solve these problems at high resolution and when matter is involved, because they are either hard to parallelize or require a large amount of computational resources. Here we present a new solver for linear and nonlinear elliptic problems that is designed to scale with resolution and to parallelize on computing clusters. To achieve this we employ a discontinuous Galerkin discretization, an iterative multigrid-Schwarz preconditioned Newton-Krylov algorithm, and a task-based parallelism paradigm. To accelerate convergence of the elliptic solver we have developed novel subdomain-preconditioning techniques. We find that our multigrid-Schwarz preconditioned elliptic solves achieve iteration counts that are independent of resolution, and our task-based parallel programs scale over 200 million degrees of freedom to at least a few thousand cores. Our new code solves a classic initial data problem for binary black holes faster than the spectral code SpEC when distributed to only eight cores, and in a fraction of the time on more cores. It is publicly accessible in the next-generation SpECTRE numerical relativity code. Our results pave the way for highly parallel elliptic solves in numerical relativity and beyond., 25 pages, 20 figures, published version. Results are reproducible with the ancillary input files

We propose two frequency-domain filters to analyze ringdown signals of binaryblack hole mergers. The first rational filter is constructed based on a set of(arbitrary) quasi-normal modes (QNMs) of the remnant black holes, whereas thesecond full filter comes from the transmissivity of the remnant black holes.The two filters can remove corresponding QNMs from original time-domainringdowns, while changing early inspiral signals in a trivial way - merely atime and phase shift. After filtering out dominant QNMs, we can visualize theexistence of various subdominant effects. For example, by applying our filtersto a GW150914-like numerical relativity (NR) waveform, we find second-ordereffects in the (l = 4, m = 4), (l = 5, m = 4) and (l = 5, m = 5) harmonics; thespherical-spheroidal mixing mode in the (l = 2,m = 2) harmonic; and a mixingmode in the (l = 2,m = 1) harmonic due to a gravitational recoil. In another NRsimulation where two component spins are anti-aligned with the orbital angularmomentum, we also find retrograde modes. The filters are sensitive to theremnant properties (i.e., mass and spin) and thus have a potential applicationto future data analyses and parameter estimations. We also investigate thestability of the full filter. Its connection to the instability of QNM spectrais discussed.

We present a new study of remnant black hole properties from 13 binary black hole systems, numerically evolved using the Spectral Einstein Code. The mass, spin, and recoil velocity of each remnant were determined quasi-locally from apparent horizon data and asymptotically from Bondi data $(h, \psi_4, \psi_3, \psi_2, \psi_1)$ computed at future null infinity using SpECTRE's Cauchy characteristic evolution. We compare these independent measurements of the remnant properties in the bulk and on the boundary of the spacetime, giving insight into how well asymptotic data are able to reproduce local properties of the remnant black hole in numerical relativity. We also discuss the theoretical framework for connecting horizon quantities to asymptotic quantities and how it relates to our results. This study recommends a simple improvement to the recoil velocities reported in the Simulating eXtreme Spacetimes waveform catalog, provides an improvement to future surrogate remnant models, and offers new analysis techniques for evaluating the physical accuracy of numerical simulations., Comment: 14 pages, 4 figures, 1 table; published Physical Review D

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE's Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during postprocessing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincar\'e charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincar\'e transformations., Comment: 18 pages, 11 figures; Published in Physical Review D

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 $\texttt{sxs}$, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws., 13 pages, 9 figures; Corrected a bug affecting the mismatch plot in Figure 7

We give full details regarding the new Cauchy-characteristic evolution (CCE)system in SpECTRE. The implementation is built to provide streamlinedflexibility for either extracting waveforms during the process of a SpECTREbinary compact object simulation, or as a standalone module for extractingwaveforms from worldtube data provided by another code base. Using our recentlypresented improved analytic formulation, the CCE system is free of pure-gaugelogarithms that would spoil the spectral convergence of the scheme. Itgracefully extracts all five Weyl scalars, in addition to the news and thestrain. The SpECTRE CCE system makes significant improvements on previousimplementations in modularity, ease of use, and speed of computation.

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 Collaboration's $\texttt{SpEC}$ code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using $\texttt{SpECTRE}$'s version of a Cauchy-characteristic extraction. We find that we can accurately resolve the strain's traditional $m=0$ memory modes and some of the $m\not=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. Lastly, we find that signal-to-noise ratios of memory for LIGO, the Einstein Telescope (ET), and the Laser Interferometer Space Antenna (LISA) with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or Minimal Waveform models., 20 pages, 11 figures; 10.1103/PhysRevD.102.104007. Corrected a minor sign error in Eqs. 27, 40, 42, 43, and 51

We present a method of calculating the strong-field gravitational lensing caused by many analytic and numerical spacetimes. We use this procedure to calculate the distortion caused by isolated black holes and by numerically evolved black hole binaries. We produce both demonstrative images illustrating details of the spatial distortion and realistic images of collections of stars taking both lensing amplification and redshift into account. On large scales the lensing from inspiraling binaries resembles that of single black holes, but on small scales the resulting images show complex and in some cases self-similar structure across different angular scales., 10 pages, 12 figures. Supplementary images and movies can be found at http://www.black-holes.org/the-science-numerical-relativity/numerical-relativity/gravitational-lensing

We construct effective-one-body waveform models suitable for data analysis with LISA for extreme-mass ratio inspirals in quasi-circular, equatorial orbits about a spinning supermassive black hole. The accuracy of our model is established through comparisons against frequency-domain, Teukolsky-based waveforms in the radiative approximation. The calibration of eight high-order post-Newtonian parameters in the energy flux suffices to obtain a phase and fractional amplitude agreement of better than 1 radian and 1 % respectively over a period between 2 and 6 months depending on the system considered. This agreement translates into matches higher than 97 % over a period between 4 and 9 months, depending on the system. Better agreements can be obtained if a larger number of calibration parameters are included. Higher-order mass ratio terms in the effective-one-body Hamiltonian and radiation-reaction introduce phase corrections of at most 30 radians in a one year evolution. These corrections are usually one order of magnitude larger than those introduced by the spin of the small object in a one year evolution. These results suggest that the effective-one-body approach for extreme mass ratio inspirals is a good compromise between accuracy and computational price for LISA data analysis purposes.

We construct effective-one-body waveform models suitable for data analysis with LISA for extreme-mass ratio inspirals in quasi-circular, equatorial orbits about a spinning supermassive black hole. The accuracy of our model is established through comparisons against frequency-domain, Teukolsky-based waveforms in the radiative approximation. The calibration of eight high-order post-Newtonian parameters in the energy flux suffices to obtain a phase and fractional amplitude agreement of better than 1 radian and 1 % respectively over a period between 2 and 6 months depending on the system considered. This agreement translates into matches higher than 97 % over a period between 4 and 9 months, depending on the system. Better agreements can be obtained if a larger number of calibration parameters are included. Higher-order mass ratio terms in the effective-one-body Hamiltonian and radiation-reaction introduce phase corrections of at most 30 radians in a one year evolution. These corrections are usually one order of magnitude larger than those introduced by the spin of the small object in a one year evolution. These results suggest that the effective-one-body approach for extreme mass ratio inspirals is a good compromise between accuracy and computational price for LISA data analysis purposes., Comment: 21 pages, 8 figures, submitted to Phys. Rev. D