Issue |
A&A
Volume 520, September-October 2010
|
|
---|---|---|
Article Number | A60 | |
Number of page(s) | 11 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201014944 | |
Published online | 04 October 2010 |
Oscillation mode lifetimes of red giants observed during the initial and first anticentre long run of CoRoT![[*]](/icons/foot_motif.png)
S. Hekker1,2 - C. Barban3 - F. Baudin4 - J. De Ridder2 - T. Kallinger5 - T. Morel6 - W. J. Chaplin1 - Y. Elsworth1
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 signal-to-noise 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
signal-to-noise ratio and possible gaps do have a non-negligible 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: late-type - stars: oscillations - methods: observational - techniques: photometric
1 Introduction
The uninterupted long timeseries of high-precision photometric data obtained with the CoRoT satellite (Baglin et al. 2006) opened many new possibilities for red-giant asteroseismology. Red giants have extended turbulent atmospheres in which solar-like 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 (,
e-folding 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

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 high-order (n) and low-degree ()
and follow the asymptotic relation Tassoul (1980), for which we use here the following equation:
with





or similar relations by Stello et al. (2009) and Mosser et al. (2010). We use this relation to check whether we indeed found


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 red-giant 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: Colour-magnitude 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 red-giant stars also observed during IR, the green triangles represent the red-giant 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 anti-centre direction (LRa01) in the
exo-planet 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 near-IR 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) near-IR 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 main-sequence and red-giant stars. In this way we found 305 red-giant 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 colour-magnitude 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 dot-dashed 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 anti-centre 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 noise-free 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 (


with Oi the observed power at frequency i and Mi 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

![]() |
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 one-to-one 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 one-to-one 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 height-to-background 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
height-to-background 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 root-mean-square 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 noise-free 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 height-to-background 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 look-up table to compute ``unbiased'' lifetimes for observed lifetimes between 5.4 and

![]() |
Figure 7: Top: power spectrum for CoRoT 102732890 (black solid line), with a fit for nine oscillation modes (blue dash-dot 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 red-giant candidates selected from near-IR 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 non-oversampled 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 root-mean-squared 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 |
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 dashed-dotted 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 red-giant stars we have data from the initial and second long run of CoRoT at our disposal with high enough signal-to-noise 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 root-mean-squared 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 non-detection 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 (dot-dashed 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
signal-to-noise 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 non-detection 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
non-radial 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) red-giant 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 11-15 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 non-radial 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 solar-like oscillations
in red-giant 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 Exo-Dat 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, AV
= 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 near-IR 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 low-mass He-burning stars in the red clump (Miglio et al. 2009). Both 102732890 and 102767771 are therefore most likely red-clump 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
H-shell and He-shell 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 H-R 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 red-giant 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 PRODEX-GAIA DPAC. This research has made use of the Exo-Dat database, operated at LAM-OAMP, 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ínez-Roger, 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 Astro-Seismology 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 SP-1306, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 279 [Google Scholar]
- Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA SP-1306, 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 [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 [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://idoc-corot.ias.u-psud.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: Colour-magnitude 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 red-giant stars also observed during IR, the green triangles represent the red-giant 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 dot-dashed 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 one-to-one 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 height-to-background 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 dash-dot 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 |
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 dashed-dotted 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 (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.