-
General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method
Authors:
Yoonsoo Kim,
Elias R. Most,
William Throwe,
Saul A. Teukolsky,
Nils Deppe
Abstract:
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical approach for simulating such force-free plasmas, which allows for high accuracy in smooth regions as well as capturing dissipation in current sheets. Using a high…
▽ More
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical approach for simulating such force-free plasmas, which allows for high accuracy in smooth regions as well as capturing dissipation in current sheets. Using a high-order accurate discontinuous Galerkin method augmented with a conservative finite-difference method, we demonstrate efficient global simulations of black hole and neutron star magnetospheres. In addition to a series of challenging test problems, we show that our approach can-depending on the physical properties of the system and the numerical implementation-be up to 10x more efficient than conventional simulations, with a speedup of 2-3x for most problems we consider in practice.
△ Less
Submitted 12 June, 2024; v1 submitted 1 April, 2024;
originally announced April 2024.
-
Striking the right tone: toward a self-consistent framework for measuring black hole ringdowns
Authors:
Teagan A. Clarke,
Maximiliano Isi,
Paul D. Lasky,
Eric Thrane,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
Keefe Mitman,
Jordan Moxon,
Kyle C. Nelli,
William Throwe,
Nils L. Vu
Abstract:
The ringdown portion of a binary black hole merger consists of a sum of modes, each containing an infinite number of tones that are exponentially damped sinusoids. In principle, these can be measured as gravitational-waves with observatories like LIGO/Virgo/KAGRA, however in practice it is unclear how many tones can be meaningfully resolved. We investigate the consistency and resolvability of the…
▽ More
The ringdown portion of a binary black hole merger consists of a sum of modes, each containing an infinite number of tones that are exponentially damped sinusoids. In principle, these can be measured as gravitational-waves with observatories like LIGO/Virgo/KAGRA, however in practice it is unclear how many tones can be meaningfully resolved. We investigate the consistency and resolvability of the overtones of the quadrupolar $\ell = m = 2$ mode by starting at late times when the gravitational waveform is expected to be well-approximated by the $\ell m n = 220$ tone alone. We present a Bayesian inference framework to measure the tones in numerical relativity data. We measure tones at different start times, checking for consistency: we classify a tone as stably recovered if and only if the 95\% credible intervals for amplitude and phase at time $t$ overlap with the credible intervals at all subsequent times. We test a set of tones including the first four overtones of the fundamental mode and the 320 tone and find that the 220 and 221 tones can be measured consistently with the inclusion of additional overtones. The 222 tone measurements can be stabilised when we include the 223 tone, but only in a narrow time window, after which it is too weak to measure. The 223 tone recovery appears to be unstable, and does not become stable with the introduction of the 224 tone. We find that $N=3$ tones can be stably recovered simultaneously. However, when analysing $N \geq 4$ tones, the amplitude of one tone is consistent with zero. Thus, within our framework, one can identify only $N=3$ tones with non-zero amplitude that are simultaneously stable.
△ Less
Submitted 11 June, 2024; v1 submitted 5 February, 2024;
originally announced February 2024.
-
Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling
Authors:
Hengrui Zhu,
Justin L. Ripley,
Frans Pretorius,
Sizheng Ma,
Keefe Mitman,
Robert Owen,
Michael Boyle,
Yitian Chen,
Nils Deppe,
Lawrence E. Kidder,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
William Throwe,
Nils L. Vu
Abstract:
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 un…
▽ More
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.
△ Less
Submitted 1 January, 2024;
originally announced January 2024.
-
Black Hole Spectroscopy for Precessing Binary Black Hole Coalescences
Authors:
Hengrui Zhu,
Harrison Siegel,
Keefe Mitman,
Maximiliano Isi,
Will M. Farr,
Michael Boyle,
Nils Deppe,
Lawrence E. Kidder,
Sizheng Ma,
Jordan Moxon,
Kyle C. Nelli,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Vijay Varma,
Nils L. Vu
Abstract:
To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulation…
▽ More
To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulations, we confirm that the relative fundamental quasinormal mode amplitudes of precessing systems are related to those of non-precessing systems by a simple rotation, and that additional structure in the spectrum 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 ${(\ell,m)=(2,\pm 2)}$ quasinormal modes. These results build upon previous works on waveform modeling, and are consistent with a recent ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521.
△ Less
Submitted 13 December, 2023;
originally announced December 2023.
-
A positivity-preserving adaptive-order finite-difference scheme for GRMHD
Authors:
Nils Deppe,
Lawrence E. Kidder,
Saul A. Teukolsky,
Marceline S. Bonilla,
François Hébert,
Yoonsoo Kim,
Mark A. Scheel,
William Throwe,
Nils L. Vu
Abstract:
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. Ou…
▽ More
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.
△ Less
Submitted 18 January, 2024; v1 submitted 7 June, 2023;
originally announced June 2023.
-
Nonlinearities in Black Hole Ringdowns
Authors:
Keefe Mitman,
Macarena Lagos,
Leo C. Stein,
Sizheng Ma,
Lam Hui,
Yanbei Chen,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Jordan Moxon,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils L. Vu
Abstract:
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…
▽ More
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.
△ Less
Submitted 22 February, 2023; v1 submitted 15 August, 2022;
originally announced August 2022.
-
High-accuracy numerical models of Brownian thermal noise in thin mirror coatings
Authors:
Nils L. Vu,
Samuel Rodriguez,
Tom Włodarczyk,
Geoffrey Lovelace,
Harald P. Pfeiffer,
Gabriel S. Bonilla,
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
Jordan Moxon,
William Throwe
Abstract:
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 c…
▽ More
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.
△ Less
Submitted 11 July, 2023; v1 submitted 12 November, 2021;
originally announced November 2021.
-
High Precision Ringdown Modeling: Multimode Fits and BMS Frames
Authors:
Lorena Magaña Zertuche,
Keefe Mitman,
Neev Khera,
Leo C. Stein,
Michael Boyle,
Nils Deppe,
François Hébert,
Dante A. B. Iozzo,
Lawrence E. Kidder,
Jordan Moxon,
Harald P. Pfeiffer,
Mark A. Scheel,
Saul A. Teukolsky,
William Throwe,
Nils Vu
Abstract:
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 im…
▽ More
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.
△ Less
Submitted 11 May, 2022; v1 submitted 29 October, 2021;
originally announced October 2021.
-
Simulating magnetized neutron stars with discontinuous Galerkin methods
Authors:
Nils Deppe,
François Hébert,
Lawrence E. Kidder,
William Throwe,
Isha Anantpurkar,
Cristóbal Armaza,
Gabriel S. Bonilla,
Michael Boyle,
Himanshu Chaudhary,
Matthew D. Duez,
Nils L. Vu,
Francois Foucart,
Matthew Giesler,
Jason S. Guo,
Yoonsoo Kim,
Prayush Kumar,
Isaac Legred,
Dongjun Li,
Geoffrey Lovelace,
Sizheng Ma,
Alexandra Macedo,
Denyz Melchor,
Marlo Morales,
Jordan Moxon,
Kyle C. Nelli
, et al. (11 additional authors not shown)
Abstract:
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 relativisti…
▽ More
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/$ΛΠ^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.
△ Less
Submitted 28 June, 2022; v1 submitted 24 September, 2021;
originally announced September 2021.
-
High-accuracy waveforms for black hole-neutron star systems with spinning black holes
Authors:
Francois Foucart,
Alexander Chernoglazov,
Michael Boyle,
Tanja Hinderer,
Max Miller,
Jordan Moxon,
Mark A. Scheel,
Nils Deppe,
Matthew D. Duez,
Francois Hebert,
Lawrence E. Kidder,
William Throwe,
Harald P. Pfeiffer
Abstract:
The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-pr…
▽ More
The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-precessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to non-spinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy Characteristic Extraction. The new waveforms have $<0.1\,{\rm rad}$ phase errors during inspiral, rising to $\sim (0.2-0.4)\,{\rm rad}$ errors at merger, and $\lesssim 1\%$ error in their amplitude. We compute the faithfulness of recent analytical models to these numerical results, and find that models specifically designed for BHNS binaries perform well ($F>0.99$) for binaries seen face-on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features.
△ Less
Submitted 27 October, 2020;
originally announced October 2020.
-
What does a binary black hole merger look like?
Authors:
Andy Bohn,
William Throwe,
François Hébert,
Katherine Henriksson,
Darius Bunandar,
Nicholas W. Taylor,
Mark A. Scheel
Abstract:
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 lensin…
▽ More
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.
△ Less
Submitted 22 April, 2015; v1 submitted 28 October, 2014;
originally announced October 2014.
