On detecting the large separation in the autocorrelation of stellar oscillation times series

The observations carried out by the space missions CoRoT and Kepler provide a large set of asteroseismic data. Their analysis requires an efficient procedure first to determine if the star is reliably showing solar-like oscillations, second to measure the so-called large separation, third to estimate the asteroseismic information that can be retrieved from the Fourier spectrum. We develop in this paper a procedure, based on the autocorrelation of the seismic Fourier spectrum. We have searched for criteria able to predict the output that one can expect from the analysis by autocorrelation of a seismic time series. First, the autocorrelation is properly scaled for taking into account the contribution of white noise. Then, we use the null hypothesis H0 test to assess the reliability of the autocorrelation analysis. Calculations based on solar and CoRoT times series are performed in order to quantify the performance as a function of the amplitude of the autocorrelation signal. We propose an automated determination of the large separation, whose reliability is quantified by the H0 test. We apply this method to analyze a large set of red giants observed by CoRoT. We estimate the expected performance for photometric time series of the Kepler mission. Finally, we demonstrate that the method makes it possible to distinguish l=0 from l=1 modes. The envelope autocorrelation function has proven to be very powerful for the determination of the large separation in noisy asteroseismic data, since it enables us to quantify the precision of the performance of different measurements: mean large separation, variation of the large separation with frequency, small separation and degree identification.Comment: A&A, in pres

On the detection of Lorentzian profiles in a power spectrum: A Bayesian approach using ignorance priors

Aims. Deriving accurate frequencies, amplitudes, and mode lifetimes from stochastically driven pulsation is challenging, more so, if one demands that realistic error estimates be given for all model fitting parameters. As has been shown by other authors, the traditional method of fitting Lorentzian profiles to the power spectrum of time-resolved photometric or spectroscopic data via the Maximum Likelihood Estimation (MLE) procedure delivers good approximations for these quantities. We, however, show that a conservative Bayesian approach allows one to treat the detection of modes with minimal assumptions (i.e., about the existence and identity of the modes). Methods. We derive a conservative Bayesian treatment for the probability of Lorentzian profiles being present in a power spectrum and describe an efficient implementation that evaluates the probability density distribution of parameters by using a Markov-Chain Monte Carlo (MCMC) technique. Results. Potentially superior to "best-fit" procedure like MLE, which only provides formal uncertainties, our method samples and approximates the actual probability distributions for all parameters involved. Moreover, it avoids shortcomings that make the MLE treatment susceptible to the built-in assumptions of a model that is fitted to the data. This is especially relevant when analyzing solar-type pulsation in stars other than the Sun where the observations are of lower quality and can be over-interpreted. As an example, we apply our technique to CoRoT observations of the solar-type pulsator HD 49933.Comment: 12 pages, 11 figures, accepted for publication in Astronomy and Astrophysic

Mode width fitting with a simple bayesian approach. Application to CoRoT targets HD 181420 and HD 49933

We investigate the asteroseismology of two solar-like targets as observed with the CoRoT satellite, with particular attention paid to the mode fitting. HD 181420 and HD 49933 are typical CoRoT solar-like targets (156 and 60-day runs). The low signal-to-noise ratio (SNR) of about 3-10 prevents us from unambiguously identifying the individual oscillation modes. In particular, convergence problems appear at the edges of the oscillation spectrum. HD 181420 and HD 49933 are typical CoRoT solar-like targets (156 and 60-day runs). The low signal-to-noise ratio (SNR) of about 3-10 prevents us from unambiguously identifying the individual oscillation modes. In particular, convergence problems appear at the edges of the oscillation spectrum. We apply a Bayesian approach to the analysis of these data. We compare the global fitting of the power spectra of this time series, obtained by the classical maximum likelihood (MLE) and the maximum a posteriori (MAP) estimators. We examine the impact of the choice of the priors upon the fitted parameters. We also propose to reduce the number of free parameters in the fitting, by replacing the individual estimate of mode height associated with each overtone by a continuous function of frequency (Gaussian profile). The MAP appears as a powerful tool to constrain the global fits, but it must be used carefully and only with reliable priors. The mode width of the stars increases with the frequency over all the oscillation spectrum.Comment: 10 pages, 9 figures, 2 table

The art of fitting p-mode spectra: Part II. Leakage and noise covariance matrices

In Part I we have developed a theory for fitting p-mode Fourier spectra assuming that these spectra have a multi-normal distribution. We showed, using Monte-Carlo simulations, how one can obtain p-mode parameters using 'Maximum Likelihood Estimators'. In this article, hereafter Part II, we show how to use the theory developed in Part I for fitting real data. We introduce 4 new diagnostics in helioseismology: the $(m,\nu)$ echelle diagramme, the cross echelle diagramme, the inter echelle diagramme, and the ratio cross spectrum. These diagnostics are extremely powerful to visualize and understand the covariance matrices of the Fourier spectra, and also to find bugs in the data analysis code. These diagrammes can also be used to derive quantitative information on the mode leakage and noise covariance matrices. Numerous examples using the LOI/SOHO and GONG data are given.Comment: 17 pages with tex and ps files, submitted to A&A, [email protected]

The art of fitting p-mode spectra: Part I. Maximum Likelihood Estimation

In this article we present our state of the art of fitting helioseismic p-mode spectra. We give a step by step recipe for fitting the spectra: statistics of the spectra both for spatially unresolved and resolved data, the use of Maximum Likelihood estimates, the statistics of the p-mode parameters, the use of Monte-Carlo simulation and the significance of fitted parameters. The recipe is applied to synthetic low-resolution data, similar to those of the LOI, using Monte-Carlo simulations. For such spatially resolved data, the statistics of the Fourier spectrum is assumed to be a multi-normal distribution; the statistics of the power spectrum is \emph{not} a $\chi^{2}$ with 2 degrees of freedom. Results for $l=1$ shows that all parameters describing the p modes can be obtained without bias and with minimum variance provided that the leakage matrix is known. Systematic errors due to an imperfect knowledge of the leakage matrix are derived for all the p-mode parameters.Comment: 13 pages, ps file gzipped. Submitted to A&

Surface Activity and Oscillation Amplitudes of Red Giants in Eclipsing Binaries

Among 19 red-giant stars belonging to eclipsing binary systems that have been identified in Kepler data, 15 display solar-like oscillations. We study whether the absence of mode detection in the remaining 4 is an observational bias or possibly evidence of mode damping that originates from tidal interactions. A careful analysis of the corresponding Kepler light curves shows that modes with amplitudes that are usually observed in red giants would have been detected if they were present. We observe that mode depletion is strongly associated with short-period systems, in which stellar radii account for 16-24 % of the semi-major axis, and where red-giant surface activity is detected. We suggest that when the rotational and orbital periods synchronize in close binaries, the red-giant component is spun up, so that a dynamo mechanism starts and generates a magnetic field, leading to observable stellar activity. Pressure modes would then be damped as acoustic waves dissipate in these fields.Comment: 11 pages, 10 figures, Accepted in Ap

The CoRoT target HD175726: an active star with weak solar-like oscillations

Context. The CoRoT short runs give us the opportunity to observe a large variety of late-type stars through their solar-like oscillations. We report observations of the star HD175726 that lasted for 27 days during the first short run of the mission. The time series reveals a high-activity signal and the power spectrum presents an excess due to solar-like oscillations with a low signal-to-noise ratio. Aims. Our aim is to identify the most efficient tools to extract as much information as possible from the power density spectrum. Methods. The most productive method appears to be the autocorrelation of the time series, calculated as the spectrum of the filtered spectrum. This method is efficient, very rapid computationally, and will be useful for the analysis of other targets, observed with CoRoT or with forthcoming missions such as Kepler and Plato. Results. The mean large separation has been measured to be 97.2+-0.5 microHz, slightly below the expected value determined from solar scaling laws.We also show strong evidence for variation of the large separation with frequency. The bolometric mode amplitude is only 1.7+-0.25 ppm for radial modes, which is 1.7 times less than expected. Due to the low signal-to-noise ratio, mode identification is not possible for the available data set of HD175726. Conclusions. This study shows the possibility of extracting a seismic signal despite a signal-to-noise ratio of only 0.37. The observation of such a target shows the efficiency of the CoRoT data, and the potential benefit of longer observing runs.Comment: 8 pages. Accepted in A&