Improving the efficiency of the detection of gravitational wave signals from inspiraling compact binaries: Chebyshev interpolation

Inspiraling compact-object binary systems are promising gravitational wave sources for ground and space-based detectors. The time-dependent signature of these sources is a well-characterized function of a relatively small number of parameters; thus, the favored analysis technique makes use of matched filtering and maximum likelihood methods. As the parameters that characterize the source model vary, so do the templates against which the detector data are compared in the matched filter. For small variations in the parameters, the filter responses are closely correlated. Current analysis methodology samples a bank of filters whose parameter values are chosen so that the correlation between successive samples from successive filters in the bank is 97%. Correspondingly, the additional information available with each successive template evaluation is, in a real sense, only 3% of that already provided by the nearby templates. The reason for such a dense coverage of parameter space is to minimize the chance that a real signal, near the detection threshold, will be missed by the parameter space sampling. Here we investigate the use of Chebyshev interpolation for reducing the number of templates that must be evaluated to obtain the same analysis sensitivity. Additionally, rather than focus on the “loss” of signal-to-noise associated with the finite number of filters in the template bank, we evaluate the receiver operating characteristic (ROC) as a measure of the effectiveness of an analysis technique. The ROC relates the false alarm probability to the false dismissal probability of an analysis, which are the quantities that bear most directly on the effectiveness of an analysis scheme. As a demonstration, we compare the present “dense sampling” analysis methodology with the “interpolation” methodology using Chebyshev polynomials, restricted to one dimension of the multidimensional analysis problem by plotting the ROC curves. We find that the interpolated search can be arranged to have the same false alarm and false dismissal probabilities as the dense sampling strategy using 25% fewer templates. Generalized to the two-dimensional space used in the computationally limited current analyses, this suggests a factor of 2 increase in computational efficiency; generalized to the full seven-dimensional parameter space that characterizes the signal associated with an eccentric binary system of spinning neutron stars or black holes, it suggests an order of magnitude increase in computational efficiency. Since the computational cost of the analysis is driven almost exclusively by the matched filter evaluations, a reduction in the number of template evaluations translates directly into an increase in computational efficiency; additionally, since the computational cost of the analysis is large, the increased efficiency translates also into an increase in the size of the parameter space that can be analyzed and, thus, the science that can be accomplished with the data.