Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
EPJ Web of Conferences 101 , 0 1 00 9 (2015) DOI: 10.1051/epjconf/ 201 5 101 0 1 0 0 9 C Owned by the authors, published by EDP Sciences, 2015 Peak Bagging of red giant stars observed by Kepler: first results with a new method based on Bayesian nested sampling Enrico Corsaro1,a and Joris De Ridder1 Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium Abstract. The peak bagging analysis, namely the fitting and identification of single oscillation modes in stars’ power spectra, coupled to the very high-quality light curves of red giant stars observed by Kepler, can play a crucial role for studying stellar oscillations of different flavor with an unprecedented level of detail. A thorough study of stellar oscillations would thus allow for deeper testing of stellar structure models and new insights in stellar evolution theory. However, peak bagging inferences are in general very challenging problems due to the large number of observed oscillation modes, hence of free parameters that can be involved in the fitting models. Efficiency and robustness in performing the analysis is what may be needed to proceed further. For this purpose, we developed a new code implementing the Nested Sampling Monte Carlo (NSMC) algorithm, a powerful statistical method well suited for Bayesian analyses of complex problems. In this talk we show the peak bagging of a sample of high signal-to-noise red giant stars by exploiting recent Kepler datasets and a new criterion for the detection of an oscillation mode based on the computation of the Bayesian evidence. Preliminary results for frequencies and lifetimes for single oscillation modes, together with acoustic glitches, are therefore presented. 1 Introduction Since the first detection of non-radial oscillations in cool giant stars [8], the asteroseismology of field red giants (RGs) has encountered a substantial increase, especially thanks to the advent of the spacebased photometric missions CoRoT (e.g. [20,19]) and Kepler (e.g. [16, 15]), where the latter allowed also for the study of RGs in open clusters (e.g. [25,7]). The discovery of so-called mixed modes [1] has led to a significant improvement in our understanding of the internal structure and evolution of RGs. In particular, the characteristic spacing in period of mixed modes frequencies provides a direct way to disentangle H-shell and He-core burning RGs [2, 22] and has been used already to classify about 13,000 targets observed in the Kepler field of view [24]. Therefore, RGs are among the most interesting and useful types of stars to refine models and test thoroughly codes of stellar structure and evolution. With the present 4-years long datasets available from Kepler for many RGs, we have the possibility to perform very detailed analyses aimed at measuring asteroseismic properties such as frequency, amplitude, and lifetime, of each individual oscillation mode observed in the Fourier spectra of the stars. This particular analysis, usually known as peak bagging plays a key role to exploit the full potential of the informative power spectra of the stars. However, due to the large number of modes populating the power spectrum of a RG (up to ∼100), the peak bagging usually turns into a computationally demanding analysis [11,17,12,5]. The complexity of the analysis is even higher considering that the mode identification of all the oscillation peaks is necessary for a proper modeling of the oscillations. This implies that also testing the significance of individual oscillation peaks becomes an important a e-mail: [email protected] 7KLVLVDQ2SHQ$FFHVVDUWLFOHGLVWULEXWHGXQGHUWKHWHUPVRIWKH&UHDWLYH&RPPRQV$WWULEXWLRQ/LFHQVHZKLFKSHUPLWV XQUHVWULFWHGXVHGLVWULEXWLRQDQGUHSURGXFWLRQLQDQ\PHGLXPSURYLGHGWKHRULJLQDOZRUNLVSURSHUO\FLWHG Article available at http://www.epj-conferences.org or http://dx.doi.org/10.1051/epjconf/201510101009 EPJ Web of Conferences Fig. 1. Resulting background fit of the star KIC 12008916, as derived by DIAMONDS. The original PSD is shown in gray, while a smoothed version with boxcar width set to Δν/5, with Δν taken from [13], is shown as a black line to guide the eye. The red thick line represents the background model without the Gaussian envelope, while the red dotted line accounts for the additional Gaussian component. The individual components of the background model are shown by blue dot-dashed lines. task within the process. For this purpose, [5] developed a new code termed DIAMONDS that implements the nested sampling Monte Carlo algorithm [23], a sampling algorithm well suited for model comparison and efficient high-dimensional inferences in a robust Bayesian perspective [5, 4]. 2 Observations and Data To select the final sample of stars for the analysis we started from the original sample of RGs observed in long cadence (LC) by Kepler and studied by [13] and by [6], consisting in a total of 1111 stars. From the original sample we only considered those stars with power spectra having an oscillation envelope with frequency of maximum power νmax 110 μHz, according to the stellar population study done by [18], with νmax values provided by [13]. We then checked for available measurements of their g-mode period spacing ΔΠ1 , as measured by [22]. This cross-match led to 21 stars observed by Kepler from Q0 till Q17.1, resulting in a frequency resolution of νbin 0.008 μHz for each power spectral density (PSD). The Kepler light curves were corrected according to [10, 9]. In this work we focus on a low-mass low-luminosity RGB star, KIC 12008916, a particular type of RG that is part of a small population of stars (about 5 % of the entire population of RGs observed by Kepler), having high signal-to-noise ratio (SNR) oscillations. Due to the clear oscillation pattern and the large number of modes, this star is well suited for testing stellar evolution theory and pulsation codes, as well as to constrain stellar formation rates in our galaxy. 3 Peak Bagging model and results The first step for the peak bagging analysis of the stars is to estimate the background signal in their PSD. We follow the probabilistic study done on a large sample of stars by [14], who found that two separate granulation components with a characteristic frequency close to the that of the region containing the oscillations, together with a third component related to long-term variations, best represent the background signal of the RGs selected. Thus, the background model that we adopt can be expressed as the sum of three super-Lorentzian profiles, a flat noise component, and a Gaussian envelope to model the power excess of the oscillations. The resulting background fit for KIC 12008916 is shown in Fig. 1. The subsequent step in analyzing the stellar PSDs is to adopt an adequate fitting model for the oscillation pattern contained in the region of the power excess. In this work we adopt a method similar to that used by [5] (see Sect. 6.3). We restrict our frequency range of analysis to the region containing the oscillations, which we identify as νmax ± 3.5σenv , with σenv the standard deviation from the Gaussian envelope used to fit the background. We thus take into account some substantial differences characterizing the oscillation pattern in RGs, as explained in the following. 01009-p.2 The Space Photometry Revolution In the oscillation pattern of all the RGs, one can observe forests of modes between consecutive quadrupole and radial modes, which are known as dipole mixed modes, namely oscillations with a mixed character between that of pure p and g modes. Since g modes have the highest mode inertia because they propagate in high-density regions, their lifetime τ is significantly longer than that of pure pressure modes (e.g. see [3]). As a consequence, even with an observing time T obs > 4 years now made available from Kepler, mixed modes with a g-dominated character appear still unresolved, while others having a more p mode-like character are partially or sometimes even fully resolved. In this work we distinguish the resolved or partially resolved peaks, for which T obs τ, from those that are unresolved, for which T obs τ. In the former case, the oscillation peak profile that we adopt is that of a Lorentzian and it is given by Pres,0 (ν) = A20 / (πΓ0 ) 2 , 0 1 + 4 ν−ν Γ0 (1) where A0 , Γ0 , ν0 are the amplitude, the linewidth and the frequency, respectively, and represent the three free parameters to be estimated during the fitting process. According to the Fourier analysis, an oscillation peak that is not resolved has a profile represented as [3]: Punres,0 (ν) = H0 sinc 2 π (ν − ν0 ) , νbin (2) where H0 and ν0 are the height in PSD units and the central frequency in μHz of the oscillation peak, respectively, and must be estimated during the fitting process, while νbin is the frequency resolution introduced in Sect. 2 and it is a fixed value. In this case fitting the height is preferred since it is an observable and we have no information about the linewidth of the peak, while fitting a Lorentzian profile would represent an inadequate model as it involves a spurious fitting parameter represented by the linewidth. Then, the final peak bagging model can be represented as a mixture of Lorentzian and sinc2 functions. Some of the preliminary results for the peak bagging analysis of KIC 12008916 are shown in Fig. 2, where one can observe a depression of the linewidths of dipole mixed modes in correspondence to the position of the frequency of maximum power νmax , and in Fig. 3, with the evidence of acoustic glitches for the quadrupole and radial modes. 3.1 Peak significance Testing the significance of an oscillation peak in the PSD of a star is a crucial aspect that has to be considered for providing a reliable set of modes that can be used for the modeling of the oscillations, hence for investigating the stellar structure and the evolution of the star. For this purpose, following [5] (Sect. 6.5) we consider a detection probability for a single peak, pB , defined as pB ≡ EB EA + EB (3) where EA and EB are the Bayesian evidences for the models excluding and including, respectively, the oscillation peak to be tested, and are computed with DIAMONDS. The detection limit is usually set to pB 0.99 as derived from a test involving a set of 1000 simulations of peak detection. In other words, we fit a chunk of PSD by considering on one side a model of the peak with the background and on the other side a model of the background alone. The best model selected with the Bayesian model comparison encompassed in Eq. (3) is then telling us whether the peak is detected or not. The method presented in this work proves to be very efficient for this type of RG and will therefore be of great usefulness for the analysis of thousands of stars observed by Kepler. 01009-p.3 EPJ Web of Conferences Fig. 2. Linewidths of resolved dipole mixed modes for KIC 12008916 as a function of the frequency in the PSD. Fig. 3. Signature of the acoustic glitches for KIC 12008916 in the large frequency separation Δν for the quadrupole and radial modes. References 1. Beck, P. G., Bedding, T. R., Mosser, B., et al., Science 332, (2011) 205 2. Bedding, T. R., Mosser, B., Huber, D., et al., Nature 471, (2011) 608 3. Christensen-Dalsgaard, J., Sol. Phys. 220, (2004) 137 4. Corsaro, E. & De Ridder, J., Garcı́a, R. A. A&A arXiv:1503.08821 [astro-ph.IM] (2015) 5. Corsaro, E. & De Ridder, J., A&A 571, (2014) A71 6. Corsaro, E., Fröhlich, H.-E., Bonanno, A., et al., MNRAS 430, (2013) 2313 7. Corsaro, E., Stello, D., Huber, D., et al., ApJ757, (2012) 190 8. De Ridder, J., Barban, C., Baudin, F., et al., Nature 459, (2009) 398 9. Garcı́a, R. A., Mathur, S., Pires, S., et al., A&A 568, (2014) A10 10. Garcı́a, R. A., Hekker, S., Stello, D., et al., MNRAS 414, (2011) L6 11. Gruberbauer, M., Kallinger, T., Weiss, W. W., & Guenther, D. B., A&A 506, (2009) 1043 12. Handberg, R. & Campante, T. L., A&A 527, (2011) A56 13. Huber, D., Bedding, T. R., Stello, D., et al., ApJ 743, (2011) 143 14. Kallinger, T., De Ridder, J., Hekker, S., et al., A&A 570, (2014) A41 15. Kallinger, T., Hekker, S., Mosser, B., et al., A&A 541, (2012) A51 16. Kallinger, T., Mosser, B., Hekker, S., et al., A&A 522, (2010a) A1 17. Kallinger, T., Weiss, W. W., Barban, C., et al., A&A 509, (2010b) A77 18. Miglio, A., Montalbán, J., Baudin, F., et al., A&A 503, (2009) L21 19. Mosser, B., Barban, C., Montalbán, J., et al., A&A 532, (2011a) A86 20. Mosser, B., Belkacem, K., Goupil, M. J., et al., A&A 525, (2011b) L9 21. Mosser, B., Goupil, M. J., Belkacem, K., et al., A&A 548, (2012a) A10 22. Mosser, B., Goupil, M. J., Belkacem, K., et al., A&A 540, (2012b) A143 23. Skilling, J., AIP Conference Proceedings 735, (2004) 395 24. Stello, D., Huber, D., Bedding, T. R., et al., ApJL 765, (2013) L41 25. Stello, D., Meibom, S., Gilliland, R. L., et al., ApJ 739, (2011) 13 01009-p.4