Hankel low-rank matrix approximation for gravitational-wave data analysis
Abstract
Next-generation gravitational-wave (GW) detectors, such as the Laser Interferometer Space Antenna (LISA), will observe vast numbers of overlapping signals. Disentangling these signals from instrumental noise and from one another constitutes a significant data analysis challenge. We explore a denoising technique based on embedding time series into Hankel matrices: a superposition of (damped) sinusoids corresponds to a matrix of rank . Thus, the problem of signal extraction is reduced to a structured low-rank approximation problem. Using synthetic data tailored to GW applications, we benchmark three Hankel-based algorithms: ESPRIT, Cadzow iterations, and iteratively reweighted least squares (IRLS). Our test scenarios include isolated and multi-component monochromatic signals, the resolution of sources with closely spaced frequencies, and the recovery of black hole quasinormal modes (QNM). All three algorithms achieve near-optimal performance consistent with Fisher matrix bounds, evidenced by an inverse-square scaling of the mismatch with the signal-to-noise ratio. Furthermore, a proof-of-concept application to numerical relativity waveforms validates the ability of these algorithms to extract QNM frequencies from ringdown signals. Hankel low-rank approximation therefore offers a transparent, computationally efficient avenue for preprocessing GW time series.
Cite
@article{arxiv.2603.15740,
title = {Hankel low-rank matrix approximation for gravitational-wave data analysis},
author = {Nicholas Geissler and Vladimir Strokov and Christian Kümmerle and Sergey Kushnarev and Emanuele Berti},
journal= {arXiv preprint arXiv:2603.15740},
year = {2026}
}
Comments
17 pages, 5 figures; for associated implementation, see https://github.com/CosmoVlad/lisa-hankel/tree/main/pipelines