Notes on the integration of numerical relativity waveforms
read the original abstract
A primary goal of numerical relativity is to provide estimates of the wave strain, $h$, from strong gravitational wave sources, to be used in detector templates. The simulations, however, typically measure waves in terms of the Weyl curvature component, $\psi_4$. Assuming Bondi gauge, transforming to the strain $h$ reduces to integration of $\psi_4$ twice in time. Integrations performed in either the time or frequency domain, however, lead to secular non-linear drifts in the resulting strain $h$. These non-linear drifts are not explained by the two unknown integration constants which can at most result in linear drifts. We identify a number of fundamental difficulties which can arise from integrating finite length, discretely sampled and noisy data streams. These issues are an artifact of post-processing data. They are independent of the characteristics of the original simulation, such as gauge or numerical method used. We suggest, however, a simple procedure for integrating numerical waveforms in the frequency domain, which is effective at strongly reducing spurious secular non-linear drifts in the resulting strain.
This paper has not been read by Pith yet.
Forward citations
Cited by 2 Pith papers
-
Exact Mass Conservation in Binary Neutron Star Merger Simulations
A rescaling algorithm for artificial atmospheres achieves exact mass and electron number conservation to round-off precision in binary neutron star merger simulations.
-
Scalar fields around black hole binaries in LIGO-Virgo-KAGRA
Semi-analytic waveform model for scalar environments around black hole binaries is validated against numerical relativity and applied to LIGO-Virgo-KAGRA data to obtain upper limits on scalar densities with tentative ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.