New technique for determining a pulsar period: Waterfall principal component analysis

Abstract

This paper describes a new technique for determining the optimal period of a pulsar and consequently its light curve. The implemented technique makes use of the Principal Component Analysis (PCA) applied to the so-called waterfall diagram, which is a bidimensional representation of the pulsar acquired data. In this context we have developed the python package pywpf to easily retrieve the period with the presented method. We applied this technique to sets of data of the brightest pulsars in visible light that we obtained with the fast photon counter Iqueye. Our results are compared with those obtained by different and more classical analyses (e.g., epoch folding), showing that the periods so determined agree within the errors, and that the errors associated to the waterfall-PCA folding technique are slightly smaller than those obtained by the ${\chi }^2$ epoch folding technique. We also simulated extremely noisy situations, showing that by means of a new merit function associated to the waterfall-PCA folding it is possible to get more confidence on the determined period with respect to the ${\chi }^2$ epoch folding technique.

Accepted for publication by Astronomy & Astrophysics (A&A)

References in corpus (7)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• PINT: A Modern Software Package for Pulsar Timing
• Weighted principal component analysis: a weighted covariance eigendecomposition approach
• Precise optical timing of PSR J1023+0038, the first millisecond pulsar detected with Aqueye+ in Asiago
• The Crab pulsar seen with Aqueye at Asiago Cima Ekar Observatory

Cited by in corpus (1)

• Millisecond pulsars phenomenology under the light of graph theory