Issue 
A&A
Volume 520, SeptemberOctober 2010



Article Number  A60  
Number of page(s)  11  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201014944  
Published online  04 October 2010 
Oscillation mode lifetimes of red giants observed during the initial and first anticentre long run of CoRoT^{}
S. Hekker^{1,2}  C. Barban^{3}  F. Baudin^{4}  J. De Ridder^{2}  T. Kallinger^{5}  T. Morel^{6}  W. J. Chaplin^{1}  Y. Elsworth^{1}
1  University of Birmingham, School of Physics and Astronomy, Edgbaston, Birmingham B15 2TT, UK
2  Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
3  LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris,
92195 Meudon Cedex, France
4  Institute d'Astrophysique Spatiale, UMR 8617, Université Paris XI, Bâtiment 121, 91405 Orsay Cedex, France
5  Department of Physics and Astronomy, University of British Colombia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada
6  Institut d'Astrophysique et de Géophysique, Université de Liège, Allée du 6 Août, 4000 Liège, Belgium
Received 5 May 2010 / Accepted 13 June 2010
Abstract
Context. Long timeseries of data increase the frequency
resolution in the power spectrum. This allows for resolving
stochastically excited modes with long mode lifetimes, as well as
features that are close together in frequency. The CoRoT fields
observed during the initial run and second long run partly overlap, and
stars in this overlapping field observed in both runs are used to
create timeseries with a longer timespan than available from the
individual runs.
Aims. We aim to measure the mode lifetimes of red giants and
compare them with theoretical predictions. We also investigate the
dependence of the mode lifetimes on frequency and the degree of the
oscillation modes.
Methods. We perform simulations to investigate the influence of
the gap in the data between the initial and second long run, the total
length of the run and the signaltonoise ratio on the measured mode
lifetime. This provides us with a correction factor to apply to the
mode lifetimes measured from a maximum likelihood fit to the
oscillation frequencies.
Results. We find that the length of the timeseries, the
signaltonoise ratio and possible gaps do have a nonnegligible effect
on the measurements of the mode lifetime of stochastically excited
oscillation modes, but we can correct for it. For the four stars for
which we can perform a fit of the oscillation frequencies, we find that
the mode lifetimes depend on frequency and on degree of the mode, which
quantitatively agrees with theoretical predictions.
Key words: stars: latetype  stars: oscillations  methods: observational  techniques: photometric
1 Introduction
The uninterupted long timeseries of highprecision photometric data obtained with the CoRoT satellite (Baglin et al. 2006) opened many new possibilities for redgiant asteroseismology. Red giants have extended turbulent atmospheres in which solarlike oscillations are stochastically excited and intrinsically damped. These oscillations have periods of a few hours, and can have mode lifetimes ranging from days to hundreds of days (Dupret et al. 2009). The mode lifetime of a stochastically excited mode is the characteristic time it takes for an excited mode to damp exponentially. This is an important parameter as it is directly related to the mode inertia and therefore the cumulated work integral. From this work integral we can derive in which part of the star the oscillations are driven and damped (see Dupret et al. 2009, for more details). Dupret et al. (2009) also predict that mode lifetimes depend on the degree of the mode and on frequency. This dependence follows similar trends for stars in different evolutionary phases, although the exact variation in the mode lifetimes with frequency changes depending on the evolutionary phase of the star. This frequency and mode dependence of the mode lifetime has not yet been observed in red giants and will be investigated here.
The mode lifetimes (,
efolding time for the amplitude) can be measured from the width (FWHM) of the frequency peaks ()
of an oscillation mode in the power spectrum as
To measure the mode lifetime it is therefore important to have enough frequency resolution to resolve the linewidth of the oscillation mode. Examples of stars with long (>50 day) and short (<20 day) lifetimes were first shown by De Ridder et al. (2009) and Carrier et al. (2010) and investigated for a larger number of stars by Baudin et al. (2010). Where Baudin et al. (2010) propose an automated procedure to fit the line width and height for the three most prominent oscillation modes in a star, with the purpose of a statistical study of the amplitudes of the modes, we adopt here an individual analysis for a smaller sample of stars, but with longer timeseries. The data we have at our disposal are described in Sect. 2. Because we do not have ``infinitely'' long timeseries and there are considerable gaps in those series, we perform simulations to identify a possible bias on in Eq. (1). This bias will be taken into account in computing the ``unbiased'' mode lifetime. These simulations are described in Sect. 3.
We performed a detailed analysis of four stars. For each star we
then computed an échelle diagram (frequency versus frequency modulo
large frequency separation, Grec et al. 1983).
From this échelle diagram we aim to identify the degree of the
oscillation modes, or at least identify modes with the same degree to
investigate the dependence of the mode lifetime on degree. This is
possible when the observed oscillation modes have highorder (n) and lowdegree ()
and follow the asymptotic relation Tassoul (1980), for which we use here the following equation:
with the oscillation frequency, the large separation between modes of same degree and consecutive orders, which is inversely proportional to the sound travel time through the star. D_{0} is sensitive to deeper layers in the star and to the surface layers. The large separation can be estimated from the relation between the frequency of maximum oscillation power ( ) and as found in an ensemble study of the first long run of CoRoT observations (Hekker et al. 2009):
or similar relations by Stello et al. (2009) and Mosser et al. (2010). We use this relation to check whether we indeed found and not a (sub)multiple of .
In general three small separations are defined: is the spacing between adjacent modes with and , is the spacing between adjacent modes with and , and is the amount by which the modes are offset from the midpoints between the mode on either side. If Eq. (2) holds then , and , although there is evidence that for the majority of redgiant stars is negative (Bedding et al. 2010; Carrier et al. 2010), i.e. the modes are located on the right side of the midpoints between the modes instead of on the left side of the midpoint as is the case for the Sun. This is most likely caused by stellar evolution and the asymptotic relation given above may not be valid for evolved stars.
2 Observations
Figure 1: Colourmagnitude diagram of the field stars observed during LRa01 using the 2MASS J and Ks photometric passbands (black dots) (Deleuil et al. 2009). The blue diamonds represent the redgiant stars also observed during IR, the green triangles represent the redgiant candidates which show clear power excess and the red crosses represent the four stars for which we could fit the individual oscillation modes. 

Open with DEXTER 
For the present analyses we use CoRoT data taken during the Initial Run (IR) and first long run in the anticentre direction (LRa01) in the exoplanet field. For information on CoRoT and technical details about the CoRoT data reduction we refer to Baglin & Chaintreuil (2006). The IR lasted 58 days from 2nd February to 31st March 2007 during which the satellite was pointed towards a field with the centre at equatorial coordinates (, ) = (102.60, 1.7). During the 131 day long LRa01, which lasted from October 2007 to March 2008, the satellite observed a field around equatorial coordinates (, ) = (102.72, 0.2). The typical timestep of the observations is 32 s, although some targets are observed at 512 s intervals.
Visual and nearIR photometry are available for all stars in the IR and LRa01 field (Deleuil et al. 2009). These colours are affected by reddening, but as shown by Bessell & Brett (1988) nearIR colours are least affected and provide a first estimate of the spectral type. A J  K versus K colour magnitude diagram is shown in Fig. 1.
We use these colours to distinguish between mainsequence and redgiant stars. In this way we found 305 redgiant candidates based on colours (blue diamonds in Fig. 1), which are observed during both the IR and LRa01. We computed a combined timeseries for these candidates, in which we ``merged'' the timeseries. Simulations to investigate the effect of merging the IR and LRa01 datasets are described in Sect. 3. Then we computed power spectra and global fits to them. These fits are composed of a power law to account for the background of the spectra due to effects such as granulation, and a Gaussian fit to the oscillation power excess, similar to the fitting described by Hekker et al. (2009). For 24 stars (green triangles in Fig. 1) we found clear power excess and we investigated the oscillations in more detail. For four stars the noise level was low enough to perform fitting of individual oscillation modes. These four stars are indicated with red crosses in the colourmagnitude diagram in Fig. 1 and the detailed analysis of these stars is described in Sect. 4.
Figure 2: Window functions for a 394 day timeseries (black solid line), the same timeseries with a large gap, i.e., no data between day 58 and day 263 (red dashed line) and for the same timeseries with the data of the first 58 and last 131 days ``merged'' (green dotdashed line). The window functions are oversampled by a factor of ten for visual purposes. 

Open with DEXTER 
3 Simulations
The large gap of 205 days between the initial run (58 days) and the first long run in the anticentre direction (131 days) has a large influence on the window function (Fig. 2) by introducing prominent sidelobes, and these have to be taken into account in the analysis of the data. To investigate how we should treat the gap and what effect the gap has on the resulting parameters obtained from a fit to the power spectrum we performed simulations using the simulator described by Chaplin et al. (1997). First we performed simulations of noisefree timeseries which cover the full timespan of 394 days with 512 s cadence. Each timeseries has one stochastically excited oscillation mode at a frequency in the frequency range of red giants observed with CoRoT, i.e., 5, 15, 25, ...., 95 Hz and for each frequency we simulated ten different mode lifetimes of 5 up to 95 days. For each combination of frequency and mode lifetime we simulated 1000 realisations.
In these simulations of the full timespan we introduced a gap by removing all data between day 58 and day 263, but not altering any of the time stamps. In addition, we also merge the timeseries data of the first 58 days and last 131 days by removing the gap and altering the timestamps (see Fig. 2 for the window functions of these scenarios). The latter is an option for red giants as these have stochastically excited modes with finite mode lifetimes. These stochastic modes do not have a coherent phase and one may therefore in principle merge data sets with different epochs.
We analysed for all three cases, i.e., ``full'', ``gap'' and ``merge'',
the individual realisations. In all instances we fitted the power
spectrum using a maximum likelihood technique (e.g., Anderson et al. 1990), where the final model M
used for comparison is a convolution of the model with the power
spectrum of the observed window function of the data set, normalised to
unit total area (see Eq. (4)). This takes the redistribution of power caused by gaps in the data into account.
In this model, each of the j oscillation peaks are fitted with the amplitude (A), the central frequency ( ) and FWHM (B) as free parameters. To maximise the possibility of obtaining the best values, i.e., the solution of the global instead of local maximum, we use different input parameters for the linewidth and height of the frequency peak. From the total of twelve fits performed for each timeseries we choose the one with the maximum likelihood. Before we used the result we also tested if we could consider the fitted function as a good approximation of the parent distribution by means of reduced figureofmerit function (MERIT , see Anderson et al. 1990, and references therein):
with O_{i} the observed power at frequency i and M_{i} the model at the same frequency and f the number of degrees of freedom, i.e., the difference between the number of independent data points in the fit and the number of free parameters in the model. A fit with MERIT less than one can be considered a good fit to the data.
Figure 3: Average of the lifetimes observed for 1000 simulated realisations as a function of their input lifetime. The ``full'', ``gap'' and ``merge'' cases are indicated with red asterisks, blue triangles and green diamonds respectively. Errors on the mean are plotted, but are in most cases smaller than the symbol size. The black dotted line is the onetoone relation. 

Open with DEXTER 
The results immediately showed that the frequency at which the oscillation mode occurred did not have any influence on the fitting and the resulting lifetime. Therefore we computed the lifetime of each realisation at one frequency (35 Hz) and computed the mean of these results. The standard error on the mean of the sample is computed as with the standard deviation of the sample and N the number of measured lifetimes over which the mean is computed. These are shown in Fig. 3. This figure clearly shows that the measured lifetimes are equal to or lower than the lifetimes that were used for the simulations and that the offset from the onetoone relations is larger for the ``gap'' and ``merge'' cases than for the ``full'' case. For the ``full'' timeseries the results become less reliable above 60 days while the ``gap'' and ``merge'' results are consistent only for lifetimes below 20 days.
The limited timespan of the observations compared to the mode lifetimes could be the cause of the discrepancy between our results and the input values. To investigate this further, we repeated the simulations and analysis for timeseries with half (``short'') and twice (``long'') the timespan of the initial simulations, i.e., 197 and 788 days and without a gap. We compared the resulting average lifetimes computed from the analysis of individual realisations. These are shown in Fig. 4. Indeed we see that for the ``long'' timeseries the results are consistent with or overestimating the input for all mode lifetimes investigated here, while for the ``short'' timeseries the discrepancy between input and observed lifetime start to emerge at lifetimes of about 25 days. This confirms that the total time span of the observations (T) influences the measured mode lifetime (), and that we can only measure `unbiased' mode lifetimes when (roughly) . We note that the ``gap'', the ``merge'' and the ``short'' timeseries, which have similar effective timespans, all show similar results.
The results for the ``gap'' and ``merge'' timeseries are approximately equal for lifetimes up to 35 days. For higher input lifetimes the mean of the observed lifetimes for the ``merge'' timeseries are systematically lower than for the ``gap'' timeseries. The reason for this increasing difference is probably due to the different frequency resolutions and might therefore vary as a function of the (fractional) length of the gap.
Figure 4: Same as Fig. 3 but here the red asterisks, blue triangles and green diamonds indicate the ``full'', ``long'' and ``short'' cases, respectively. 

Open with DEXTER 
Figure 5: Same as Fig. 3 for the ``gap'' (triangles) and ``merge'' (diamonds) cases with 1.5, 1, 0.5, 0.25 and 0.125 times the gap in the observations. These are indicated in red, black, green, blue and yellow, respectively. The results for the full timeseries are indicated with the light blue asterisks. 

Open with DEXTER 
To test the influence of the length of the gap on the resulting lifetimes, we repeated the analyses as described above for timeseries with a total length of 394 days, and different gaps of 1.5, 0.5, 0.25 and 0.125 times the length of the gap in the observed data. The results are shown in Fig. 5. From these simulations it becomes clear that the deviations of the results shown for the ``gap'' and ``merge'' timeseries increase with increasing length of the gap, i.e. they depend on the effective timespan of the data. We see that for timeseries with a shorter effective timespan we can measure only shorter ``unbiased'' mode lifetimes, which is in agreement with the results from the ``short'' and ``long'' timeseries. Furthermore, the length of the gap in the observed timeseries (205 days) is longer than the lifetimes considered here. This changes when we investigate the same lifetimes for data with shorter gaps. Nevertheless, we find smooth trends between the observed and input lifetimes for all gaps and mode lifetimes considered here.
We also investigated the influence of the heighttobackground ratio of the oscillation modes on the measurement of the mode lifetime. We did this with simulations of ``merge'' timeseries, similar to the observed timeseries we analysed, on which we added noise such that the heighttobackground ratio (H/B) of the oscillation mode is of the order of 10, 25 and 50. In these simulations the noise has been added to the simulated timeseries, while the height of the oscillation mode (H) was derived from the rootmeansquare of the flux (A) and the input mode lifetime () as . Results for these simulations with different H/B are shown in Fig. 6. These results showed that the results are to a large extent insensitive to H/B, although the inclusion of noise in the simulations provided different results compared to results from noisefree simulations (see Fig. 3). These differences are most likely due to the noise and due to the addition of an extra free parameter (background level) in the fitting.
Figure 6: Top: same as Fig. 3 for the ``merge'' case with heighttobackground ratio (H/B) of the order of 10, 25 and 50 indicated in red, green and blue, respectively. The black solid line indicates the polynomial used to obtain the ``unbiased'' lifetime. Bottom: the difference between the ``observed'' lifetimes and the fit. 

Open with DEXTER 
From these simulations we conclude that we need to apply a correction
to the ``observed'' mode lifetime to obtain the ``unbiased'' lifetime
whether we leave the gap in the data or we merge the timeseries of the
initial and long run. To avoid selection of frequencies which have high
power due to the more severe sidelobes in the window function of the
data with a gap, we have performed the analysis on observed power
spectra of merged timeseries. The lifetimes of the observed modes then
need to be corrected. For this correction we use a third order
polynomial fit to the lifetimes measured from the simulations including
noise:
with y the observed lifetime and x the input lifetime as shown in Fig. 6. This shows that the fit is not exact, but that the deviation from the ``observed'' lifetimes is less than five days for all lifetimes and H/B. In the analysis of the observed stars we will use this correction as a lookup table to compute ``unbiased'' lifetimes for observed lifetimes between 5.4 and 35 days, i.e. the range of observed lifetimes over which the correction is computed. For any observed lifetimes outside this range we can infer only lower or upper limits.
Figure 7: Top: power spectrum for CoRoT 102732890 (black solid line), with a fit for nine oscillation modes (blue dashdot line) and a fit with ten oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit with nine frequencies and ten frequencies (same colours and linestyles as in the top panel). The value of s for this frequency interval is 8.1. See text in Sect. 4 for a description of s. 

Open with DEXTER 
Table 1: Parameters for the fit to CoRoT 102732890.
4 Individual stars
From the redgiant candidates selected from nearIR colours and power excess as described in Sect. 2 we have four candidates for which we can fit the individual oscillations. The fitting is performed to the nonoversampled power spectra in a similar way as for the simulations (Sect. 3). The fitting is performed in two steps. First the background, for which we use the power law fit to account for effects such as activity and granulation as described in Sect. 2 is fitted. Then, we compute a global maximum likelihood fit to all oscillation modes. In this fit the previously determined background is kept fixed, while the frequency, height and line width of each oscillation mode and the noise level are free parameters. No dependence of parameters on the degree of the mode and thus on the azimuthal order is included in the fitting. This is both because we do not know the degree of the modes and because we assume that the stars are slow rotators and the rotational splitting is negligible.
The selection of the oscillation modes is based on a statistical test of the binned power spectrum. For this test we bin the power spectrum over intervals of three frequency bins. Then we compute the probability of the power in the binned power spectrum to be due to noise. This probability is computed using a distribution in which we take the width and number of the bins into account in the degrees of freedom. Frequencies at which the probability of the power not being due to noise is larger than 80% are selected as candidate oscillation frequencies. After performing the fit, we also verify that in the ratio of the observed to the fitted spectra no prominent peaks are left. Therefore we compute the relative height (s) for the investigated frequency range for which the probability of observing at least one spike with this height due to noise is 10%, following the formulation by Chaplin et al. (2002) and references therein.
From the resulting height (H) and width () for each oscillation frequency we also compute the rootmeansquared amplitude ( ). This parameter is a measure of the total power of a mode and reflects the balance between the damping and the excitation of the mode. H and are not independent and this is taken into account in the computation of the error in A using a correlation coefficient of 0.9 following the method described by Chaplin et al. (2000).
4.1 CoRoT 102732890
CoRoT 102732890 is a star with frequencies showing maximum oscillation power at about 26 Hz. Our statistical test indicated nine frequency intervals to have more than 80% probability to be due to signal. We fit for these frequencies and the results are shown in Fig. 7 together with the ratio of the observed power spectrum to the fit. Indeed, in the fitted region no clear oscillation modes are present anymore. The MERIT is also below one (0.79), which indicates a good fit. However, one could argue that the power around 22 Hz is not due to one oscillation mode but due to two modes. Therefore, we also perform a fit with ten frequencies. According to the MERIT this is also a good fit (0.79). The maximum likelihood is (slightly) higher for the fit with ten frequencies compared to the fit with nine frequencies. We perform a likelihood ratio test (e.g., Appourchaux et al. 1995) to see if the increase in the maximum likelihood justifies the addition of three free parameters. In this particular case we see that there is a 35% probability to get this increase without the addition of an oscillation mode. From this probability we cannot draw any firm conclusions and therefore we continue the analysis using both nine and ten frequencies. The values used for the fits are listed in Table 1 together with the computed amplitude and ``unbiased'' mode lifetime.
From the observed and Eq. (3), we expect a large separation of about 3.2 Hz. Indeed when we plot the ten frequencies in an échelle diagram with a folding frequency of 3.15 Hz, the different mode sequences become visible (see Fig. 8).
Inspection of Fig. 8 shows that the modes are located in two roughly vertical regions. The interpretation of this is not unambiguous. One could argue that the broad ridge consists of two ridges, namely the and ridge, which seems to be the case for the fit with ten frequencies (black symbols in Fig. 8). We could also interpret the broad ridge as the ridge where we see more scatter due to less efficient trapping of the modes, (see e.g., Bedding et al. (2010), and Dupret et al. (2009) for more details on mode trapping). The latter might be the case for the fit with nine frequencies as the two bottom black symbols of the broad ridge are then replaced by the red point.
Figure 8: Echelle diagram of CoRoT 102732890, using ten frequencies in black symbols, with a size proportional to the height of the fitted oscillations. When the two modes at 22 Hz are replaced by one mode then the two lower left black symbols are replaced by the red symbol. 

Open with DEXTER 
4.2 CoRoT 102788308
CoRoT 102788308 has a strikingly clear and regular pattern around a frequency of about 50 Hz. As seen in Figs. 9 and 10 for the power spectrum and échelle diagram, respectively, and Table 2 for the fitted parameters. The MERIT of the global fit to the power spectrum is 1.2, which indicates that this is a reasonable fit to the data. When looking at the ratio of the observations to the fit we indeed still see some structure, with a main feature at about 46.3 Hz. At this frequency we see for many CoRoT stars the fourth harmonic of a one day feature. Therefore we do not fit this structure in this case, although we cannot exclude that part of this signal is due to stellar oscillations. In addition, the peaks at frequencies about 35, 40, 65 and 68 Hz might be due to oscillations, but the height to background ratio is not large enough to confirm that and therefore these are not taken into account in the fit.
Figure 9: Top: power spectrum for CoRoT 102788308 (black solid line), with a fit for nine oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The dashed horizontal line indicates s (see text for more details). 

Open with DEXTER 
From the of about 50 Hz and Eq. (3), we expect a of about 5.4 Hz. Indeed we find = 5.3 Hz by folding the frequencies. The échelle diagram shows two very clear ridges, i.e., the and ridges (although the actual degree of the ridges is not identified) with one point outside the central ridge in Fig. 10. This mode has most likely been shifted due to an avoided crossing and in that case the centre ridge is the = 1 ridge.
4.3 CoRoT 102762620
Table 2: Parameters for the fit to CoRoT 102788308.
CoRoT 102762620 shows oscillation frequencies at slightly lower frequencies than CoRoT 102732890. With the statistical test of the smoothed power spectrum we find oscillations roughly between 17 and 27 Hz (see Fig. 11 and Table 3 for the fitting parameters).
From the frequency of maximum oscillation power of about 22 Hz we predict a of about 2.8 Hz. From the observations we find a consistent value of 2.7 Hz for this star. This results in one narrow ridge and a second ridge which consist either of two ridges relatively close together or one ridge with larger scatter.
4.4 CoRoT 102767771
CoRoT 102767771 shows oscillations at about the same frequencies as 102732890 (see Fig. 13 and Table 4 for the fitting parameters), although the power of the oscillations in this star is less than for 102732890. We find a similar value for of 3.1 Hz (Fig. 14) as for 102732890, which is in line with the predictions. The échelle diagram shows two clear ridges and only the mode with lowest frequency seems to be located outside the right ridge. From the current analysis we cannot say whether this ridge is curved and this mode is located on the ridge, whether the frequency is altered due to avoided crossing, or whether this is induced by noise.
5 Discussion and conclusions
5.1 Simulations
Simulations of long (a few hundred days) timeseries of data, with different timespans and different length of gaps reveal that for short mode lifetimes we can recover the input lifetime of the mode, irrespective of the presence of a gap in the data, or whether we merge the data. For longer lifetimes we see a general under estimation of the ``observed'' mode lifetimes compared to the input lifetimes, in which the ``observed'' mode lifetime seems to reach a maximum value. The investigations of longer and shorter timeseries reveal that for roughly the observed lifetimes are consistent with the ``unbiased'' lifetime. For longer lifetimes a correction has to be applied, although the uncertainty in the resulting values increases with increasing lifetimes. We also investigated the trends in the observed lifetimes as a function of the length of the gap and whether merging the data has a large influence. As expected the shorter the gap the better. Furthermore, the results for the data with a gap and merged together differ for longer lifetimes, but as a correction is necessary in both cases, we concluded that we could merge the data of the initial run and second long run of CoRoT observations. This has the advantage that sidelobes of the window function are avoided.
Figure 10: Echelle diagram of CoRoT 102788308, with symbol sizes proportional to the height of the fitted oscillations. The dashed lines indicate the same ridge with the dasheddotted line the midpoint between them. 

Open with DEXTER 
Figure 11: Top: power spectrum for CoRoT 102762620 (black solid line), with a fit for nine oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.0 (horizontal dashed line). 

Open with DEXTER 
5.2 Observations
For four redgiant stars we have data from the initial and second long run of CoRoT at our disposal with high enough signaltonoise ratio to perform a global fit to the individual oscillation frequencies. The oscillation modes were identified using a statistical test of the binned power spectrum. In this global fit the frequency, linewidth and height are free parameters for each identified oscillation mode. From the parameters obtained from the fitting we also compute the rootmeansquared amplitude and the mode lifetime. The latter includes a correction using the third order polynomial obtained from the averaged mode lifetimes of the simulations including noise of merged data.
5.2.1 = 2 modes
Interestingly for none of the four stars did we find a clear = 2 ridge in the échelle diagram. For 102732890 and 102762620 one ridge is broad and could be interpreted as an and ridge but for 102788308 and 102767771 only two ridges, separated by 0.5 are present. The nondetection of the modes could possibly be due to a lack of frequency resolution, in which case the and modes are interpreted as one oscillation or due to the lack of height in the power spectrum. We investigate both possibilities.
Table 3: Same as Table 2 but for parameters of the fit to CoRoT 102762620.
From the individual oscillations observed for 102788308, we can obtain the differences between consecutive frequencies with reasonable accuracy. We find here two typical differences: 2.4 and 2.8 Hz. This means that the amount of offset of the modes from the midpoints between the = 0 modes ( ) is of the order of 0.2 Hz. If we interpret the centre ridge as ridge (Fig. 10) then this ridge is mostly located on the left side of the midpoint (dotdashed vertical line in Fig. 10) of the modes (dashed vertical lines in Fig. 10). The position of the ridge on the left side of the midpoint between modes might indicate that this star is in an earlier evolutionary state and that it seems to be justified to use the prediction of the small separation from the asymptotic approximation (Eq. (2)). When we use this approximation and the estimate of of the observed frequencies, we find that should be of the order 0.6 Hz. The frequency at 57.2 Hz lies by this amount outside the centre ridge but no other convincing features are present, although the frequency resolution of the dataset is sufficient to resolve this separation.
Figure 12: Echelle diagram of CoRoT 102762620, with symbol sizes proportional to the height of the fitted oscillations. 

Open with DEXTER 
Figure 13: Top: power spectrum for CoRoT 102767771 (black solid line), with a fit for eight oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.1. 

Open with DEXTER 
Table 4: Same as Table 2 but for parameters of the fit to CoRoT 102767771.
Also for 102767771 we find evidence for only two ridges and we seem to miss the ridge. In this case this might be caused by a combination of a lower signaltonoise ratio and the low frequencies of the oscillations. At these low frequencies we expect also a smaller separation between the and modes compared to the 0.6 Hz estimated for 102788308 and the frequency resolution of the spectrum is therefore more critical to detect them, although the frequency resolution of 0.06 Hz should be sufficient. If the modes are very wide, i.e., short mode lifetimes or unresolved mode pairs, a single observed power excess hump might actually consist of two modes. Fitting the two modes as one will overestimate the line width considerably and the correction to compute the ``unbiased'' mode lifetime (Eq. (6)) might also no longer be valid, as this is determined for a single oscillation mode.
The theoretical computations by Dupret et al. (2009) predict that the trapping of modes is more efficient than the trapping of the modes due to the increased size of the evanescent region for modes of higher degree. Due to this more efficient trapping the modes are predicted to have similar heights in the power spectrum as radial modes. In this prediction the visibility effects are not taken into account. For the Sun it is known that modes with degrees , 1, 2 have relative heights 1.0, 1.5 and 0.5. These relative heights depend on limb darkening and might be different for red giants compared to the Sun. For the stars investigated here we computed the mean of the fitted heights per ridge and the ratio of these mean values. The results are listed in Table 5 and cover a considerable range (note that the inverse of the ratio should also be considered as the degree of the ridges ``A'' and ``B'' is unknown). This might indicate that the relative visibilities of modes with different degrees in red giants are different from that of the Sun and this might cause the nondetection of modes, although this would be in contradiction with red giants for which modes have been observed (e.g., De Ridder et al. 2009; Carrier et al. 2010).
It might also be possible that the lifetimes of the modes are long and that these modes are not resolved. Unresolved nonradial modes with long lifetimes have smaller heights than the closest radial mode (Dupret et al. 2009). With longer observing runs, these modes can be resolved, which will improve the chances to detect them.
Despite the difficulties to detect modes encountered for the stars investigated here, these modes have been observed in a considerable number of (but not all) redgiant stars using both CoRoT and Kepler data with shorter timespans (Bedding et al. 2010; De Ridder et al. 2009; Huber et al. 2010; Kallinger et al. 2010a; Carrier et al. 2010). The reason why the modes are not observable in the stars investigated here is not entirely clear. The fact that for some stars modes are observed and that for other stars, such as the ones investigated here, we cannot observe these modes, might also contain information about their internal structures. The small number statistics of only four stars does not allow us to draw any firm conclusions at this stage.
5.2.2 Mode lifetimes
Theory predicts that mode lifetimes depend on frequency and on the degree of the mode. In general the mode lifetimes decrease for modes with increasing frequencies and are longer for modes with higher degrees (Dupret et al. 2009). We test this prediction using the four stars we analysed in detail, by plotting the linewidth of each mode as a function of frequency and perform a linear fit (see Fig. 15). For three of the four stars there is a clear increase of the linewidth with frequency present, i.e., the mode lifetimes decrease with increasing frequency.
Also a considerable number of lifetimes observed in the four stars are in agreement with the scaling prediction presented by Chaplin et al. (2009). For the effective temperatures mentioned in Table 5 lifetimes of the order 1115 days can be expected. We indeed see in Fig. 16 that many results are present in this regime.
Furthermore, we investigate the dependence of lifetimes on the degree of the modes (see Fig. 16). The most obvious example of this dependence is present in 102767771 where all modes in the left (``A'') ridge in the échelle diagram in Fig. 14 have lifetimes shorter than 20 days, while all modes in the right (``B'') ridge have lifetimes longer than 20 days. Although we could not identify the degree of the ridges ``A'' and ``B'', this result is consistent with predictions by Dupret et al. (2009), namely that modes with different degrees can have systematically different mode lifetimes. For the other stars the dependence of the mode lifetime on degree is not so clear, although the two longest mode lifetimes for 102788308 are both present in the ridge indicated with the dashed vertical line in Fig. 10 (ridge ``A''). The three modes with the same lifetime of 58 days (102762620) are also located in the same ridge.
Table 5: Global parameters for the four stars.
From the échelle diagrams we were not yet able to identify the degrees of the modes. Because Dupret et al. (2009) predict that nonradial modes have longer lifetimes than radial modes, we use that as additional information to identify the degree of the modes. For 102767771 this would indicate that the right ridge in Fig. 14 (ridge ``B'') is the ridge and the left ridge (``A'') the = 0 ridge.
The investigation of mode lifetimes as performed here on four CoRoT stars will benefit from longer uninterrupted timeseries on a larger number of stars. This will be possible in the near future. The Kepler satellite will be taking data of a preselected set of stars during its entire lifetime of nominally 3.5 years. This increased timespan with only short interruptions and the lower noise level of the photometric timeseries compared to CoRoT or any other facility currently in operation, will provide data which will be even better suited for investigating the mode lifetimes of solarlike oscillations in redgiant stars. Analysis of the currently available Kepler data with a timespan of 230 days is in progress.
Figure 14: Echelle diagram of CoRoT 102767771, with symbol sizes proportional to the height of the fitted oscillations. 

Open with DEXTER 
Figure 15: Linewidth of the oscillation modes of the four stars studied in detail as a function of frequency. The dashed lines indicated linear fits to the results. 

Open with DEXTER 
Figure 16: Mode lifetimes for ridge ``A'' and ``B'' in the échelle diagrams of the four stars separately. The degree of the ridges are not known and the assignment of ``A'' and ``B'' to a ridge is arbitrary. The symbols are the same as in Fig. 15. The oscillation modes for which we could only find a lower limit are plotted with a lifetime around 100 days, but with a small offset. An offset is also introduced in two of the three modes with a lifetime of 58.3 days (102762620). 

Open with DEXTER 
5.2.3 Stellar parameters
For the four stars we investigated in detail we also compute the asteroseismic mass and radius using , and , as described by Kallinger et al. (2010b). For this determination we need the effective temperatures. The effective temperatures have been derived using ExoDat photometric data (Deleuil et al. 2009) using the relations of Alonso et al. (1999), as described by Baudin et al. (2010). Because of large, systematic discrepancies between the reddening values along the various lines of sight estimated from the extinction maps of Dobashi et al. (2005) and Rowles & Froebrich (2009), a representative value, A_{V} = 0.6 mag, was adopted for all stars. In view of the quite uncertain reddening, and although we note that there is a good agreement between the values derived from optical and nearIR colour indices, only the temperatures based on 2MASS data are considered for the determination of the mass and radius. The values of these stellar parameters are listed in Table 5.
5.2.4 Comparison with models of Dupret et al. (2009)
The frequency of maximum oscillation power provides an indication of the evolutionary state of the stars. The stars 102732890 and 102767771 have a of about 26 Hz. Stars with roughly between 25 and 40 Hz are commonly observed (Hekker et al. 2009) and are to a large extent lowmass Heburning stars in the red clump (Miglio et al. 2009). Both 102732890 and 102767771 are therefore most likely redclump stars. The oscillations in 102762620 are centred around 22 Hz, which is comparable with model E of Dupret et al. (2009). This is a model of a 2 star high in the red giant branch. For 102788308 we find oscillations with Hz. These frequencies are comparable to models B, D and E of Dupret et al. (2009), which represent a 2 star ascending the giant branch, and 3 Hshell and Heshell burning models respectively. The models all have slightly higher masses than the asteroseismic masses derived from the observations. The exact positions of the stars in the HR diagram is as yet unknown, due to the difference in mass between the models and the observations, although we can still conclude that these four stars are likely to be in different evolutionary phases on the redgiant branch.
The frequency pattern of 102788308 is strikingly regular and the linewidths and trend of the linewidth with frequency for this star are different from that of the others. This is most likely due to the different evolutionary state of this star. This will be further investigated using detailed modelling of these stars, which will be presented in subsequent publications.
AcknowledgementsS.H., W.J.C. and Y.E. acknowledge financial support from the UK Science and Technology Facilities Council (STFC). T.M. acknowledges financial support from Belspo for contract PRODEXGAIA DPAC. This research has made use of the ExoDat database, operated at LAMOAMP, Marseilles, France, on behalf of the CoRoT/Exoplanet programme. We would like to thank the referee R. Gilliland for his useful report that improved the paper.
References
 Alonso, A., Arribas, S., & MartínezRoger, C. 1999, A&AS, 140, 261 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Anderson, E. R., Duvall, Jr., T. L., & Jefferies, S. M. 1990, ApJ, 364, 699 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Toutain, T., Gough, D. O., & Kosovichev, A. 1995, in GONG 1994. Helio and AstroSeismology from the Earth and Space, ed. R. K. Ulrich, E. J. Rhodes Jr., & W. Dappen, ASP Conf. Ser., 76, 314 [Google Scholar]
 Baglin, A., & Chaintreuil, S. 2006, in ESA SP1306, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 279 [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA SP1306, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33 [Google Scholar]
 Baudin, F., Barban, C., Belkacem, K., et al. 2010, A&A, submitted [Google Scholar]
 Bedding, T. R., Huber, D., Stello, D., et al. 2010, ApJ, 713, L176 [NASA ADS] [CrossRef] [Google Scholar]
 Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 [NASA ADS] [CrossRef] [Google Scholar]
 Carrier, F., De Ridder, J., Baudin, F., et al. 2010, A&A, 509, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chaplin, W. J., Elsworth, Y., Howe, R., et al. 1997, MNRAS, 287, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Elsworth, Y., Isaak, G. R., Miller, B. A., & New, R. 2000, MNRAS, 313, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Elsworth, Y., Isaak, G. R., et al. 2002, MNRAS, 336, 979 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Houdek, G., Karoff, C., Elsworth, Y., & New, R. 2009, A&A, 500, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 De Ridder, J., Barban, C., Baudin, F., et al. 2009, Nature, 459, 398 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Deleuil, M., Meunier, J. C., Moutou, C., et al. 2009, AJ, 138, 649 [NASA ADS] [CrossRef] [Google Scholar]
 Dobashi, K., Uehara, H., Kandori, R., et al. 2005, PASJ, 57, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Dupret, M., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Grec, G., Fossat, E., & Pomerantz, M. A. 1983, Sol. Phys., 82, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Hekker, S., Kallinger, T., Baudin, F., et al. 2009, A&A, 506, 465 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huber, D., Bedding, T., Stello, D., et al. 2010, ApJ, submitted [Google Scholar]
 Kallinger, T., Mosser, B., Hekker, S., et al. 2010a, A&A, submitted [Google Scholar]
 Kallinger, T., Weiss, W. W., Barban, C., et al. 2010b, A&A, 509, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miglio, A., Montalbán, J., Baudin, F., et al. 2009, A&A, 503, L21 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Mosser, B., Belkacem, K., Goupil, M.J., et al. 2010, A&A, 517, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rowles, J., & Froebrich, D. 2009, MNRAS, 395, 1640 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Chaplin, W. J., Basu, S., Elsworth, Y., & Bedding, T. R. 2009, MNRAS, 400, L80 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
 ... CoRoT^{}
 The CoRoT space mission which was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. Light curves can be retrieved from the CoRoT archive: http://idoccorot.ias.upsud.fr.
All Tables
Table 1: Parameters for the fit to CoRoT 102732890.
Table 2: Parameters for the fit to CoRoT 102788308.
Table 3: Same as Table 2 but for parameters of the fit to CoRoT 102762620.
Table 4: Same as Table 2 but for parameters of the fit to CoRoT 102767771.
Table 5: Global parameters for the four stars.
All Figures
Figure 1: Colourmagnitude diagram of the field stars observed during LRa01 using the 2MASS J and Ks photometric passbands (black dots) (Deleuil et al. 2009). The blue diamonds represent the redgiant stars also observed during IR, the green triangles represent the redgiant candidates which show clear power excess and the red crosses represent the four stars for which we could fit the individual oscillation modes. 

Open with DEXTER  
In the text 
Figure 2: Window functions for a 394 day timeseries (black solid line), the same timeseries with a large gap, i.e., no data between day 58 and day 263 (red dashed line) and for the same timeseries with the data of the first 58 and last 131 days ``merged'' (green dotdashed line). The window functions are oversampled by a factor of ten for visual purposes. 

Open with DEXTER  
In the text 
Figure 3: Average of the lifetimes observed for 1000 simulated realisations as a function of their input lifetime. The ``full'', ``gap'' and ``merge'' cases are indicated with red asterisks, blue triangles and green diamonds respectively. Errors on the mean are plotted, but are in most cases smaller than the symbol size. The black dotted line is the onetoone relation. 

Open with DEXTER  
In the text 
Figure 4: Same as Fig. 3 but here the red asterisks, blue triangles and green diamonds indicate the ``full'', ``long'' and ``short'' cases, respectively. 

Open with DEXTER  
In the text 
Figure 5: Same as Fig. 3 for the ``gap'' (triangles) and ``merge'' (diamonds) cases with 1.5, 1, 0.5, 0.25 and 0.125 times the gap in the observations. These are indicated in red, black, green, blue and yellow, respectively. The results for the full timeseries are indicated with the light blue asterisks. 

Open with DEXTER  
In the text 
Figure 6: Top: same as Fig. 3 for the ``merge'' case with heighttobackground ratio (H/B) of the order of 10, 25 and 50 indicated in red, green and blue, respectively. The black solid line indicates the polynomial used to obtain the ``unbiased'' lifetime. Bottom: the difference between the ``observed'' lifetimes and the fit. 

Open with DEXTER  
In the text 
Figure 7: Top: power spectrum for CoRoT 102732890 (black solid line), with a fit for nine oscillation modes (blue dashdot line) and a fit with ten oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit with nine frequencies and ten frequencies (same colours and linestyles as in the top panel). The value of s for this frequency interval is 8.1. See text in Sect. 4 for a description of s. 

Open with DEXTER  
In the text 
Figure 8: Echelle diagram of CoRoT 102732890, using ten frequencies in black symbols, with a size proportional to the height of the fitted oscillations. When the two modes at 22 Hz are replaced by one mode then the two lower left black symbols are replaced by the red symbol. 

Open with DEXTER  
In the text 
Figure 9: Top: power spectrum for CoRoT 102788308 (black solid line), with a fit for nine oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The dashed horizontal line indicates s (see text for more details). 

Open with DEXTER  
In the text 
Figure 10: Echelle diagram of CoRoT 102788308, with symbol sizes proportional to the height of the fitted oscillations. The dashed lines indicate the same ridge with the dasheddotted line the midpoint between them. 

Open with DEXTER  
In the text 
Figure 11: Top: power spectrum for CoRoT 102762620 (black solid line), with a fit for nine oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.0 (horizontal dashed line). 

Open with DEXTER  
In the text 
Figure 12: Echelle diagram of CoRoT 102762620, with symbol sizes proportional to the height of the fitted oscillations. 

Open with DEXTER  
In the text 
Figure 13: Top: power spectrum for CoRoT 102767771 (black solid line), with a fit for eight oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.1. 

Open with DEXTER  
In the text 
Figure 14: Echelle diagram of CoRoT 102767771, with symbol sizes proportional to the height of the fitted oscillations. 

Open with DEXTER  
In the text 
Figure 15: Linewidth of the oscillation modes of the four stars studied in detail as a function of frequency. The dashed lines indicated linear fits to the results. 

Open with DEXTER  
In the text 
Figure 16: Mode lifetimes for ridge ``A'' and ``B'' in the échelle diagrams of the four stars separately. The degree of the ridges are not known and the assignment of ``A'' and ``B'' to a ridge is arbitrary. The symbols are the same as in Fig. 15. The oscillation modes for which we could only find a lower limit are plotted with a lifetime around 100 days, but with a small offset. An offset is also introduced in two of the three modes with a lifetime of 58.3 days (102762620). 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.