Abstract:Astrophysical time series often contain periodic signals. The large and growing volume of time series data from photometric surveys demands computationally efficient methods for detecting and characterizing such signals. The most efficient algorithms available for this purpose are those that exploit the $\mathcal{O}(N\log N)$ scaling of the Fast Fourier Transform (FFT). However, these methods are not optimal for non-sinusoidal signal shapes. Template fits (or periodic matched filters) optimize sensitivity for a priori known signal shapes but at a significant computational cost. Current implementations of template periodograms scale as $\mathcal{O}(N_f N_{obs})$, where $N_f$ is the number of trial frequencies and $N_{obs}$ is the number of lightcurve observations, and due to non-convexity, they do not guarantee the best fit at each trial frequency, which can lead to spurious results. In this work, we present a non-linear extension of the Lomb-Scargle periodogram to obtain a template-fitting algorithm that is both accurate (globally optimal solutions are obtained except in pathological cases) and computationally efficient (scaling as $\mathcal{O}(N_f\log N_f)$ for a given template). The non-linear optimization of the template fit at each frequency is recast as a polynomial zero-finding problem, where the coefficients of the polynomial can be computed efficiently with the non-equispaced fast Fourier transform. We show that our method, which uses truncated Fourier series to approximate templates, is an order of magnitude faster than existing algorithms for small problems ($N\lesssim 10$ observations) and 2 orders of magnitude faster for long base-line time series with $N_{obs} \gtrsim 10^4$ observations. An open-source implementation of the fast template periodogram is available at <a class="link-external link-https" href="https://www.github.com/PrincetonUniversity/FastTemplatePeriodogram" rel="external noopener nofollow">this https URL</a>.

An application of Saddlepoint approximation for period detection of stellar light observations

Automated period detection from variable stars' time series database

Improving Power Spectrum Estimation Using Multitapering: Efficient Asteroseismic Analyses for Understanding Stars, the Milky Way, and Beyond

Investigating Light Curve Modulation via Kernel Smoothing. I. Application to 53 fundamental mode and first-overtone Cepheids in the LMC

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

GRAPE: Genetic Routine for Astronomical Period Estimation

Modelling stochastic and quasi-periodic behaviour in stellar time-series: Gaussian process regression versus power-spectrum fitting

Optimal frequency-domain analysis for spacecraft time series: Introducing the missing-data multitaper power spectrum estimator

Inflection point in the power spectrum of stellar brightness variations. I. The model

Measuring periods in aperiodic light curves -- Applying the GPS method to infer rotation periods of solar-like stars

Automatic detection of plateau phases in light curves of variable stars

On the Detection and Characterization of Quasiperiodic Oscillations in Astronomical Time Series: Gamma-Ray Burst X-Ray Light Curves as a Test Case

Long rotation period main-sequence stars from Kepler SAP light curves

Testing whether a signal is strictly periodic. Application to disentangling planets and stellar activity in radial velocities

A statistical primer on classical period-finding techniques in astronomy

A new representation of the light curve and its power density spectrum

FINKER: Frequency Identification through Nonparametric KErnel Regression in astronomical time series