Optimal design of interpolation methods for time-delay interferometry

TL;DR

This paper develops a new interpolation kernel design for time-delay interferometry in space-based gravitational wave detectors, improving accuracy, reducing computational cost, and avoiding spectral artifacts.

Contribution

It introduces a novel method for designing smooth interpolation kernels with fewer coefficients, tailored for high-precision TDI applications.

Findings

A 22-coefficient kernel meets LISA's picometre precision requirements.

The new kernel reduces computational cost compared to Lagrange interpolation.

Enhanced robustness against data artifacts and spectral glitches.

Abstract

Time-delay interferometry (TDI) suppresses laser frequency noise by forming linear combinations of time-shifted interferometric measurements. The time-shift operation is implemented by interpolating discretely sampled data. To enable in-band laser noise reduction by eight to nine orders of magnitude, interpolation has to be performed with high accuracy. Optimizing the design of those interpolation methods is the focus of this work. Previous research that studied constant time-shifts suggested Lagrange interpolation as the interpolation method for TDI. Its transfer function performs well at low frequency but requires a high number of coefficients. Furthermore, when applied in TDI we observed prominent time-domain features when a time-varying shift scanned over a pure integer sample shift. To limit this effect we identify an additional requirement for the interpolation kernel: when considering time-varying shifts the interpolation kernel must be sufficiently smooth to avoid unwanted time-domain transitions that produce glitch-like features in power spectral density estimates. The Lagrange interpolation kernel exhibits a discontinuous first derivative by construction, which is insufficient for the application to LISA or other space-based GW observatories. As a solution we propose a novel design method for interpolation kernels that respect a predefined requirement on in-band interpolation residuals and that possess continuous derivatives up to a prescribed order. Using this method we show that an interpolation kernel with 22 coefficients is sufficient to respect LISA's picometre-requirement and to allow for a continuous first derivative which suppresses the magnitude of the time-domain transition adequately. The reduction from 42 (Lagrange interpolation) to 22 coefficients enables us to save computational cost and increases robustness against artefacts in the data.