Numerical convergence of model Cauchy-characteristic extraction and matching

Abstract

Gravitational waves provide a powerful enhancement to our understanding of fundamental physics. To make the most of their detection we need to accurately model the entire process of their emission and propagation toward interferometers. Cauchy-characteristic extraction and matching are methods to compute gravitational waves at null infinity, a mathematical idealization of detector location, from numerical relativity simulations. Both methods can in principle contribute to modeling by providing highly accurate gravitational waveforms. An underappreciated subtlety in realizing this potential is posed by the (mere) weak hyperbolicity of the particular PDE systems solved in the characteristic formulation of the Einstein field equations. This shortcoming results from the popular choice of Bondi-like coordinates. So motivated, we construct toy models that capture that PDE structure and study Cauchy-characteristic extraction and matching with them. Where possible we provide energy estimates for their solutions and perform careful numerical norm convergence tests to demonstrate the effect of weak hyperbolicity on Cauchy-characteristic extraction and matching. Our findings strongly indicate that, as currently formulated, Cauchy-characteristic matching for the Einstein field equations would provide solutions that are, at best, convergent at an order lower than expected for the numerical method, and may be unstable. In contrast, under certain conditions, the extraction method can provide properly convergent solutions. Establishing however that these conditions hold for the aforementioned characteristic formulations is still an open problem.

20 pages, 11 figures, 2 tables, data and code can be found at http://dx.doi.org/10.5281/zenodo.7981429 and https://github.com/ThanasisGiannakopoulos/model_CCE_CCM_public, updated to match published version

References in corpus (19)

• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Advanced LIGO
• A Gravitational Wave Detector with Cosmological Reach
• Extraction of Gravitational Waves in Numerical Relativity
• Adding Gravitational Memory to Waveform Catalogs using BMS Balance Laws
• Regularity of the Einstein Equations at Future Null Infinity
• Strategies for the Characteristic Extraction of Gravitational Waveforms
• Tetrad formalism for numerical relativity on conformally compactified constant mean curvature hypersurfaces
• Extending Gravitational Wave Extraction Using Weyl Characteristic Fields
• A framework for large-scale relativistic simulations in the characteristic approach
• Spectral Characteristic Evolution: A New Algorithm for Gravitational Wave Propagation
• Comparing Remnant Properties from Horizon Data and Asymptotic Data in Numerical Relativity
• Gauge Drivers for the Generalized Harmonic Einstein Equations
• High-accuracy waveforms for black hole-neutron star systems with spinning black holes
• Global simulations of Minkowski space-time including space-like infinity
• The affine-null formulation of the gravitational equations: spherical case
• Gauge structure of the Einstein field equations in Bondi-like coordinates
• The Bondi problem revisited: a spectral domain decomposition code
• A symmetric hyperbolic formulation of the vacuum Einstein equations in affine-null coordinates

Cited by in corpus (8)

• Fully relativistic three-dimensional Cauchy-characteristic matching for physical degrees of freedom
• Merging black holes with Cauchy-characteristic matching: Computation of late-time tails
• A waveform model for the missing quadrupole mode from black hole coalescence: memory effect and ringdown of the $(\ell =2,m=0)$ spherical harmonic
• Einstein-Klein-Gordon system via Cauchy-characteristic evolution: Computation of memory and ringdown tail
• Simulations of gravitational collapse in null coordinates: III. Hyperbolicity
• Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions
• Quasi-Normal Mode Ringing of Binary Black Hole Mergers in Scalar-Gauss-Bonnet Gravity
• Ringdown mode amplitudes of charged binary black holes