Measurements of Stellar Properties through Asteroseismology: A Tool for Planet Transit Studies
Oscillations occur in stars of most masses and essentially all stages of evolution. Asteroseismology is the study of the frequencies and other properties of stellar oscillations, from which we can extract fundamental parameters such as density, mass, radius, age and rotation period. We present an overview of asteroseismic analysis methods, focusing on how this technique may be used as a tool to measure stellar properties relevant to planet transit studies. We also discuss details of the Kepler Asteroseismic Investigation – the use of asteroseismology on the Kepler mission in order to measure basic stellar parameters. We estimate that applying asteroseismology to stars observed by Kepler will allow the determination of stellar mean densities to an accuracy of 1%, radii to 2–3%, masses to 5%, and ages to 5–10% of the main-sequence lifetime. For rotating stars, the angle of inclination can also be determined.
Asteroseismology, a tool for transit studies] Measurements of Stellar Properties through
A Tool for Planet Transit Studies Kjeldsen, Bedding & Christensen-Dalsgaard] Hans Kjeldsen, Timothy R. Bedding and
Jørgen Christensen-Dalsgaard 2008 \volume253 \pagerange1 – 10 \jnameTransiting Planets \editors
Asteroseismology – the study of stellar oscillations – is a relatively new and growing research field in astrophysics. The analysis of frequencies and other properties of stellar oscillations allows us to constrain fundamental parameters of stars such as density, mass, radius, age, rotation period and chemical composition.
Oscillations are found in stars of most masses and essentially all stages of evolution. The amplitudes and phases are controlled by the energetics and dynamics of the near-surface layers. The frequencies are determined by the internal sound-speed and density structure of the star, as well as rotation and (in some cases) magnetic fields. Observationally, the frequencies can be determined with exceedingly high accuracy compared to any other quantity relevant to the internal properties of the stars. Analysis of the observed frequencies, including comparison with computed stellar models, allows determination of the properties of the stellar interiors and tests of the physics used in the model computation, applied under extreme conditions that cannot be matched in terrestrial laboratories.
Rotation induces fine structure in the frequency spectrum, in the form of rotational splitting. The observed frequencies are determined by averages over the stellar interior of the rotation rate, which in general varies with position within the star. By comparing with independent determinations of the surface rotation rate, or from rotational splittings for a sufficient broad variety of modes, information about this variation can be obtained.
Some of the basic parameters that may be obtained from asteroseismology are crucial for understanding the fundamental properties of exoplanet systems. Hence, asteroseismology will potentially be an excellent tool for characterizing planet transit systems and, specifically, it will be able to measure stellar radii to a relative accuracy of 2–3%.
For a recent review on the present state of asteroseismology, see [Aerts \etal (2008), Aerts \etal (2008)] and [Bedding & Kjeldsen (2008), Bedding & Kjeldsen (2008)]. In Figure 1 we show examples of power spectra for four different stars. In each case, we can clearly see the excess energy arising from stellar oscillations. The detailed properties of those oscillations form the observational basis for the asteroseismic analysis.
2 The relation between the stellar properties and the frequencies
To obtain information about stellar properties and the underlying physics, the observed oscillation frequencies must be compared with those computed from models. Major uncertainties affect the treatment of hydrodynamical processes in stellar interiors. These include convective motions, with possible overshoot into the surrounding convectively stable regions, and flows and instabilities related to internal rotation. Furthermore, the evolution of rotation as a star evolves, with internal redistribution and possibly surface loss of angular momentum, has potentially substantial – and highly uncertain – effects on the evolution of the star. Despite those uncertainties in the knowledge of some of the detailed physics, we are able to obtain quite detailed basic properties for a number of general stellar parameters such as density and the amount of hydrogen in the core. The frequencies of stellar p-mode oscillations are related to the sound travel time across the star. Since the sound speed, , is given by
(where the approximation assumes an ideal gas), we are basically measuring the average ratio between pressure and density in the stellar interior.
The exact frequencies of stellar oscillations depend on the detailed structure of the star and on the physical properties of the gas. However, one may use asymptotic theory to derive a simple relation for the mode frequencies (see, for example, [Christensen-Dalsgaard (2004), Christensen-Dalsgaard 2004]). Each oscillation mode is characterized by three integers: the radial order , the angular degree and the azimuthal order . The results of the asymptotic analysis for a non-rotating star give the mode frequencies as
Here, (the so-called large separation) depends on the average stellar density, is sensitive to the sound speed near the core and is sensitive to the surface layers. It is conventional to define , the so-called small separation, as the frequency spacing between adjacent modes with and . These separations are indicated for the Sun in Figure 2, together with the similar quantity . We can further define to be the amount by which modes are offset from the midpoint between the modes on either side. If the asymptotic relation holds exactly, then it follows that .
In practice, the asymptotic relation does not hold exactly, even for the Sun. For example, the large separation depends on . However, we can define average values for and , and then can calculate how those parameters depend on the stellar properties and the stellar evolutionary stage. An example of such a calculation is shown by the lines in Figure 3. In this figure we also show the measured position of the Sun, together with that of Centauri A ([Bedding et al.(2004)Bedding, Kjeldsen, Butler, et al., Bedding \etal 2004]) and Centauri B ([Kjeldsen \etal (2005), Kjeldsen \etal 2005]). This diagram confirms that Cen A has a higher mass than the Sun and is more evolved, while Cen B has a lower mass than the Sun and is less evolved.
The exact positions of the tracks shown in Figure 3 will depend on the detailed physical properties, including chemical composition. However, the basic stellar parameters can still be estimated with some robustness. We have used simulations to estimate the accuracy with which various parameters can be measured, assuming the stellar temperature is already known to within 150–200 K and the heavy-element abundance is known to within a factor of two. We find that
stellar mean densities can be measured to 1%,
stellar radii can be measured to 2–3%,
stellar masses can be measured to 5% and
stellar ages can be measured to 5–10% of the main-sequence lifetime.
3 Stellar Rotation
Equation 2 applies to a non-rotating star, for which the frequency of each mode depends only on the radial order and the angular degree . For a rotating star, the frequencies also depend on the azimuthal degree, (which takes on values of ). Provided the rotation is slow, the frequencies are well approximated by
Here, is related to the inverse of average of the internal rotation period and we should therefore be able to infer stellar rotation periods directly from the oscillation frequencies.
In principle, the inclination of the rotation axis can also be determined because it affects the relative amplitudes of the different azimuthal degrees. This effect is discussed in detail by [Gizon & Solanki (2003), Gizon & Solanki (2003)]. For example, they show that for intensity observations (oscillations detected in photometry) the dipole multiplets (, , 0, 1) will have relative mode powers given by
In Figures 4, 5 and 6 we show examples of rotational splitting for the Sun. The splitting of the = 1 and 2 modes is clearly seen. The relative strengths of the different components is a result of the observed inclination of the solar rotation axis. Since we observe the Sun in the equatorial plane, we have almost no sensitivity to the modes with =1, =0 and those with =2, =. The structure of the power spectrum therefore indicates that the solar interior is rotating in the surface equatorial plane of the Sun, approximately corresponding to the plane of the ecliptic.
4 The Kepler Asteroseismic Investigation
The NASA Kepler mission ([Borucki \etal (2008), Borucki \etal 2008]) will fly a wide-field Schmidt camera with 0.95m aperture, staring at a single field continuously for at least 3.5 years. Although the mission’s principal aim is to locate transiting extrasolar planets, it will provide an unprecedented opportunity to make asteroseismic observations of a wide variety of stars. This will give the opportunity to measure global properties and internal structure for a large number of stars across a broad range of different types. Plans are now being developed to exploit this opportunity to the fullest. In particular, asteroseismic analysis of Kepler data will provide an accurate determination of the radius for a large fraction of the stars found to host planetary systems, as determined from the transit analysis of the Kepler data, as well as estimates of the ages of the systems. Through investigation of a broad range of stars, the asteroseismic investigation will also substantially improve our understanding of general stellar evolution, and hence strengthen the use of such modelling to further constrain the properties and evolution of the stars and systems investigated in the Kepler extra-solar planet programme. The purpose of the Kepler Asteroseismic Investigation (e.g., [Christensen-Dalsgaard \etal (2007), Christensen-Dalsgaard \etal 2007]) is to ensure that full use is made of this potential to benefit the Kepler investigations of extra-solar planetary systems.
Kepler will give the possibility of observing at two cadences (1 minute and 30 minutes). Due to the high frequency of the stellar oscillations Kepler will not be able to make asteroseismic observations of all the 170,000 stars observed in long cadence mode (30 minute cadence). We will therefore need a target selection programme in order to optimize the asteroseismic part of the Kepler programme. Based on simulations, we have estimated that seismology can be done on unevolved main-sequence stars down to magnitude , and on evolved main-sequence stars (which have higher amplitudes) down to . For subgiant stars we may reach even fainter. It should be noted that Kepler will be the first mission where asteroseismology will be applied to a large number of planet-hosting stars.
Finally, we note that accurate measurements of variations in the light travel times for planets orbiting classical pulsating stars (the pulsations are here used as accurate clocks) may provide a unique tool for the Kepler mission to detect a number of additional planets. The technique requires accurate measurements of phase variations of the pulsations (better than seconds) and the method will be most sensitive to long-period planets (1–2 years). For details on this technique we refer to [Silvotti \etal (2007), Silvotti \etal (2007)].
- Aerts, C., Christensen-Dalsgaard, J., Cunha, M. & Kurtz, D. W., 2008, Solar Phys., in the press [arXiv:0803.3527]
- Arentoft, T., Kjeldsen, H., Bedding, T. R., \etal 2008, ApJ, in the press
- Bedding, T. R. & Kjeldsen, H. 2008, in: G. van Belle (ed.), 14th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ASP Conf. Ser., Vol. 384 (San Francisco: ASP), 21
- Bedding, T. R., Kjeldsen, H., Butler, R. P., et al. 2004, ApJ 614, 380
- Borucki, W. \etal 2008, these Proceedings
- Christensen-Dalsgaard, J. 1993, in: T. M. Brown (ed.), Proc. GONG 1992: Seismic investigation of the Sun and stars, ASP Conf. Ser., Vol. 42 (San Francisco: ASP), 347
- Christensen-Dalsgaard, J. 2004, Solar Phys. 220, 137
- Christensen-Dalsgaard, J., Arentoft, T., Brown, T. M., Gilliland, R. L., Kjeldsen, H., Borucki, W. J. & Koch, D. 2007, Comm. in Asteroseismology 150, 350
- Frandsen, S., Carrier, F., Aerts, C., \etal 2008, A&A 394, L5
- García, R. A., Turck-Chièze, S., Boumier, P., \etal 2005, A&A 442, 385
- Gizon, L. & Solanki, S. K. 2003, ApJ 589, 1009
- Kjeldsen, H., Bedding, T. R., Butler, R. P., \etal 2005, ApJ 635, 1281
- Silvotti, R., Schuh, S., Janulis, R., \etal Nature 449, 189
- Ulrich, R. K., García, R. A., Robillot, J.-M. \etal 2000, A&A 364, 799