Notes on the integration of numerical relativity waveforms
Abstract
A primary goal of numerical relativity is to provide estimates of the wave strain, , from strong gravitational wave sources, to be used in detector templates. The simulations, however, typically measure waves in terms of the Weyl curvature component, . Assuming Bondi gauge, transforming to the strain reduces to integration of twice in time. Integrations performed in either the time or frequency domain, however, lead to secular non-linear drifts in the resulting strain . 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.
Cite
@article{arxiv.1006.1632,
title = {Notes on the integration of numerical relativity waveforms},
author = {Christian Reisswig and Denis Pollney},
journal= {arXiv preprint arXiv:1006.1632},
year = {2011}
}
Comments
23 pages, 10 figures, matches final published version