Beyond prewhitening: detection of gravity modes and their period spacings in slowly pulsating B stars using the multitaper F-test

Gravity modes in main-sequence stars have traditionally been studied using a prewhitening approach, which iteratively identifies modes in the Fourier domain and subsequently tunes their frequencies, amplitudes, and phases through time-domain regression. While effective, this method becomes inefficient when analysing large volumes of long time-series data and often relies on subjective stopping criteria to determine the number of iterations. We aim to perform frequency extraction of gravity modes in slowly pulsating B (SPB) stars using a statistically robust, data-driven approach based on advanced power spectrum and harmonic analysis techniques. Our approach employs the multitaper non-uniform fast Fourier transform, mtNUFFT, a power spectrum estimator that addresses several statistical limitations of traditional methods such as the Lomb-Scargle periodogram. We apply its extension, the multitaper F-test, to extract coherent gravity modes from 4-year Kepler light curves of SPB stars and to search for period spacing patterns among the extracted modes. The multitaper F-test enables fast and accurate extraction of the properties of gravity modes with quasi-infinite lifetimes, preferentially selecting modes that exhibit purely periodic behaviour. Although the method typically extracts fewer frequencies than conventional prewhitening, it recovers most known modes and, in some cases, reveals new ones. We also find evidence for gravity modes with long but finite lifetimes, and detect more than one period spacing pattern in some of the studied SPB stars. Overall, the multitaper F-test offers a more objective and statistically sound alternative to prewhitening. It scales efficiently to large datasets containing thousands of pulsators, and has the potential to facilitate mode identification and to distinguish between the different excitation mechanisms operating in SPB stars.