Numerical convergence of model Cauchy-characteristic extraction and matching

TL;DR

This paper investigates the numerical convergence of Cauchy-characteristic extraction and matching methods in gravitational wave modeling, revealing potential stability issues and conditions for proper convergence due to weak hyperbolicity in the equations.

Contribution

It constructs toy models to analyze the impact of weak hyperbolicity on these methods, providing energy estimates and convergence tests to identify stability and convergence limitations.

Findings

Cauchy-characteristic matching may be unstable or converge at lower than expected order.

Cauchy-characteristic extraction can achieve proper convergence under certain conditions.

Weak hyperbolicity significantly affects the numerical stability of these methods.

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.