Search for: All records

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

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

Abstract Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes (BBHs) can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for subdominant modes. A system of higher order absorbing boundary conditions which greatly reduces this problem was introduced in earlier work (Buchman and Sarbach 2006Class. Quantum Grav.236709). In this paper, we devise two new implementations of this boundary condition system in the Spectral Einstein Code (SpEC), and test them in both linear multipolar gravitational wave and inspiralling mass ratio 7:1 BBH simulations. One of our implementations in particular is shown to be extremely robust and to produce accuracy superior to the standard freezing-Ψ₀boundary condition usually used bySpEC. 

Abstract Gravitational memory effects are predictions of general relativity that are characterized by an observable effect that persists after the passage of gravitational waves. In recent years, they have garnered particular interest, both due to their connection to asymptotic symmetries and soft theorems and because their observation would serve as a unique test of the nonlinear nature of general relativity. Apart from the more commonly known displacement and spin memories, however, there are other memory effects predicted by Einstein’s equations that are associated with more subleading terms in the asymptotic expansion of the Bondi-Sachs metric. In this paper, we write explicit expressions for these higher memory effects in terms of their charge and flux contributions. Further, by using a numerical relativity simulation of a binary black hole merger, we compute the magnitude and morphology of these terms and compare them to those of the displacement and spin memory. We find that, although these terms are interesting from a theoretical perspective, due to their small magnitude they will be particularly challenging to observe with current and future detectors. 

One of the most promising avenues to perform numerical evolutions in theories beyond general relativity is the approach, a proposal in which new “driver” equations are added to the evolution equations in a way that allows for stable numerical evolutions. In this direction, we extend the numerical relativity code p to evolve a “fixed” version of scalar Gauss-Bonnet theory in the decoupling limit, a phenomenologically interesting theory that allows for hairy black hole solutions in vacuum. We focus on isolated black hole systems both with and without linear and angular momentum, and propose a new driver equation to improve the recovery of such stationary solutions. We demonstrate the effectiveness of the latter by numerically evolving black holes that undergo spontaneous scalarization using different driver equations. Finally, we evaluate the accuracy of the obtained solutions by comparing with the original unaltered theory. Published by the American Physical Society2024 

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the lowT/∣W∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.