Improving Power Spectral Estimation using Multitapering: Precise asteroseismic modeling of stars, exoplanets, and beyond

by   Aarya A. Patil, et al.

Asteroseismic time-series data have imprints of stellar oscillation modes, whose detection and characterization through time-series analysis allows us to probe stellar interiors physics. Such analyses usually occur in the Fourier domain by computing the Lomb-Scargle (LS) periodogram, an estimator of the power spectrum underlying unevenly-sampled time-series data. However, the LS periodogram suffers from the statistical problems of (1) inconsistency (or noise) and (2) bias due to high spectral leakage. In addition, it is designed to detect strictly periodic signals but is unsuitable for non-sinusoidal periodic or quasi-periodic signals. Here, we develop a multitaper spectral estimation method that tackles the inconsistency and bias problems of the LS periodogram. We combine this multitaper method with the Non-Uniform Fast Fourier Transform () to more precisely estimate the frequencies of asteroseismic signals that are non-sinusoidal periodic (e.g., exoplanet transits) or quasi-periodic (e.g., pressure modes). We illustrate this using a simulated and the Kepler-91 red giant light curve. Particularly, we detect the Kepler-91b exoplanet and precisely estimate its period, 6.246 ± 0.002 days, in the frequency domain using the multitaper F-test alone. We also integrate into the package to obtain a Kepler-91 age estimate of 3.96 ± 0.48 Gyr. This 36% improvement in age precision relative to the 4.27 ± 0.75 Gyr APOKASC-2 (uncorrected) estimate illustrates that has promising implications for Galactic archaeology, in addition to stellar interiors and exoplanet studies. Our frequency analysis method generally applies to time-domain astronomy and is implemented in the public Python package , available at <>.


page 4

page 6

page 8

page 10

page 13

page 17

page 20


A Fast Template Periodogram for Detecting Non-sinusoidal Fixed-shape Signals in Irregularly Sampled Time Series

Astrophysical time series often contain periodic signals. The large and ...

The Sliding Window Discrete Fourier Transform

This paper introduces a new tool for time-series analysis: the Sliding W...

Magnitude-squared coherence: A powerful tool for disentangling Doppler planet discoveries from stellar activity

If Doppler searches for earth-mass, habitable planets are to succeed, ob...

Enhancing Representation Learning for Periodic Time Series with Floss: A Frequency Domain Regularization Approach

Time series analysis is a fundamental task in various application domain...

A new segmentation method for the homogenisation of GNSS-derived IWV time-series

Homogenization is an important and crucial step to improve the usage of ...

The Permutation-Spectrum Test: Identifying Periodic Signals using the Maximum Fourier Intensity

This paper examines the problem of testing whether a discrete time-serie...

Please sign up or login with your details

Forgot password? Click here to reset