Issue |
A&A
Volume 525, January 2011
|
|
---|---|---|
Article Number | L9 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201015440 | |
Published online | 07 December 2010 |
Letters to the Editor
The universal red-giant oscillation pattern
An automated determination with CoRoT data
1
LESIA, CNRS, Université Pierre et Marie Curie, Université Denis Diderot,
Observatoire de Paris,
92195
Meudon Cedex,
France
e-mail: benoit.mosser@obspm.fr
2
Institut d’Astrophysique et de Géophysique, Université de Liège,
Allée du 6 Août 17,
4000
Liège,
Belgium
3
School of Physics and Astronomy, University of Birmingham,
Edgbaston,
Birmingham
B15 2TT,
UK
4
Institute for Astronomy (IfA), University of Vienna,
Türkenschanzstrasse
17, 1180
Vienna,
Austria
5
Department of Physics and Astronomy, University of British
Columbia, Vancouver, BC
V6T 1Z1,
Canada
6
Instituut voor Sterrenkunde, K. U. Leuven, Celestijnenlaan 200D,
3001
Leuven,
Belgium
7
Institut d’Astrophysique Spatiale, UMR 8617, Université Paris XI,
Bâtiment 121,
91405
Orsay Cedex,
France
8
Centro de Física Computacional, Department of Physics, University
of Coimbra, 3004-516
Coimbra,
Portugal
Received:
20
July
2010
Accepted:
9
November
2010
Aims. The CoRoT and Kepler satellites have provided thousands of red-giant oscillation spectra. The analysis of these spectra requires efficient methods of identifying all eigenmode parameters.
Methods. The assumption of new scaling laws allowed us to construct a theoretical oscillation pattern. We then obtained a highly precise determination of the large separation by correlating the observed patterns with this reference.
Results. We demonstrate that this pattern is universal and are able to unambiguously assign the eigenmode radial orders and angular degrees. This solves one of the remaining problems of asteroseismology, hence allowing precise theoretical investigation of red-giant interiors.
Key words: stars: oscillations / stars: interiors / methods: data analysis / methods: analytical
© ESO, 2010
1. Introduction
Red giants are evolved stars that have depleted the hydrogen in their cores and are no longer able to generate energy from core-hydrogen burning. The physical processes taking place in their interiors are currently fairly poorly understood. Observations with the space-borne mission CoRoT have revealed the oscillation pattern (De Ridder et al. 2009) of many of these stars, which is a crucial step toward probing their internal structure. Before the advent of the CoRoT data, complex oscillation patterns were explained by short-lived modes of oscillation (Stello et al. 2004; Barban et al. 2007). The new era of the space-borne missions CoRoT and Kepler has dramatically increased the amount and quality of the available asteroseismic data of red giants (Hekker et al. 2009; Bedding et al. 2010a; Mosser et al. 2010; Huber et al. 2010). Analysis of the oscillation eigenmodes now allows seismic inferences to be drawn about the internal structure.
Identifying the angular degree and radial order of the eigenmodes represents a first and crucial step in asteroseismic analysis. The values of the eigenfrequencies can be related to the order and degree of a mode with the commonly used asymptotic equation: (1)where νn,ℓ is the eigenfrequency of a mode with radial order n and angular degree ℓ; Δν is the mean value of the large separation (Δν ≃ νn + 1,ℓ − νn,ℓ), and δν0ℓ a second-order term, or small separation, which depends on the mode degree. This form, which is similar to the original expression developed for the Sun 30 years ago (Tassoul 1980), is useful for analysing the observations and performing the complete mode identification. It assumes that δν00 equals 0. The parameter ε comprises two parts: the offset due to the mode propagation in the uppermost layers of the star, and the second-order term of the asymptotic approximation, which is sensitive to the gradient of sound speed in the stellar interior.
In the absence of accurate determinations of the individual mode frequencies, the global seismic parameters used in the asymptotic expression above are important indicators of the physical parameters of the star. The large separation gives a measure of the mean stellar density; the small separation δν0ℓ describes the stratification of the central regions. Unfortunately, the methods currently used to determine the global oscillation parameters suffer from various sources of uncertainty (Hekker et al. 2009; Huber et al. 2009; Mathur et al. 2010). First, the stochastic excitation of the modes gives rise to variability in the amplitudes, resulting in an apparently irregular comb structure; second, the finite mode lifetime blurs the estimates of the eigenfrequencies; third, estimates are affected by the stellar noise and granulation signal superimposed on the oscillations (Mosser & Appourchaux 2009). In fact, simulations have shown that the impact of realization noise on the measurement of the large separation Δν, can be much stronger than the background noise for red giants (Hekker et al. 2010).
In an analysis of a subsample of Kepler red giants, Huber et al. (2010) have shown the regularity of the oscillation spectra of these stars. In this paper, we show that all red giants have a regular pattern, as modelled recently by Montalbán et al. (2010). We propose a method that allows us to tag all the modes with their appropriate radial order and angular degree, regardless of whether there are any of the perturbing effects described above.
2. Method
The method of mitigating the effects of realization noise uses Eq. (1) in a dimensionless form: (2)In a departure from the previous practice, we have assumed that ε obeys a scaling law ε = A + Blog Δν, as derived from the observation of thousands of CoRoT targets (Mosser et al. 2010) and as observed by Huber et al. (2010). This is justified by the observation that scaling laws apparently govern all global asteroseismic parameters (Hekker et al. 2009; Stello et al. 2009; Bedding et al. 2010a; Mosser et al. 2010) and is equivalent to assuming that the underlying physics of ε varies with the global stellar parameters. As the mixed nature of dipole modes (ℓ = 1) is more pronounced, we did not include them in the template, but only those doublets that correspond to the eigenmodes with even degrees (νn − 1,2 and νn,0), with equal amplitudes. As a first guess, we set the small separation d02 at −0.14 and then allowed it to vary with the value of Δν according to the same relationship as given for ε. For constructing the peaks of the template, we also used the scaling laws of the Gaussian excess power derived by Mosser et al. (2010). We assumed that the mode lifetime varies as Δν-1 and used mode widths equal to about a few percent of Δν. Finally, we stress that no background model is needed.
![]() |
Fig. 1 CoRoT red giant power spectra stacked into an image after sorting on the large separation. One line corresponds to one star. a) The first sorting is based on the large separation before any correction. b) The blurred aspect disappears once the correction has been performed and reveals the clear comb structure common to all red giants. The vertical lines at 23.2 μHz are the signature of the low-Earth orbit (Auvergne et al. 2009). |
![]() |
Fig. 2 Spectra presented in Fig. 1b rearranged to have the dimensionless frequency on the abscissa and the mean large separation on the ordinate. The hyperbolic branch is here the signature of the low-Earth orbit. |
The measurement of the large separation is performed in two steps. First, an initial-guess value Δνguess of the large separation is computed by an automated pipeline (Mosser & Appourchaux 2009). This is used to form the initial synthetic template to correlate with the real spectrum. The best correlation between the observed and synthetic spectra then provides the corrected value of the large separation. The template was iteratively adjusted by varying its parameters to maximize the correlation.
In Fig. 1, we show the results with all high signal-to-noise CoRoT data (Mosser et al. 2010). In both cases the graphs show the spectra arranged in strips with the colour representing the strength of the signal, as in Gilliland et al. (2010). The spectra are sorted by increasing large separation with the smallest large separation at the top of the plots: in the upper plot we use the output from a conventional pipeline and in the lower one the corrected value. The remarkably regular structure within the oscillation spectra in the lower plot reveals the signature of comb-like structure of the asymptotic relationship in Eq. (2) already reported in De Ridder et al. (2009), Carrier et al. (2010), Bedding et al. (2010a), and Huber et al. (2010). Furthermore, it validates the scaling law in ε included in the reference template. The global agreement of all high signal-to-noise spectra of bright targets with the synthetic pattern (Fig. 2) shows that these oscillation patterns are homologous and that the red-giant oscillation pattern is universal.
We further found that the template is significantly improved if it takes account of the linear dependence of the large separation in frequency, expressed by the degree-dependent gradient αℓ = (dlog Δν / dn)ℓ: (3)with νmax the frequency of maximum oscillation amplitude. The corrected values of Δν are derived from this template. The values of the 12 free parameters that account for the variations in frequency of the parameter ε, of the small separations d0ℓ, and of the gradients αℓ as derived from the fits to more than 6000 eigenmodes (Fig. 3) are given in Table 1. The value of ε, defined modulo 1, is fixed thanks to the extrapolation to the solar case. As noted by Carrier et al. (2010) and Huber et al. (2010), the small separation d01 is negative. All these fits are consistent with Kepler results.
The average value of the correction from Δνguess to Δν is of the order of 2.5%. For the largest giants, with the lowest values of Δν and the lowest ratio νmax / Δν, the correction can be as high as 6%. The absolute difference |Δνguess − Δν| can be 10 times the estimate of the stellar noise contribution. At low frequency, with an observing run not much longer than the mode lifetimes (Baudin et al. 2010), the realization noise dominates the background noise and the mean accuracy of the determination of Δν is uniform, at about 0.015 μHz. We are aware that any bias in the ε(Δν) input relation will induce a bias in the results. We therefore took care to insure that the iterative process is unbiased. Analysis of synthetic data indicates that the precision gained with this method is about 10 times better than obtained with conventional methods (Hekker et al. 2010).
Fits of the ridges.
![]() |
Fig. 3 Complete identification of the ridges, automatically derived from the eigenfrequencies extracted with a height-to-background ratio greater than 3. Each colour corresponds to a different mode degree (radial modes in red, dipole modes in dark blue, ℓ = 2 modes in green, ℓ = 3 modes in light blue). The solid grey lines superimposed on the ridges indicate the fits of ε for each radial order (with indication of the radial order n). The fits of d01, d02 and d03 are superimposed on the respective ridges (respectively dash-dot, dot, and dash lines for ℓ = 1, 2, and 3). The dark dashed lines, derived from the scaling law dealing with the oscillation excess power, delineate the region where the modes have noticeable amplitudes (Mosser et al. 2010). |
3. Discussion
This new method based on a simple hypothesis and an automated procedure removes any ambiguity in the identification of the modes (Fig. 3), despite the complexity induced by mixed modes. Mode identification is derived by looking at the closest ridge. In particular, we provide a straightforward determination of the mode radial orders, which were previously unknown. Radial eigenfrequencies are located at (4) Ridges have already been shown in previous works. While De Ridder et al. (2009) and Carrier et al. (2010) looked at single stars separately, Bedding et al. (2010a) and Huber et al. (2010) used manual fine-tuning of the large separation to align the radial modes of a large sample of stars. However, the radial modes were identified in only one third of the spectra by Huber et al. (2010), but they also showed the ridges with varying ε in the folded and collapsed power spectrum.
In most regions of the oscillation spectra, we observe the presence of both radial and non-radial modes. Realization noise causes the height of the individual modes to show considerable variability, but on average, the ratio between the dipole and radial mode height is almost independent of νmax, although there is some reduction in the strength of the dipole mode at very low νmax. We also make it clear that the larger spread of the ridges corresponding to dipole modes (Fig. 3) comes from the presence of many mixed modes, as already noticed (Dupret et al. 2009; Bedding et al. 2010a). The universal pattern makes it easier to identify them, opening up the possibility of exploring the conditions in the inner layers of the red giants.
Despite their low amplitudes and the resulting poor signal, ℓ = 3 modes have been detected in Kepler data on red giants (Bedding et al. 2010a; Huber et al. 2010). Our results represent the first such detection in CoRoT data. This identification gives access to the fine structure of the oscillation spectra, because modes of different angular degree probe different depths within the star. Their detection and complete characterization will first be derived from the universal pattern, then the small differences to this pattern will be exploited to characterize a given object in detail (Miglio et al. 2010).
More than 75% of the red-giant candidates with brighter magnitudes than mR = 13 observed with CoRoT show solar-like oscillations. In the remainder, we observe a high proportion of classical pulsators or of giants with such a large radius that the oscillations occur at too low a frequency for positive detection. In a very limited number of cases at very low frequency, the possible confusion between radial and dipole modes is not clearly solved. This confusion increases toward dimmer targets with lower quality time series. Among the positive detections of bright stars, we did not observe any outliers when performing the correlation with the universal pattern. For this procedure to be effective, we require that the modes have a significant height-to-background ratio. Hence, for all high signal-to-noise targets (Mosser et al. 2010) we are able to derive corrected values for the large frequency spacing.
It is recognized that the majority of the red giants in the CoRoT field of view are in their post-flash helium-burning phase (Miglio et al. 2009). In terms of stellar evolution, demonstration of the universal regular pattern of red giants proves that these red giants have similar and homologous interior structures. On the other hand, despite the agreement of the fit in ε with the solar value, we have verified that the method does not work with subgiants or main-sequence stars (Barban et al. 2009; Benomar et al. 2009; Bedding et al. 2010b; Deheuvels et al. 2010). We explain this by the wide range of evolutionary phases that pertain outside the red clump, which certainly will cause a complex dependence of ε on more parameters than just the large separation. This disagreement reinforces the homogenous properties of red-giant stars. The clearly observed pattern confirms that, on average, the linewidths of the modes are significantly smaller than the separation of the even mode pairs and hence makes short mode lifetimes unlikely (Stello et al. 2004). From the quasi-uniform width of the ridges (Fig. 3), we can see that the lifetimes of the modes increase significantly with decreasing large separation, contrary to Huber et al. (2010). It also suggests that very complex previously observed oscillation spectra may be an artefact of noise.
![]() |
Fig. 4 Variation in stellar mass and radius as a function of the large separation Δν. The vertical line corresponds to the error box in large separation and radius. |
Finally, the better determination of the large separation Δν helps us to enhance the accuracy of the estimates of the stellar mass and radius as done in Mosser et al. (2010), with typical error bars of 12 and 5%, respectively, instead of 20 and 8% (Fig. 4). This accuracy, achieved without the need of stellar modelling (Kallinger et al. 2010), constitutes an important progress compared to the current photometric determination and demonstrates the power of asteroseismology.
4. Conclusion
We have shown with CoRoT observations that the red-giant oscillation spectrum is very regular and can be described by its underlying universal pattern. This was modelled in parallel by Montalbán et al. (2010). As a consequence, the precise measures of the large separation and the scaling relation of the parameter ε allow us to provide an unambiguous detection of the radial orders and angular degrees of the modes. Since the method is able to mitigate the realization noise, we consider it to give the most precise determination of the large separation available.
It remains important to interpret the physical meaning of the scaling law for the term ε in the Tassoul expression. We will have to disentangle the contributions of the surface and of the inner region. This will require investigating the term ε in the context of the second-order corrections of the Tassoul development. Very long-duration observations with Kepler will help for this task.
Despite the uniform aspect of the oscillation spectra, many differences invisible in the global approach are revealed by a detailed analysis of each individual spectrum, as, for example, the modulation with frequency of the large separation (Mosser et al. 2010). Study of this variation will give access to the most accurate analysis of the red-giant interior structure available so far. This summarizes the power of asteroseismology: first, the regular pattern provides the identification of the individual modes; second, the difference to this regular pattern unveils the detailed interior structure. We are confident that the shifts to the regular pattern will be explained by mass, age, and/or metallicity effects.
Similar analysis can be performed for oscillations in subgiants and solar-like stars (Michel et al. 2008; Appourchaux et al. 2008). Thanks to the wide variety of evolutionary stages among those stars, we expect the Tassoul parameter ε to depend on more than the large separation.
Acknowledgments
This work was supported by the Centre National d’Études Spatiales (CNES). It is based on observations with CoRoT. The research has made use of the Exo-Dat database, operated at LAM-OAMP, Marseille, France, on behalf of the CoRoT/Exoplanet program. K.B. acknowledges financial support through a postdoctoral fellowship from the Subside fédéral pour la recherche 2010, University of Liège. F.J.G.P. acknowledges FCT’s grant SFRH/BPD/37491/2007. Y.E. and S.H. acknowledge financial support from the UK Science and Technology Facilities Council (STFC). The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007–2013)/ERC grant agreement No. 227224 (PROSPERITY), as well as from the Research Council of K.U. Leuven grant agreement GOA/2008/04.
References
- Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barban, C., Matthews, J. M., De Ridder, J., et al. 2007, A&A, 468, 1033 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barban, C., Deheuvels, S., Baudin, F., et al. 2009, A&A, 506, 51 [NASA ADS] [CrossRef] [EDP Sciences] [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. 2010a, ApJ, 713, L176 [NASA ADS] [CrossRef] [Google Scholar]
- Bedding, T. R., Kjeldsen, H., Campante, T. L., et al. 2010b, ApJ, 713, 935 [NASA ADS] [CrossRef] [Google Scholar]
- Benomar, O., Baudin, F., Campante, T. L., et al. 2009, A&A, 507, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carrier, F., De Ridder, J., Baudin, F., et al. 2010, A&A, 509, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- De Ridder, J., Barban, C., Baudin, F., et al. 2009, Nature, 459, 398 [Google Scholar]
- Deheuvels, S., Bruntt, H., Michel, E., et al. 2010, A&A, 515, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dupret, M., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131 [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]
- Hekker, S., Elsworth, Y., De Ridder, J., et al. 2010, A&A, accepted [arXiv:1008.2959] [Google Scholar]
- Huber, D., Stello, D., Bedding, T. R., et al. 2009, Commun. Asteroseismol., 160, 74 [NASA ADS] [CrossRef] [Google Scholar]
- Huber, D., Bedding, T. R., Stello, D., et al. 2010, ApJ, 723, 1607 [NASA ADS] [CrossRef] [Google Scholar]
- Kallinger, T., Mosser, B., Hekker, S., et al. 2010, A&A, 522, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mathur, S., García, R. A., Régulo, C., et al. 2010, A&A, 511, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Miglio, A., Montalbán, J., Baudin, F., et al. 2009, A&A, 503, L21 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Miglio, A., Montalbán, J., Carrier, F., et al. 2010, A&A, 520, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Montalbán, J., Miglio, A., Noels, A., Scuflaire, R., & Ventura, P. 2010, ApJ, 721, L182 [NASA ADS] [CrossRef] [Google Scholar]
- Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mosser, B., Belkacem, K., Goupil, M., et al. 2010, A&A, 517, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stello, D., Kjeldsen, H., Bedding, T. R., et al. 2004, Sol. Phys., 220, 207 [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]
All Tables
All Figures
![]() |
Fig. 1 CoRoT red giant power spectra stacked into an image after sorting on the large separation. One line corresponds to one star. a) The first sorting is based on the large separation before any correction. b) The blurred aspect disappears once the correction has been performed and reveals the clear comb structure common to all red giants. The vertical lines at 23.2 μHz are the signature of the low-Earth orbit (Auvergne et al. 2009). |
In the text |
![]() |
Fig. 2 Spectra presented in Fig. 1b rearranged to have the dimensionless frequency on the abscissa and the mean large separation on the ordinate. The hyperbolic branch is here the signature of the low-Earth orbit. |
In the text |
![]() |
Fig. 3 Complete identification of the ridges, automatically derived from the eigenfrequencies extracted with a height-to-background ratio greater than 3. Each colour corresponds to a different mode degree (radial modes in red, dipole modes in dark blue, ℓ = 2 modes in green, ℓ = 3 modes in light blue). The solid grey lines superimposed on the ridges indicate the fits of ε for each radial order (with indication of the radial order n). The fits of d01, d02 and d03 are superimposed on the respective ridges (respectively dash-dot, dot, and dash lines for ℓ = 1, 2, and 3). The dark dashed lines, derived from the scaling law dealing with the oscillation excess power, delineate the region where the modes have noticeable amplitudes (Mosser et al. 2010). |
In the text |
![]() |
Fig. 4 Variation in stellar mass and radius as a function of the large separation Δν. The vertical line corresponds to the error box in large separation and radius. |
In the text |
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.