Fourier Analysis of Gapped Time Series: Improved Estimates of Solar and Stellar Oscillation Parameters

Quantitative helioseismology and asteroseismology require very precise measurements of the frequencies, amplitudes, and lifetimes of the global modes of stellar oscillation. The precision of these measurements depends on the total length (T), quality, and completeness of the observations. Except in a few simple cases, the effect of gaps in the data on measurement precision is poorly understood, in particular in Fourier space where the convolution of the observable with the observation window introduces correlations between different frequencies. Here we describe and implement a rather general method to retrieve maximum likelihood estimates of the oscillation parameters, taking into account the proper statistics of the observations. Our fitting method applies in complex Fourier space and exploits the phase information. We consider both solar-like stochastic oscillations and long-lived harmonic oscillations, plus random noise. Using numerical simulations, we demonstrate the existence of cases for which our improved fitting method is less biased and has a greater precision than when the frequency correlations are ignored. This is especially true of low signal-to-noise solar-like oscillations. For example, we discuss a case where the precision of the mode frequency estimate is increased by a factor of five, for a duty cycle of 15%. In the case of long-lived sinusoidal oscillations, a proper treatment of the frequency correlations does not provide any significant improvement; nevertheless, we confirm that the mode frequency can be measured from gapped data with a much better precision than the 1/T Rayleigh resolution.

[1]  Maurice G. Kendall,et al.  The Advanced Theory of Statistics, Vol. 2: Inference and Relationship , 1979 .

[2]  J. Harvey,et al.  Asymmetries of solar oscillation line profiles , 1993 .

[3]  S. Jefferies,et al.  Modeling of solar oscillation power spectra , 1990 .

[4]  Hans Kjeldsen,et al.  Observations of solar-like oscillations , 2006, astro-ph/0609770.

[5]  P. Goldreich,et al.  Solar seismology. II - The stochastic excitation of the solar p-modes by turbulent convection , 1977 .

[6]  L. Gizon,et al.  The art of fitting p-mode spectra , 1998 .

[7]  W. E. Larson,et al.  Nato Advanced Research Workshop , 1989 .

[8]  Robert F. Stein,et al.  Excitation of Radial P-Modes in the Sun and Stars , 2004 .

[9]  L. Gizon,et al.  The art of fitting p-mode spectra - I. Maximum likelihood estimation , 1997 .

[10]  T. Sekii,et al.  On measuring low-degree p-mode frequency splitting with full-disc integrated data , 2000 .

[11]  Brian Warner,et al.  Asteroseismology of the DOV star PG 1159 - 035 with the Whole Earth Telescope , 1991 .

[12]  R. Paul Butler,et al.  Ultra-High-Precision Velocity Measurements of Oscillations in α Centauri A , 2004 .

[13]  Kenneth G. Libbrecht,et al.  On the Ultimate Accuracy of Solar Oscillation Frequency Measurements , 1992 .

[14]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[15]  Laurent Gizon,et al.  Determining the Inclination of the Rotation Axis of a Sun-like Star , 2003 .