We produce the first astrophysically relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with ...parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasinormal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio ≳ 180 – 240 , with the precise value depending on the dimension of the GR waveform family used in data analysis.
In classical electrodynamics, an accelerating charged body emits radiation and experiences a corresponding radiation-reaction force, or self-force. We extend to higher order in the total charge a ...previous rigorous derivation of the electromagnetic self-force in flat spacetime by Gralla, Harte, and Wald. The method introduced by Gralla, Harte, and Wald computes the self-force from the Maxwell field equations and conservation of stress-energy in a limit where the charge, size, and mass of the body go to zero, and it does not require regularization of a singular self-field. For our higher-order computation, an adjustment of the definition of the mass of the body is necessary to avoid including self-energy from the electromagnetic field sourced by the body in the distant past. We derive the evolution equations for the mass, spin, and center-of-mass position of the body through second order. We derive, for the first time, the second-order acceleration dependence of the evolution of the spin (self-torque), as well as a mixing between the extended body effects and the acceleration-dependent effects on the overall body motion.
We study the possibility that particle production during inflation could source observable gravity waves on scales relevant for cosmic microwave background experiments. A crucial constraint on such ...scenarios arises because particle production can also source inflaton perturbations and might ruin the usual predictions for a nearly scale-invariant spectrum of nearly Gaussian curvature fluctuations. To minimize this effect, we consider two models of particle production in a sector that is only gravitationally coupled to the inflaton. For a single instantaneous burst of massive particle production, we find that localized features in the scalar spectrum and bispectrum might be observable, but gravitational wave signatures are unlikely to be detectable (due to the suppressed quadrupole moment of nonrelativistic quanta) without invoking some additional effects. We also consider a model with a rolling pseudoscalar that leads to a continuous production of relativistic gauge field fluctuations during inflation. Here we find that gravitational waves from particle production can actually exceed the usual inflationary vacuum fluctuations in a regime where non-Gaussianity is consistent with observational limits. In this model observable B-mode polarization can be obtained for any choice of inflaton potential, and the amplitude of the signal is not necessarily correlated with the scale of inflation.
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 Cauchy-characteristic 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 signal-to-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.
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-van der Burg-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. Last, 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.
We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or ...perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in postprocessing. These methods have been used to produce Ψ4 and h waveforms for the Simulating eXtreme Spacetimes (SXS) waveform catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy characteristic extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors.
We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at I+ from numerical relativity simulations. ...Cauchy-characteristic evolution combines an interior solution of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to I+. The foundation of our improved algorithm is a comprehensive method of handling the gauge transformations between the arbitrarily specified coordinates of the interior Cauchy evolution and the unique (up to Bondi-Metzner-Sachs group transformations) Bondi-Sachs coordinate system of the exterior characteristic evolution. We present a reformulated set of characteristic evolution equations better adapted to numerical implementation. In addition, we develop a method to ensure that the angular coordinates used in the volume during the characteristic evolution are asymptotically inertial. This provides a direct route to an expanded set of waveform outputs and is guaranteed to avoid pure-gauge logarithmic dependence that has caused trouble for previous spectral implementations of the characteristic evolution equations. We construct a set of Weyl scalars compatible with the Bondi-like coordinate systems used in characteristic evolution and determine simple, easily implemented forms for the asymptotic Weyl scalars in our suggested set of coordinates.