Abstract
The advent of space-based missions like Kepler has revolutionized the study of solar-type stars, particularly through the measurement and modeling of their resonant modes of oscillation. Here we analyze a sample of 66 Kepler main-sequence stars showing solar-like oscillations as part of the Kepler seismic LEGACY project. We use Kepler short-cadence data, of which each star has at least 12 months, to create frequency-power spectra optimized for asteroseismology. For each star, we identify its modes of oscillation and extract parameters such as frequency, amplitude, and line width using a Bayesian Markov chain Monte Carlo "peak-bagging" approach. We report the extracted mode parameters for all 66 stars, as well as derived quantities such as frequency difference ratios, the large and small separations and the behavior of line widths with frequency and line widths at with , for which we derive parametrizations; and behavior of mode visibilities. These average properties can be applied in future peak-bagging exercises to better constrain the parameters of the stellar oscillation spectra. The frequencies and frequency ratios can tightly constrain the fundamental parameters of these solar-type stars, and mode line widths and amplitudes can test models of mode damping and excitation.
Export citation and abstract BibTeX RIS
1. Introduction
The study of stars and extrasolar planets via the properties of their host stars has experienced a revolution in recent years (Chaplin & Miglio 2013; Christensen-Dalsgaard 2016). This largely arose from the successful application of asteroseismology using observations from the CoRoT (Baglin et al. 2009) and Kepler missions (Gilliland et al. 2010). This application has been made possible by extracting high-precision parameters from the stellar frequency-power spectra owing to the long time-baseline and photometric quality of these space missions.
Asteroseismology allows the determination of fundamental stellar parameters such as mass, radius, and age through modeling of individual mode frequencies or frequency-difference ratios. The Kepler mission has already provided stellar parameters for a number of stars, including planetary hosts, using average seismic parameters (Chaplin et al. 2011c; Silva Aguirre et al. 2012; Huber et al. 2013; Chaplin et al. 2014), individual frequencies (Basu et al. 2010; Christensen-Dalsgaard et al. 2010; Howell et al. 2012; Metcalfe et al. 2012, 2014, 2015; Lund et al. 2014b; Van Eylen et al. 2014; Campante et al. 2015; Silva Aguirre et al. 2015), and frequency-difference ratios (Silva Aguirre et al. 2013, 2015; Lebreton & Goupil 2014).
The high precision of extracted mode frequencies further allows the study of ionization zones and the convective envelope boundary from acoustic glitches (Houdek & Gough 2007; Silva Aguirre et al. 2011; Mazumdar et al. 2014; Verma et al. 2014), and one may also learn about the physics of the excitation and damping of the oscillation modes from measured mode line widths, amplitudes, and visibilities (Houdek et al. 1999; Samadi et al. 2005, 2007; Houdek 2006; Belkacem et al. 2012).
Frequencies and line widths have been reported for several solar-like and subgiant stars observed by Kepler by Appourchaux et al. (2012a, 2012b, 2014), and planet-hosting stars by Davies et al. (2016). In this paper, we analyze a sample of 66 main-sequence (MS) solar-like stars observed for at least 12 months by the Kepler mission. We extracted mode parameters by "peak-bagging"13 the frequency-power spectra of the stars using a Markov chain Monte Carlo (MCMC) method (Handberg & Campante 2011) and used the Bayesian quality control presented by Davies et al. (2016). For each star, we report values for the mode frequencies, amplitudes, line widths, and visibilities. Additionally, we provide summary descriptions for each of the above quantities, such as average seismic parameters derived from the frequencies and prescriptions of the mode line widths against frequency. The frequencies reported here will be modeled in the accompanying paper by Silva Aguirre et al. 2017, hereafter Paper II. The lessons learned from the presented analysis will be useful for the study of MS solar-like oscillators with the TESS (Ricker et al. 2014) and PLATO (Rauer et al. 2013) missions, and the continued analysis of these data from K2 (Chaplin et al. 2015; Lund et al. 2016a, 2016b).
The paper is structured as follows. Section 2 describes the target sample, including the preparation of Kepler data and spectroscopic properties. Section 3 is devoted to the parameter estimation from the MCMC peak-bagging, including a description of the fitting strategy, the adopted Bayesian quality assurance, and the derivation of frequency difference ratios and their correlations. In Section 4, we present our results from the peak-bagging for the mode frequencies, focusing specifically on frequency errors and average seismic parameters in Section 4.1; amplitudes in Section 4.2; line widths in Section 4.3; and visibilities in Section 4.4. In Section 5, we give an example of the output generated for each of the analyzed stars. We conclude in Section 6.
2. Target Sample
Our sample consists of 66 solar-type oscillators observed by the Kepler satellite, all part of the KASC (Kjeldsen et al. 2010) working group 1 (WG1) sample of solar-like p-mode oscillators. All stars have short-cadence (SC; ) observations with an observing baseline of at least 12 months, and represent some of the highest signal-to-noise solar-like oscillators observed by Kepler. The sample consists only of main-sequence (MS) and slightly more evolved subgiant stars. These have frequency structures corresponding to the "Simple" or "F-type" categories by Appourchaux et al. (2012b), i.e., none of the stars show obvious bumped dipole modes. The sample was peak-bagged as part of the Kepler dwarf seismic "LEGACY" project, with the asteroseismic modeling of extracted parameters presented in Paper II. In Figure 1, the sample is shown in a Kiel-diagram ( versus ), with parameters adopted from Paper II; for additional details on the sample see Table 1. We note that all targets from the Appourchaux et al. (2014) study of oscillation mode line widths (which included data up to Quarter 12) are part of our sample, with the exception of KIC 3424541, 3733735, 10355856, and 10909629. These four stars were classified as F-type by Appourchaux et al. (2014), but were omitted from our sample because of possible mixed-mode structures.
Table 1. Parameters for the Targets in the Studied Sample
KIC | Popular | Kp | Number of | Category | Braketing | Missing | [Fe/H] | LOS | |||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Name | (mag) | (μHz) | (μHz) | Modes | Quarters | Quarters | (K) | (dex) | (cgs; dex) | () | () | ||
1435467 | ⋯ | 8.88 | 1407 | 70.4 | 46 | F-like | 5.1–17.2 | ⋯ | 6326 ± 77 | 0.01 ± 0.10 | 4.100 | −66.52 ± 0.10 | 11.90 ± 0.50 |
2837475 | ⋯ | 8.48 | 1558 | 75.7 | 54 | F-like | 5.1–17.2 | ⋯ | 6614 ± 77 | 0.01 ± 0.10 | 4.163 | −13.28 ± 0.10 | 23.30 ± 0.50 |
3427720 | ⋯ | 9.11 | 2737 | 120.1 | 36 | Simple | 5.1–17.2 | ⋯ | 6045 ± 77 | −0.06 ± 0.10 | 4.387 | −22.91 ± 0.10 | 2.90 ± 0.50 |
3456181 | ⋯ | 9.66 | 970 | 52.3 | 45 | F-like | 5.1–11.3 | 6, 10 | 6384 ± 77 | −0.15 ± 0.10 | 3.950 | −50.12 ± 0.10 | 8.50 ± 0.50 |
3632418 | Cassie | 8.22 | 1167 | 60.7 | 54 | F-like | 5.1–17.2 | ⋯ | 6193 ± 77 | −0.12 ± 0.10 | 4.024 | −19.11 ± 0.10 | 8.50 ± 0.50 |
3656476 | Java | 9.52 | 1925 | 93.2 | 38 | Simple | 5.1–17.2 | 6, 10, 14 | 5668 ± 77 | 0.25 ± 0.10 | 4.225 | −13.29 ± 0.10 | 2.30 ± 0.50 |
3735871 | ⋯ | 9.71 | 2863 | 123.0 | 34 | Simple | 5.1–17.2 | ⋯ | 6107 ± 77 | −0.04 ± 0.10 | 4.396 | 6.90 ± 0.10 | 4.80 ± 0.50 |
4914923 | Vitto | 9.46 | 1817 | 88.5 | 38 | Simple | 5.1–17.2 | 6 | 5805 ± 77 | 0.08 ± 0.10 | 4.197 | −39.16 ± 0.10 | 3.40 ± 0.50 |
5184732 | Kitty | 8.16 | 2089 | 95.5 | 49 | Simple | 7.1–17.2 | ⋯ | 5846 ± 77 | 0.36 ± 0.10 | 4.255 | 15.41 ± 0.10 | 4.00 ± 0.50 |
5773345 | ⋯ | 9.16 | 1101 | 57.3 | 45 | F-like | 6.1–11.3 | 9 | 6130 ± 84(5) | 0.21 ± 0.09(5) | 3.993 | ⋯ | ⋯ |
5950854 | ⋯ | 10.96 | 1927 | 96.6 | 26 | Simple | 5.1–10.3 | 6, 7.2 | 5853 ± 77 | −0.23 ± 0.10 | 4.238 | −42.49 ± 0.10 | 3.00 ± 0.50 |
6106415 | Perky | 7.18 | 2249 | 104.1 | 49 | Simple | 6.1–16.3 | 9, 13 | 6037 ± 77 | −0.04 ± 0.10 | 4.295 | −14.80 ± 0.10 | 4.90 ± 0.50 |
6116048 | Nunny | 8.42 | 2127 | 100.8 | 49 | Simple | 5.1–17.2 | ⋯ | 6033 ± 77 | −0.23 ± 0.10 | 4.254 | −53.26 ± 0.10 | 4.00 ± 0.50 |
6225718 | Saxo2 | 7.50 | 2364 | 105.7 | 59 | Simple | 6.1–17.2 | ⋯ | 6313 ± 77 | −0.07 ± 0.10 | 4.319 | −1.32 ± 0.10 | 5.50 ± 0.50 |
6508366 | Baloo | 8.97 | 958 | 51.6 | 50 | F-like | 5.1–17.2 | ⋯ | 6331 ± 77 | −0.05 ± 0.10 | 3.942 | 2.62 ± 0.10 | 22.50 ± 0.50 |
6603624 | Saxo | 9.09 | 2384 | 110.1 | 44 | Simple | 5.1–17.2 | ⋯ | 5674 ± 77 | 0.28 ± 0.10 | 4.320 | −58.82 ± 0.10 | 0.70 ± 0.50 |
6679371 | ⋯ | 8.73 | 942 | 50.6 | 55 | F-like | 5.1–17.2 | ⋯ | 6479 ± 77 | 0.01 ± 0.10 | 3.934 | −23.58 ± 0.10 | 17.30 ± 0.50 |
6933899 | Fred | 9.62 | 1390 | 72.1 | 39 | Simple | 5.1–17.2 | ⋯ | 5832 ± 77 | −0.01 ± 0.10 | 4.079 | −6.97 ± 0.10 | 3.60 ± 0.50 |
7103006 | ⋯ | 8.86 | 1168 | 59.7 | 54 | F-like | 5.1–17.2 | ⋯ | 6344 ± 77 | 0.02 ± 0.10 | 4.015 | −22.36 ± 0.10 | 12.10 ± 0.50 |
7106245 | ⋯ | 10.79 | 2398 | 111.4 | 24 | Simple | 5.1–15.3 | ⋯ | 6068 ± 102(3) | −0.99 ± 0.19(3) | 4.310 | ⋯ | ⋯ |
7206837 | Bagheera | 9.77 | 1653 | 79.1 | 45 | F-like | 5.1–17.2 | ⋯ | 6305 ± 77 | 0.10 ± 0.10 | 4.163 | −18.54 ± 0.10 | 9.30 ± 0.50 |
7296438 | ⋯ | 10.09 | 1848 | 88.7 | 32 | Simple | 7.1–11.3 | ⋯ | 5775 ± 77 | 0.19 ± 0.10 | 4.201 | 3.36 ± 0.10 | 1.80 ± 0.50 |
7510397 | ⋯ | 7.77 | 1189 | 62.2 | 47 | Simple | 7.1–17.2 | 16 | 6171 ± 77 | −0.21 ± 0.10 | 4.036 | −34.10 ± 0.10 | 6.40 ± 0.50 |
7680114 | Simba | 10.07 | 1709 | 85.1 | 41 | Simple | 5.1–17.2 | 6, 7.2, 10 | 5811 ± 77 | 0.05 ± 0.10 | 4.172 | −58.93 ± 0.10 | 3.00 ± 0.50 |
7771282 | ⋯ | 10.77 | 1465 | 72.5 | 32 | F-like | 5.1–11.3 | 6 | 6248 ± 77 | −0.02 ± 0.10 | 4.112 | −0.38 ± 0.10 | 8.30 ± 0.50 |
7871531 | ⋯ | 9.25 | 3456 | 151.3 | 35 | Simple | 5.1–17.2 | ⋯ | 5501 ± 77 | −0.26 ± 0.10 | 4.478 | −20.65 ± 0.10 | 0.90 ± 0.50 |
7940546 | Akela | 7.40 | 1117 | 58.8 | 58 | F-like | 7.1–17.2 | ⋯ | 6235 ± 77 | −0.20 ± 0.10 | 4.000 | −3.03 ± 0.10 | 9.10 ± 0.50 |
7970740 | ⋯ | 7.78 | 4197 | 173.5 | 46 | Simple | 6.1–17.2 | ⋯ | 5309 ± 77 | −0.54 ± 0.10 | 4.539 | −60.24 ± 0.10 | 0.00 ± 0.50 |
8006161 | Doris | 7.36 | 3575 | 149.4 | 54 | Simple | 5.1–17.2 | ⋯ | 5488 ± 77 | 0.34 ± 0.10 | 4.494 | −45.56 ± 0.10 | 0.70 ± 0.50 |
8150065 | ⋯ | 10.74 | 1877 | 89.3 | 24 | Simple | 5.1–10.3 | 6, 7.2 | 6173 ± 101(3) | −0.13 ± 0.15(3) | 4.220 | ⋯ | ⋯ |
8179536 | ⋯ | 9.46 | 2075 | 95.1 | 39 | Simple | 5.1–11.3 | 6 | 6343 ± 77 | −0.03 ± 0.10 | 4.255 | −31.40 ± 0.10 | 9.90 ± 0.50 |
8228742 | Horace | 9.37 | 1190 | 62.1 | 44 | Simple | 5.1–17.2 | ⋯ | 6122 ± 77 | −0.08 ± 0.10 | 4.032 | 10.71 ± 0.10 | 6.10 ± 0.50 |
8379927 | Arthur | 6.96 | 2795 | 120.3 | 49 | Simple | 2.1–17.2 | 2.2, 2.3, 3, 4 | 6067 ± 120(1) | −0.10 ± 0.15(1) | ⋯ | ⋯ | |
8394589 | ⋯ | 9.52 | 2397 | 109.5 | 44 | Simple | 5.1–17.2 | ⋯ | 6143 ± 77 | −0.29 ± 0.10 | 4.322 | 22.58 ± 0.10 | 6.40 ± 0.50 |
8424992 | ⋯ | 10.32 | 2534 | 120.6 | 22 | Simple | 7.1–10.3 | ⋯ | 5719 ± 77 | −0.12 ± 0.10 | 4.359 | −87.63 ± 0.10 | 1.30 ± 0.50 |
8694723 | ⋯ | 8.88 | 1471 | 75.1 | 53 | Simple | 5.1–17.2 | ⋯ | 6246 ± 77 | −0.42 ± 0.10 | 4.113 | 15.88 ± 0.10 | 7.10 ± 0.50 |
8760414 | Pucky | 9.62 | 2455 | 117.2 | 44 | Simple | 5.1–17.2 | ⋯ | 5873 ± 77 | −0.92 ± 0.10 | 4.320 | −115.64 ± 0.10 | 2.50 ± 0.50 |
8938364 | Java2 | 10.11 | 1675 | 85.7 | 41 | Simple | 6.1–17.2 | ⋯ | 5677 ± 77 | −0.13 ± 0.10 | 4.173 | −68.12 ± 0.10 | 2.40 ± 0.50 |
9025370 | ⋯ | 8.85 | 2989 | 132.6 | 28 | Simple | 5.1–17.2 | ⋯ | 5270 ± 180(2) | −0.12 ± 0.18(2) | 4.423 | ⋯ | ⋯ |
9098294 | ⋯ | 9.76 | 2315 | 108.9 | 34 | Simple | 5.1–17.2 | ⋯ | 5852 ± 77 | −0.18 ± 0.10 | 4.308 | −71.72 ± 0.10 | 3.00 ± 0.50 |
9139151 | Carlsberg | 9.18 | 2690 | 117.3 | 35 | Simple | 5.1–17.2 | ⋯ | 6302 ± 77 | 0.10 ± 0.10 | 4.382 | −29.06 ± 0.10 | 5.50 ± 0.50 |
9139163 | Punto | 8.33 | 1730 | 81.2 | 57 | F-like | 5.1–17.2 | ⋯ | 6400 ± 84(5) | 0.15 ± 0.09(5) | 4.200 | ⋯ | ⋯ |
9206432 | ⋯ | 9.08 | 1866 | 84.9 | 49 | F-like | 5.1–12.3 | 7 | 6538 ± 77 | 0.16 ± 0.10 | 4.220 | −1.73 ± 0.10 | 6.70 ± 0.50 |
9353712 | ⋯ | 10.84 | 934 | 51.5 | 41 | F-like | 5.1–12.3 | 6, 7.2 | 6278 ± 77 | −0.05 ± 0.10 | 3.943 | −46.67 ± 0.10 | 6.80 ± 0.50 |
9410862 | ⋯ | 10.71 | 2279 | 107.4 | 33 | Simple | 5.1–15.3 | ⋯ | 6047 ± 77 | −0.31 ± 0.10 | 4.300 | −56.84 ± 0.10 | 3.80 ± 0.50 |
9414417 | ⋯ | 9.58 | 1155 | 60.1 | 54 | F-like | 6.1–17.2 | 7 | 6253 ± 75(6) | −0.13 ± 0.10(6) | 4.016 | ⋯ | ⋯ |
9812850 | ⋯ | 9.47 | 1255 | 64.7 | 49 | F-like | 5.1–17.2 | ⋯ | 6321 ± 77 | −0.07 ± 0.10 | 4.053 | 31.18 ± 0.10 | 12.50 ± 0.50 |
9955598 | ⋯ | 9.44 | 3617 | 153.3 | 31 | Simple | 5.1–17.2 | ⋯ | 5457 ± 77 | 0.05 ± 0.10 | 4.497 | −28.48 ± 0.10 | 1.10 ± 0.50 |
9965715 | ⋯ | 9.34 | 2079 | 97.2 | 40 | Simple | 5.1–13.3 | 7, 11 | 5860 ± 180(2) | −0.44 ± 0.18(2) | 4.272 | ⋯ | ⋯ |
10068307 | ⋯ | 8.18 | 995 | 53.9 | 49 | Simple | 7.1–17.2 | ⋯ | 6132 ± 77 | −0.23 ± 0.10 | 3.967 | −14.78 ± 0.10 | 6.40 ± 0.50 |
10079226 | ⋯ | 10.07 | 2653 | 116.3 | 31 | Simple | 7.1–10.3 | ⋯ | 5949 ± 77 | 0.11 ± 0.10 | 4.366 | −37.15 ± 0.10 | 4.00 ± 0.50 |
10162436 | ⋯ | 8.61 | 1052 | 55.7 | 51 | F-like | 5.1–17.2 | 7, 10, 11, 15 | 6146 ± 77 | −0.16 ± 0.10 | 3.981 | −52.92 ± 0.10 | 6.40 ± 0.50 |
10454113 | Pinocha | 8.62 | 2357 | 105.1 | 54 | Simple | 5.1–17.2 | ⋯ | 6177 ± 77 | −0.07 ± 0.10 | 4.314 | −21.22 ± 0.10 | 6.10 ± 0.50 |
10516096 | Manon | 9.46 | 1690 | 84.4 | 40 | Simple | 5.1–17.2 | 6, 10.1 | 5964 ± 77 | −0.11 ± 0.10 | 4.169 | 1.28 ± 0.10 | 4.60 ± 0.50 |
10644253 | Mowgli | 9.16 | 2900 | 123.1 | 34 | Simple | 5.1–17.2 | ⋯ | 6045 ± 77 | 0.06 ± 0.10 | 4.396 | −18.91 ± 0.10 | 3.20 ± 0.50 |
10730618 | ⋯ | 10.45 | 1282 | 66.3 | 39 | F-like | 0–11.3 | 6, 7.2 | 6150 ± 180(2) | −0.11 ± 0.18(2) | 4.062 | ⋯ | ⋯ |
10963065 | Rudy | 8.77 | 2204 | 103.2 | 42 | Simple | 2.3–17.2 | 3, 4, 8, 9, 12, 16 | 6140 ± 77 | −0.19 ± 0.10 | 4.277 | −54.95 ± 0.10 | 4.50 ± 0.50 |
11081729 | ⋯ | 9.03 | 1968 | 90.1 | 40 | F-like | 5.1–17.2 | ⋯ | 6548 ± 83 | 0.11 ± 0.10 | 4.245 | 0.27 ± 0.10 | 24.10 ± 0.50 |
11253226 | Tinky | 8.44 | 1591 | 76.9 | 58 | F-like | 5.1–17.2 | ⋯ | 6642 ± 77 | −0.08 ± 0.10 | 4.173 | 10.65 ± 0.10 | 14.40 ± 0.50 |
11772920 | ⋯ | 9.66 | 3675 | 157.7 | 27 | Simple | 5.1–17.2 | ⋯ | 5180 ± 180(2) | −0.09 ± 0.18(2) | 4.500 | ⋯ | ⋯ |
12009504 | Dushera | 9.32 | 1866 | 88.2 | 43 | Simple | 5.1–17.2 | ⋯ | 6179 ± 77 | −0.08 ± 0.10 | 4.211 | 12.82 ± 0.10 | 7.70 ± 0.50 |
12069127 | ⋯ | 10.70 | 885 | 48.4 | 39 | F-like | 5.1–11.3 | 6 | 6276 ± 77 | 0.08 ± 0.10 | 3.912 | −25.33 ± 0.10 | 6.00 ± 0.50 |
12069424 | 16 Cyg A | 5.86 | 2188 | 103.3 | 53 | Simple | 6.1–17.2 | ⋯ | 5825 ± 50(4) | 0.10 ± 0.03(4) | 4.287 | −27.35 ± 0.10 | 2.80 ± 0.50 |
12069449 | 16 Cyg B | 6.09 | 2561 | 116.9 | 52 | Simple | 6.1–17.2 | ⋯ | 5750 ± 50(4) | 0.05 ± 0.02(4) | 4.353 | −27.82 ± 0.10 | 2.10 ± 0.50 |
12258514 | Barney | 8.08 | 1513 | 74.8 | 45 | Simple | 5.1–17.2 | ⋯ | 5964 ± 77 | −0.00 ± 0.10 | 4.126 | −18.98 ± 0.10 | 3.90 ± 0.50 |
12317678 | ⋯ | 8.74 | 1212 | 63.5 | 57 | F-like | 5.1–17.2 | ⋯ | 6580 ± 77 | −0.28 ± 0.10 | 4.048 | −58.14 ± 0.10 | 8.40 ± 0.50 |
Note: and [Fe/H] from (1) Pinsonneault et al. (2012); (2) Pinsonneault et al. (2014); (3) the SAGA project (Casagrande et al. 2014, see http://www.mso.anu.edu.au/saga/saga_home.html); (4) Ramírez et al. (2009); (5) Chaplin et al. (2014); or (6) Huber et al. (2013). Metallicities, temperatures, line of sight (LOS) velocities, and values are adopted from the SPC analysis of the targets unless otherwise indicated (see table notes). Systematic uncertainties of ±59 K () and ±0.062 dex ([Fe/H]) have been added in quadrature as suggested by Torres et al. (2012). Values for are adopted from Paper II. The table lists the KIC value and other popular names given, if any; Kepler magnitude (Kp); frequency of maximum amplitude and large separation ()—for uncertainties in and see Table 2; number of peak-bagged modes; the category according to Appourchaux et al. (2012b); first-last quartes during which the targets were observed in SC, and which quarters were missing in between.
2.1. Data Preparation
For most targets, data were taken continuously from Quarter 5 (Q5) through Q17. To minimize gaps in the time series, data from the initial short quarters (Q0 or Q1) were omitted unless continuous with the subsequent data. Table 1 lists the quarters used for each target. Light curves were constructed from pixel data downloaded from the KASOC database,14 using the procedure developed by S. Bloemen (2016, private communication) to define pixel masks for aperture photometry. The light curves were then corrected using the KASOC filter (see Handberg & Lund 2014). Briefly, the light curves were first corrected for jumps and concatenated. They were then median filtered using two filters of different widths—one long, one short—with the final filter being a weighted sum of the two filters based on the variability in the light curve. For the four Kepler objects of interest (KOIs) in the sample (KICs 3632418, 9414417, 9955598, and 10963065) an iterative removal of the planetary transits was performed based on the planetary phase-curve (see Handberg & Lund 2014 for further details).
The power density spectrum (PDS) returned from the KASOC filter is made from a weighted least-squares sine-wave fitting, single-sided calibrated, normalized to Parseval's theorem, and converted to power density by dividing by the integral of the spectral window (Kjeldsen 1992; Kjeldsen & Frandsen 1992).
2.2. Atmospheric and Stellar Parameters
We have obtained atmospheric parameters from the Stellar Parameters Classification tool (SPC; see Buchhave et al. 2012), with data from the Tillinghast Reflector Echelle Spectrograph (TRES; Szentgyorgyi & Furész 2007; Fürész 2008) on the 1.5-m Tillinghast telescope at the F. L. Whipple Observatory. Information from the SPC analysis is available on the Kepler Community Follow-up Observing Program (CFOP) website.15 In the SPC derivation of parameters, values were fixed to the asteroseismic values given in Chaplin et al. (2014) to decrease the impact on uncertainties from correlations between , log g, and . We added in quadrature to the derived uncertainties on and systematic uncertainties of ±59 K and ±0.062 dex, as suggested by Torres et al. (2012). For a subset of targets, spectroscopic values were taken from the literature (Table 1). We also list in Table 1 the line-of-sight (LOS) velocities derived from the SPC analysis, which should be used in any modeling efforts using individual frequencies to account for the Doppler shift of the frequencies (Davies et al. 2014b). In Figure 2, we show the values of these Doppler frequency shifts, which in some cases exceed the uncertainties on the individual frequencies. Even if the frequency shift is small compared to the uncertainties on the mode frequencies, it is systematic and should therefore always be corrected to avoid biases in the stellar modeling. The SPC LOS values have been corrected by to put the velocities onto the IAU system. This correction is primarily accounting for the fact that the CfA library of synthetic spectra does not include the solar gravitational redshift. Stellar parameters used in this paper, such as masses and radii, are adopted from the modeling effort presented in Paper II.
Download figure:
Standard image High-resolution image2.3. Sun-as-a-star Data
As part of the project, the Sun was fitted in the same manner as the sample targets (see Section 3). This was done primarily to test the modeling efforts presented in Paper II against a known reference, and at the same time to assess the returns from the peak-bagging. The power spectrum was produced from data from VIRGO16 (Fröhlich 2009) on board the SoHO17 spacecraft (Fröhlich et al. 1995; Frohlich et al. 1997). Specifically, a time series was created from a weighted sum of the green and red channels of the VIRGO Sun photometers (SPM) with central wavelengths of 500 nm (green), and 862 nm (red). Weights were selected such that the response-function weighted centroid wavelength from the two SPM channels matched that from the Kepler response function (641.7 nm). The two-component light curves were filtered individually using a 30 day median filter and then summed in relative flux units with the appropriate weights (green: 0.785; red: 0.215). The solar time series had a length of 1150 days (corresponding to ∼3.15 years, or the approximate duration of 13 Quarters). This is the typical time series length for targets in the sample. To find the level to which the spectrum should be degraded, the magnitude distribution was computed for the sample, including also stars that have a mixed-mode character. The median magnitude of closely matches that of KIC 9139151 and so noise was added to the solar time series to match the level of this star.
The solar data set will primarily be used for estimates relating to frequencies, such as and but not for analysis of line widths, amplitudes, or visibilities. This is because one cannot, with the simple weighting of relatively narrow band filters done here, assume that the measurements of amplitudes and visibilities adhere strictly to what would be observed with Kepler.
3. Parameter Estimation
3.1. Oscillation Spectrum Model
To model the power spectrum, we described each oscillation mode as a Lorentzian function (Lnlm), which corresponds to the shape of a stochastically excited and intrinsically damped mode (Batchelor 1953; Kumar et al. 1988):
Each mode is characterized by the frequency of the zonal (m = 0) component, a height , an FWHM mode width , and a rotational splitting (assumed constant with frequency; see Lund et al. 2014c, for a discussion of the impact of differential rotation on the constancy of ). In , is the geometrical factor that sets the relative visibilities between the 2l + 1 (azimuthal) m-components as a function of the stellar inclination i⋆ (see, Dziembowski 1977; Gizon & Solanki 2003); denotes the squared visibility (power units) of a non-radial mode relative to a radial mode at the same frequency, i.e.,
from the spatial filtering resulting from integrating the intensity for a mode of a given degree over the stellar surface; then denotes the height of the radial mode of the order of n. The use of assumes equipartition of energy between modes of different angular degrees, thus only with a dependence on frequency. This is a good assumption for stochastically excited low degree high order acoustic modes, observed for many lifetimes (see, e.g., Christensen-Dalsgaard & Gough 1982).
The full model to fit to the power spectrum is then given by a series of the Lorentzian functions in Equation (1) as
Here denotes the adopted background function; W is a constant white shot-noise component; describes the apodization of the signal power at frequency ν from the ∼1-minute sampling of the temporal signal (see, e.g., Chaplin et al. 2011b; Kallinger et al. 2014) and is given by
where Δt gives the integration time for the observations.18 For the background, we used the function (Harvey et al. 1993; Andersen et al. 1994)
which characterizes a temporal signal from granulation having an exponentially decaying autocovariance, with a power of the temporal decay rate as gives the characteristic timescale of the ith background component; the corresponding root-mean-square (rms) variation of the component in the time domain. The normalization constants are such that the integral (for positive frequencies) of the background component equals , in accordance with the Parseval-Plancherel theorem (see, e.g., Michel et al. 2009; Karoff et al. 2013; Kallinger et al. 2014).
In fitting Equation (3) to the power spectrum, we varied the mode amplitude (square-root of integrated mode power) rather than the mode height to decrease the correlation with (Toutain & Appourchaux 1994). To obtain the height (S) in power density units from the varied amplitude (A), we used the relation (Fletcher et al. 2006; Chaplin et al. 2008b):
This is a valid approximation for a single-sided power spectrum when the modes are well resolved, i.e., when the observing duration greatly exceeds the mode lifetime . We note that Anl and were varied for radial modes only (l = 0), and then linearly interpolated to the frequencies of the non-radial modes. The fitting of the power spectrum then finally involved estimating the parameters .
3.2. Fitting Strategy
Parameters were estimated in a Bayesian manner from a global peak-bagging fit to the power spectrum including all parameters (see, e.g., Handberg & Campante 2011). This was done by mapping the posterior probability of the parameters given the data D and any prior information I, which from Bayes' theorem is given as
Here is the prior probability assigned to the parameters given I, and is the likelihood of the observed data D given the parameters . , known as the evidence, is given by the integral of the numerator over the full parameter space, and thus acts as a normalization. The evidence is unnecessary in the mapping of the relative posterior distribution, so we end up mapping:
where logarithmic units are adopted for numerical stability, and C is a constant. Assuming a 2-dof statistic for the power spectrum relative to the limit spectrum in Equation (3) (Gabriel 1994), the logarithm of the likelihood for a given observed power, Oj, relates to the limit spectrum, , as (see Duvall & Harvey 1986; Anderson et al. 1990; Toutain & Appourchaux 1994):
Mapping of Equation (7) was performed using an affine invariant ensemble Markov chain Monte Carlo (MCMC) sampler (see Goodman & Weare 2010), specifically via the Python implementation emcee by Foreman-Mackey et al. (2013). For a given fit, we employed 500 so-called walkers that were initiated by sampling from the priors of the model parameters (see Section 3.2.2). Each of these was run for at least 2000 steps. We further adopted parallel tempering using five temperatures, with tempering parameters determined according to Benomar et al. (2009), and a thinning of the MCMC chains by a factor of 10. As part of the post-processing, the appropriate burn-in for a given target and whether sufficient mixing had been achieved was determined by (1) visual inspection of the chain traces, (2) using the Geweke diagnostic (Geweke 1992), and (3) by assessing the length of the chain compared to the autocorrelation time (giving the number of independent draws from the target distribution).
Final parameter estimates were obtained from the median (frequencies) or mode (amplitudes, line widths, and visibilities) of the marginalized posterior probability density functions (PDFs)—with the MCMC sampling, the marginalization is obtained naturally and the PDF for a given parameter is simply given by the normalized distribution of the samples of the parameter. A measure for the parameter uncertainty is given by the credible interval as the interval spanning the 68.27% highest probability density (HPD) of the PDF.
3.2.1. Mode Identification and Initial Guesses
Before the peak-bagging can commence, initial guesses must be defined for the mode-frequencies to include in Equation (3), and the modes must further be identified in terms of their angular degree l. For acoustic modes of high radial order n and low angular degree l the frequencies may be approximated by the asymptotic relation (Tassoul 1980; Scherrer et al. 1983):
Here is the large separation, given by the average frequency spacing between consecutive overtones n for modes of a given l; is a dimensionless offset sensitive to the surface layers (see, e.g., Gough 1986; Pérez Hernández & Christensen-Dalsgaard 1998; Roxburgh 2016); D0 is sensitive to the sound-speed gradient near the stellar core (Scherrer et al. 1983; Christensen-Dalsgaard 1993). Mode identification was then, by and large, achieved via visual inspection of échelle diagrams (Grec et al. 1983; Bedding 2011). Here, modes of a given l will form vertical ridges for the correct average large separation. The identification of l and radial order n was checked against the relation for as a function of (Figure 3), where can be found from the échelle diagram (Figure 4) by the vertical position of the radial degree (l = 0) ridge (see White et al. 2011a, 2011b).
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageFor this study, consisting of high S/N oscillation signals, the identification was relatively simple. Initial guesses for mode frequencies were primarily defined by hand from smoothed versions of the power density spectra. These were checked against frequencies returned from applying the pseudo-global fitting method of Fletcher et al. (2009).
The power spectrum was fitted in the range of to , where and denote the minimum and maximum mode frequency included in the peak-bagging. Before the peak-bagging fit, a background-only fit was performed in the range from 5 μHz to the SC Nyquist frequency (∼8496 μHz). In this fit, the power from solar-like oscillations was accounted for by a Gaussian envelope centered at Using the posterior distributions from the background-only fit as priors in the peak-bagging allowed us to constrain the background in the relatively narrow frequency range included.
3.2.2. Prior Functions
For mode frequencies, we adopted 14 μHz wide top-hat priors centered on the initial guesses of the frequencies. Top-hat priors were also adopted for the rotational frequency splitting and inclination i⋆, with the inclination sampled from the range of −90° to 180°. The reason for the extended range in inclination is that if the solution is close to either or any sharp truncation from a prior at these values will make it difficult to properly sample these extreme values. The final posterior on the inclination was then obtained from folding the samples onto the range from 0° to 90°.
For the amplitudes and line widths, we adopted a modified Jeffrey's prior, given as (see, e.g., Handberg & Campante 2011)
which behaves as a uniform prior when and a standard scale invariant Jeffrey's prior when . The maximum of the prior occurs at
For mode visibilities, we adopted truncated Gaussian functions as priors, defined as
with D given as:
Here, erf denotes the error function, and σ give the chosen mode value and width of the Gaussian, and and give the lower and upper truncations of the Gaussian. We specifically adopted for l = 1, for l = 2, and for l = 3. Similarly, we adopted truncated Gaussian priors for the parameters of the background. Here we used as the Gaussian mode value () the median of the posteriors from the background-only fit (see Section 3.2.1), and we adopted , , and .
To ensure that the mode identification did not swap for neighboring l = 0 and 2 modes, which is a risk especially at high frequency where the small frequency separation is small compared to the mode line width, we added the prior constraint that on that it must be positive. In principle, could be negative in the event of bumped l = 2 modes; however, because the stars were screened for bumped l = 1 modes and the strength of an avoided l = 2 crossing is expected to be lower than that of a l = 1 mode due to the larger evanescent region (see, e.g., Aerts et al. 2010; Deheuvels & Michel 2011), we do not expect values of for these stars.
3.2.3. Quality Assurance
For each fitted mode, we computed a metric for the quality of the fit in the same manner as detailed in Davies et al. (2016), see also Appourchaux (2004) and Appourchaux et al. (2012b). Briefly, we first ran a fast null hypothesis (H0) test to identify which modes had unambiguous detections, and for which the probability of detection conditioned on the data D needed to be explicitly determined and evaluated. This was done because the explicit determination of is computationally expensive. In the fast H0 test it was assessed whether the S/N in the background corrected power spectrum (D) for a given proposed mode, with the power binned across a number of frequencies to account for the spread in power from the mode line width, was consistent with a pure noise spectrum or whether the H0 hypothesis could be rejected at the level (Appourchaux 2003a, 2004; Lund et al. 2012). When the high S/N modes had been identified in this manner the probability of detection was computed for the remainder low-S/N modes.
In the computation of , both the probability of D assuming H0, , and the probability of the alternative hypothesis H1 of a detected mode, , need to be estimated. The latter was assessed by integrating the probability of measuring the data given a model over a range of mode parameters —this integration was achieved by marginalizing over , with the parameter space sampled over the posteriors from the peak-bagging using emcee. Specifically, a mixture-model was used in which both and were optimized simultaneously to give , the probability of observing the data given the model of a given set of modes with parameters :
Here, the parameter pa, ranging between 0 and 1, then gives the probability of the detection of the given set of modes. This probability was kept free in the emcee run, and in the end was assessed from the posterior distribution of pa (Hogg et al. 2010; Farr et al. 2015). We finally report the Bayes factor K, given as the median of the posterior probability distributions of the natural logarithm of the ratio of over , as
The value of can then be assessed qualitatively on the Kass & Raftery (1995) scale as follows.
For a detailed account of the quality control we refer to Davies et al. (2016).
3.3. Derived Quantities and Correlations
Besides the parameters included in the model of the power spectrum, we computed parameters for derived quantities, such as frequency difference ratios. First, we derived the frequency ratios defined as (Roxburgh & Vorontsov 2003)
Here, d01 and d10 are the smooth five-point small frequency separations defined as
and the large separation is
These ratios are useful for model fitting, where they can be used instead of individual frequencies (Lebreton & Goupil 2014; Silva Aguirre et al. 2015) because they are largely insensitive to the stellar surface layers (Kjeldsen et al. 2008b; Ball & Gizon 2014; Roxburgh 2015; Ball et al. 2016). For further details on the use of these ratios, we refer to Roxburgh & Vorontsov (2003, 2013), Roxburgh (2005), Otí Floranes et al. (2005), and Silva Aguirre et al. (2011, 2013).
Second, we calculated the second differences:
which are useful for studying acoustic glitches from the base of the convection zone and the position of the second helium ionization zone (see, e.g., Basu et al. 1994, 2004; Houdek & Gough 2007; Mazumdar et al. 2014). Figure 5 gives an example of the second differences for KIC 6225718 (Saxo2), together with the best-fitting glitch model from G. Houdek et al. (2017, in preparation). The second differences shown in Figure 5 were computed using Equation (20) on the frequencies given in Table 6, and frequency uncertainties were taken as the average of the asymmetric uncertainties.
Download figure:
Standard image High-resolution imageIn computing these derived quantities we used the full posterior probability distributions (PPDs) of the individual frequencies entering in the descriptions, rather than simply using the median value for the PPD of a given frequency. This ensured that any asymmetries, and deviations from a Gaussian shape, in general, that might be for the PPDs of the individual frequencies were properly propagated to the description of the derived quantity. The final value and credible interval were then computed from the distribution of the quantity in the same manner as for the parameters describing the model power spectrum. Using the full distribution also allowed us to easily compute the correlations between the above frequency differences and ratios, such that these might be included as a covariance matrix in any fit to the quantities.
The parameter correlations were calculated in a robust way using the median absolute deviation (MAD) correlation coefficient (see Pasman & Shevlyakov 1987; Shevlyakov & Smirnov 2011). The MAD estimator is given by the median of the absolute deviation around the median. We opted for instead of the standard Pearson product-moment correlation coefficient, because the latter would be very susceptible to even a single walker in the MCMC optimization straying away from the stationary solution. The between two parameters x and y was calculated as follows.
where u and v are the robust principle variables for x and y:
Figure 6 gives an example of the correlation matrix between the frequency ratios for a given star in the form of a Hinton diagram (Hinton et al. 1986).
Download figure:
Standard image High-resolution image4. Results
Below we present some of the conclusions that can be drawn on the different parameters extracted from the peak-bagging. Results on the rotational splittings and inclination angles will be presented in a separate paper. All results will be made available on the KASOC database (see footnote 14).
4.1. Mode Frequencies
4.1.1. Frequency Uncertainties
A proper understanding of the uncertainties on mode frequencies is important, because they will ultimately limit the precision with which stellar parameters can be estimated from modeling the individual frequencies. We may compare the frequency uncertainties from the peak-bagging with expectations from an analytical maximum likelihood (ML) approach (see, e.g., Libbrecht 1992; Toutain & Appourchaux 1994; Ballot et al. 2008). It should be noted that the ML estimator (if unbiased) reaches the minimum variance bound, in accordance with the Cramér-Rao theorem (Cramér 1946; Rao 1945). Thus one should expect uncertainties at least as large as those from the ML estimator (MLE). The standard way of obtaining uncertainties on an ML estimator is from inverting the negative Hessian matrix, and therefore the standard parameters come from the diagonal elements of the resulting variance-covariance matrix. The Hessian itself is obtained from the matrix of second derivatives of the log-likelihood function with respect to the parameters. Assuming an isolated mode as described by Equation (1) and a likelihood function as given in Equation (9), one may follow Libbrecht (1992) and Toutain & Appourchaux (1994) in defining a theoretical Hessian corresponding to the average of a large number of realizations. From this, the predicted frequency uncertainty for an isolated mode of a given l is given as (Ballot et al. 2008)
where (the signal-to-noise ratio) is the level of the background divided by the mode height; i⋆ gives the stellar inclination; xs is the reduced splitting, given by and is the observing duration. For radial modes (l = 0), the factor fl only depends on β and is given by (Libbrecht 1992)
For the general version of , see Ballot et al. (2008).
In Figure 7, we show the individual frequency uncertainties obtained from the peak-bagging as a function of , i.e., the signal-to-noise ratio of the mode. The left panel shows the behavior for the radial (l = 0) modes. We see that the uncertainties for a given value of depend on β as expected from Equation (25), and for a given β a clear dependence is seen as a function of . Given the relatively small spread in , the vertical gradient largely depicts the gradient in Γ and, as expected, the largest values of are seen toward low values because of the correlation between Γ and mode heights. Comparing with Equation (24), we find that the median uncertainties from our MCMC peak-bagging are ∼1.2 times larger than those predicted from the ML estimator.
Download figure:
Standard image High-resolution imageThe right panel of Figure 7 shows the uncertainties for all extracted modes. The uncertainties for each angular degree again follow the expected trend against β. Some extra scatter is expected from i⋆ and xs, but overall the l > 0 modes are seen to follow the trend of the l = 0 ones. The kernel density estimates of the uncertainties for a given l, obtained by representing each sample with a Gaussian kernel, show that the uncertainties of dipole (l = 1) modes are overall the lowest, followed by l = 0 and l = 2. This is expected because (Equation (2)), hence l = 1 modes will typically have the highest S/N. The l = 3 are seen not to follow this trend, probably because only the very highest amplitude l = 3 were selected for fitting.
It is reassuring to see that the measured uncertainties follow expectations to this level, also considering that Equation (24) is constructed for an isolated mode without factoring in potential contributions from close-by neighboring modes. We note that the demonstration of the agreement is useful for predicting uncertainty yields for future missions like TESS (Ricker et al. 2014) and PLATO (Rauer et al. 2013), and could potentially be used to adjust uncertainties derived from a much faster MLE fitting (see, e.g., Appourchaux et al. 2012b) in cases where the sheer number of stars investigated would render the relatively slow MCMC approach impractical.
4.1.2. Average Seismic Parameters
For each star, we computed the average seismic parameters, including and , from the mode frequencies (see Table 2). The value of was obtained by fitting a Gaussian function to the extracted amplitudes as a function of frequency. The value found for the solar of 3078 ± 13 μHz is in agreement with the value of by Huber et al. (2011). We obtained and from a fit to the extracted mode frequencies with an extended version of the asymptotic relation, as in Lund et al. (2014a; see also Mosser et al. 2011, 2013). This fit was made in a Bayesian manner using emcee, with a likelihood function assuming Gaussian frequency errors, and with final parameter values and uncertainties given by the posterior medians and HPDs. We initially adopted the following formula:
Here gives the value of the large separation at The latter corresponds to the radial order which may take on a non-integer value. Specifically, is obtained for a given l from interpolating the measured mode frequencies () against radial order (nl) to the frequency of This description gives a direct estimate of, for instance, the small frequency separations , given as the amount by which l = 1 modes are offset from the neighboring l = 0 modes (), and (see, e.g., Bedding 2011). The value of is a good probe of the evolutionary state of the star because it is sensitive to the sound-speed gradient of the core, which in turn varies with the composition. The frequency dependence of these frequency separations is captured by the change in D0, assumed to change linearly with n (or frequency), as found, for instance, for the Sun by Elsworth et al. (1990; see also Anguera Gubau et al. 1992; Toutain & Fröehlich 1992). Lastly, the change in the large separation, assumed quadratic in n, is included and mimics the overall curvature of the ridges in the échelle diagram.
Table 2. Values from the Fit of Equation (27) to the Mode Frequencies (see Figure 8)
KIC | ||||||||
---|---|---|---|---|---|---|---|---|
(μHz) | (μHz) | (μHz) | (μHz) | |||||
1435467 | 1406.7 | 70.369 | 0.223 | 1.114 | 2.907 | −0.126 | 5.682 | −0.116 |
2837475 | 1557.6 | 75.729 | 0.179 | 0.911 | 2.235 | −0.017 | 6.417 | −0.126 |
3427720 | 2737.0 | 120.068 | 0.287 | 1.356 | 3.487 | −0.053 | 10.186 | −0.124 |
3456181 | 970.0 | 52.264 | 0.216 | 0.988 | 3.551 | −0.075 | 4.370 | −0.030 |
3632418 | 1166.8 | 60.704 | 0.232 | 1.114 | 3.697 | −0.041 | 4.189 | −0.053 |
3656476 | 1925.0 | 93.194 | 0.207 | 1.445 | 4.608 | −0.095 | 4.554 | −0.157 |
3735871 | 2862.6 | 123.049 | 0.260 | 1.325 | 3.823 | −0.004 | 11.028 | −0.077 |
4914923 | 1817.0 | 88.531 | 0.233 | 1.377 | 4.534 | −0.119 | 5.378 | −0.118 |
5184732 | 2089.3 | 95.545 | 0.216 | 1.374 | 2.772 | −0.104 | 5.863 | −0.097 |
5773345 | 1101.2 | 57.303 | 0.257 | 1.077 | 1.754 | −0.093 | 4.131 | −0.047 |
5950854 | 1926.7 | 96.629 | 0.208 | 1.431 | 4.988 | −0.212 | 4.713 | −0.232 |
6106415 | 2248.6 | 104.074 | 0.254 | 1.343 | 3.422 | −0.114 | 6.881 | −0.118 |
6116048 | 2126.9 | 100.754 | 0.258 | 1.336 | 3.687 | −0.143 | 6.034 | −0.155 |
6225718 | 2364.2 | 105.695 | 0.274 | 1.225 | 3.207 | −0.053 | 8.741 | −0.062 |
6508366 | 958.3 | 51.553 | 0.223 | 1.006 | 2.880 | −0.088 | 2.535 | −0.108 |
6603624 | 2384.0 | 110.128 | 0.201 | 1.492 | 2.801 | −0.159 | 4.944 | −0.201 |
6679371 | 941.8 | 50.601 | 0.181 | 0.880 | 2.861 | −0.016 | 3.143 | 0.044 |
6933899 | 1389.9 | 72.135 | 0.255 | 1.319 | 5.314 | 0.013 | 4.910 | −0.063 |
7103006 | 1167.9 | 59.658 | 0.211 | 0.978 | 2.504 | −0.045 | 4.471 | 0.031 |
7106245 | 2397.9 | 111.376 | 0.246 | 1.392 | 3.489 | −0.113 | 6.529 | −0.265 |
7206837 | 1652.5 | 79.131 | 0.253 | 1.054 | 2.106 | −0.054 | 6.619 | −0.094 |
7296438 | 1847.8 | 88.698 | 0.242 | 1.358 | 4.505 | −0.055 | 5.079 | −0.135 |
7510397 | 1189.1 | 62.249 | 0.258 | 1.112 | 3.971 | 0.003 | 4.370 | −0.016 |
7680114 | 1709.1 | 85.145 | 0.238 | 1.368 | 5.039 | −0.043 | 4.980 | −0.108 |
7771282 | 1465.1 | 72.463 | 0.368 | 1.117 | 3.527 | −0.118 | 5.058 | −0.020 |
7871531 | 3455.9 | 151.329 | 0.285 | 1.504 | 2.142 | −0.149 | 7.350 | −0.143 |
7940546 | 1116.6 | 58.762 | 0.217 | 1.075 | 3.985 | −0.002 | 4.346 | 0.012 |
7970740 | 4197.4 | 173.541 | 0.272 | 1.455 | 2.356 | −0.097 | 7.901 | −0.268 |
8006161 | 3574.7 | 149.427 | 0.195 | 1.547 | 3.061 | −0.084 | 9.680 | −0.150 |
8150065 | 1876.9 | 89.264 | 0.403 | 1.163 | 3.027 | 0.036 | 6.357 | 0.012 |
8179536 | 2074.9 | 95.090 | 0.277 | 1.153 | 3.137 | −0.082 | 8.245 | −0.041 |
8228742 | 1190.5 | 62.071 | 0.244 | 1.158 | 4.371 | 0.007 | 4.517 | −0.048 |
8379927 | 2795.3 | 120.288 | 0.232 | 1.311 | 3.676 | −0.058 | 10.932 | −0.062 |
8394589 | 2396.7 | 109.488 | 0.234 | 1.267 | 3.382 | −0.061 | 7.979 | −0.106 |
8424992 | 2533.7 | 120.584 | 0.120 | 1.517 | 2.678 | −0.179 | 5.190 | −0.282 |
8694723 | 1470.5 | 75.112 | 0.296 | 1.113 | 5.339 | 0.005 | 5.879 | 0.012 |
8760414 | 2455.3 | 117.230 | 0.295 | 1.400 | 4.403 | −0.280 | 5.132 | −0.291 |
8938364 | 1675.1 | 85.684 | 0.235 | 1.444 | 6.491 | −0.046 | 5.184 | −0.119 |
9025370 | 2988.6 | 132.628 | 0.205 | 1.475 | 3.099 | −0.066 | 9.141 | −0.126 |
9098294 | 2314.7 | 108.894 | 0.251 | 1.439 | 3.331 | −0.185 | 5.265 | −0.228 |
9139151 | 2690.4 | 117.294 | 0.240 | 1.337 | 3.557 | −0.019 | 10.050 | −0.093 |
9139163 | 1729.8 | 81.170 | 0.241 | 1.007 | 2.079 | −0.024 | 6.213 | 0.040 |
9206432 | 1866.4 | 84.926 | 0.135 | 0.958 | 3.235 | 0.001 | 7.115 | 0.004 |
9353712 | 934.3 | 51.467 | 0.254 | 1.095 | 3.536 | −0.038 | 3.907 | −0.090 |
9410862 | 2278.8 | 107.390 | 0.223 | 1.343 | 3.625 | −0.181 | 6.098 | −0.239 |
9414417 | 1155.3 | 60.115 | 0.237 | 1.045 | 3.572 | −0.029 | 4.648 | −0.006 |
9812850 | 1255.2 | 64.746 | 0.240 | 1.067 | 2.654 | −0.134 | 4.418 | −0.032 |
9955598 | 3616.8 | 153.283 | 0.195 | 1.529 | 2.796 | −0.095 | 8.941 | −0.172 |
9965715 | 2079.3 | 97.236 | 0.373 | 1.139 | 3.685 | −0.102 | 7.958 | −0.105 |
10068307 | 995.1 | 53.945 | 0.247 | 1.131 | 4.220 | 0.014 | 3.799 | −0.041 |
10079226 | 2653.0 | 116.345 | 0.264 | 1.350 | 3.236 | −0.014 | 9.387 | 0.098 |
10162436 | 1052.0 | 55.725 | 0.242 | 1.106 | 3.746 | −0.033 | 3.706 | −0.031 |
10454113 | 2357.2 | 105.063 | 0.283 | 1.206 | 3.059 | −0.063 | 9.426 | −0.079 |
10516096 | 1689.8 | 84.424 | 0.257 | 1.318 | 5.001 | −0.058 | 5.248 | −0.089 |
10644253 | 2899.7 | 123.080 | 0.250 | 1.313 | 3.863 | 0.004 | 11.378 | −0.138 |
10730618 | 1282.1 | 66.333 | 0.239 | 1.032 | 2.651 | 0.019 | 4.556 | 0.113 |
10963065 | 2203.7 | 103.179 | 0.297 | 1.275 | 3.567 | −0.081 | 7.083 | −0.058 |
11081729 | 1968.3 | 90.116 | 0.242 | 1.020 | 3.056 | −0.103 | 6.602 | 0.010 |
11253226 | 1590.6 | 76.858 | 0.183 | 0.920 | 1.748 | −0.132 | 6.973 | 0.039 |
11772920 | 3674.7 | 157.746 | 0.238 | 1.516 | 2.366 | −0.082 | 7.849 | −0.131 |
12009504 | 1865.6 | 88.217 | 0.289 | 1.200 | 3.533 | −0.072 | 6.117 | −0.066 |
12069127 | 884.7 | 48.400 | 0.204 | 1.061 | 3.399 | −0.037 | 3.650 | −0.102 |
12069424 | 2188.5 | 103.277 | 0.246 | 1.437 | 3.392 | −0.152 | 5.274 | −0.181 |
12069449 | 2561.3 | 116.929 | 0.201 | 1.461 | 2.707 | −0.117 | 6.045 | −0.133 |
12258514 | 1512.7 | 74.799 | 0.209 | 1.281 | 4.227 | −0.061 | 4.827 | −0.056 |
12317678 | 1212.4 | 63.464 | 0.231 | 0.928 | 3.883 | −0.032 | 5.273 | −0.061 |
Note. is obtained from the fit to the extracted modes. Plots of the fitted parameters are shown in Figures 9–12; the values for are shown in Figure 3.
We found this description to perform poorly across the range of stars in the sample, which is to be expected for more evolved stars (see, e.g., Gabriel 1989; Christensen-Dalsgaard 1991). To that end, we modified the formula as follows.
where the term has been replaced by , which takes on independent values for l > 0. Thereby, we optimize independently for the separations and . In the fit, we also included with a Gaussian prior from the fit to the mode amplitudes; the values of , and correlate strongly, but we include to properly marginalize over the uncertainty of the pivoting . Figure 8 gives an example of the fit of Equation (27) to KIC 8228742 (Horace). All parameters from the fits are listed in Table 2.
Download figure:
Standard image High-resolution imageFigure 9 shows the estimated values for versus together with the empirical relation by Huber et al. (2011), with which we find an excellent agreement. We can even see the expected mass gradient in the residuals, with higher mass stars having slightly lower for a given Also shown is the change in the large separation, which is found to be near-constant at and thus always with the same concavity sign, corresponding to a positive gradient of the large separation with frequency. Here one should remember that a constant value of against would correspond to a linear change in with frequency.
Download figure:
Standard image High-resolution imageFor the Sun, we obtain a value of of 134.91 ± 0.02 μHz; comparing this to the value of by Huber et al. (2011), we would expect a slightly smaller value, because the curvature is incorporated in our estimation. However, the value of by Toutain & Fröehlich (1992), who also used a second-order version of the asymptotic relation, compares very well to our estimate.
The estimated values of have already been shown in Figure 3. For the estimation of , it should be noted that this value is strongly anti-correlated with , and a small change in can thus induce a significant change in (White et al. 2011b). This will especially take effect if the estimate of is off, because defines the pivoting point for the change in , , and the curvatures. Like White et al. (2011b), we see an offset between the estimated and those obtained from model tracks, which is ascribed to the effects of the incorrect modeling of the stellar surface layers. In Figure 3, we show the -tracks adopted from White et al. (2011a), which were computed from evolutionary tracks from the Aarhus STellar Evolution Code (ASTEC; Christensen-Dalsgaard 2008), neglecting diffusion and core overshoot and with Z0 = 0.017. It should also be noted that the from models were derived in a slightly different manner than used here, as described by White et al. (2011a). We therefore also estimated the values for and using the White et al. (2011a) method, namely from a weighted fit of the asymptotic function in Equation (10) to the radial mode frequencies as a function of radial order n, and with weights given by a Gaussian with an FWHM of 0.25 Comparing the values for and from the fit of Equation (10) versus Equation (27), we only find minor differences in estimated values. For the maximum absolute difference was ∼0.36 μHz, with no systematic differences; for a maximum difference of ∼0.14 was found, and again with no systematic differences.
In Figures 10 and 11, we show the estimated values for the separations and , together with their gradients in n (or frequency); the results for in Figure 11 are given in the form of a modified C–D diagram (see Christensen-Dalsgaard 1993; White et al. 2011a, 2011b, 2012). Because all values for the changes in and are either zero or negative, we see that the small separations are virtually all decreasing functions of frequency. In Figure 11, we show again the ASTEC tracks from White et al. (2011a). It is interesting to see the degree to which the asymptotic relations to first (Equation (10)) and second (Equation (26)) order are satisfied. From these, one would expect a ratio of from assuming —in Figure 12, we show this ratio for the sample. As seen, only a few stars, including the Sun, actually adhere to the expectation from the asymptotic relation. This was similarly found by Christensen-Dalsgaard & Houdek (2010)19 for α Cen A (Bedding et al. 2004), while α Cen B (Kjeldsen et al. 2005) was found to fulfill the asymptotic relation. As explained by Christensen-Dalsgaard & Houdek (2010), this reflects a rapid variation of the sound speed in the core of α Cen A from its more evolved state due to its higher mass compared to the B component. Similarly, Verma et al. (2014) found an increasing departure from the asymptotic description of the oscillation frequencies while fitting the signatures of the acoustic glitches. The departure becomes noticeable when a peak in the Brunt–Väisälä frequency N develops just outside the stellar core toward the end of the MS and becomes comparable or higher than the lowest frequency fitted. To check this evolutionary explanation, we plot in the right panel of Figure 12 the against the central hydrogen content Xc, because this serves as a good probe for the evolutionary state. As seen, stars with high Xc, i.e., stars that are less evolved, adhere better to the asymptotic regime than the evolved stars with low Xc. The values of Xc were obtained as part of the modeling with the BASTA pipeline; we refer to Paper II for further details. While a difference from a ratio of can be explained from an evolutionary viewpoint, it is interesting to note that the majority of the stars analyzed indeed do not follow the asymptotic relation. This is a cautionary note to the use of Equation (10) and/or Equation (26) for extracting average seismic parameters as here, or if used for predicting the location of oscillation modes frequencies in the power spectrum from measured in an independent manner. For details on the physics behind the curvatures, we refer to Tassoul (1980), Houdek & Gough (2007), Cunha & Metcalfe (2007), and Mosser et al. (2013) and references therein.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution image4.2. Mode Amplitudes
Mode amplitudes were measured in the peak-bagging as described in Section 3.1. Figure 13 displays the radial degree five-point Epanechnikov (Epanechnikov 1969; Hastie et al. 2009) smoothed amplitudes for all the stars in the sample against frequency (left) and normalized against distance from (right). One clearly sees the overall decrease in amplitude with increasing (Kjeldsen et al. 2005; Arentoft et al. 2008), and thus decreasing (Kjeldsen & Bedding 1995, 2011; Appourchaux et al. 2014). The value of was estimated from the amplitudes as the frequency where the modes of oscillation show their maximum amplitude,
Download figure:
Standard image High-resolution imageFigure 14 gives the measured values for against (Table 3). We see both the expected change with in addition to a mass gradient across the overall decrease in with (Huber et al. 2011). In Figure 14, we also give a comparison to the amplitudes estimated from the scaling relations of Huber et al. (2011):
with t = 1.32 ± 0.02, s = 0.838 ± 0.002, and r = 2 and Corsaro et al. (2013):
with , t = 1.27 ± 0.04, r = 3.47 ± 0.09, and . In converting the bolometric amplitudes from these relations to the Kepler bandpass, we used the root-mean-square value of from Michel et al. (2009) and the temperature dependent bolometric correction by Ballot et al. (2011), which is specific to the Kepler spectral response function. Overall, we see a reasonable agreement with a scatter within approximately 25% as also found by Huber et al. (2011). At low amplitudes, the Huber et al. (2011) relation is seen to provide the best agreement, while the Corsaro et al. (2013) relation has an overall lower scatter across the amplitude range. We note that the relations tested here were calibrated against values obtained using the method of Kjeldsen et al. (2005, 2008a) and the high range occupied by the stars in the current analysis is sparsely covered in the calibrations of the scaling relations, so a better relationship than observed cannot readily be expected (see more below). In addition, amplitudes will have natural scatter due to the impact of activity (Chaplin et al. 2011a).
Download figure:
Standard image High-resolution imageTable 3. Values for the Line Width at (Figure 20) from the Fit of Equation (30)
KIC | α | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
1435467 | 5.18 ± 0.17 | 4.37 | 4.30 | 4.25 | 0.54 | 5.19 ± 0.18 | ⋯ | ⋯ | ⋯ | ⋯ |
2837475 | 6.39 ± 0.2 | 3.87 | 3.67 | 3.71 | 0.91 | 6.42 ± 0.21 | ⋯ | ⋯ | ⋯ | ⋯ |
3427720 | 1.88 ± 0.15 | 2.62 | 2.55 | 2.38 | 2.49 | 5.0 ± 2.96 | 0.36 | 4428 | 2695 | 1309 ± 711 |
3456181 | 4.07 ± 0.19 | 5.60 | 5.92 | 5.48 | 1.11 | 4.06 ± 0.2 | ⋯ | ⋯ | ⋯ | ⋯ |
3632418 | 2.72 ± 0.16 | 5.04 | 5.04 | 4.76 | 1.76 | 4.67 ± 0.71 | 0.56 | 1698 | 1188 | 448 ± 183 |
3656476 | 0.69 ± 0.01 | 4.70 | 4.13 | 4.12 | 2.50 | 0.69 ± 0.01 | ⋯ | ⋯ | ⋯ | ⋯ |
3735871 | 2.25 ± 0.25 | 2.52 | 2.38 | 2.30 | 2.79 | 7.34 ± 4.29 | 0.28 | 4717 | 2825 | 1426 ± 595 |
4914923 | 1.16 ± 0.08 | 4.62 | 4.19 | 4.15 | 3.52 | 6.2 ± 2.4 | 0.18 | 2600 | 1847 | 666 ± 170 |
5184732 | 1.36 ± 0.06 | 4.08 | 3.84 | 3.88 | 3.52 | 7.65 ± 2.61 | 0.17 | 3247 | 2117 | 941 ± 184 |
5773345 | 3.37 ± 0.13 | 5.69 | 5.17 | 5.52 | 0.58 | 3.37 ± 0.14 | ⋯ | ⋯ | ⋯ | ⋯ |
5950854 | 0.83 ± 0.01 | 4.12 | 3.77 | 3.87 | 0.56 | 0.83 ± 0.01 | ⋯ | ⋯ | ⋯ | ⋯ |
6106415 | 1.64 ± 0.07 | 3.86 | 3.60 | 3.63 | 3.49 | 5.53 ± 1.08 | 0.29 | 3242 | 2274 | 837 ± 166 |
6116048 | 1.62 ± 0.09 | 4.02 | 3.75 | 3.67 | 3.26 | 5.74 ± 1.35 | 0.28 | 3102 | 2098 | 796 ± 182 |
6225718 | 2.58 ± 0.12 | 3.53 | 3.39 | 3.29 | 2.32 | 5.77 ± 0.67 | 0.43 | 3076 | 2316 | 611 ± 121 |
6508366 | 5.22 ± 0.17 | 5.02 | 5.21 | 4.75 | 1.95 | 5.21 ± 0.15 | ⋯ | ⋯ | ⋯ | ⋯ |
6603624 | 0.56 ± 0.02 | 4.15 | 3.55 | 3.58 | 4.42 | 0.57 ± 0.02 | ⋯ | ⋯ | ⋯ | ⋯ |
6679371 | 4.53 ± 0.18 | 5.56 | 5.33 | 5.28 | 0.45 | 5.51 ± 0.74 | 0.69 | 1339 | 833 | 295 ± 190 |
6933899 | 1.3 ± 0.07 | 5.55 | 5.07 | 4.94 | 3.15 | 6.0 ± 2.75 | 0.21 | 2218 | 1422 | 671 ± 178 |
7103006 | 5.09 ± 0.16 | 4.74 | 4.29 | 4.58 | 1.17 | 5.09 ± 0.17 | ⋯ | ⋯ | ⋯ | ⋯ |
7106245 | 1.65 ± 0.14 | 3.39 | 2.96 | 2.91 | 1.71 | 1.64 ± 0.14 | ⋯ | ⋯ | ⋯ | ⋯ |
7206837 | 4.34 ± 0.32 | 4.32 | 3.85 | 4.18 | 0.75 | 8.25 ± 2.83 | 0.50 | 2854 | 1720 | 952 ± 433 |
7296438 | 1.22 ± 0.12 | 4.59 | 4.24 | 4.11 | 3.87 | 8.39 ± 4.56 | 0.14 | 2794 | 1835 | 764 ± 205 |
7510397 | 2.42 ± 0.15 | 3.62 | 3.52 | 3.40 | 1.90 | 4.18 ± 0.71 | 0.55 | 1663 | 1226 | 413 ± 176 |
7680114 | 1.16 ± 0.09 | 4.71 | 4.40 | 4.40 | 3.70 | 4.68 ± 1.95 | 0.24 | 2438 | 1736 | 620 ± 191 |
7771282 | 3.29 ± 0.28 | 4.20 | 4.14 | 3.96 | 0.71 | 3.29 ± 0.29 | ⋯ | ⋯ | ⋯ | ⋯ |
7871531 | 1.21 ± 0.12 | 1.91 | 1.71 | 1.77 | 1.94 | 4.83 ± 3.18 | 0.20 | 6026 | 3161 | 1780 ± 653 |
7940546 | 2.92 ± 0.14 | 5.29 | 5.24 | 5.51 | 1.56 | 5.15 ± 0.69 | 0.53 | 1578 | 1177 | 409 ± 126 |
7970740 | 1.99 ± 0.15 | 1.65 | 1.35 | 1.51 | 4.25 | 6.85 ± 3.6 | 0.26 | 6875 | 3910 | 1949 ± 735 |
8006161 | 1.17 ± 0.06 | 1.94 | 1.69 | 1.73 | 4.75 | 6.0 ± 3.27 | 0.19 | 5581 | 3507 | 1575 ± 478 |
8150065 | 2.37 ± 0.3 | 3.55 | 3.38 | 3.23 | 2.01 | 2.36 ± 0.28 | ⋯ | ⋯ | ⋯ | ⋯ |
8179536 | 3.51 ± 0.26 | 3.42 | 3.35 | 3.25 | 1.88 | 7.05 ± 3.12 | 0.47 | 3940 | 2112 | 1378 ± 486 |
8228742 | 2.05 ± 0.13 | 5.30 | 5.20 | 4.94 | 2.39 | 4.19 ± 0.68 | 0.48 | 1640 | 1189 | 383 ± 117 |
8379927 | 2.43 ± 0.11 | 2.21 | 2.09 | 2.05 | 2.45 | 6.12 ± 1.22 | 0.39 | 4091 | 2731 | 1046 ± 258 |
8394589 | 2.15 ± 0.16 | 3.49 | 3.33 | 3.24 | 1.71 | 5.79 ± 1.69 | 0.28 | 3185 | 2233 | 637 ± 208 |
8424992 | 1.14 ± 0.1 | 2.97 | 2.77 | 2.67 | 2.99 | 1.14 ± 0.1 | ⋯ | ⋯ | ⋯ | ⋯ |
8694723 | 3.08 ± 0.16 | 5.38 | 5.08 | 5.07 | 1.98 | 5.42 ± 0.79 | 0.56 | 2129 | 1472 | 548 ± 192 |
8760414 | 1.25 ± 0.08 | 3.96 | 3.49 | 3.52 | 2.97 | 5.06 ± 2.23 | 0.23 | 3659 | 2349 | 943 ± 277 |
8938364 | 0.8 ± 0.01 | 5.22 | 4.61 | 4.43 | 2.49 | 0.8 ± 0.01 | ⋯ | ⋯ | ⋯ | ⋯ |
9025370 | 1.39 ± 0.08 | 1.64 | 1.45 | 1.45 | 1.60 | 1.38 ± 0.09 | ⋯ | ⋯ | ⋯ | ⋯ |
9098294 | 1.3 ± 0.08 | 3.55 | 3.21 | 3.16 | 2.81 | 7.81 ± 4.74 | 0.16 | 3986 | 2317 | 1275 ± 355 |
9139151 | 1.99 ± 0.14 | 2.88 | 2.72 | 2.66 | 3.03 | 7.78 ± 3.84 | 0.25 | 4136 | 2697 | 1169 ± 399 |
9139163 | 5.28 ± 0.13 | 3.76 | 3.64 | 3.81 | 1.79 | 5.28 ± 0.13 | ⋯ | ⋯ | ⋯ | ⋯ |
9206432 | 5.87 ± 0.27 | 3.54 | 3.47 | 3.48 | 0.88 | 5.88 ± 0.28 | ⋯ | ⋯ | ⋯ | ⋯ |
9353712 | 3.22 ± 0.22 | 5.54 | 5.68 | 5.46 | 1.20 | 3.2 ± 0.21 | ⋯ | ⋯ | ⋯ | ⋯ |
9410862 | 1.89 ± 0.25 | 3.75 | 3.60 | 3.67 | 2.40 | 7.2 ± 3.8 | 0.24 | 3337 | 2204 | 846 ± 414 |
9414417 | 3.47 ± 0.23 | 5.55 | 5.18 | 5.29 | 1.04 | 5.81 ± 1.31 | 0.57 | 1910 | 1191 | 605 ± 299 |
9812850 | 5.68 ± 0.19 | 4.61 | 4.51 | 4.85 | 0.69 | 5.67 ± 0.2 | ⋯ | ⋯ | ⋯ | ⋯ |
9955598 | 0.77 ± 0.03 | 2.07 | 1.79 | 1.94 | 4.98 | 0.77 ± 0.01 | ⋯ | ⋯ | ⋯ | ⋯ |
9965715 | 3.14 ± 0.34 | 3.89 | 3.54 | 3.59 | 1.82 | 5.92 ± 1.9 | 0.47 | 3024 | 2103 | 792 ± 528 |
10068307 | 2.15 ± 0.13 | 5.98 | 5.86 | 5.54 | 2.15 | 3.91 ± 0.45 | 0.53 | 1388 | 1015 | 340 ± 102 |
10079226 | 2.05 ± 0.35 | 2.81 | 2.63 | 2.63 | 3.36 | 6.73 ± 4.24 | 0.27 | 4248 | 2601 | 1236 ± 585 |
10162436 | 3.08 ± 0.08 | 5.38 | 5.40 | 5.42 | 1.82 | 3.08 ± 0.09 | ⋯ | ⋯ | ⋯ | ⋯ |
10454113 | 4.12 ± 0.15 | 3.19 | 2.87 | 3.01 | 1.60 | 4.13 ± 0.15 | ⋯ | ⋯ | ⋯ | ⋯ |
10516096 | 1.56 ± 0.1 | 4.74 | 4.40 | 4.26 | 3.44 | 5.33 ± 2.47 | 0.28 | 2629 | 1723 | 768 ± 298 |
10644253 | 2.21 ± 0.22 | 2.22 | 2.15 | 2.06 | 3.16 | 11.06 ± 4.51 | 0.19 | 4253 | 2852 | 1099 ± 282 |
10730618 | 5.07 ± 0.25 | 4.64 | 4.24 | 4.49 | 0.41 | 5.09 ± 0.27 | ⋯ | ⋯ | ⋯ | ⋯ |
10963065 | 2.16 ± 0.11 | 3.72 | 3.62 | 3.45 | 3.36 | 8.44 ± 3.9 | 0.25 | 3921 | 2158 | 1261 ± 383 |
11081729 | 6.19 ± 0.3 | 3.11 | 2.99 | 3.01 | 1.79 | 6.19 ± 0.3 | ⋯ | ⋯ | ⋯ | ⋯ |
11253226 | 5.8 ± 0.17 | 4.20 | 3.62 | 3.98 | 0.40 | 5.8 ± 0.17 | ⋯ | ⋯ | ⋯ | ⋯ |
11772920 | 0.77 ± 0.03 | 1.83 | 1.49 | 1.55 | 2.69 | 0.77 ± 0.01 | ⋯ | ⋯ | ⋯ | ⋯ |
12009504 | 2.38 ± 0.14 | 3.94 | 3.84 | 3.65 | 2.31 | 6.06 ± 2.12 | 0.38 | 2915 | 1821 | 819 ± 334 |
12069127 | 3.6 ± 0.23 | 6.00 | 6.00 | 5.79 | 0.93 | 3.61 ± 0.23 | ⋯ | ⋯ | ⋯ | ⋯ |
12069424 | 0.98 ± 0.04 | 3.96 | 3.55 | 3.56 | 3.33 | 5.47 ± 1.11 | 0.18 | 3187 | 2181 | 825 ± 108 |
12069449 | 0.91 ± 0.05 | 3.49 | 3.05 | 3.08 | 4.10 | 6.06 ± 1.56 | 0.15 | 3716 | 2579 | 966 ± 148 |
12258514 | 1.69 ± 0.09 | 4.74 | 4.52 | 4.35 | 2.92 | 4.01 ± 0.41 | 0.42 | 1939 | 1510 | 376 ± 79 |
12317678 | 5.83 ± 0.16 | 4.75 | 4.98 | 4.57 | 1.10 | 5.82 ± 0.15 | ⋯ | ⋯ | ⋯ |
Note. Amax gives the mode amplitude at from a fit to the individual mode amplitudes against frequency; and give the amplitudes from the smoothing method by Kjeldsen et al. (2008a; Figure 15), with different values for the effective number of modes per radial order c; α, Γα, ΔΓdip, Wdip, and give the obtained parameters from the fit of Equation (30). The values of and are only given if the inclusion of the lorentzian in Equation (30) gave a better fit than the power-law-only fit. gives the full-width-half-maximum of the lorentzian line width dip.
4.2.1. Amplitudes from Smoothed Amplitude Spectra
With amplitudes measured from individual modes, it is interesting to see how these compare to those obtained from the often adopted method by Kjeldsen et al. (2005, 2008a). This is especially worthwhile because it is often amplitudes from this method that are extracted by automated analysis pipelines (see, e.g., Huber et al. 2009; Hekker et al. 2010; Mathur et al. 2010) and thus used in calibrating scaling relations (see, e.g., Huber et al. 2011; Stello et al. 2011; Corsaro et al. 2013).
Amplitudes were estimated following the prescription of Kjeldsen et al. (2008a). Here the power density spectrum is first convolved with a Gaussian filter with an FWHM of , to produce a spectrum with a single power hump from the oscillations. A noise background is then fitted to the smoothed spectrum and subtracted, after which the spectrum is multiplied by /c and the square-root is taken to convert to amplitude. The multiplication of /c converts to power per radial mode, with c giving the effective number of radial modes per order, i.e., (see Equation (2)). Concerning the order of the different steps in the method, we subtract instead first the fitted background function from the peak-bagging, and then apply the smoothing—this removes the potential bias on both the amplitude and central frequency of the smoothed power hump from mixing background and oscillation power, and then try to fit the now smoothed underlying background with a non-smoothed background function.
The value of c is often kept fixed rather than estimated for a given star following, for instance, the results by Ballot et al. (2011). Bedding et al. (2010) estimated a value of c = 3.04 for the Sun as seen by Kepler using the method of Bedding et al. (1996) and Kjeldsen et al. (2008a) and adopting a mean observing wavelength of —this value has since been used by, e.g., Huber et al. (2011), Stello et al. (2011), and Corsaro et al. (2013). In our estimation of amplitudes, we tried both the value of c = 3.04 and the one obtained from the visibilities measured in the peak-bagging .
Figure 15 shows the result from this exercise; we find that the values from the above method are systematically offset (fixed bias) from those estimated from individual frequencies, with a median relative offset of −6%; no proportional bias is seen. This offset could in part originate from the somewhat arbitrary choice of for the smoothing window—however, if this was the only contributor, one might expect a proportional offset rather than a constant one. The identified offset fully corroborates the results by Verner et al. (2011) who found a systematic offset between −15 and −2% from pipeline analysis of simulated data. We observe the same median offset for the different values of c, but the scatter is slightly lower for the values using c from the measured visibilities. Considering that and thus c depends on various stellar parameters, such as , , and log g, it should be expected that adopting a constant c adds to the scatter. The values of from the smoothing method, with both values of c are given in Table 3.
Download figure:
Standard image High-resolution image4.3. Mode Line Widths
Following Appourchaux et al. (2014), we adopt the following relation for the line widths Γ against mode frequency:
where the first part describes a power-law trend with an exponent α and a value of at The second part describes a Lorentzian dip in the line widths in , centered on with a width and an depth . For the fit, the values of were calculated from the Γ posterior distributions; the -distribution is typically better approximated by a normal distribution (Toutain & Appourchaux 1994). We note that while Equation (30) matches that reported by Appourchaux et al. (2014), these authors did in fact use a formula with the Lorentzian subtracted rather than added as in Equation (30) (Appourchaux et al. 2016). Both formulations can be used, but the reported parameters will naturally differ; in Equation (30) is for instance constrained to a value between 0 and 1 if a dip in Γ is to be produced, whereas it can take on any value above 1 if the Lorentzian is instead subtracted. We chose the formulation where the Lorentzian is added in log-space because we found that this provided better constrained fits from a lower correlation between the , , and parameters. We note that the parameter has two solutions, one higher than and one lower—in our fits, we chose the solution with a prior on . All stars were fitted using both the full version of Equation (30) and using only the first power-law component—Figure 16 shows for all stars in the sample the measured line widths against a proxy for the radial order. The fit was made using emcee (Foreman-Mackey et al. 2013), with the parameter values and uncertainties given by the posterior distribution medians and 68% HPD intervals. We evaluate which of the two types of fit is the best based on a visual inspection of the fits.
Download figure:
Standard image High-resolution imageFigure 17 shows three examples where the full version of Equation (30) was deemed better—as seen, the line widths show a clear depression around We omitted fitting the solar line widths from the degraded VIRGO data, because in computing relations for the fit parameters we are interested in as accurate and precise data as possible.
Download figure:
Standard image High-resolution imageThe parameters from the fits are given in Table 3 and shown in Figure 19 as a function of and . For the dependencies of the FWHM of the Lorentzian line width dip in frequency units, we used the transformation (Appourchaux et al. 2016). The amplitude of the Lorentzian dip may be calculated using the transformation .20 Like Appourchaux et al. (2014) for their sample of 22 Kepler stars, we find a clear correlation between the Lorentzian width ( or ) and while the power-law exponent (α) and the dip amplitude () are found to correlate most strongly with . Only for the fits without the Lorentzian component does the value of correspond to the line width at which is seen to correlate with . We return to the overall behavior of the line width at against below. Our estimates of agree well with updated values from Appourchaux et al. (2016). We make linear fits to the different parameters (P) against as and as . Specifically, we perform a Deming regression (Deming 1943), where uncertainties in both the dependent () and independent variables () are considered via the following merit function (Press et al. 1993):
which serves as our log-likelihood function in the emcee optimization. The coefficients from the different fits are given in Table 4. For the fit of against we include only the values from the fit including the Lorentzian in Equation (30). Some of the scatter in the relations will likely be due to activity; the dip in line widths for the Sun has for instance been found to depend on the solar cycle (see, e.g., Komm et al. 2000a, 2000b). We note that the fit coefficients in Table 4 may well be used to define prior functions and initial guesses in future peak-bagging exercises.
Table 4. Parameters from Linear Fits to the Fitted Values of Equation (30) Against and
b | a | b | a | |
---|---|---|---|---|
α | 2.95 ± 0.16 | 0.39 ± 0.08 | −25.5 ± 1.4 | 29.1 ± 1.5 |
3.08 ± 0.98 | 3.32 ± 0.50 | 6.3 ± 6.6 | −1.8 ± 6.9 | |
−0.47 ± 0.06 | 0.62 ± 0.04 | 3.5 ± 0.4 | −3.3 ± 0.4 | |
4637 ± 237 | −141 ± 138 | −28021 ± 2964 | 31971 ± 3105 | |
2984 ± 31 | 60 ± 18 | −23818 ± 1909 | 26785 ± 1983 | |
1253 ± 162 | −85 ± 96 | −5649 ± 1093 | 6550 ± 1151 |
Note. The linear fits are indicated as red lines in Figure 19. For the fits of only values from the full fit of Equation (30) were included.
Download table as: ASCIITypeset image
For a star with a temperature of the Sun, we find that it should have an exponent of —this comes close to that measured by Komm et al. (2000b) of α ≈ 3.31 for frequencies below 2450 μHz, which is about the lower limit of what one would observe for a similar star with Kepler. In general, the exponent for the Sun has been measured below the dip where higher values have been found, e.g., α ≈ 5 by Libbrecht & Woodard (1991), α ≈ 7 by Chaplin et al. (1997), and α ≈ 8 by Davies et al. (2014a)—it is therefore not surprising that a lower value of α is found when the measurement is done across the dip. Like Appourchaux et al. (2014), we find that the amplitude of the line width dip decreases with . This behavior is contrary to theoretical damping-rate estimates (Balmforth 1992; Houdek 1996), which assume solar-calibrated convection parameters, such as the mixing-length and anisotropy of the turbulent velocity field, for damping-rate calculations in other stars. Updated calculations of mode damping rates (G. Houdek et al. 2017, in preparation) are, however, able to capture the overall behavior of the line widths for stars in the current study spanning a large range in —these new calculations adopt, in addition to the standard ingredients described in Houdek et al. (1999) and Houdek (2006), turbulent pressure profiles and T − τ relations calibrated to 3D hydrodynamical simulation results by Trampedach et al. (2014). An example of such an improved damping-rate computation provides Figure 18, which compares estimated mode line widths with observations for KIC 6933899 (Fred), using stellar parameters from the best-fitting ASTFIT model (Paper II).
Download figure:
Standard image High-resolution imageThe estimate of Γ at was obtained by a Monte Carlo sampling from the posteriors of the fit parameters in Equation (30) and the estimate of from the mode amplitudes. The values obtained are plotted against in Figure 20 and given in Table 3. We have fitted two relations for the line width at namely, the power-law relation by Appourchaux et al. (2012a):
and the exponential relation used, for instance, by Corsaro et al. (2012):
In the fits of Equations (32) and (33), we complement our set of line widths with those from peak-bagging of 42 giants in NGC 6819 by R. Handberg et al. (2017, in preparation), where we adopt an uncertainty of 100 K on for all stars, and 19 red giants by Corsaro et al. (2015), where we adopt values from Pinsonneault et al. (2012; from the infra-red flux method). The line widths from R. Handberg et al. (2017, in preparation) are given by the average of the radial modes; those from Corsaro et al. (2015) are given as the average over three radial modes centered in the mode with the highest amplitude. We fit both relations to the line widths using an orthogonal distance regression (ODR; Boggs & Rogers 1990) in order to take into account both the uncertainties on Γ and . This procedure is employed in a Monte Carlo (MC) run, where we in each iteration draw at random 50 stars to include in the fit. For Equation (32), we obtain from the MC distributions the following values for the fitting parameters: , , and β ≈ 15.3 ± 1.9; for Equation (33) we obtain and . If we fit only the red giants by a constant value for the line width, we find a value of . For the parameter uncertainties, we have added in quadrature the median from the internal uncertainties from the ODR and the standardized MAD of the best-fit ODR values from the MC run. We note here that neglecting the uncertainties in from using, for instance, an ordinary least-squares (OLS) approach will affect the fit and the parameters in Equation (32). Both relations perform reasonably well, but the fit of Equation (32) gives a slightly lower value. The resulting relations are displayed in Figure 20. Also shown in Figure 20 are the fits obtained by Appourchaux et al. (2012a) (parameters from their Table 3) and Corsaro et al. (2012).
We see that the cooler MS targets appear to be outliers to the overall relation between line width and temperature. From Figure 19, we can see that the value of has a proportional bias with respect to hence the higher the the further below does the dip in line widths appear. For Figure 20, the line widths were estimated at and thus away from the line width dip for the highest stars, it is therefore not surprising that these deviate from the overall relation. An offset could potentially also be caused by frequency shifts from stellar activity cycles (Chaplin et al. 2007, 2008a), the evidence of which is currently being studied by A. R. G. Santos et al. (2017, in preparation). The line width for KIC 7970740 was omitted in the fits because it appears to be a particularly strong outlier at K.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageWe find that our estimates agree well with those from Appourchaux et al. (2012a) with the main difference seen at low temperatures, where the fit is the least constrained if no line widths from giants are included. The estimates from Corsaro et al. (2012) appear to be higher at low temperatures than the line widths from R. Handberg et al. (2017, in preparation)—this could be due to the method of estimating Γ in Corsaro et al. (2012) from collapsed échelle diagrams for a group of stars. Some scatter is to be expected, because the line widths will depend on parameters besides —Belkacem et al. (2012, 2013) specifically suggest a scaling with and log g. As noted by Appourchaux et al. (2012a), the uncertainty on currently sets the limitations for obtaining a more well constrained prediction for Γ at which may then be confronted with theoretical calculations. As for the dip amplitude, activity and stellar cycles will play a role in the scatter around the mean relation for Γ at (Komm et al. 2000b; Chaplin et al. 2008a). For other analysis of the variation of Γ at we refer to Goldreich & Kumar (1991), Houdek et al. (1999), Chaplin et al. (2009), Baudin et al. (2011), and Belkacem et al. (2012, 2013).
4.4. Mode Visibilities
Mode visibilities are estimated as part of the peak-bagging (see Equation (2)); these are given in Table 5. Figure 21 shows the obtained visibilities as a function of , and mass. The total visibilities are computed from the combination of the MCMC chains for the individual visibilities. We see no strong correlation with , (not shown), or (not shown), and some scatter is observed for the visibilities about their median values. From median binned values, there are slight indications of structure in the visibilities against the plotted parameters; the visibilities do, for instance, seem to peak around a temperature of ∼6100 K and have a depression around a mass of . None of the dependencies are, however, very strong, but they do appear similar in shape for different angular degrees.
Download figure:
Standard image High-resolution imageTable 5. Extracted Mode Visibilities for Angular Degrees
KIC | ||||
---|---|---|---|---|
1435467 | ⋯ | |||
2837475 | ⋯ | |||
3427720 | ⋯ | |||
3456181 | ⋯ | |||
3632418 | ⋯ | |||
3656476 | ||||
3735871 | ⋯ | |||
4914923 | ||||
5184732 | ||||
5773345 | ⋯ | |||
5950854 | ⋯ | |||
6106415 | ||||
6116048 | ||||
6225718 | ||||
6508366 | ⋯ | |||
6603624 | ||||
6679371 | ⋯ | |||
6933899 | ⋯ | |||
7103006 | ⋯ | |||
7106245 | ⋯ | |||
7206837 | ⋯ | |||
7296438 | ⋯ | |||
7510397 | ⋯ | |||
7680114 | ||||
7771282 | ⋯ | |||
7871531 | ⋯ | |||
7940546 | ⋯ | |||
7970740 | ⋯ | |||
8006161 | ||||
8150065 | ⋯ | |||
8179536 | ⋯ | |||
8228742 | ⋯ | |||
8379927 | ⋯ | |||
8394589 | ⋯ | |||
8424992 | ⋯ | |||
8694723 | ⋯ | |||
8760414 | ||||
8938364 | ||||
9025370 | ⋯ | |||
9098294 | ⋯ | |||
9139151 | ⋯ | |||
9139163 | ⋯ | |||
9206432 | ⋯ | |||
9353712 | ⋯ | |||
9410862 | ⋯ | |||
9414417 | ⋯ | |||
9812850 | ⋯ | |||
9955598 | ⋯ | |||
9965715 | ⋯ | |||
10068307 | ⋯ | |||
10079226 | ⋯ | |||
10162436 | ⋯ | |||
10454113 | ||||
10516096 | ⋯ | |||
10644253 | ⋯ | |||
10730618 | ⋯ | |||
10963065 | ⋯ | |||
11081729 | ⋯ | |||
11253226 | ⋯ | |||
11772920 | ⋯ | |||
12009504 | ⋯ | |||
12069127 | ⋯ | |||
12069424 | ||||
12069449 | ||||
12258514 | ⋯ | |||
12317678 | ⋯ |
Note. Total visibilities have been constructed from combining the MCMC chains of the individual visibilities. Uncertainties are obtained from the 68% HPD intervals of the posterior PDFs.
In comparing to the theoretically estimated values by Ballot et al. (2011), which are calculated considering the Kepler bandpass, we see that in median l = 1 visibilities are lower than expected, those for l = 2 and l = 3 are larger than expected, and the total visibilities (up to l = 3) are only slightly lower than expected. The difference is most pronounced for the l = 3 modes. It is important to note that because l = 3 visibilities have only been measured for a subset of stars, it might be that the remainder of the stars (with an S/N too low for a visual detection of l = 3 modes) have visibilities in agreement with theory. The comparisons do, however, qualitatively match those obtained by Mosser et al. (2012) for giants observed by Kepler; the l = 3 modes disagree most with theory in their results too. Similar discrepancies with theory have also been observed for individual CoRoT and Kepler targets analyzed by Deheuvels et al. (2010), Mathur et al. (2013), and Lund et al. (2014a), and the Sun by Salabert et al. (2011). Given these discrepancies, we discourage the adoption of fixed mode visibilities in peak-bagging exercises. To assess which other parameters in the peak-bagging contribute the most to the uncertainties on the visibilities, and thus which parameters will be most affected by adopting fixed visibilities, we show in Figure 22 with a box-plot the correlations between the visibilities and the remainder of the fitted parameters. Correlations have been estimated using Spearman's rank correlation (Spearman 1904) from the MCMC chains, because this better catches dependencies that are monotonic but not necessarily linear (as done in Pearson's correlation; Pearson 1895). As seen, the correlations are overall quite small with median values within ρ = ±0.2 in all cases. The visibilities are seen to primarily correlate with the parameters of the noise background, as might be expected with an increasing correlation with increasing angular degree l; a small change in, say, N has a relatively larger impact on the visibility of l = 3 modes compared to that of l = 1 modes and one should therefore expect a larger correlation. This suggests that fixing the visibilities might bias the fit of the background and vice versa for fixed backgrounds. One should, however, also be cautious with the extracted visibilities, because an inappropriate model for the background might bias the measured values.
Download figure:
Standard image High-resolution imageConsidering the total visibility, we note that Ballot et al. (2011) include modes up to l = 4, but this should not significantly affect any observed discrepancy. The median value obtained for is 3.07, which is close to the value of c = 3.04 (Bedding et al. 2010) often adopted in estimating radial mode amplitudes via the method of Kjeldsen et al. (2005, 2008a). In Section 4.2.1, we did indeed also find that the amplitudes from the smoothing method were equal in median from using c = 3.04 and , with only a slightly reduced scatter from using the measured total visibility. The value of c = 3.04 (or alternatively c = 3.07) should therefore serve as a reasonably good choice for analysis of amplitudes for a large sample of stars, but the small systematic offset found in Section 4.2.1 should be remembered when comparing with theory. We note also that Mosser et al. (2012) find a mean value of c = 3.06 from Kepler giants with between approximately 4000 and 5100 K, but with a larger scatter than for our values. This does, however, indicate that the visibilities do not increase with decreasing temperature as suggested by theory. If a trend exists with and/or , which should be the main parameters determining the visibilities, they cannot be clearly discerned from the observations.
Some of the discrepancies, and scatter in observed values, can likely be explained by some of the simplifications adopted in the calculations by Ballot et al. (2011) and in general. These include, for instance, the neglect of non-adiabatic effects and a height dependence on mode amplitudes in the stellar atmosphere (see, e.g., Baldner & Schou 2012; Schou 2014, 2015). Furthermore, phenomena such as spots and other local surface features will influence the measured visibility. An effect will also come from the way in which the stellar limb-darkening (LD) is described; the calculation of LD coefficients (LDCs) for parametrized LD laws will, e.g., depend on the description of convection in the adopted atmosphere models, the method used for integrating the specific flux, and the resolution used when fitting a parametrized LD law to these flux values. Ballot et al. (2011) compared their visibilities with those obtained using LDC from Sing (2010; specific to the Kepler bandpass) and found that their values were generally higher. Similarly, Claret & Bloemen (2011) compared their Kepler LDCs to those of Sing (2010) and found differences. This indicates that some systematic uncertainty should be added to the theoretically derived visibilities. In any case, Sing (2010) showed from fits to planetary transit curves that model computations of LDC generally disagree with those derived empirically. The treatment of the LD is, however, likely a secondary effect—as shown in Lund et al. (2014a) the specific LD law adopted (neglecting the linear one) only has minor effects on the visibilities, and the shape has to be changed by a large amount away from the limb in order to take effect. This, and the fact that discrepancies between measured and modeled visibilities are found for the Sun (Salabert et al. 2011), where the LD is well known as a function of wavelength (Neckel & Labs 1994), indicates that the simplified assumptions concerning the mode physics likely are the main contributor to the discrepancies. We note, however, that modes of l = 3 would be relatively more affected by details of the LD, because of the considerably stronger cancellation for l = 3 modes (total in the absence of LD), compared to modes of l = 1, 2 (Lund et al. 2014a). It could also be questioned if the assumption of equipartition of energy between modes of different angular degree holds true.
4.4.1. Detection of Modes
For solar-like oscillators observed by Kepler the highest angular degree of modes is typically l = 2. Only for the highest S/N cases can higher degree modes be identified—in the current sample, octupole l = 3 modes were identified and included in the peak-bagging in 14 such targets. It is, however, possible to obtain information on the combined signal for l = 3 modes in a given star. Such a signal could, for instance, be used in estimating the small separation , which can contribute a constraint in modeling efforts (see, e.g., Bellinger et al. 2016).
To optimize the combined signal from modes of a given angular degree, we applied the method outlined in Lund et al. (2014a). Briefly, the échelle diagram is collapsed along the vertical direction after having first stretched the frequency scale such that modes of a given l form a straight ridge—this ensures that all mode power from the particular l-values is co-added with as little spread as possible. Individual frequencies from best-fitting ASTFIT models were used to stretch the frequency scale (see Paper II). Before the power is co-added, the power spectrum is smoothed to account also for the mode line widths and a potential rotational splitting. In Figure 23, we show in a contour plot the collapsed échelle diagrams for all stars, optimized to increase the signal of l = 3 modes. To further increase the visibility of these modes, we first divided out the model of the power spectrum from the peak-bagging, including also fitted l = 3 modes. We see that for most stars a clearly detectable signal from l = 3 is present. It is beyond the scope of this paper to estimate for the stars and test the statistical significance of the signals, but Figure 23 should indicate that obtaining this information would be possible. Some of the signal seen around l < 3 modes in the residual spectra is likely due contamination from higher degree modes (l ≥ 5) or deviations in the mode shape from the assumed simple Lorentzian—asymmetry of the modes would, for instance, leave some unaccounted for power in the residual spectrum.
Download figure:
Standard image High-resolution imageFor the highest S/N targets, one may further look for indications of hexadecapole l = 4 modes. In Figure 24, we show the collapsed échelle diagrams optimized for the detection of l = 4 for three targets with the strongest apparent l = 4 signals, namely, KIC 6603624 (Saxo), KIC 7510397 (middle), and KIC 12069449 (16 Cyg B). As seen, an excess in collapsed power in each of these cases falls close to the expected position at from the ASTFIT model. The signal seen in 16 Cyg A (not shown) and B corresponds well to the signal found in Lund et al. (2014a) from a shorter Kepler data set. Curiously, the peak around ∼0.85 seen in Saxo and 16 Cyg B appears to coincide with the expected position of dotriacontapole l = 5 modes—whether this signal truly is from l = 5 modes requires further investigation. In any case, it is clear that for a star like KIC 7510397 the l = 0,2 modes will be polluted by l = 5 modes; similarly will l = 6,9 modes pollute the l = 1 signal, etc. (Appourchaux & Virgo Team 1998; Lund et al. 2014a).
Download figure:
Standard image High-resolution image5. Example Output
For each of the 66 stars in the sample, we have provided a set of outputs from the peak-bagging. Tables and plots for all stars will be available in the online appendices. The extracted parameters may in addition be obtained from the results section of the KASOC database.21
The outputs include, first of all, a table with the mode information from the MCMC peak-bagging, like the one given in Table 6 for KIC 6225718 (Saxo2). The table gives for a given mode the angular degree, radial order, frequency, amplitude, line width, and the natural logarithm of the Bayes factor K from the quality control in Section 3.2.3. The uncertainties on the mode parameters are obtained from the 68% HPD interval of the posterior probability distributions. As amplitudes and line widths are only fitted to radial l = 0 modes, we only give these values for these modes. The visibilities for a given star can be found in Table 5. We note that the radial orders given are obtained from matching to the expected value as a function of (see Section 3.2.1)—we suggest that this be used with some caution in modeling efforts and checked independently. A table is also included with the derived frequency difference ratios and second differences (see Section 3.3), examples of these are given in Tables 7 and 8.
Table 6. Extracted Mode Parameters and Quality Control (Equation (3.2.3)) for KIC 6225718 (Saxo2).
n | l | Line width | ln K | ||
---|---|---|---|---|---|
11 | 1 | 1351.15 | ⋯ | ⋯ | 2.2 |
12 | 0 | 1407.23 | 0.79 | 2.50 | >6 |
12 | 1 | 1454.25 | ⋯ | ⋯ | >6 |
13 | 0 | 1510.10 | 0.99 | 2.50 | >6 |
13 | 1 | 1558.45 | ⋯ | ⋯ | >6 |
13 | 2 | 1605.68 | ⋯ | ⋯ | 1.45 |
14 | 0 | 1615.12 | 1.16 | 2.60 | >6 |
14 | 1 | 1664.09 | ⋯ | ⋯ | >6 |
14 | 2 | 1711.40 | ⋯ | ⋯ | 3.14 |
15 | 0 | 1720.35 | 1.46 | 2.30 | >6 |
15 | 1 | 1769.65 | ⋯ | ⋯ | >6 |
15 | 2 | 1816.19 | ⋯ | ⋯ | >6 |
16 | 0 | 1825.41 | 1.89 | 2.19 | >6 |
16 | 1 | 1873.88 | ⋯ | ⋯ | >6 |
16 | 2 | 1919.97 | ⋯ | ⋯ | >6 |
17 | 0 | 1929.05 | 2.29 | 2.84 | >6 |
17 | 1 | 1977.35 | ⋯ | ⋯ | >6 |
17 | 2 | 2023.80 | ⋯ | ⋯ | >6 |
18 | 0 | 2032.68 | 2.77 | 2.67 | >6 |
18 | 1 | 2081.57 | ⋯ | ⋯ | >6 |
18 | 2 | 2128.62 | ⋯ | ⋯ | >6 |
19 | 0 | 2137.59 | 3.14 | 2.50 | >6 |
19 | 1 | 2186.89 | ⋯ | ⋯ | >6 |
19 | 2 | 2234.70 | ⋯ | ⋯ | >6 |
20 | 0 | 2243.42 | 3.47 | 2.22 | >6 |
19 | 3 | 2281.61 | ⋯ | ⋯ | 3.01 |
20 | 1 | 2293.05 | ⋯ | ⋯ | >6 |
20 | 2 | 2340.63 | ⋯ | ⋯ | >6 |
21 | 0 | 2349.64 | 3.46 | 2.61 | >6 |
20 | 3 | 2385.57 | ⋯ | ⋯ | 3.94 |
21 | 1 | 2399.39 | ⋯ | ⋯ | >6 |
21 | 2 | 2446.71 | ⋯ | ⋯ | >6 |
22 | 0 | 2455.69 | 3.46 | 3.03 | >6 |
21 | 3 | 2493.08 | ⋯ | ⋯ | 3.66 |
22 | 1 | 2505.34 | ⋯ | ⋯ | >6 |
22 | 2 | 2552.85 | ⋯ | ⋯ | >6 |
23 | 0 | 2561.29 | 3.18 | 4.01 | >6 |
22 | 3 | 2598.55 | ⋯ | ⋯ | 3.12 |
23 | 1 | 2611.20 | ⋯ | ⋯ | >6 |
23 | 2 | 2658.63 | ⋯ | ⋯ | >6 |
24 | 0 | 2666.49 | 2.72 | 5.23 | >6 |
24 | 1 | 2717.47 | ⋯ | ⋯ | >6 |
24 | 2 | 2765.05 | ⋯ | ⋯ | >6 |
25 | 0 | 2773.06 | 2.41 | 6.75 | >6 |
25 | 1 | 2824.15 | ⋯ | ⋯ | >6 |
25 | 2 | 2872.28 | ⋯ | ⋯ | >6 |
26 | 0 | 2879.34 | 1.96 | 7.60 | >6 |
26 | 1 | 2931.24 | ⋯ | ⋯ | >6 |
26 | 2 | 2978.49 | ⋯ | ⋯ | 3.78 |
27 | 0 | 2987.15 | 1.70 | 7.53 | >6 |
27 | 1 | 3038.67 | ⋯ | ⋯ | >6 |
27 | 2 | 3084.55 | ⋯ | ⋯ | 1.45 |
28 | 0 | 3092.80 | 1.46 | 8.85 | >6 |
28 | 1 | 3145.65 | ⋯ | ⋯ | 4.92 |
28 | 2 | 3194.64 | ⋯ | ⋯ | 2.18 |
29 | 0 | 3204.41 | 0.98 | 5.87 | 3.81 |
29 | 1 | 3251.96 | ⋯ | ⋯ | 4.2 |
29 | 2 | 3302.59 | ⋯ | ⋯ | 1.22 |
30 | 0 | 3314.17 | 1.12 | 11.64 | 3.81 |
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
Table 7. Example of Calculated Mode Frequency Difference Ratios (Equation (16)) for KIC 6225718 (Saxo2)
n | ||
---|---|---|
r02 | 14 | |
r02 | 15 | |
r02 | 16 | |
r02 | 17 | |
r02 | 18 | |
r02 | 19 | |
r02 | 20 | |
r02 | 21 | |
r02 | 22 | |
r02 | 23 | |
r02 | 24 | |
r02 | 25 | |
r02 | 26 | |
r02 | 27 | |
r02 | 28 | |
r02 | 29 | |
r10 | 12 | |
r01 | 13 | |
r10 | 13 | |
r01 | 14 | |
r10 | 14 | |
r01 | 15 | |
r10 | 15 | |
r01 | 16 | |
r10 | 16 | |
r01 | 17 | |
r10 | 17 | |
r01 | 18 | |
r10 | 18 | |
r01 | 19 | |
r10 | 19 | |
r01 | 20 | |
r10 | 20 | |
r01 | 21 | |
r10 | 21 | |
r01 | 22 | |
r10 | 22 | |
r01 | 23 | |
r10 | 23 | |
r01 | 24 | |
r10 | 24 | |
r01 | 25 | |
r10 | 25 | |
r01 | 26 | |
r10 | 26 | |
r01 | 27 | |
r10 | 27 | |
r01 | 28 | |
r10 | 28 | |
r01 | 29 |
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
Table 8. Example of Calculated Second Differences (Equation (20)) for KIC 6225718 (Saxo2)
n | l | |
---|---|---|
12 | 1 | 1.49 |
13 | 0 | 1.64 |
13 | 1 | 1.20 |
14 | 0 | 0.37 |
14 | 1 | 0.13 |
14 | 2 | −1.08 |
15 | 0 | −0.23 |
15 | 1 | −1.27 |
15 | 2 | −0.96 |
16 | 1 | −0.76 |
16 | 2 | 0.14 |
16 | 0 | −1.35 |
17 | 0 | 0.02 |
17 | 2 | 0.97 |
17 | 1 | 0.69 |
18 | 2 | 1.25 |
18 | 1 | 1.13 |
18 | 0 | 1.30 |
19 | 2 | −0.15 |
19 | 1 | 0.82 |
19 | 0 | 0.89 |
20 | 2 | 0.13 |
20 | 3 | −1.08 |
20 | 1 | 0.15 |
20 | 0 | 0.42 |
21 | 3 | −0.96 |
21 | 1 | −0.37 |
21 | 2 | 0.16 |
21 | 0 | −0.16 |
22 | 1 | −0.09 |
22 | 2 | −0.42 |
22 | 0 | −0.45 |
23 | 2 | 0.43 |
23 | 0 | −0.44 |
23 | 1 | 0.38 |
24 | 0 | 1.38 |
24 | 2 | 0.72 |
24 | 1 | 0.40 |
25 | 0 | −0.20 |
25 | 2 | −1.27 |
25 | 1 | 0.47 |
26 | 2 | −0.15 |
26 | 0 | 1.51 |
26 | 1 | 0.36 |
27 | 0 | −1.98 |
27 | 1 | −0.35 |
27 | 2 | 4.80 |
28 | 2 | −4.34 |
28 | 1 | −0.51 |
28 | 0 | 6.76 |
29 | 0 | −1.46 |
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
For each star, a number of plots are also prepared. These include (1) a visualization of the obtained "best fit" from the peak-bagging as in Figure 25, given by a plot of the power spectrum overlain with the best-fit model and with an indication of the extracted modes; (2) an échelle diagram overlain with the extracted frequencies as in the left panel of Figure 26; (3) a plot of the derived frequency difference ratios as shown in the right panel of Figure 26; (4) a plot of the derived second differences as shown in Figure 5 for KIC 6225718 (Saxo2), but without model fits; (5) a plot of the extracted mode amplitudes and line widths as shown in Figure 27.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution image6. Conclusions
In this paper, we have presented the mode parameters for a sample of 66 MS solar-like oscillators, the frequencies of which are modeled in Paper II. In addition to the individual mode frequencies, we have constructed frequency difference ratios and their correlations for the use in modeling efforts. We also report for each star the values for the mode line widths, amplitudes, and visibilities. For each of these quantities, we have derived summary parameters and descriptions, such as the average seismic parameters and and the behavior of mode line widths against frequency.
The reported parameters were derived through peak-bagging of the power spectra using an MCMC optimization scheme. This resulted in posterior probability distributions for each of the included parameters, from which credible intervals and correlations were directly obtained. We found that the derived frequency uncertainties adhere to the expectations from theory in terms of S/N, observing time, and line widths, and are only higher by a factor of ∼1.2 compared to estimates from MLE. This observation will be useful for predicting seismic modeling yields for future missions such as TESS and PLATO. As a quality control on the detection of the reported modes, we performed a Bayesian hypothesis testing that for each mode gave the probability of detection (Appourchaux et al. 2012b; Davies et al. 2016). Our main conclusions are as follows.
- 1.The derived values for and agree with empirically derived relations from the Kepler mission. We derived parameters from an extended version of the standard asymptotic frequency relations, including mode ridge curvatures and variations of and small frequency separations. From the small frequency separations, we further found that most stars deviate from the asymptotic description by an amount that correlates with the evolutionary state of the star, i.e., the central hydrogen content.
- 2.The measured amplitudes at for our sample largely follow the expected trend from empirical relations from Kepler. We also identified a systematic offset of approximately −6% between the maximum amplitudes obtained from the modes and those obtained from the smoothing method by Kjeldsen et al. (2008a). This corroborates the findings by Verner et al. (2011). This systematic offset should be corrected for whenever the two methods are compared.
- 3.For the line widths, we adopted the frequency dependence of Appourchaux et al. (2014) given by an overall power-law dependence and a Lorentzian dip near We fitted this relation for all stars and were able to derive simple relations between the parameters of the fit and and These were found to confirm the results by Appourchaux et al. (2014). Such relations will be useful for future simulations of solar-like oscillators and may be compared to theory. We also obtained a fit for the dependence of the line width at complementing our values with line widths from 42 giants in NGC 6819 (R. Handberg et al. 2017, in preparation). The obtained dependence largely agreed with that found by Appourchaux et al. (2012a), except for stars with low .
- 4.Concerning the estimated mode visibilities, we found that those for l = 1 were slightly lower than expected from the theoretical calculations by Ballot et al. (2011), whereas those for l = 2 and 3 were larger than expected, especially for the l = 3 modes. We found no overall dependence on , which is also evident from the fact that Mosser et al. (2012) found a mean value of the total visibility at nearly the same level as here for Kepler giants, which have values lower by about 1000–2000 K. Some structure was observed in the visibilities against and mass, but it was not possible to say directly if this is simply due to scatter from the measurements or if they have some underlying physical explanation. Applying the method of Lund et al. (2014a) enabled us to identify power from l = 3 modes in most stars, and for some high S/N targets even for l = 4.
Funding for this Discovery mission is provided by NASA's Science Mission Directorate. The authors acknowledge the dedicated team behind the Kepler and K2 missions, without whom this work would not have been possible. We thank Thierry Appourchaux for useful discussions on the properties of the line width fits, Ian Roxburgh for discussions on mode frequencies, and Enrico Corsaro for supplying line widths for a set of 19 red giants. We thank Jens Jessen-Hansen for inspiration to the paper title. Funding for the Stellar Astrophysics Centre (SAC) is provided by The Danish National Research Foundation (Grant DNRF106). The research was supported by the ASTERISK project (ASTERoseismic Investigations with SONG and Kepler) funded by the European Research Council (Grant agreement no.: 267864). M.N.L. acknowledges the support of The Danish Council for Independent Research Natural Science (Grant DFF-4181-00415). V.S.A. and T.R.W. acknowledges support from VILLUM FONDEN (research grant 10118). G.R.D. and W.J.C. acknowledge the support of the UK Science and Technology Facilities Council (STFC). S.B. is partially supported by NSF grant AST-1514676 and NASA grant NNX16A109G. D.H. acknowledges support by the Australian Research Council's Discovery Projects funding scheme (project number DE140101364) and support by the National Aeronautics and Space Administration under Grant NNX14AB92G issued through the Kepler Participating Scientist Program. W.H.B. acknowledges research funding by Deutsche Forschungsgemeinschaft (DFG) under grant SFB 963/1 "Astrophysical flow instabilities and turbulence" (Project A18).
Footnotes
- 13
First coined by R. Ulrich circa (1983, private communication) from the analogy to hill climbing, where it refers to reaching the summits of a collection of peaks, the term was later re-introduced by Appourchaux (2003b).
- 14
- 15
- 16
Variability of Solar Irradiance and Gravity Oscillations.
- 17
Solar and Heliospheric Observatory.
- 18
In Kepler, Δt equals the sampling time wherefore x sometimes is given as , where is the Nyquist frequency—this is, however, an imprecise definition, because relates to the sampling time whereas the apodization relates to the integration time.
- 19
Expressed in the units and .
- 20
- 21