A publishing partnership

Standing on the Shoulders of Dwarfs: the Kepler Asteroseismic LEGACY Sample. I. Oscillation Mode Parameters

, , , , , , , , , , , , , , , , , , , and

Published 2017 January 30 © 2017. The American Astronomical Society. All rights reserved.
, , Citation Mikkel N. Lund et al 2017 ApJ 835 172 DOI 10.3847/1538-4357/835/2/172

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

This article is corrected by 2017 ApJ 850 110

0004-637X/835/2/172

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 ${\rm{\Delta }}\nu $ and $\delta {\nu }_{02};$ the behavior of line widths with frequency and line widths at ${\nu }_{\max }$ with ${T}_{\mathrm{eff}}$, 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; ${\rm{\Delta }}t=58.89\,{\rm{s}}$) 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 (${T}_{\mathrm{eff}}$ versus $\mathrm{log}g$ ), 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.

Figure 1. Refer to the following caption and surrounding text.

Figure 1. Kiel-diagram of sample stars, with ${T}_{\mathrm{eff}}$ and $[\mathrm{Fe}/{\rm{H}}]$ from the spectroscopic input (Table 1) and $\mathrm{log}g$ from the modeling in Paper II. Stellar evolutionary tracks have been calculated using the Garching Stellar Evolution Code (GARSTEC; Weiss & Schlattl 2008) adopting [Fe/H] = 0. The marker type indicates whether the star is considered to be of "Simple" or "F-type" character according to Appourchaux et al. (2012b). The dotted lines show lines of constant ${\nu }_{\max },$ according to the simple scaling as ${\nu }_{\max }\propto g\sqrt{{T}_{\mathrm{eff}}\,}$, in steps of 250 μHz.

Standard image High-resolution image

Table 1.  Parameters for the Targets in the Studied Sample

KIC Popular Kp ${\nu }_{\max }$ ${\rm{\Delta }}\nu $ Number of Category Braketing Missing ${T}_{\mathrm{eff}}$ [Fe/H] $\mathrm{log}g$ LOS $v\sin i$
  Name (mag) (μHz) (μHz) Modes   Quarters Quarters (K) (dex) (cgs; dex) ($\mathrm{km}\,{{\rm{s}}}^{-1}$) ($\mathrm{km}\,{{\rm{s}}}^{-1}$)
1435467 8.88 1407 70.4 46 F-like 5.1–17.2 6326 ± 77 0.01 ± 0.10 4.100${}_{-0.009}^{+0.009}$ −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${}_{-0.007}^{+0.007}$ −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${}_{-0.005}^{+0.004}$ −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${}_{-0.007}^{+0.005}$ −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${}_{-0.007}^{+0.005}$ −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${}_{-0.008}^{+0.010}$ −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${}_{-0.007}^{+0.007}$ 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${}_{-0.008}^{+0.010}$ −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${}_{-0.010}^{+0.008}$ 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${}_{-0.008}^{+0.007}$
5950854 10.96 1927 96.6 26 Simple 5.1–10.3 6, 7.2 5853 ± 77 −0.23 ± 0.10 4.238${}_{-0.007}^{+0.007}$ −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${}_{-0.009}^{+0.009}$ −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${}_{-0.009}^{+0.012}$ −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${}_{-0.005}^{+0.007}$ −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${}_{-0.005}^{+0.007}$ 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${}_{-0.004}^{+0.005}$ −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${}_{-0.007}^{+0.008}$ −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${}_{-0.008}^{+0.009}$ −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${}_{-0.007}^{+0.007}$ −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${}_{-0.010}^{+0.008}$
7206837 Bagheera 9.77 1653 79.1 45 F-like 5.1–17.2 6305 ± 77 0.10 ± 0.10 4.163${}_{-0.008}^{+0.007}$ −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${}_{-0.009}^{+0.010}$ 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${}_{-0.007}^{+0.004}$ −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${}_{-0.008}^{+0.010}$ −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.007}^{+0.007}$ −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${}_{-0.007}^{+0.005}$ −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${}_{-0.002}^{+0.002}$ −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${}_{-0.004}^{+0.005}$ −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${}_{-0.007}^{+0.007}$ −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${}_{-0.008}^{+0.008}$
8179536 9.46 2075 95.1 39 Simple 5.1–11.3 6 6343 ± 77 −0.03 ± 0.10 4.255${}_{-0.010}^{+0.010}$ −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${}_{-0.005}^{+0.004}$ 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) ${4.388}_{-0.008}^{+0.007}$
8394589 9.52 2397 109.5 44 Simple 5.1–17.2 6143 ± 77 −0.29 ± 0.10 4.322${}_{-0.008}^{+0.008}$ 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${}_{-0.007}^{+0.007}$ −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${}_{-0.007}^{+0.009}$ 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${}_{-0.007}^{+0.003}$ −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${}_{-0.007}^{+0.002}$ −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${}_{-0.007}^{+0.004}$
9098294 9.76 2315 108.9 34 Simple 5.1–17.2 5852 ± 77 −0.18 ± 0.10 4.308${}_{-0.005}^{+0.007}$ −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${}_{-0.008}^{+0.008}$ −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${}_{-0.009}^{+0.008}$
9206432 9.08 1866 84.9 49 F-like 5.1–12.3 7 6538 ± 77 0.16 ± 0.10 4.220${}_{-0.005}^{+0.007}$ −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${}_{-0.007}^{+0.005}$ −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${}_{-0.009}^{+0.008}$ −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${}_{-0.005}^{+0.005}$
9812850 9.47 1255 64.7 49 F-like 5.1–17.2 6321 ± 77 −0.07 ± 0.10 4.053${}_{-0.008}^{+0.009}$ 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${}_{-0.007}^{+0.005}$ −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${}_{-0.008}^{+0.009}$
10068307 8.18 995 53.9 49 Simple 7.1–17.2 6132 ± 77 −0.23 ± 0.10 3.967${}_{-0.004}^{+0.004}$ −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${}_{-0.005}^{+0.005}$ −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${}_{-0.005}^{+0.005}$ −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${}_{-0.005}^{+0.005}$ −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${}_{-0.010}^{+0.011}$ 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${}_{-0.007}^{+0.008}$ −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${}_{-0.008}^{+0.007}$
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${}_{-0.011}^{+0.011}$ −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.010}^{+0.009}$ 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${}_{-0.005}^{+0.004}$ 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${}_{-0.005}^{+0.008}$
12009504 Dushera 9.32 1866 88.2 43 Simple 5.1–17.2 6179 ± 77 −0.08 ± 0.10 4.211${}_{-0.007}^{+0.005}$ 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${}_{-0.005}^{+0.004}$ −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${}_{-0.007}^{+0.007}$ −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${}_{-0.007}^{+0.005}$ −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${}_{-0.004}^{+0.003}$ −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${}_{-0.008}^{+0.009}$ −58.14 ± 0.10 8.40 ± 0.50

Note: ${T}_{\mathrm{eff}}$ 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 $v\sin i$ values are adopted from the SPC analysis of the targets unless otherwise indicated (see table notes). Systematic uncertainties of ±59 K (${T}_{\mathrm{eff}}$) and ±0.062 dex ([Fe/H]) have been added in quadrature as suggested by Torres et al. (2012). Values for $\mathrm{log}g$ 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 $({\nu }_{\max })$ and large separation (${\rm{\Delta }}\nu $)—for uncertainties in ${\nu }_{\max }$ and ${\rm{\Delta }}\nu $ 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.

Download table as:  ASCIITypeset images: 1 2

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, $\mathrm{log}g$ values were fixed to the asteroseismic values given in Chaplin et al. (2014) to decrease the impact on uncertainties from correlations between ${T}_{\mathrm{eff}}$, log g, and $[\mathrm{Fe}/{\rm{H}}]$ . We added in quadrature to the derived uncertainties on ${T}_{\mathrm{eff}}$ and $[\mathrm{Fe}/{\rm{H}}]$ 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 $-0.61\,\mathrm{km}\,{{\rm{s}}}^{-1}$ 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.

Figure 2. Refer to the following caption and surrounding text.

Figure 2. Frequency shift at ${\nu }_{\max }$ from line-of-sight velocities of the stars in the sample (see Table 1). The uncertainty indicates for a given target the minimum frequency uncertainty of the five radial modes nearest ${\nu }_{\max }.$ The frequency uncertainty on individual modes is in several cases lower than the line-of-sight frequency shift. Large or small, the systematic shift can thus cause bias in the modeling if left uncorrected.

Standard image High-resolution image

2.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 $\mathrm{Kp}\,\approx 9.17$ 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 ${\rm{\Delta }}\nu $ and ${\nu }_{\max },$ 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):

Equation (1)

Each mode is characterized by the frequency ${\nu }_{{nl}}$ of the zonal (m = 0) component, a height ${H}_{{nlm}}={{ \mathcal E }}_{{lm}}({i}_{\star }){\tilde{V}}_{l}^{2}{S}_{n0}$, an FWHM mode width ${{\rm{\Gamma }}}_{{nl}}$, and a rotational splitting ${\nu }_{s}$ (assumed constant with frequency; see Lund et al. 2014c, for a discussion of the impact of differential rotation on the constancy of ${\nu }_{s}$). In ${H}_{{nlm}}$, ${{ \mathcal E }}_{{lm}}({i}_{\star })$ 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); ${\tilde{V}}_{l}^{2}$ denotes the squared visibility (power units) of a non-radial mode relative to a radial mode at the same frequency, i.e.,

Equation (2)

from the spatial filtering resulting from integrating the intensity for a mode of a given degree over the stellar surface; ${S}_{n0}$ then denotes the height of the radial mode of the order of n. The use of ${\tilde{V}}_{l}^{2}$ 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

Equation (3)

Here $N(\nu )$ denotes the adopted background function; W is a constant white shot-noise component; ${\eta }^{2}(\nu )$ 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

Equation (4)

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)

Equation (5)

which characterizes a temporal signal from granulation having an exponentially decaying autocovariance, with a power of the temporal decay rate as $-2/{\alpha }_{i};$ ${\tau }_{i}$ gives the characteristic timescale of the ith background component; ${\sigma }_{i}$ the corresponding root-mean-square (rms) variation of the component in the time domain. The normalization constants ${\xi }_{i}$ are such that the integral (for positive frequencies) of the background component equals ${\sigma }_{i}^{2}$, 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 ${{\rm{\Gamma }}}_{{nl}}$ (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):

Equation (6)

This is a valid approximation for a single-sided power spectrum when the modes are well resolved, i.e., when the observing duration ${T}_{\mathrm{obs}}$ greatly exceeds the mode lifetime $2/\pi {\rm{\Gamma }}$. We note that Anl and ${{\rm{\Gamma }}}_{{nl}}$ 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 ${\boldsymbol{\Theta }}=\{{\nu }_{s},{i}_{\star },W,{\tau }_{i},{\sigma }_{i},{\alpha }_{i},{A}_{n0},{{\rm{\Gamma }}}_{n0},{\nu }_{{nl}},{\tilde{V}}_{l}\}$.

3.2. Fitting Strategy

Parameters were estimated in a Bayesian manner from a global peak-bagging fit to the power spectrum including all parameters ${\boldsymbol{\Theta }}$ (see, e.g., Handberg & Campante 2011). This was done by mapping the posterior probability of the parameters ${\boldsymbol{\Theta }}$ given the data D and any prior information I, which from Bayes' theorem is given as

Equation (7)

Here $p({\boldsymbol{\Theta }}| I)$ is the prior probability assigned to the parameters ${\boldsymbol{\Theta }}$ given I, and $p(D| {\boldsymbol{\Theta }},I)$ is the likelihood of the observed data D given the parameters ${\boldsymbol{\Theta }}$. $p(D| I)$, 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:

Equation (8)

where logarithmic units are adopted for numerical stability, and C is a constant. Assuming a ${\chi }^{2}$ 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, ${ \mathcal P }({\nu }_{j};{\boldsymbol{\Theta }})$, as (see Duvall & Harvey 1986; Anderson et al. 1990; Toutain & Appourchaux 1994):

Equation (9)

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):

Equation (10)

Here ${\rm{\Delta }}\nu $ is the large separation, given by the average frequency spacing between consecutive overtones n for modes of a given l; epsilon 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 epsilon as a function of ${T}_{\mathrm{eff}}$ (Figure 3), where epsilon 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).

Figure 3. Refer to the following caption and surrounding text.

Figure 3. Measured values of epsilon against ${T}_{\mathrm{eff}}$. The color indicates the modeled mass of the stars using the results from the BASTA pipeline (Paper II). Shown are also epsilon-evolutionary tracks from White et al. (2011a), calculated from ASTEC evolutionary tracks with Z0 = 0.017 (Christensen-Dalsgaard 2008), for masses going from $0.8\,{M}_{\odot }$ to $1.6\,{M}_{\odot }$ in steps of $0.1\,{M}_{\odot }$.

Standard image High-resolution image
Figure 4. Refer to the following caption and surrounding text.

Figure 4. Échelle diagrams for three sample stars with different values of ${\nu }_{\max };$ KIC 1435467 (left), KIC 6106415 (middle), and KIC 7970740 (right). The color scale goes from white (low power) to black (high power), and the power spectra have been background corrected and smoothed by a 2 μHz Epanechnikov filter (Epanechnikov 1969; Hastie et al. 2009). The échelle spectra have been shifted along the abscissa for a better rendering (see individual labels), with the l = 1 ridges being the right-most in each case and the l = 0, 2 ridges the left-most.

Standard image High-resolution image

For 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 ${f}_{\min }-5{\rm{\Delta }}\nu $ to ${f}_{\max }+5{\rm{\Delta }}\nu $, where ${f}_{\min }$ and ${f}_{\max }$ 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 ${\nu }_{\mathrm{nq}}$ (∼8496 μHz). In this fit, the power from solar-like oscillations was accounted for by a Gaussian envelope centered at ${\nu }_{\max }.$ 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 ${\nu }_{s}$ 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 ${i}_{\star }=0^\circ $ or ${i}_{\star }=90^\circ $ 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)

Equation (11)

which behaves as a uniform prior when $\theta \ll {\theta }_{\mathrm{uni}}$ and a standard scale invariant Jeffrey's prior when $\theta \gg {\theta }_{\mathrm{uni}}$. The maximum of the prior occurs at ${\theta }_{\max }.$

For mode visibilities, we adopted truncated Gaussian functions ${ \mathcal N }({\theta }_{0},\sigma ,{\theta }_{\min },{\theta }_{\max })$ as priors, defined as

Equation (12)

with D given as:

Equation (13)

Here, erf denotes the error function, ${\theta }_{0}$ and σ give the chosen mode value and width of the Gaussian, and ${\theta }_{\min }$ and ${\theta }_{\max }$ give the lower and upper truncations of the Gaussian. We specifically adopted ${ \mathcal N }(1.5,1.5,0,3)$ for l = 1, ${ \mathcal N }(0.5,0.5,0,1)$ for l = 2, and ${ \mathcal N }(0.05,0.05,0,0.5)$ for l = 3. Similarly, we adopted truncated Gaussian priors for the parameters of the background. Here we used as the Gaussian mode value (${\theta }_{0}$) the median of the posteriors from the background-only fit (see Section 3.2.1), and we adopted $\sigma =0.1\,{\theta }_{0}$, ${\theta }_{\min }=0.1\,{\theta }_{0}$, and ${\theta }_{\max }=10\,{\theta }_{0}$.

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 $\delta {\nu }_{02}={\nu }_{n,0}-{\nu }_{n-\mathrm{1,2}}$ is small compared to the mode line width, we added the prior constraint that on $\delta {\nu }_{02}$ that it must be positive. In principle, $\delta {\nu }_{02}$ 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 $\delta {\nu }_{02}\lt 0$ 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 $p({\mathrm{Det}}_{n,l}| D)$ conditioned on the data D needed to be explicitly determined and evaluated. This was done because the explicit determination of $p({\mathrm{Det}}_{n,l}| D)$ 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 $p(D| {H}_{0})=0.001$ level (Appourchaux 2003a, 2004; Lund et al. 2012). When the high S/N modes had been identified in this manner the probability of detection $p({\mathrm{Det}}_{n,l}| D)$ was computed for the remainder low-S/N modes.

In the computation of $p({\mathrm{Det}}_{n,l}| D)$, both the probability of D assuming H0, $p(D| {H}_{0})$, and the probability of the alternative hypothesis H1 of a detected mode, $p(D| {H}_{1})$, 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 ${\boldsymbol{\theta }}$—this integration was achieved by marginalizing over $p(D| {H}_{1},{\boldsymbol{\theta }})$, with the parameter space sampled over the posteriors from the peak-bagging using emcee. Specifically, a mixture-model was used in which both $p(D| {H}_{0})$ and $p(D| {H}_{1})$ were optimized simultaneously to give $p(D| {\boldsymbol{\theta }},{p}_{a})$, the probability of observing the data given the model of a given set of modes with parameters ${\boldsymbol{\theta }}$:

Equation (14)

Here, the parameter pa, ranging between 0 and 1, then gives the probability of the detection $p(D| {\mathrm{Det}}_{l})$ 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 $p(D| {\mathrm{Det}}_{l})$ over $p(D| {H}_{0})$, as

Equation (15)

The value of $\mathrm{ln}K$ 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)

Equation (16)

Here, d01 and d10 are the smooth five-point small frequency separations defined as

Equation (17)

Equation (18)

and the large separation is

Equation (19)

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:

Equation (20)

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.

Figure 5. Refer to the following caption and surrounding text.

Figure 5.

Top: example of calculated second differences ${{\rm{\Delta }}}_{2}\nu (n,l)$ from Equation (20), for the star KIC 6225718 (Saxo2). The frequencies used to calculate ${{\rm{\Delta }}}_{2}\nu (n,l)$ are given in Table 6, and frequency uncertainties were taken as the average of the asymmetric uncertainties. A clear oscillation is seen in ${{\rm{\Delta }}}_{2}\nu $, indicating acoustic glitches. The full line gives the acoustic glitch fit to ${{\rm{\Delta }}}_{2}\nu (n,l)$ by G. Houdek et al. (2017, in preparation). Only the values of ${{\rm{\Delta }}}_{2}\nu $ used in the fitting by G. Houdek et al. (2017, in preparation) are shown. Bottom: individual components from the acoustic glitches to ${{\rm{\Delta }}}_{2}\nu (n,l)$ (full line in top panel), showing the contributions from the first ($\mathrm{HeI}$) and second ($\mathrm{HeII}$) stages of helium ionization, and the base of the convective zone ($\mathrm{BCZ}$). (The complete figure set (67 images) is available.)

Standard image High-resolution image

In 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 ${r}_{\mathrm{MAD}}$ (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 ${r}_{\mathrm{MAD}}$ 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 ${r}_{\mathrm{MAD}}$ between two parameters x and y was calculated as follows.

Equation (21)

where u and v are the robust principle variables for x and y:

Equation (22)

Equation (23)

Figure 6 gives an example of the correlation matrix between the ${r}_{\mathrm{01,10}}$ frequency ratios for a given star in the form of a Hinton diagram (Hinton et al. 1986).

Figure 6. Refer to the following caption and surrounding text.

Figure 6. Hinton diagram showing the correlations between frequency difference ratios ${r}_{\mathrm{01,10}}(n)$ of different radial order n for the star KIC 6106415 (Perky). The size and color are proportional to the correlation ${r}_{\mathrm{MAD}}$.

Standard image High-resolution image

4. 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)

Equation (24)

where $\beta =B/H$ (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 ${x}_{s}=2{\nu }_{s}/{\rm{\Gamma }};$ and ${T}_{\mathrm{obs}}$ is the observing duration. For radial modes (l = 0), the factor fl only depends on β and is given by (Libbrecht 1992)

Equation (25)

For the general version of ${f}_{l}(\beta ,{x}_{s},{i}_{\star })$, see Ballot et al. (2008).

In Figure 7, we show the individual frequency uncertainties obtained from the peak-bagging as a function of ${\beta }^{-1}$, 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 $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$ depend on β as expected from Equation (25), and for a given β a clear dependence is seen as a function of $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$. Given the relatively small spread in ${T}_{\mathrm{obs}}$, the vertical gradient largely depicts the gradient in Γ and, as expected, the largest values of $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$ are seen toward low ${\beta }^{-1}$ 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.

Figure 7. Refer to the following caption and surrounding text.

Figure 7. Left: frequency uncertainties as a function of the height-to-background ratio (${\beta }^{-1}$) for the radial modes (l = 0). The color scale indicates the value of $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$ which, for a given value of β, is expected to define the spread in the uncertainties (Equation (24)). Dashed red lines show the behavior of the β-dependent factor $\sqrt{{f}_{0}(\beta )}$ (Equation (25)) for the values 0.009, 0.018, and 0.027 of $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$. The $\sqrt{{f}_{0}(\beta )}$ lines nicely follow the points of a given $\sqrt{{\rm{\Gamma }}/{T}_{\mathrm{obs}}}$, showing that the uncertainties behave as expected from an analytical approach (Section 4.1.1). Right: uncertainties for all mode frequencies extracted from the peak-bagging. The color indicates the angular degree l of the modes. The small right panel gives the Gaussian kernel density estimates of the uncertainties for each l.

Standard image High-resolution image

The 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 ${\tilde{V}}_{1}^{2}\gt {\tilde{V}}_{0}^{2}\gt {\tilde{V}}_{2}^{2}$ (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 ${\nu }_{\max }$ and ${\rm{\Delta }}\nu $, from the mode frequencies (see Table 2). The value of ${\nu }_{\max }$ was obtained by fitting a Gaussian function to the extracted amplitudes as a function of frequency. The value found for the solar ${\nu }_{\max }$ of 3078 ± 13 μHz is in agreement with the value of ${\nu }_{\max ,\odot }=3090\pm 30\,\mu \mathrm{Hz}$ by Huber et al. (2011). We obtained ${\rm{\Delta }}\nu $ and epsilon 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:

Equation (26)

Here ${\rm{\Delta }}{\nu }_{0}$ gives the value of the large separation at ${\nu }_{\max }.$ The latter corresponds to the radial order ${n}_{{\nu }_{\max }},$ which may take on a non-integer value. Specifically, ${n}_{{\nu }_{\max },l}$ is obtained for a given l from interpolating the measured mode frequencies (${\nu }_{l}$) against radial order (nl) to the frequency of ${\nu }_{\max }.$ This description gives a direct estimate of, for instance, the small frequency separations $\delta {\nu }_{01}$, given as the amount by which l = 1 modes are offset from the neighboring l = 0 modes ($\delta {\nu }_{01}=({\nu }_{n,0}\ -{\nu }_{n+\mathrm{1,0}})/2-{\nu }_{n,1}=2{D}_{0}$), and $\delta {\nu }_{02}={\nu }_{n,l=0}-{\nu }_{n-1,l=2}=6{D}_{0}$ (see, e.g., Bedding 2011). The value of $\delta {\nu }_{02}$ 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 ${\nu }_{\max }$ ${\rm{\Delta }}\nu $ $d{\rm{\Delta }}\nu /{dn}$ epsilon $\delta {\nu }_{01}$ $d\delta {\nu }_{01}/{dn}$ $\delta {\nu }_{02}$ $d\delta {\nu }_{02}/{dn}$
  (μHz) (μHz)     (μHz)   (μHz)  
1435467 1406.7${}_{-8.4}^{+6.3}$ 70.369${}_{-0.033}^{+0.034}$ 0.223${}_{-0.009}^{+0.010}$ 1.114${}_{-0.009}^{+0.009}$ 2.907${}_{-0.125}^{+0.125}$ −0.126${}_{-0.021}^{+0.020}$ 5.682${}_{-0.233}^{+0.252}$ −0.116${}_{-0.040}^{+0.040}$
2837475 1557.6${}_{-9.2}^{+8.2}$ 75.729${}_{-0.042}^{+0.041}$ 0.179${}_{-0.012}^{+0.012}$ 0.911${}_{-0.011}^{+0.011}$ 2.235${}_{-0.188}^{+0.177}$ −0.017${}_{-0.028}^{+0.028}$ 6.417${}_{-0.392}^{+0.408}$ −0.126${}_{-0.048}^{+0.049}$
3427720 2737.0${}_{-17.7}^{+10.7}$ 120.068${}_{-0.032}^{+0.031}$ 0.287${}_{-0.008}^{+0.010}$ 1.356${}_{-0.006}^{+0.006}$ 3.487${}_{-0.080}^{+0.079}$ −0.053${}_{-0.014}^{+0.014}$ 10.186${}_{-0.112}^{+0.100}$ −0.124${}_{-0.022}^{+0.022}$
3456181 970.0${}_{-5.9}^{+8.3}$ 52.264${}_{-0.039}^{+0.041}$ 0.216${}_{-0.011}^{+0.011}$ 0.988${}_{-0.014}^{+0.013}$ 3.551${}_{-0.167}^{+0.158}$ −0.075${}_{-0.023}^{+0.024}$ 4.370${}_{-0.248}^{+0.251}$ −0.030${}_{-0.043}^{+0.043}$
3632418 1166.8${}_{-3.8}^{+3.0}$ 60.704${}_{-0.018}^{+0.019}$ 0.232${}_{-0.004}^{+0.004}$ 1.114${}_{-0.006}^{+0.005}$ 3.697${}_{-0.062}^{+0.064}$ −0.041${}_{-0.010}^{+0.009}$ 4.189${}_{-0.115}^{+0.102}$ −0.053${}_{-0.019}^{+0.017}$
3656476 1925.0${}_{-6.3}^{+7.0}$ 93.194${}_{-0.020}^{+0.018}$ 0.207${}_{-0.008}^{+0.009}$ 1.445${}_{-0.004}^{+0.004}$ 4.608${}_{-0.042}^{+0.041}$ −0.095${}_{-0.011}^{+0.010}$ 4.554${}_{-0.051}^{+0.055}$ −0.157${}_{-0.014}^{+0.016}$
3735871 2862.6${}_{-26.5}^{+16.6}$ 123.049${}_{-0.046}^{+0.047}$ 0.260${}_{-0.019}^{+0.019}$ 1.325${}_{-0.009}^{+0.008}$ 3.823${}_{-0.124}^{+0.125}$ −0.004${}_{-0.029}^{+0.027}$ 11.028${}_{-0.208}^{+0.202}$ −0.077${}_{-0.053}^{+0.047}$
4914923 1817.0${}_{-5.2}^{+6.3}$ 88.531${}_{-0.019}^{+0.019}$ 0.233${}_{-0.008}^{+0.007}$ 1.377${}_{-0.004}^{+0.004}$ 4.534${}_{-0.042}^{+0.048}$ −0.119${}_{-0.010}^{+0.009}$ 5.378${}_{-0.074}^{+0.068}$ −0.118${}_{-0.021}^{+0.022}$
5184732 2089.3${}_{-4.1}^{+4.4}$ 95.545${}_{-0.023}^{+0.024}$ 0.216${}_{-0.006}^{+0.005}$ 1.374${}_{-0.005}^{+0.005}$ 2.772${}_{-0.047}^{+0.044}$ −0.104${}_{-0.009}^{+0.009}$ 5.863${}_{-0.065}^{+0.063}$ −0.097${}_{-0.015}^{+0.015}$
5773345 1101.2${}_{-6.6}^{+5.7}$ 57.303${}_{-0.027}^{+0.030}$ 0.257${}_{-0.007}^{+0.007}$ 1.077${}_{-0.010}^{+0.009}$ 1.754${}_{-0.101}^{+0.122}$ −0.093${}_{-0.016}^{+0.016}$ 4.131${}_{-0.261}^{+0.268}$ −0.047${}_{-0.046}^{+0.047}$
5950854 1926.7${}_{-20.4}^{+21.9}$ 96.629${}_{-0.107}^{+0.102}$ 0.208${}_{-0.029}^{+0.032}$ 1.431${}_{-0.020}^{+0.023}$ 4.988${}_{-0.178}^{+0.215}$ −0.212${}_{-0.038}^{+0.038}$ 4.713${}_{-0.235}^{+0.273}$ −0.232${}_{-0.074}^{+0.078}$
6106415 2248.6${}_{-3.9}^{+4.6}$ 104.074${}_{-0.026}^{+0.023}$ 0.254${}_{-0.005}^{+0.005}$ 1.343${}_{-0.005}^{+0.005}$ 3.422${}_{-0.047}^{+0.049}$ −0.114${}_{-0.008}^{+0.008}$ 6.881${}_{-0.070}^{+0.066}$ −0.118${}_{-0.013}^{+0.014}$
6116048 2126.9${}_{-5.0}^{+5.5}$ 100.754${}_{-0.017}^{+0.017}$ 0.258${}_{-0.005}^{+0.005}$ 1.336${}_{-0.003}^{+0.003}$ 3.687${}_{-0.052}^{+0.048}$ −0.143${}_{-0.009}^{+0.009}$ 6.034${}_{-0.070}^{+0.068}$ −0.155${}_{-0.016}^{+0.014}$
6225718 2364.2${}_{-4.6}^{+4.9}$ 105.695${}_{-0.017}^{+0.018}$ 0.274${}_{-0.005}^{+0.005}$ 1.225${}_{-0.004}^{+0.004}$ 3.207${}_{-0.060}^{+0.061}$ −0.053${}_{-0.010}^{+0.011}$ 8.741${}_{-0.088}^{+0.085}$ −0.062${}_{-0.017}^{+0.016}$
6508366 958.3${}_{-3.6}^{+4.6}$ 51.553${}_{-0.047}^{+0.046}$ 0.223${}_{-0.009}^{+0.009}$ 1.006${}_{-0.017}^{+0.017}$ 2.880${}_{-0.138}^{+0.132}$ −0.088${}_{-0.019}^{+0.019}$ 2.535${}_{-0.199}^{+0.205}$ −0.108${}_{-0.028}^{+0.030}$
6603624 2384.0${}_{-5.6}^{+5.4}$ 110.128${}_{-0.012}^{+0.012}$ 0.201${}_{-0.004}^{+0.004}$ 1.492${}_{-0.002}^{+0.002}$ 2.801${}_{-0.027}^{+0.029}$ −0.159${}_{-0.006}^{+0.006}$ 4.944${}_{-0.034}^{+0.031}$ −0.201${}_{-0.008}^{+0.008}$
6679371 941.8${}_{-5.0}^{+5.1}$ 50.601${}_{-0.029}^{+0.029}$ 0.181${}_{-0.007}^{+0.008}$ 0.880${}_{-0.010}^{+0.011}$ 2.861${}_{-0.116}^{+0.123}$ −0.016${}_{-0.016}^{+0.017}$ 3.143${}_{-0.266}^{+0.265}$ 0.044${}_{-0.034}^{+0.034}$
6933899 1389.9${}_{-3.6}^{+3.9}$ 72.135${}_{-0.018}^{+0.018}$ 0.255${}_{-0.005}^{+0.005}$ 1.319${}_{-0.005}^{+0.004}$ 5.314${}_{-0.041}^{+0.042}$ 0.013${}_{-0.009}^{+0.008}$ 4.910${}_{-0.054}^{+0.054}$ −0.063${}_{-0.015}^{+0.014}$
7103006 1167.9${}_{-6.9}^{+7.2}$ 59.658${}_{-0.030}^{+0.029}$ 0.211${}_{-0.008}^{+0.007}$ 0.978${}_{-0.009}^{+0.010}$ 2.504${}_{-0.107}^{+0.140}$ −0.045${}_{-0.018}^{+0.019}$ 4.471${}_{-0.295}^{+0.354}$ 0.031${}_{-0.048}^{+0.045}$
7106245 2397.9${}_{-28.7}^{+24.0}$ 111.376${}_{-0.061}^{+0.063}$ 0.246${}_{-0.032}^{+0.025}$ 1.392${}_{-0.012}^{+0.011}$ 3.489${}_{-0.118}^{+0.110}$ −0.113${}_{-0.031}^{+0.032}$ 6.529${}_{-0.167}^{+0.189}$ −0.265${}_{-0.068}^{+0.069}$
7206837 1652.5${}_{-11.7}^{+10.6}$ 79.131${}_{-0.039}^{+0.037}$ 0.253${}_{-0.011}^{+0.010}$ 1.054${}_{-0.010}^{+0.011}$ 2.106${}_{-0.147}^{+0.140}$ −0.054${}_{-0.023}^{+0.023}$ 6.619${}_{-0.417}^{+0.419}$ −0.094${}_{-0.073}^{+0.079}$
7296438 1847.8${}_{-12.6}^{+8.5}$ 88.698${}_{-0.036}^{+0.040}$ 0.242${}_{-0.015}^{+0.015}$ 1.358${}_{-0.008}^{+0.009}$ 4.505${}_{-0.081}^{+0.073}$ −0.055${}_{-0.022}^{+0.020}$ 5.079${}_{-0.098}^{+0.088}$ −0.135${}_{-0.028}^{+0.029}$
7510397 1189.1${}_{-4.4}^{+3.4}$ 62.249${}_{-0.020}^{+0.020}$ 0.258${}_{-0.004}^{+0.004}$ 1.112${}_{-0.006}^{+0.006}$ 3.971${}_{-0.059}^{+0.072}$ 0.003${}_{-0.009}^{+0.009}$ 4.370${}_{-0.086}^{+0.086}$ −0.016${}_{-0.017}^{+0.015}$
7680114 1709.1${}_{-6.5}^{+7.1}$ 85.145${}_{-0.043}^{+0.039}$ 0.238${}_{-0.007}^{+0.007}$ 1.368${}_{-0.009}^{+0.010}$ 5.039${}_{-0.054}^{+0.050}$ −0.043${}_{-0.011}^{+0.011}$ 4.980${}_{-0.072}^{+0.074}$ −0.108${}_{-0.016}^{+0.018}$
7771282 1465.1${}_{-18.7}^{+27.0}$ 72.463${}_{-0.079}^{+0.069}$ 0.368${}_{-0.039}^{+0.042}$ 1.117${}_{-0.019}^{+0.021}$ 3.527${}_{-0.238}^{+0.241}$ −0.118${}_{-0.050}^{+0.049}$ 5.058${}_{-0.445}^{+0.484}$ −0.020${}_{-0.104}^{+0.089}$
7871531 3455.9${}_{-26.5}^{+19.3}$ 151.329${}_{-0.023}^{+0.025}$ 0.285${}_{-0.008}^{+0.009}$ 1.504${}_{-0.004}^{+0.003}$ 2.142${}_{-0.068}^{+0.071}$ −0.149${}_{-0.014}^{+0.013}$ 7.350${}_{-0.169}^{+0.155}$ −0.143${}_{-0.049}^{+0.044}$
7940546 1116.6${}_{-3.6}^{+3.3}$ 58.762${}_{-0.029}^{+0.029}$ 0.217${}_{-0.005}^{+0.004}$ 1.075${}_{-0.009}^{+0.009}$ 3.985${}_{-0.071}^{+0.072}$ −0.002${}_{-0.011}^{+0.011}$ 4.346${}_{-0.123}^{+0.133}$ 0.012${}_{-0.019}^{+0.019}$
7970740 4197.4${}_{-18.4}^{+21.2}$ 173.541${}_{-0.068}^{+0.060}$ 0.272${}_{-0.005}^{+0.005}$ 1.455${}_{-0.008}^{+0.010}$ 2.356${}_{-0.084}^{+0.083}$ −0.097${}_{-0.010}^{+0.011}$ 7.901${}_{-0.165}^{+0.169}$ −0.268${}_{-0.026}^{+0.025}$
8006161 3574.7${}_{-10.5}^{+11.4}$ 149.427${}_{-0.014}^{+0.015}$ 0.195${}_{-0.005}^{+0.005}$ 1.547${}_{-0.002}^{+0.002}$ 3.061${}_{-0.046}^{+0.041}$ −0.084${}_{-0.007}^{+0.008}$ 9.680${}_{-0.063}^{+0.070}$ −0.150${}_{-0.012}^{+0.012}$
8150065 1876.9${}_{-32.4}^{+38.1}$ 89.264${}_{-0.121}^{+0.134}$ 0.403${}_{-0.047}^{+0.048}$ 1.163${}_{-0.030}^{+0.029}$ 3.027${}_{-0.203}^{+0.198}$ 0.036${}_{-0.063}^{+0.079}$ 6.357${}_{-0.348}^{+0.342}$ 0.012${}_{-0.128}^{+0.135}$
8179536 2074.9${}_{-12.0}^{+13.8}$ 95.090${}_{-0.054}^{+0.058}$ 0.277${}_{-0.019}^{+0.019}$ 1.153${}_{-0.013}^{+0.012}$ 3.137${}_{-0.164}^{+0.171}$ −0.082${}_{-0.032}^{+0.034}$ 8.245${}_{-0.315}^{+0.352}$ −0.041${}_{-0.070}^{+0.073}$
8228742 1190.5${}_{-3.7}^{+3.4}$ 62.071${}_{-0.021}^{+0.022}$ 0.244${}_{-0.005}^{+0.005}$ 1.158${}_{-0.007}^{+0.006}$ 4.371${}_{-0.061}^{+0.057}$ 0.007${}_{-0.010}^{+0.010}$ 4.517${}_{-0.087}^{+0.082}$ −0.048${}_{-0.020}^{+0.019}$
8379927 2795.3${}_{-5.7}^{+6.0}$ 120.288${}_{-0.018}^{+0.017}$ 0.232${}_{-0.005}^{+0.005}$ 1.311${}_{-0.003}^{+0.003}$ 3.676${}_{-0.062}^{+0.062}$ −0.058${}_{-0.011}^{+0.011}$ 10.932${}_{-0.083}^{+0.096}$ −0.062${}_{-0.020}^{+0.022}$
8394589 2396.7${}_{-9.4}^{+10.5}$ 109.488${}_{-0.035}^{+0.034}$ 0.234${}_{-0.011}^{+0.011}$ 1.267${}_{-0.007}^{+0.006}$ 3.382${}_{-0.091}^{+0.089}$ −0.061${}_{-0.021}^{+0.022}$ 7.979${}_{-0.161}^{+0.164}$ −0.106${}_{-0.045}^{+0.045}$
8424992 2533.7${}_{-28.1}^{+27.0}$ 120.584${}_{-0.064}^{+0.062}$ 0.120${}_{-0.038}^{+0.035}$ 1.517${}_{-0.010}^{+0.012}$ 2.678${}_{-0.103}^{+0.100}$ −0.179${}_{-0.032}^{+0.033}$ 5.190${}_{-0.143}^{+0.122}$ −0.282${}_{-0.055}^{+0.046}$
8694723 1470.5${}_{-4.1}^{+3.7}$ 75.112${}_{-0.021}^{+0.019}$ 0.296${}_{-0.005}^{+0.005}$ 1.113${}_{-0.005}^{+0.005}$ 5.339${}_{-0.067}^{+0.068}$ 0.005${}_{-0.009}^{+0.010}$ 5.879${}_{-0.108}^{+0.111}$ 0.012${}_{-0.021}^{+0.020}$
8760414 2455.3${}_{-8.3}^{+9.1}$ 117.230${}_{-0.018}^{+0.022}$ 0.295${}_{-0.007}^{+0.007}$ 1.400${}_{-0.004}^{+0.003}$ 4.403${}_{-0.053}^{+0.045}$ −0.280${}_{-0.012}^{+0.009}$ 5.132${}_{-0.059}^{+0.063}$ −0.291${}_{-0.015}^{+0.014}$
8938364 1675.1${}_{-5.8}^{+5.2}$ 85.684${}_{-0.020}^{+0.018}$ 0.235${}_{-0.006}^{+0.007}$ 1.444${}_{-0.004}^{+0.004}$ 6.491${}_{-0.035}^{+0.044}$ −0.046${}_{-0.009}^{+0.010}$ 5.184${}_{-0.052}^{+0.048}$ −0.119${}_{-0.014}^{+0.013}$
9025370 2988.6${}_{-16.9}^{+20.0}$ 132.628${}_{-0.024}^{+0.030}$ 0.205${}_{-0.015}^{+0.016}$ 1.475${}_{-0.005}^{+0.004}$ 3.099${}_{-0.062}^{+0.065}$ −0.066${}_{-0.017}^{+0.018}$ 9.141${}_{-0.119}^{+0.113}$ −0.126${}_{-0.030}^{+0.037}$
9098294 2314.7${}_{-10.4}^{+9.2}$ 108.894${}_{-0.022}^{+0.023}$ 0.251${}_{-0.008}^{+0.009}$ 1.439${}_{-0.005}^{+0.004}$ 3.331${}_{-0.057}^{+0.053}$ −0.185${}_{-0.011}^{+0.012}$ 5.265${}_{-0.088}^{+0.086}$ −0.228${}_{-0.025}^{+0.024}$
9139151 2690.4${}_{-9.0}^{+14.5}$ 117.294${}_{-0.032}^{+0.031}$ 0.240${}_{-0.010}^{+0.011}$ 1.337${}_{-0.006}^{+0.006}$ 3.557${}_{-0.077}^{+0.083}$ −0.019${}_{-0.020}^{+0.017}$ 10.050${}_{-0.158}^{+0.162}$ −0.093${}_{-0.057}^{+0.055}$
9139163 1729.8${}_{-5.9}^{+6.2}$ 81.170${}_{-0.036}^{+0.042}$ 0.241${}_{-0.005}^{+0.005}$ 1.007${}_{-0.010}^{+0.010}$ 2.079${}_{-0.098}^{+0.109}$ −0.024${}_{-0.010}^{+0.012}$ 6.213${}_{-0.215}^{+0.218}$ 0.040${}_{-0.028}^{+0.026}$
9206432 1866.4${}_{-14.9}^{+10.3}$ 84.926${}_{-0.051}^{+0.046}$ 0.135${}_{-0.013}^{+0.013}$ 0.958${}_{-0.012}^{+0.012}$ 3.235${}_{-0.210}^{+0.207}$ 0.001${}_{-0.031}^{+0.027}$ 7.115${}_{-0.411}^{+0.388}$ 0.004${}_{-0.055}^{+0.056}$
9353712 934.3${}_{-8.3}^{+11.1}$ 51.467${}_{-0.104}^{+0.091}$ 0.254${}_{-0.011}^{+0.012}$ 1.095${}_{-0.031}^{+0.038}$ 3.536${}_{-0.159}^{+0.146}$ −0.038${}_{-0.018}^{+0.018}$ 3.907${}_{-0.250}^{+0.236}$ −0.090${}_{-0.031}^{+0.031}$
9410862 2278.8${}_{-16.6}^{+31.2}$ 107.390${}_{-0.053}^{+0.050}$ 0.223${}_{-0.021}^{+0.020}$ 1.343${}_{-0.010}^{+0.009}$ 3.625${}_{-0.121}^{+0.116}$ −0.181${}_{-0.032}^{+0.027}$ 6.098${}_{-0.204}^{+0.201}$ −0.239${}_{-0.086}^{+0.081}$
9414417 1155.3${}_{-4.6}^{+6.1}$ 60.115${}_{-0.024}^{+0.024}$ 0.237${}_{-0.007}^{+0.006}$ 1.045${}_{-0.008}^{+0.008}$ 3.572${}_{-0.092}^{+0.106}$ −0.029${}_{-0.014}^{+0.015}$ 4.648${}_{-0.200}^{+0.211}$ −0.006${}_{-0.030}^{+0.031}$
9812850 1255.2${}_{-7.0}^{+9.1}$ 64.746${}_{-0.068}^{+0.067}$ 0.240${}_{-0.011}^{+0.012}$ 1.067${}_{-0.020}^{+0.021}$ 2.654${}_{-0.166}^{+0.156}$ −0.134${}_{-0.025}^{+0.025}$ 4.418${}_{-0.307}^{+0.297}$ −0.032${}_{-0.045}^{+0.048}$
9955598 3616.8${}_{-29.6}^{+21.2}$ 153.283${}_{-0.032}^{+0.029}$ 0.195${}_{-0.011}^{+0.011}$ 1.529${}_{-0.004}^{+0.005}$ 2.796${}_{-0.076}^{+0.073}$ −0.095${}_{-0.017}^{+0.016}$ 8.941${}_{-0.126}^{+0.143}$ −0.172${}_{-0.030}^{+0.029}$
9965715 2079.3${}_{-10.4}^{+9.2}$ 97.236${}_{-0.042}^{+0.041}$ 0.373${}_{-0.015}^{+0.016}$ 1.139${}_{-0.009}^{+0.009}$ 3.685${}_{-0.129}^{+0.132}$ −0.102${}_{-0.027}^{+0.024}$ 7.958${}_{-0.269}^{+0.254}$ −0.105${}_{-0.058}^{+0.063}$
10068307 995.1${}_{-2.7}^{+2.8}$ 53.945${}_{-0.020}^{+0.019}$ 0.247${}_{-0.003}^{+0.004}$ 1.131${}_{-0.006}^{+0.007}$ 4.220${}_{-0.058}^{+0.055}$ 0.014${}_{-0.008}^{+0.008}$ 3.799${}_{-0.088}^{+0.083}$ −0.041${}_{-0.014}^{+0.014}$
10079226 2653.0${}_{-44.3}^{+47.7}$ 116.345${}_{-0.052}^{+0.059}$ 0.264${}_{-0.028}^{+0.027}$ 1.350${}_{-0.011}^{+0.010}$ 3.236${}_{-0.199}^{+0.205}$ −0.014${}_{-0.046}^{+0.043}$ 9.387${}_{-0.371}^{+0.401}$ 0.098${}_{-0.098}^{+0.098}$
10162436 1052.0${}_{-4.2}^{+4.0}$ 55.725${}_{-0.039}^{+0.035}$ 0.242${}_{-0.004}^{+0.004}$ 1.106${}_{-0.012}^{+0.013}$ 3.746${}_{-0.080}^{+0.078}$ −0.033${}_{-0.009}^{+0.009}$ 3.706${}_{-0.117}^{+0.115}$ −0.031${}_{-0.018}^{+0.017}$
10454113 2357.2${}_{-9.1}^{+8.2}$ 105.063${}_{-0.033}^{+0.031}$ 0.283${}_{-0.009}^{+0.009}$ 1.206${}_{-0.007}^{+0.007}$ 3.059${}_{-0.110}^{+0.095}$ −0.063${}_{-0.020}^{+0.016}$ 9.426${}_{-0.165}^{+0.178}$ −0.079${}_{-0.033}^{+0.033}$
10516096 1689.8${}_{-5.8}^{+4.6}$ 84.424${}_{-0.025}^{+0.022}$ 0.257${}_{-0.008}^{+0.009}$ 1.318${}_{-0.006}^{+0.005}$ 5.001${}_{-0.061}^{+0.054}$ −0.058${}_{-0.012}^{+0.013}$ 5.248${}_{-0.083}^{+0.084}$ −0.089${}_{-0.022}^{+0.023}$
10644253 2899.7${}_{-22.8}^{+21.3}$ 123.080${}_{-0.055}^{+0.056}$ 0.250${}_{-0.016}^{+0.016}$ 1.313${}_{-0.011}^{+0.010}$ 3.863${}_{-0.139}^{+0.140}$ 0.004${}_{-0.030}^{+0.029}$ 11.378${}_{-0.155}^{+0.192}$ −0.138${}_{-0.049}^{+0.049}$
10730618 1282.1${}_{-12.7}^{+14.6}$ 66.333${}_{-0.064}^{+0.061}$ 0.239${}_{-0.017}^{+0.016}$ 1.032${}_{-0.018}^{+0.018}$ 2.651${}_{-0.218}^{+0.199}$ 0.019${}_{-0.037}^{+0.040}$ 4.556${}_{-0.454}^{+0.428}$ 0.113${}_{-0.080}^{+0.073}$
10963065 2203.7${}_{-6.3}^{+6.7}$ 103.179${}_{-0.027}^{+0.027}$ 0.297${}_{-0.008}^{+0.008}$ 1.275${}_{-0.005}^{+0.005}$ 3.567${}_{-0.071}^{+0.072}$ −0.081${}_{-0.014}^{+0.014}$ 7.083${}_{-0.096}^{+0.103}$ −0.058${}_{-0.020}^{+0.022}$
11081729 1968.3${}_{-12.6}^{+11.0}$ 90.116${}_{-0.047}^{+0.048}$ 0.242${}_{-0.017}^{+0.016}$ 1.020${}_{-0.011}^{+0.011}$ 3.056${}_{-0.264}^{+0.251}$ −0.103${}_{-0.026}^{+0.031}$ 6.602${}_{-0.664}^{+0.605}$ 0.010${}_{-0.084}^{+0.087}$
11253226 1590.6${}_{-6.8}^{+10.6}$ 76.858${}_{-0.030}^{+0.026}$ 0.183${}_{-0.008}^{+0.008}$ 0.920${}_{-0.008}^{+0.008}$ 1.748${}_{-0.136}^{+0.155}$ −0.132${}_{-0.018}^{+0.020}$ 6.973${}_{-0.396}^{+0.435}$ 0.039${}_{-0.049}^{+0.048}$
11772920 3674.7${}_{-36.1}^{+55.1}$ 157.746${}_{-0.033}^{+0.032}$ 0.238${}_{-0.015}^{+0.015}$ 1.516${}_{-0.005}^{+0.004}$ 2.366${}_{-0.084}^{+0.078}$ −0.082${}_{-0.019}^{+0.017}$ 7.849${}_{-0.209}^{+0.200}$ −0.131${}_{-0.063}^{+0.059}$
12009504 1865.6${}_{-6.2}^{+7.7}$ 88.217${}_{-0.025}^{+0.026}$ 0.289${}_{-0.007}^{+0.008}$ 1.200${}_{-0.006}^{+0.006}$ 3.533${}_{-0.079}^{+0.078}$ −0.072${}_{-0.014}^{+0.014}$ 6.117${}_{-0.134}^{+0.133}$ −0.066${}_{-0.034}^{+0.032}$
12069127 884.7${}_{-8.0}^{+10.1}$ 48.400${}_{-0.048}^{+0.048}$ 0.204${}_{-0.014}^{+0.012}$ 1.061${}_{-0.018}^{+0.018}$ 3.399${}_{-0.173}^{+0.153}$ −0.037${}_{-0.025}^{+0.023}$ 3.650${}_{-0.291}^{+0.244}$ −0.102${}_{-0.039}^{+0.039}$
12069424 2188.5${}_{-3.0}^{+4.6}$ 103.277${}_{-0.020}^{+0.021}$ 0.246${}_{-0.004}^{+0.004}$ 1.437${}_{-0.004}^{+0.004}$ 3.392${}_{-0.039}^{+0.039}$ −0.152${}_{-0.006}^{+0.007}$ 5.274${}_{-0.047}^{+0.051}$ −0.181${}_{-0.010}^{+0.010}$
12069449 2561.3${}_{-5.6}^{+5.0}$ 116.929${}_{-0.013}^{+0.012}$ 0.201${}_{-0.004}^{+0.004}$ 1.461${}_{-0.002}^{+0.002}$ 2.707${}_{-0.031}^{+0.032}$ −0.117${}_{-0.006}^{+0.006}$ 6.045${}_{-0.044}^{+0.037}$ −0.133${}_{-0.010}^{+0.009}$
12258514 1512.7${}_{-2.9}^{+3.3}$ 74.799${}_{-0.015}^{+0.016}$ 0.209${}_{-0.004}^{+0.004}$ 1.281${}_{-0.004}^{+0.004}$ 4.227${}_{-0.043}^{+0.042}$ −0.061${}_{-0.008}^{+0.007}$ 4.827${}_{-0.052}^{+0.061}$ −0.056${}_{-0.011}^{+0.011}$
12317678 1212.4${}_{-4.9}^{+5.5}$ 63.464${}_{-0.024}^{+0.025}$ 0.231${}_{-0.005}^{+0.005}$ 0.928${}_{-0.008}^{+0.006}$ 3.883${}_{-0.115}^{+0.112}$ −0.032${}_{-0.013}^{+0.013}$ 5.273${}_{-0.194}^{+0.188}$ −0.061${}_{-0.025}^{+0.023}$

Note.  ${\nu }_{\max }$ is obtained from the fit to the extracted modes. Plots of the fitted parameters are shown in Figures 912; the values for epsilon are shown in Figure 3.

Download table as:  ASCIITypeset images: 1 2

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.

Equation (27)

where the term $l(l+1){D}_{0}$ has been replaced by $\delta {\nu }_{0l}$, which takes on independent values for l > 0. Thereby, we optimize independently for the separations $\delta {\nu }_{02}$ and $\delta {\nu }_{01}$. In the fit, we also included ${\nu }_{\max }$ with a Gaussian prior from the fit to the mode amplitudes; the values of ${\nu }_{\max },$ ${\rm{\Delta }}\nu $, and epsilon correlate strongly, but we include ${\nu }_{\max }$ to properly marginalize over the uncertainty of the pivoting ${n}_{{\nu }_{\max },l}$. 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.

Figure 8. Refer to the following caption and surrounding text.

Figure 8. Example of a fit of Equation (27) to individual frequencies of KIC 8228742 (Horace), displayed in échelle diagram format. The angular degrees of the ridges are indicated in the top of the plot; the numbers on the right side of the plot give the radial orders n of the radial l = 0 modes. The full lines give the solution from the median of the posterior distributions of the fit; the dark-gray lines give 500 solutions with parameters drawn from the posterior distributions. The small separation $\delta {\nu }_{02}$ is given by the difference between the l = 0 and 2 ridges at ${\nu }_{\max }.$ The dashed line gives the l = 0 ridge offset by ${\rm{\Delta }}\nu /2$, hence $\delta {\nu }_{01}$ is given by the difference between this line and the l = 1 ridge at ${\nu }_{\max }.$ The dotted line gives the expected position of the l = 1 ridge from the asymptotic relation (Equation (10)), where $\delta {\nu }_{01}=\delta {\nu }_{02}/3$. A clear oscillatory behavior from acoustic glitches is seen for the frequencies around the median solution.

Standard image High-resolution image

Figure 9 shows the estimated values for ${\rm{\Delta }}\nu $ versus ${\nu }_{\max },$ 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 ${\rm{\Delta }}\nu $ for a given ${\nu }_{\max }.$ Also shown is the change in the large separation, which is found to be near-constant at $d{\rm{\Delta }}\nu /{dn}\,\approx \,0.25\,\mu \mathrm{Hz}$ 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 $d{\rm{\Delta }}\nu /{dn}$ against ${\rm{\Delta }}{\nu }_{0}$ would correspond to a linear change in ${\rm{\Delta }}\nu $ with frequency.

Figure 9. Refer to the following caption and surrounding text.

Figure 9. Left: measured values of ${\rm{\Delta }}\nu $ against ${\nu }_{\max }$ from fitting of Equation (27), with the color indicating the modeled mass. The dashed line shows the empirical relation by Huber et al. (2011), with the 1 and $2\sigma $ uncertainties on the relation given by the dark and light blue regions. From our values we calculate a correlation of ${\rho }_{\alpha ,\beta }=-0.99913$ between then α and β parameters of the Huber et al. (2011) relation—this was included to estimate the uncertainty regions. The bottom panel shows the residuals of the measured values to the relation. Measured values are seen to overall agree well with expectations, with a modest mass gradient across the residuals. Right: measured gradient of ${\rm{\Delta }}\nu $ with radial order, $d{\rm{\Delta }}\nu /{dn}$, defining the overall curvature of the ridges in the échelle diagram against ${\rm{\Delta }}\nu $.

Standard image High-resolution image

For the Sun, we obtain a value of ${\rm{\Delta }}{\nu }_{\odot }$ of 134.91 ± 0.02 μHz; comparing this to the value of ${\rm{\Delta }}{\nu }_{\odot }=135.1\pm 0.1\,\mu \mathrm{Hz}$ by Huber et al. (2011), we would expect a slightly smaller value, because the curvature is incorporated in our estimation. However, the value of ${\rm{\Delta }}{\nu }_{\odot }=134.92\pm 0.02\,\mu \mathrm{Hz}$ 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 epsilon have already been shown in Figure 3. For the estimation of epsilon, it should be noted that this value is strongly anti-correlated with ${\rm{\Delta }}\nu $, and a small change in ${\rm{\Delta }}\nu $ can thus induce a significant change in epsilon (White et al. 2011b). This will especially take effect if the estimate of ${\nu }_{\max }$ is off, because ${\nu }_{\max }$ defines the pivoting point for the change in $\delta {\nu }_{01}$, $\delta {\nu }_{02}$, and the curvatures. Like White et al. (2011b), we see an offset between the estimated epsilon 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 epsilon-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 epsilon 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 ${\rm{\Delta }}\nu $ and epsilon 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 ${\nu }_{\max }.$ Comparing the values for ${\rm{\Delta }}\nu $ and epsilon from the fit of Equation (10) versus Equation (27), we only find minor differences in estimated values. For ${\rm{\Delta }}\nu $ the maximum absolute difference was ∼0.36 μHz, with no systematic differences; for epsilon 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 $\delta {\nu }_{01}$ and $\delta {\nu }_{02}$, together with their gradients in n (or frequency); the results for $\delta {\nu }_{02}$ 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 $\delta {\nu }_{01}$ and $\delta {\nu }_{02}$ 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 $\delta {\nu }_{02}/\delta {\nu }_{01}=3$ from assuming $\delta {\nu }_{0l}=l(l+1){D}_{0}$—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 $\delta {\nu }_{02}/\delta {\nu }_{01}$ 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 $\delta {\nu }_{02}/\delta {\nu }_{01}=3$ 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 ${\rm{\Delta }}{\nu }_{0}$ 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.

Figure 10. Refer to the following caption and surrounding text.

Figure 10. Left: measured values of $\delta {\nu }_{01}$ against ${\rm{\Delta }}\nu $ from fitting of Equation (27), with the color indicating the modeled mass. Right: measured change in $\delta {\nu }_{01}$ with radial order, $d\delta {\nu }_{01}/{dn}$ against ${\rm{\Delta }}\nu $. All slopes in $\delta {\nu }_{01}$ are seen to be either negative or consistent with zero.

Standard image High-resolution image
Figure 11. Refer to the following caption and surrounding text.

Figure 11. Left: measured values of $\delta {\nu }_{02}$ against ${\rm{\Delta }}\nu $ from fitting of Equation (27), with the color indicating the modeled mass. Shown are also the $\delta {\nu }_{02}$-evolutionary tracks from White et al. (2011a), calculated from ASTEC evolutionary tracks with Z0 = 0.017 (Christensen-Dalsgaard 2008), for masses going from $0.8\,{M}_{\odot }$ to $1.6\,{M}_{\odot }$ in steps of $0.1\,{M}_{\odot }$. For the $M=1\,{M}_{\odot }$ track we have indicated the effect of changing the metallicity to ${Z}_{0}=0.011$ (dashed) and Z0 = 0.028 (dotted). Right: measured change in $\delta {\nu }_{02}$ with radial order, $d\delta {\nu }_{02}/{dn}$ against ${\rm{\Delta }}\nu $. All slopes in $\delta {\nu }_{02}$ are seen to be either negative or consistent with zero.

Standard image High-resolution image
Figure 12. Refer to the following caption and surrounding text.

Figure 12. Left: ratio between measured values of $\delta {\nu }_{02}$ and $\delta {\nu }_{01}$ against ${\rm{\Delta }}\nu $, with the color indicating the modeled mass. The dashed line indicates the expected value of 3 from the asymptotic relation Equation (10) where $\delta {\nu }_{0l}=l(l+1){D}_{0}$. Right: ratio plotted against the modeled value for central hydrogen content Xc, which is a good indicator for the evolutionary state. More evolved stars with low Xc are seen to deviate more from the asymptotic expectation than lesser evolved stars.

Standard image High-resolution image

4.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 ${\nu }_{\max }$(right). One clearly sees the overall decrease in amplitude with increasing ${\nu }_{\max }$(Kjeldsen et al. 2005; Arentoft et al. 2008), and thus decreasing ${T}_{\mathrm{eff}}$ (Kjeldsen & Bedding 1995, 2011; Appourchaux et al. 2014). The value of ${\nu }_{\max }$ was estimated from the amplitudes as the frequency where the modes of oscillation show their maximum amplitude, ${A}_{\max }.$

Figure 13. Refer to the following caption and surrounding text.

Figure 13. Left: radial mode amplitude envelopes against frequency, with the color indicating the ${T}_{\mathrm{eff}}$. For a better visualization, the amplitudes have been smoothed with the five-point Epanechnikov filter. The dashed envelope gives the results obtained for the Sun. Right: amplitudes centered on ${\nu }_{\max }$ and plotted against a proxy for the radial order.

Standard image High-resolution image

Figure 14 gives the measured values for ${A}_{\max }$ against ${\nu }_{\max }$(Table 3). We see both the expected change with ${\nu }_{\max },$ in addition to a mass gradient across the overall decrease in ${A}_{\max }$ with ${\nu }_{\max }$(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):

Equation (28)

with t = 1.32 ± 0.02, s = 0.838 ± 0.002, and r = 2 and Corsaro et al. (2013):

Equation (29)

with $s=0.748\pm 0.015$, t = 1.27 ± 0.04, r = 3.47 ± 0.09, and $\mathrm{ln}\beta =0.321\pm 0.020$. In converting the bolometric amplitudes from these relations to the Kepler bandpass, we used the root-mean-square value of ${A}_{\mathrm{bol},\odot }=2.53\pm 0.11\,\mathrm{ppm}$ from Michel et al. (2009) and the temperature dependent bolometric correction ${c}_{{\rm{K}}}({T}_{\mathrm{eff}})$ 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 ${A}_{\max }$ values obtained using the method of Kjeldsen et al. (2005, 2008a) and the high ${\nu }_{\max }$ 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).

Figure 14. Refer to the following caption and surrounding text.

Figure 14. Left: estimates of radial mode amplitudes at ${\nu }_{\max }$ as a function of ${\nu }_{\max };$ the color indicates the modeled mass. The ${A}_{\max }$ values are obtained from a fit of the individual mode amplitudes against frequency. Right: fractional differences between observed values of ${A}_{\max }$ and those estimated by the relations of Huber et al. (2011) and Corsaro et al. (2013).

Standard image High-resolution image

Table 3.  Values for the Line Width at ${\nu }_{\max }$(Figure 20) from the Fit of Equation (30)

KIC ${\rm{\Gamma }}\,@\,{\nu }_{\max }$ ${A}_{\max }$ ${A}_{\max ,\mathrm{smo}}(c=3.04)$ ${A}_{\max ,\mathrm{smo}}(c={\tilde{V}}_{\mathrm{tot}}^{2})$ α ${{\rm{\Gamma }}}_{\alpha }$ ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$ ${W}_{\mathrm{dip}}$ ${\nu }_{\mathrm{dip}}$ ${\mathrm{FWHM}}_{\mathrm{dip}}$
  $(\mu \mathrm{Hz})$ $(\mathrm{ppm})$ $(\mathrm{ppm})$ $(\mathrm{ppm})$   $(\mu \mathrm{Hz})$ $(\mu \mathrm{Hz})$ $(\mu \mathrm{Hz})$ $(\mu \mathrm{Hz})$ $(\mu \mathrm{Hz})$
1435467 5.18 ± 0.17 4.37${}_{-0.07}^{+0.10}$ 4.30 4.25 0.54${}_{-0.22}^{+0.21}$ 5.19 ± 0.18
2837475 6.39 ± 0.2 3.87${}_{-0.05}^{+0.06}$ 3.67 3.71 0.91${}_{-0.19}^{+0.19}$ 6.42 ± 0.21
3427720 1.88 ± 0.15 2.62${}_{-0.07}^{+0.04}$ 2.55 2.38 2.49${}_{-0.96}^{+0.80}$ 5.0 ± 2.96 0.36${}_{-0.26}^{+0.13}$ 4428${}_{-1660}^{+605}$ 2695${}_{-100}^{+126}$ 1309 ± 711
3456181 4.07 ± 0.19 5.60${}_{-0.14}^{+0.14}$ 5.92 5.48 1.11${}_{-0.26}^{+0.27}$ 4.06 ± 0.2
3632418 2.72 ± 0.16 5.04${}_{-0.07}^{+0.06}$ 5.04 4.76 1.76${}_{-0.22}^{+0.18}$ 4.67 ± 0.71 0.56${}_{-0.07}^{+0.09}$ 1698${}_{-304}^{+210}$ 1188${}_{-34}^{+26}$ 448 ± 183
3656476 0.69 ± 0.01 4.70${}_{-0.13}^{+0.13}$ 4.13 4.12 2.50${}_{-0.10}^{+0.13}$ 0.69 ± 0.01
3735871 2.25 ± 0.25 2.52${}_{-0.07}^{+0.08}$ 2.38 2.30 2.79${}_{-1.24}^{+1.18}$ 7.34 ± 4.29 0.28${}_{-0.18}^{+0.10}$ 4717${}_{-1110}^{+798}$ 2825${}_{-167}^{+148}$ 1426 ± 595
4914923 1.16 ± 0.08 4.62${}_{-0.08}^{+0.12}$ 4.19 4.15 3.52${}_{-0.57}^{+0.48}$ 6.2 ± 2.4 0.18${}_{-0.07}^{+0.07}$ 2600${}_{-270}^{+201}$ 1847${}_{-25}^{+24}$ 666 ± 170
5184732 1.36 ± 0.06 4.08${}_{-0.06}^{+0.06}$ 3.84 3.88 3.52${}_{-0.41}^{+0.38}$ 7.65 ± 2.61 0.17${}_{-0.06}^{+0.06}$ 3247${}_{-305}^{+243}$ 2117${}_{-28}^{+26}$ 941 ± 184
5773345 3.37 ± 0.13 5.69${}_{-0.11}^{+0.13}$ 5.17 5.52 0.58${}_{-0.23}^{+0.26}$ 3.37 ± 0.14
5950854 0.83 ± 0.01 4.12${}_{-0.27}^{+0.37}$ 3.77 3.87 0.56${}_{-0.03}^{+0.03}$ 0.83 ± 0.01
6106415 1.64 ± 0.07 3.86${}_{-0.04}^{+0.06}$ 3.60 3.63 3.49${}_{-0.25}^{+0.25}$ 5.53 ± 1.08 0.29${}_{-0.05}^{+0.06}$ 3242${}_{-265}^{+201}$ 2274${}_{-26}^{+23}$ 837 ± 166
6116048 1.62 ± 0.09 4.02${}_{-0.06}^{+0.06}$ 3.75 3.67 3.26${}_{-0.28}^{+0.32}$ 5.74 ± 1.35 0.28${}_{-0.06}^{+0.06}$ 3102${}_{-319}^{+207}$ 2098${}_{-28}^{+27}$ 796 ± 182
6225718 2.58 ± 0.12 3.53${}_{-0.03}^{+0.04}$ 3.39 3.29 2.32${}_{-0.22}^{+0.24}$ 5.77 ± 0.67 0.43${}_{-0.05}^{+0.05}$ 3076${}_{-180}^{+139}$ 2316${}_{-23}^{+26}$ 611 ± 121
6508366 5.22 ± 0.17 5.02${}_{-0.07}^{+0.09}$ 5.21 4.75 1.95${}_{-0.12}^{+0.13}$ 5.21 ± 0.15
6603624 0.56 ± 0.02 4.15${}_{-0.07}^{+0.10}$ 3.55 3.58 4.42${}_{-0.24}^{+0.19}$ 0.57 ± 0.02
6679371 4.53 ± 0.18 5.56${}_{-0.08}^{+0.07}$ 5.33 5.28 0.45${}_{-0.45}^{+0.19}$ 5.51 ± 0.74 0.69${}_{-0.16}^{+0.14}$ 1339${}_{-389}^{+209}$ 833${}_{-48}^{+44}$ 295 ± 190
6933899 1.3 ± 0.07 5.55${}_{-0.10}^{+0.09}$ 5.07 4.94 3.15${}_{-0.69}^{+0.56}$ 6.0 ± 2.75 0.21${}_{-0.11}^{+0.08}$ 2218${}_{-304}^{+233}$ 1422${}_{-33}^{+30}$ 671 ± 178
7103006 5.09 ± 0.16 4.74${}_{-0.09}^{+0.11}$ 4.29 4.58 1.17${}_{-0.18}^{+0.18}$ 5.09 ± 0.17
7106245 1.65 ± 0.14 3.39${}_{-0.15}^{+0.11}$ 2.96 2.91 1.71${}_{-0.97}^{+0.93}$ 1.64 ± 0.14
7206837 4.34 ± 0.32 4.32${}_{-0.06}^{+0.10}$ 3.85 4.18 0.75${}_{-0.62}^{+0.36}$ 8.25 ± 2.83 0.50${}_{-0.21}^{+0.16}$ 2854${}_{-830}^{+545}$ 1720${}_{-101}^{+92}$ 952 ± 433
7296438 1.22 ± 0.12 4.59${}_{-0.22}^{+0.20}$ 4.24 4.11 3.87${}_{-1.11}^{+1.10}$ 8.39 ± 4.56 0.14${}_{-0.08}^{+0.05}$ 2794${}_{-336}^{+272}$ 1835${}_{-48}^{+49}$ 764 ± 205
7510397 2.42 ± 0.15 3.62${}_{-0.04}^{+0.05}$ 3.52 3.40 1.90${}_{-0.31}^{+0.25}$ 4.18 ± 0.71 0.55${}_{-0.08}^{+0.11}$ 1663${}_{-270}^{+199}$ 1226${}_{-36}^{+30}$ 413 ± 176
7680114 1.16 ± 0.09 4.71${}_{-0.10}^{+0.10}$ 4.40 4.40 3.70${}_{-0.77}^{+0.66}$ 4.68 ± 1.95 0.24${}_{-0.11}^{+0.09}$ 2438${}_{-311}^{+215}$ 1736${}_{-34}^{+32}$ 620 ± 191
7771282 3.29 ± 0.28 4.20${}_{-0.21}^{+0.12}$ 4.14 3.96 0.71${}_{-0.71}^{+0.66}$ 3.29 ± 0.29
7871531 1.21 ± 0.12 1.91${}_{-0.07}^{+0.09}$ 1.71 1.77 1.94${}_{-1.53}^{+0.91}$ 4.83 ± 3.18 0.20${}_{-0.15}^{+0.09}$ 6026${}_{-1405}^{+967}$ 3161${}_{-178}^{+145}$ 1780 ± 653
7940546 2.92 ± 0.14 5.29${}_{-0.07}^{+0.07}$ 5.24 5.51 1.56${}_{-0.19}^{+0.19}$ 5.15 ± 0.69 0.53${}_{-0.06}^{+0.08}$ 1578${}_{-194}^{+137}$ 1177${}_{-23}^{+23}$ 409 ± 126
7970740 1.99 ± 0.15 1.65${}_{-0.03}^{+0.03}$ 1.35 1.51 4.25${}_{-0.72}^{+0.92}$ 6.85 ± 3.6 0.26${}_{-0.17}^{+0.10}$ 6875${}_{-1384}^{+1110}$ 3910${}_{-151}^{+177}$ 1949 ± 735
8006161 1.17 ± 0.06 1.94${}_{-0.04}^{+0.03}$ 1.69 1.73 4.75${}_{-0.62}^{+0.68}$ 6.0 ± 3.27 0.19${}_{-0.13}^{+0.06}$ 5581${}_{-861}^{+621}$ 3507${}_{-71}^{+73}$ 1575 ± 478
8150065 2.37 ± 0.3 3.55${}_{-0.17}^{+0.25}$ 3.38 3.23 2.01${}_{-1.53}^{+1.47}$ 2.36 ± 0.28
8179536 3.51 ± 0.26 3.42${}_{-0.07}^{+0.09}$ 3.35 3.25 1.88${}_{-0.89}^{+0.78}$ 7.05 ± 3.12 0.47${}_{-0.30}^{+0.14}$ 3940${}_{-535}^{+1129}$ 2112${}_{-194}^{+195}$ 1378 ± 486
8228742 2.05 ± 0.13 5.30${}_{-0.06}^{+0.08}$ 5.20 4.94 2.39${}_{-0.30}^{+0.26}$ 4.19 ± 0.68 0.48${}_{-0.07}^{+0.09}$ 1640${}_{-191}^{+127}$ 1189${}_{-26}^{+21}$ 383 ± 117
8379927 2.43 ± 0.11 2.21${}_{-0.03}^{+0.03}$ 2.09 2.05 2.45${}_{-0.27}^{+0.31}$ 6.12 ± 1.22 0.39${}_{-0.07}^{+0.08}$ 4091${}_{-435}^{+324}$ 2731${}_{-43}^{+49}$ 1046 ± 258
8394589 2.15 ± 0.16 3.49${}_{-0.05}^{+0.08}$ 3.33 3.24 1.71${}_{-0.60}^{+0.59}$ 5.79 ± 1.69 0.28${}_{-0.08}^{+0.09}$ 3185${}_{-356}^{+231}$ 2233${}_{-34}^{+41}$ 637 ± 208
8424992 1.14 ± 0.1 2.97${}_{-0.19}^{+0.18}$ 2.77 2.67 2.99${}_{-1.42}^{+1.46}$ 1.14 ± 0.1
8694723 3.08 ± 0.16 5.38${}_{-0.05}^{+0.09}$ 5.08 5.07 1.98${}_{-0.20}^{+0.21}$ 5.42 ± 0.79 0.56${}_{-0.07}^{+0.09}$ 2129${}_{-323}^{+220}$ 1472${}_{-34}^{+34}$ 548 ± 192
8760414 1.25 ± 0.08 3.96${}_{-0.06}^{+0.11}$ 3.49 3.52 2.97${}_{-0.61}^{+0.77}$ 5.06 ± 2.23 0.23${}_{-0.12}^{+0.09}$ 3659${}_{-486}^{+359}$ 2349${}_{-45}^{+52}$ 943 ± 277
8938364 0.8 ± 0.01 5.22${}_{-0.14}^{+0.12}$ 4.61 4.43 2.49${}_{-0.13}^{+0.10}$ 0.8 ± 0.01
9025370 1.39 ± 0.08 1.64${}_{-0.05}^{+0.04}$ 1.45 1.45 1.60${}_{-0.64}^{+0.62}$ 1.38 ± 0.09
9098294 1.3 ± 0.08 3.55${}_{-0.09}^{+0.09}$ 3.21 3.16 2.81${}_{-1.17}^{+1.09}$ 7.81 ± 4.74 0.16${}_{-0.10}^{+0.06}$ 3986${}_{-650}^{+517}$ 2317${}_{-95}^{+86}$ 1275 ± 355
9139151 1.99 ± 0.14 2.88${}_{-0.06}^{+0.05}$ 2.72 2.66 3.03${}_{-0.82}^{+0.71}$ 7.78 ± 3.84 0.25${}_{-0.15}^{+0.07}$ 4136${}_{-697}^{+494}$ 2697${}_{-70}^{+72}$ 1169 ± 399
9139163 5.28 ± 0.13 3.76${}_{-0.04}^{+0.05}$ 3.64 3.81 1.79${}_{-0.14}^{+0.13}$ 5.28 ± 0.13
9206432 5.87 ± 0.27 3.54${}_{-0.09}^{+0.12}$ 3.47 3.48 0.88${}_{-0.29}^{+0.31}$ 5.88 ± 0.28
9353712 3.22 ± 0.22 5.54${}_{-0.16}^{+0.18}$ 5.68 5.46 1.20${}_{-0.36}^{+0.38}$ 3.2 ± 0.21
9410862 1.89 ± 0.25 3.75${}_{-0.10}^{+0.13}$ 3.60 3.67 2.40${}_{-0.81}^{+0.82}$ 7.2 ± 3.8 0.24${}_{-0.14}^{+0.07}$ 3337${}_{-769}^{+460}$ 2204${}_{-63}^{+65}$ 846 ± 414
9414417 3.47 ± 0.23 5.55${}_{-0.10}^{+0.08}$ 5.18 5.29 1.04${}_{-0.35}^{+0.30}$ 5.81 ± 1.31 0.57${}_{-0.13}^{+0.14}$ 1910${}_{-552}^{+374}$ 1191${}_{-65}^{+56}$ 605 ± 299
9812850 5.68 ± 0.19 4.61${}_{-0.08}^{+0.11}$ 4.51 4.85 0.69${}_{-0.18}^{+0.19}$ 5.67 ± 0.2
9955598 0.77 ± 0.03 2.07${}_{-0.05}^{+0.06}$ 1.79 1.94 4.98${}_{-0.21}^{+0.26}$ 0.77 ± 0.01
9965715 3.14 ± 0.34 3.89${}_{-0.07}^{+0.07}$ 3.54 3.59 1.82${}_{-0.57}^{+0.47}$ 5.92 ± 1.9 0.47${}_{-0.16}^{+0.15}$ 3024${}_{-802}^{+687}$ 2103${}_{-84}^{+78}$ 792 ± 528
10068307 2.15 ± 0.13 5.98${}_{-0.09}^{+0.09}$ 5.86 5.54 2.15${}_{-0.18}^{+0.16}$ 3.91 ± 0.45 0.53${}_{-0.05}^{+0.07}$ 1388${}_{-165}^{+110}$ 1015${}_{-22}^{+19}$ 340 ± 102
10079226 2.05 ± 0.35 2.81${}_{-0.15}^{+0.14}$ 2.63 2.63 3.36${}_{-1.44}^{+1.49}$ 6.73 ± 4.24 0.27${}_{-0.19}^{+0.12}$ 4248${}_{-1099}^{+743}$ 2601${}_{-165}^{+156}$ 1236 ± 585
10162436 3.08 ± 0.08 5.38${}_{-0.09}^{+0.11}$ 5.40 5.42 1.82${}_{-0.15}^{+0.16}$ 3.08 ± 0.09
10454113 4.12 ± 0.15 3.19${}_{-0.04}^{+0.05}$ 2.87 3.01 1.60${}_{-0.26}^{+0.26}$ 4.13 ± 0.15
10516096 1.56 ± 0.1 4.74${}_{-0.11}^{+0.08}$ 4.40 4.26 3.44${}_{-0.54}^{+0.48}$ 5.33 ± 2.47 0.28${}_{-0.13}^{+0.12}$ 2629${}_{-539}^{+345}$ 1723${}_{-39}^{+38}$ 768 ± 298
10644253 2.21 ± 0.22 2.22${}_{-0.05}^{+0.06}$ 2.15 2.06 3.16${}_{-1.06}^{+0.99}$ 11.06 ± 4.51 0.19${}_{-0.09}^{+0.05}$ 4253${}_{-456}^{+364}$ 2852${}_{-69}^{+74}$ 1099 ± 282
10730618 5.07 ± 0.25 4.64${}_{-0.13}^{+0.14}$ 4.24 4.49 0.41${}_{-0.40}^{+0.39}$ 5.09 ± 0.27
10963065 2.16 ± 0.11 3.72${}_{-0.05}^{+0.09}$ 3.62 3.45 3.36${}_{-0.59}^{+0.54}$ 8.44 ± 3.9 0.25${}_{-0.14}^{+0.07}$ 3921${}_{-761}^{+568}$ 2158${}_{-71}^{+64}$ 1261 ± 383
11081729 6.19 ± 0.3 3.11${}_{-0.06}^{+0.07}$ 2.99 3.01 1.79${}_{-0.31}^{+0.30}$ 6.19 ± 0.3
11253226 5.8 ± 0.17 4.20${}_{-0.05}^{+0.07}$ 3.62 3.98 0.40${}_{-0.17}^{+0.17}$ 5.8 ± 0.17
11772920 0.77 ± 0.03 1.83${}_{-0.07}^{+0.07}$ 1.49 1.55 2.69${}_{-0.11}^{+0.12}$ 0.77 ± 0.01
12009504 2.38 ± 0.14 3.94${}_{-0.06}^{+0.08}$ 3.84 3.65 2.31${}_{-0.49}^{+0.46}$ 6.06 ± 2.12 0.38${}_{-0.13}^{+0.14}$ 2915${}_{-621}^{+420}$ 1821${}_{-57}^{+49}$ 819 ± 334
12069127 3.6 ± 0.23 6.00${}_{-0.17}^{+0.18}$ 6.00 5.79 0.93${}_{-0.36}^{+0.38}$ 3.61 ± 0.23
12069424 0.98 ± 0.04 3.96${}_{-0.06}^{+0.06}$ 3.55 3.56 3.33${}_{-0.29}^{+0.28}$ 5.47 ± 1.11 0.18${}_{-0.03}^{+0.04}$ 3187${}_{-171}^{+139}$ 2181${}_{-18}^{+20}$ 825 ± 108
12069449 0.91 ± 0.05 3.49${}_{-0.06}^{+0.07}$ 3.05 3.08 4.10${}_{-0.34}^{+0.33}$ 6.06 ± 1.56 0.15${}_{-0.04}^{+0.04}$ 3716${}_{-233}^{+185}$ 2579${}_{-22}^{+23}$ 966 ± 148
12258514 1.69 ± 0.09 4.74${}_{-0.06}^{+0.07}$ 4.52 4.35 2.92${}_{-0.18}^{+0.17}$ 4.01 ± 0.41 0.42${}_{-0.04}^{+0.04}$ 1939${}_{-116}^{+87}$ 1510${}_{-13}^{+12}$ 376 ± 79
12317678 5.83 ± 0.16 4.75${}_{-0.06}^{+0.08}$ 4.98 4.57 1.10${}_{-0.12}^{+0.14}$ 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 ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$ and ${W}_{\mathrm{dip}}$ are only given if the inclusion of the lorentzian in Equation (30) gave a better fit than the power-law-only fit. ${\mathrm{FWHM}}_{\mathrm{dip}}$ gives the full-width-half-maximum of the lorentzian line width dip.

Download table as:  ASCIITypeset images: 1 2

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 $4{\rm{\Delta }}\nu $, 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 ${\rm{\Delta }}\nu $/c and the square-root is taken to convert to amplitude. The multiplication of ${\rm{\Delta }}\nu $/c converts to power per radial mode, with c giving the effective number of radial modes per order, i.e., $c\approx {\tilde{V}}_{\mathrm{tot}}^{2}$ (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 $\lambda =650\,\mathrm{nm}$—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 $c={\tilde{V}}_{\mathrm{tot}}^{2}$.

Figure 15 shows the result from this exercise; we find that the ${A}_{\max }$ 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 $4{\rm{\Delta }}\nu $ 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 ${\tilde{V}}_{\mathrm{tot}}^{2}$ and thus c depends on various stellar parameters, such as ${T}_{\mathrm{eff}}$, $[\mathrm{Fe}/{\rm{H}}]$ , and log g, it should be expected that adopting a constant c adds to the scatter. The values of ${A}_{\max }$ from the smoothing method, with both values of c are given in Table 3.

Figure 15. Refer to the following caption and surrounding text.

Figure 15. Top: comparison between ${A}_{\max }$ obtained from the smoothed amplitude spectrum following the procedure by Kjeldsen et al. (2008a) with those obtained from the individual mode amplitudes. The different markers indicate the value of c used to convert the maximum smoothed amplitude to amplitude per radial degree mode. The dashed line gives the 1 : 1 relation. Bottom: relative offset between the two amplitude measures, with the zero offset given by the black dashed line. The blue dashed–dotted line indicates the median −6% offset between the amplitudes, which is found to be the same for both values of c.

Standard image High-resolution image

4.3. Mode Line Widths

Following Appourchaux et al. (2014), we adopt the following relation for the line widths Γ against mode frequency:

Equation (30)

where the first part describes a power-law trend with an exponent α and a value of ${{\rm{\Gamma }}}_{\alpha }$ at ${\nu }_{\max }.$ The second part describes a Lorentzian dip in the line widths in $\mathrm{ln}\nu $, centered on ${\nu }_{\mathrm{dip}}$ with a width ${W}_{\mathrm{dip}}$ and an depth ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$. For the fit, the values of $\mathrm{ln}({\rm{\Gamma }})$ were calculated from the Γ posterior distributions; the $\mathrm{ln}({\rm{\Gamma }})$-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; ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$ 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 ${{\rm{\Gamma }}}_{\alpha }$, ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$, and ${W}_{\mathrm{dip}}$ parameters. We note that the parameter ${W}_{\mathrm{dip}}$ has two solutions, one higher than ${\nu }_{\max },$ and one lower—in our fits, we chose the ${W}_{\mathrm{dip}}\gt {\nu }_{\max }$ solution with a prior on ${W}_{\mathrm{dip}}$. 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.

Figure 16. Refer to the following caption and surrounding text.

Figure 16. Radial mode line widths against a proxy for the radial order, with the color indicating the ${T}_{\mathrm{eff}}$. For a better visualization, the line widths have been smoothed with a five-point Epanechnikov filter.

Standard image High-resolution image

Figure 17 shows three examples where the full version of Equation (30) was deemed better—as seen, the line widths show a clear depression around ${\nu }_{\max }.$ 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.

Figure 17. Refer to the following caption and surrounding text.

Figure 17. Radial mode line widths as a function of frequency for KIC 12258514 (top), KIC 6116048 (middle), and KIC 12069449 (bottom). For all three stars, the line widths are best fit with the Lorentzian included in Equation (30), this fit is indicated by the full red lines; the dotted line gives the power-law component of the fit. The shaded dark and light blue regions indicate the 1 and 2σ credible regions of the fits. The dashed vertical lines give the ${\nu }_{\max }$ values.

Standard image High-resolution image

The parameters from the fits are given in Table 3 and shown in Figure 19 as a function of ${\nu }_{\max }$ and ${T}_{\mathrm{eff}}$. For the dependencies of the FWHM of the Lorentzian line width dip in frequency units, we used the transformation ${\mathrm{FWHM}}_{\mathrm{dip}}\ ={\nu }_{\mathrm{dip}}(| \sqrt{{{\rm{W}}}_{\mathrm{dip}}/{\nu }_{\max }}-\sqrt{{\nu }_{\max }/{{\rm{W}}}_{\mathrm{dip}}}| )$ (Appourchaux et al. 2016). The amplitude of the Lorentzian dip may be calculated using the transformation ${A}_{\mathrm{dip}}=\exp (| \mathrm{ln}{\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}| )$.20 Like Appourchaux et al. (2014) for their sample of 22 Kepler stars, we find a clear correlation between the Lorentzian width (${W}_{\mathrm{dip}}$ or ${\mathrm{FWHM}}_{\mathrm{dip}}$) and ${\nu }_{\max },$ while the power-law exponent (α) and the dip amplitude (${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$) are found to correlate most strongly with ${T}_{\mathrm{eff}}$. Only for the fits without the Lorentzian component does the value of ${{\rm{\Gamma }}}_{\alpha }$ correspond to the line width at ${\nu }_{\max },$ which is seen to correlate with ${T}_{\mathrm{eff}}$. We return to the overall behavior of the line width at ${\nu }_{\max }$ against ${T}_{\mathrm{eff}}$ below. Our estimates of ${\mathrm{FWHM}}_{\mathrm{dip}}$ agree well with updated values from Appourchaux et al. (2016). We make linear fits to the different parameters (P) against ${\nu }_{\max }$ as $P=a\,({\nu }_{\max }/3090)+b$ and ${T}_{\mathrm{eff}}$ as $P=a\,({T}_{\mathrm{eff}}\,/5777)+b$. Specifically, we perform a Deming regression (Deming 1943), where uncertainties in both the dependent (${\sigma }_{x}$) and independent variables (${\sigma }_{y}$) are considered via the following merit function (Press et al. 1993):

Equation (31)

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 ${{\rm{\Gamma }}}_{\alpha }$ against ${T}_{\mathrm{eff}}$ 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 ${\nu }_{\max }$ and ${T}_{\mathrm{eff}}$

  $P=a\,({\nu }_{\max }/3090)+b$ $P=a\,({T}_{\mathrm{eff}}\,/5777)+b$
$P\,(\mu \mathrm{Hz})$ b a b a
α 2.95 ± 0.16 0.39 ± 0.08 −25.5 ± 1.4 29.1 ± 1.5
${{\rm{\Gamma }}}_{\alpha }$ 3.08 ± 0.98 3.32 ± 0.50 6.3 ± 6.6 −1.8 ± 6.9
${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$ −0.47 ± 0.06 0.62 ± 0.04 3.5 ± 0.4 −3.3 ± 0.4
${W}_{\mathrm{dip}}$ 4637 ± 237 −141 ± 138 −28021 ± 2964 31971 ± 3105
${\nu }_{\mathrm{dip}}$ 2984 ± 31 60 ± 18 −23818 ± 1909 26785 ± 1983
${\mathrm{FWHM}}_{\mathrm{dip}}$ 1253 ± 162 −85 ± 96 −5649 ± 1093 6550 ± 1151

Note.  The linear fits are indicated as red lines in Figure 19. For the fits of ${{\rm{\Gamma }}}_{\alpha }$ 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 $\alpha \approx 3.6$—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 ${T}_{\mathrm{eff}}$. 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 ${T}_{\mathrm{eff}}$ —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).

Figure 18. Refer to the following caption and surrounding text.

Figure 18. Line widths (left axis) for KIC 6933899 (Fred) and model calculated damping rates (η) by G. Houdek et al. (2017, in preparation; right axis), multiplied by 2 to correspond to the mode FWHM. The full and dashed lines give the results for two model calculations of mode damping rates. Model 1 assumes a constant value for the velocity anisotropy Φ; Model 2 has a depth-dependent velocity anisotropy, guided by 3D simulation results from Trampedach et al. (2014). Stellar parameters from the best-fitting ASTFIT model (Paper II) were used in the damping-rate calculations. The vertical dotted line indicates the value of ${\nu }_{\max }.$

Standard image High-resolution image

The estimate of Γ at ${\nu }_{\max }$ was obtained by a Monte Carlo sampling from the posteriors of the fit parameters in Equation (30) and the estimate of ${\nu }_{\max }$ from the mode amplitudes. The values obtained are plotted against ${T}_{\mathrm{eff}}$ in Figure 20 and given in Table 3. We have fitted two relations for the line width at ${\nu }_{\max },$ namely, the power-law relation by Appourchaux et al. (2012a):

Equation (32)

and the exponential relation used, for instance, by Corsaro et al. (2012):

Equation (33)

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 ${T}_{\mathrm{eff}}$ for all stars, and 19 red giants by Corsaro et al. (2015), where we adopt ${T}_{\mathrm{eff}}$ 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 ${T}_{\mathrm{eff}}$. 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: ${{\rm{\Gamma }}}_{0}\approx 0.07\,\pm 0.03\,\mu \mathrm{Hz}$, $\alpha \approx 0.91\pm 0.13\,\mu \mathrm{Hz}$, and β ≈ 15.3 ± 1.9; for Equation (33) we obtain ${{\rm{\Gamma }}}_{0}\approx 1.02\pm 0.07\,\mu \mathrm{Hz}$ and ${T}_{0}\approx 436\pm 24\,{\rm{K}}$. If we fit only the red giants by a constant value for the line width, we find a value of ${{\rm{\Gamma }}}_{0}\approx 0.12\,\pm 0.01\,\mu \mathrm{Hz}$. 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 ${T}_{\mathrm{eff}}$ 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 ${\chi }^{2}$ 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 ${\nu }_{\mathrm{dip}}$ has a proportional bias with respect to ${\nu }_{\max },$ hence the higher the ${\nu }_{\max }$ the further below ${\nu }_{\max }$ does the dip in line widths appear. For Figure 20, the line widths were estimated at ${\nu }_{\max }$ and thus away from the line width dip for the highest ${\nu }_{\max }$ 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 ${T}_{\mathrm{eff}}\,\approx 5300$ K.

Figure 19. Refer to the following caption and surrounding text.

Figure 19. Parameters from the fit of Equation (30) to the radial mode line widths of the sample stars against ${\nu }_{\max }$(left panels) and ${T}_{\mathrm{eff}}$ (right panels). The color indicates the modeled mass. Stars fitted using only the first part of Equation (30), i.e., with parameters α and ${{\rm{\Gamma }}}_{\alpha }$, are given by circles (◦); stars where also the Lorentzian component was included, i.e., using also ${\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$, ${W}_{\mathrm{dip}}$, and ${\nu }_{\mathrm{dip}}$, are plotted with squares (□). The coefficient for the fitted linear relations (solid lines) are given in Table 4. The dashed line in the ${{\rm{\Gamma }}}_{\alpha }$ against ${T}_{\mathrm{eff}}$ panel gives the relation for Γ at ${\nu }_{\max }$ by Appourchaux et al. (2012a).

Standard image High-resolution image
Figure 20. Refer to the following caption and surrounding text.

Figure 20. Line width at ${\nu }_{\max }$ against ${T}_{\mathrm{eff}}$ with the color indicating the value of log g. See Table 3 for the plotted values, where the line widths are estimated from the fit of Equation (30). The red upward triangles indicate the line widths from the peak-bagging of 42 giants in NGC 6819 by R. Handberg et al. (2017, in preparation); the blue downward triangles give the line widths from Corsaro et al. (2015) for 19 red giants. The full red line gives the fit of Equation (32) with the shaded dark and light blue regions indicating the 1 and 2σ intervals of the fit; the dashed–dotted line gives the corresponding fit of Equation (32) by Appourchaux et al. (2012a; their Table 3); the dashed line gives the fit of Equation (33); the dotted line gives the corresponding fit of Equation (33) by Corsaro et al. (2012); the dashed–dotted–dotted line gives the constant fit to only the red giants. The markers indicate if the full fit of Equation (30) is preferred (□), or only the power-law component (◦).

Standard image High-resolution image

We 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 ${T}_{\mathrm{eff}}$ —Belkacem et al. (2012, 2013) specifically suggest a scaling with ${T}_{\mathrm{eff}}$ and log g. As noted by Appourchaux et al. (2012a), the uncertainty on ${T}_{\mathrm{eff}}$ currently sets the limitations for obtaining a more well constrained prediction for Γ at ${\nu }_{\max },$ 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 ${\nu }_{\max }$(Komm et al. 2000b; Chaplin et al. 2008a). For other analysis of the variation of Γ at ${\nu }_{\max },$ 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 ${T}_{\mathrm{eff}}$, ${\nu }_{\max },$ and mass. The total visibilities are computed from the combination of the MCMC chains for the individual visibilities. We see no strong correlation with ${T}_{\mathrm{eff}}$, $[\mathrm{Fe}/{\rm{H}}]$ (not shown), or $\mathrm{log}g$ (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 $\sim 1.4\,{M}_{\odot }$. None of the dependencies are, however, very strong, but they do appear similar in shape for different angular degrees.

Figure 21. Refer to the following caption and surrounding text.

Figure 21. Visibilities as a function of ${T}_{\mathrm{eff}}$ (left), ${\nu }_{\max }$(middle), and mass (right). The color indicates the metallicity; the markers indicate the angular degree (see the legend), and shown are also the total visibilities ${\tilde{V}}_{\mathrm{tot}}^{\,2}$—the bottom panels show the l = 3 visibilities on an expanded ordinate scale. The red dashed horizontal lines give the medians for the respective visibilities, with the values indicated in red on the ordinate; the shaded red regions give the standardized median-absolute-deviation (MAD) around the median values, given as 1.4826 times the MAD. The shaded gray regions in the left panel indicate the expected theoretical values from Ballot et al. (2011) for $\mathrm{log}g\,=4.0$ and $[\mathrm{Fe}/{\rm{H}}]$ in the range of ±1. The continuous black lines give the median binned values where the span of the parameters in the different panels have been divided into 10 bins.

Standard image High-resolution image

Table 5.  Extracted Mode Visibilities for Angular Degrees $l=1-3$

KIC ${\tilde{V}}_{1}^{2}$ ${\tilde{V}}_{2}^{2}$ ${\tilde{V}}_{3}^{2}$ ${\tilde{V}}_{\mathrm{tot}}^{2}$
1435467 ${1.52}_{-0.09}^{+0.07}$ ${0.58}_{-0.05}^{+0.09}$ ${3.08}_{-0.12}^{+0.17}$
2837475 ${1.46}_{-0.08}^{+0.06}$ ${0.53}_{-0.07}^{+0.05}$ ${3.00}_{-0.14}^{+0.10}$
3427720 ${1.75}_{-0.08}^{+0.11}$ ${0.74}_{-0.05}^{+0.05}$ ${3.50}_{-0.13}^{+0.13}$
3456181 ${1.74}_{-0.12}^{+0.12}$ ${0.80}_{-0.08}^{+0.13}$ ${3.56}_{-0.21}^{+0.23}$
3632418 ${1.65}_{-0.06}^{+0.06}$ ${0.73}_{-0.04}^{+0.07}$ ${3.41}_{-0.12}^{+0.10}$
3656476 ${1.35}_{-0.06}^{+0.09}$ ${0.62}_{-0.05}^{+0.04}$ ${0.07}_{-0.01}^{+0.02}$ ${3.08}_{-0.15}^{+0.08}$
3735871 ${1.55}_{-0.09}^{+0.11}$ ${0.69}_{-0.06}^{+0.06}$ ${3.21}_{-0.11}^{+0.17}$
4914923 ${1.43}_{-0.06}^{+0.10}$ ${0.58}_{-0.05}^{+0.02}$ ${0.08}_{-0.02}^{+0.01}$ ${3.11}_{-0.14}^{+0.10}$
5184732 ${1.57}_{-0.04}^{+0.05}$ ${0.68}_{-0.03}^{+0.03}$ ${0.08}_{-0.01}^{+0.01}$ ${3.32}_{-0.06}^{+0.09}$
5773345 ${1.27}_{-0.06}^{+0.09}$ ${0.38}_{-0.06}^{+0.07}$ ${2.65}_{-0.13}^{+0.13}$
5950854 ${1.32}_{-0.17}^{+0.17}$ ${0.62}_{-0.10}^{+0.11}$ ${2.92}_{-0.20}^{+0.29}$
6106415 ${1.50}_{-0.05}^{+0.03}$ ${0.63}_{-0.02}^{+0.03}$ ${0.10}_{-0.01}^{+0.01}$ ${3.21}_{-0.04}^{+0.08}$
6116048 ${1.45}_{-0.06}^{+0.04}$ ${0.62}_{-0.02}^{+0.03}$ ${0.09}_{-0.01}^{+0.01}$ ${3.14}_{-0.06}^{+0.09}$
6225718 ${1.54}_{-0.04}^{+0.03}$ ${0.62}_{-0.02}^{+0.02}$ ${0.07}_{-0.01}^{+0.01}$ ${3.24}_{-0.06}^{+0.05}$
6508366 ${1.80}_{-0.07}^{+0.04}$ ${0.99}_{-0.05}^{+0.01}$ ${3.77}_{-0.10}^{+0.07}$
6603624 ${1.33}_{-0.05}^{+0.06}$ ${0.59}_{-0.03}^{+0.04}$ ${0.06}_{-0.01}^{+0.01}$ ${2.98}_{-0.08}^{+0.08}$
6679371 ${1.44}_{-0.07}^{+0.10}$ ${0.64}_{-0.07}^{+0.08}$ ${3.07}_{-0.12}^{+0.19}$
6933899 ${1.55}_{-0.06}^{+0.07}$ ${0.64}_{-0.04}^{+0.03}$ ${3.20}_{-0.09}^{+0.10}$
7103006 ${1.27}_{-0.07}^{+0.09}$ ${0.37}_{-0.06}^{+0.08}$ ${2.66}_{-0.15}^{+0.14}$
7106245 ${1.48}_{-0.13}^{+0.12}$ ${0.66}_{-0.08}^{+0.09}$ ${3.14}_{-0.18}^{+0.19}$
7206837 ${1.25}_{-0.04}^{+0.09}$ ${0.29}_{-0.05}^{+0.05}$ ${2.58}_{-0.12}^{+0.08}$
7296438 ${1.56}_{-0.12}^{+0.10}$ ${0.67}_{-0.06}^{+0.06}$ ${3.24}_{-0.16}^{+0.15}$
7510397 ${1.58}_{-0.05}^{+0.11}$ ${0.64}_{-0.04}^{+0.06}$ ${3.23}_{-0.10}^{+0.15}$
7680114 ${1.56}_{-0.08}^{+0.07}$ ${0.65}_{-0.04}^{+0.04}$ ${0.09}_{-0.02}^{+0.03}$ ${3.32}_{-0.13}^{+0.09}$
7771282 ${1.48}_{-0.13}^{+0.15}$ ${0.98}_{-0.21}^{+0.02}$ ${3.35}_{-0.23}^{+0.24}$
7871531 ${1.45}_{-0.11}^{+0.10}$ ${0.40}_{-0.05}^{+0.05}$ ${2.86}_{-0.14}^{+0.11}$
7940546 ${1.63}_{-0.08}^{+0.05}$ ${0.74}_{-0.06}^{+0.06}$ ${3.33}_{-0.11}^{+0.15}$
7970740 ${1.11}_{-0.04}^{+0.05}$ ${0.32}_{-0.03}^{+0.03}$ ${2.42}_{-0.05}^{+0.09}$
8006161 ${1.37}_{-0.04}^{+0.05}$ ${0.45}_{-0.02}^{+0.02}$ ${0.07}_{-0.01}^{+0.01}$ ${2.88}_{-0.05}^{+0.07}$
8150065 ${1.65}_{-0.24}^{+0.17}$ ${0.72}_{-0.14}^{+0.11}$ ${3.38}_{-0.36}^{+0.24}$
8179536 ${1.54}_{-0.08}^{+0.10}$ ${0.68}_{-0.07}^{+0.06}$ ${3.23}_{-0.14}^{+0.14}$
8228742 ${1.63}_{-0.05}^{+0.06}$ ${0.72}_{-0.04}^{+0.04}$ ${3.35}_{-0.09}^{+0.10}$
8379927 ${1.52}_{-0.03}^{+0.04}$ ${0.65}_{-0.03}^{+0.02}$ ${3.17}_{-0.05}^{+0.05}$
8394589 ${1.57}_{-0.06}^{+0.08}$ ${0.62}_{-0.04}^{+0.04}$ ${3.18}_{-0.08}^{+0.12}$
8424992 ${1.57}_{-0.22}^{+0.17}$ ${0.70}_{-0.12}^{+0.11}$ ${3.22}_{-0.25}^{+0.29}$
8694723 ${1.48}_{-0.04}^{+0.06}$ ${0.57}_{-0.02}^{+0.05}$ ${3.07}_{-0.09}^{+0.08}$
8760414 ${1.38}_{-0.06}^{+0.05}$ ${0.55}_{-0.03}^{+0.03}$ ${0.08}_{-0.02}^{+0.01}$ ${3.01}_{-0.10}^{+0.08}$
8938364 ${1.57}_{-0.07}^{+0.09}$ ${0.64}_{-0.04}^{+0.04}$ ${0.05}_{-0.02}^{+0.02}$ ${3.29}_{-0.12}^{+0.09}$
9025370 ${1.41}_{-0.07}^{+0.11}$ ${0.61}_{-0.06}^{+0.05}$ ${3.03}_{-0.12}^{+0.12}$
9098294 ${1.51}_{-0.08}^{+0.07}$ ${0.62}_{-0.04}^{+0.04}$ ${3.12}_{-0.10}^{+0.11}$
9139151 ${1.60}_{-0.08}^{+0.05}$ ${0.59}_{-0.03}^{+0.04}$ ${3.18}_{-0.08}^{+0.10}$
9139163 ${1.46}_{-0.05}^{+0.06}$ ${0.52}_{-0.05}^{+0.06}$ ${2.99}_{-0.12}^{+0.10}$
9206432 ${1.44}_{-0.08}^{+0.12}$ ${0.56}_{-0.09}^{+0.10}$ ${2.99}_{-0.15}^{+0.24}$
9353712 ${1.62}_{-0.14}^{+0.16}$ ${0.74}_{-0.14}^{+0.15}$ ${3.37}_{-0.27}^{+0.29}$
9410862 ${1.34}_{-0.10}^{+0.09}$ ${0.57}_{-0.06}^{+0.07}$ ${2.95}_{-0.16}^{+0.11}$
9414417 ${1.36}_{-0.08}^{+0.09}$ ${0.52}_{-0.06}^{+0.08}$ ${2.85}_{-0.11}^{+0.19}$
9812850 ${1.42}_{-0.10}^{+0.09}$ ${0.58}_{-0.10}^{+0.10}$ ${2.99}_{-0.18}^{+0.20}$
9955598 ${1.20}_{-0.09}^{+0.06}$ ${0.41}_{-0.04}^{+0.04}$ ${2.58}_{-0.08}^{+0.11}$
9965715 ${1.38}_{-0.10}^{+0.06}$ ${0.59}_{-0.06}^{+0.05}$ ${2.97}_{-0.16}^{+0.11}$
10068307 ${1.68}_{-0.08}^{+0.08}$ ${0.72}_{-0.07}^{+0.05}$ ${3.40}_{-0.14}^{+0.13}$
10079226 ${1.52}_{-0.20}^{+0.15}$ ${0.52}_{-0.08}^{+0.10}$ ${3.01}_{-0.21}^{+0.25}$
10162436 ${1.69}_{-0.09}^{+0.09}$ ${0.73}_{-0.08}^{+0.08}$ ${3.39}_{-0.15}^{+0.18}$
10454113 ${1.21}_{-0.04}^{+0.05}$ ${0.49}_{-0.02}^{+0.03}$ ${0.06}_{-0.01}^{+0.02}$ ${2.77}_{-0.06}^{+0.07}$
10516096 ${1.58}_{-0.06}^{+0.06}$ ${0.66}_{-0.04}^{+0.04}$ ${3.23}_{-0.09}^{+0.09}$
10644253 ${1.59}_{-0.09}^{+0.11}$ ${0.73}_{-0.05}^{+0.06}$ ${3.32}_{-0.12}^{+0.15}$
10730618 ${1.29}_{-0.09}^{+0.13}$ ${0.38}_{-0.09}^{+0.11}$ ${2.68}_{-0.19}^{+0.21}$
10963065 ${1.62}_{-0.07}^{+0.06}$ ${0.72}_{-0.03}^{+0.04}$ ${3.31}_{-0.07}^{+0.11}$
11081729 ${1.43}_{-0.11}^{+0.04}$ ${0.60}_{-0.10}^{+0.05}$ ${2.98}_{-0.15}^{+0.13}$
11253226 ${1.18}_{-0.04}^{+0.06}$ ${0.31}_{-0.03}^{+0.06}$ ${2.48}_{-0.05}^{+0.12}$
11772920 ${1.42}_{-0.12}^{+0.14}$ ${0.39}_{-0.07}^{+0.06}$ ${2.80}_{-0.15}^{+0.18}$
12009504 ${1.63}_{-0.05}^{+0.07}$ ${0.72}_{-0.05}^{+0.04}$ ${3.34}_{-0.09}^{+0.11}$
12069127 ${1.53}_{-0.12}^{+0.14}$ ${0.69}_{-0.10}^{+0.16}$ ${3.27}_{-0.27}^{+0.24}$
12069424 ${1.38}_{-0.03}^{+0.05}$ ${0.59}_{-0.02}^{+0.02}$ ${0.08}_{-0.01}^{+0.00}$ ${3.07}_{-0.07}^{+0.05}$
12069449 ${1.37}_{-0.05}^{+0.05}$ ${0.53}_{-0.02}^{+0.03}$ ${0.07}_{-0.00}^{+0.01}$ ${2.98}_{-0.08}^{+0.07}$
12258514 ${1.60}_{-0.05}^{+0.04}$ ${0.68}_{-0.03}^{+0.03}$ ${3.29}_{-0.08}^{+0.06}$
12317678 ${1.77}_{-0.07}^{+0.12}$ ${0.81}_{-0.08}^{+0.10}$ ${3.53}_{-0.17}^{+0.20}$

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.

Download table as:  ASCIITypeset images: 1 2

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.

Figure 22. Refer to the following caption and surrounding text.

Figure 22. Box-plot of Spearman's rank correlation ρ between the parameters optimized in the peak-bagging and the mode visibilities. The box whiskers give the 15th and 85th percentiles.

Standard image High-resolution image

Considering 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 ${\tilde{V}}_{\mathrm{tot}}^{\,2}$ 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 $c={\tilde{V}}_{\mathrm{tot}}^{\,2}$, 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 ${T}_{\mathrm{eff}}$ 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 ${T}_{\mathrm{eff}}$ and/or $[\mathrm{Fe}/{\rm{H}}]$ , 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 $l\geqslant 3$ 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 $\delta {\nu }_{\mathrm{1,3}}$, 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 $\delta {\nu }_{\mathrm{1,3}}$ 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.

Figure 23. Refer to the following caption and surrounding text.

Figure 23. Collapsed échelle diagrams optimized for l = 3 using the method of Lund et al. (2014). The residual power spectra were used where the peak-bagging model has been divided out. Each horizontal line in the top part of the plot corresponds to a given star, ordered in terms of ${\nu }_{\max }.$ The scale goes white (low power) to black (high power) below the 90th percentile, and from blue (low power) to white (high power) above it. The lower panel shows the average collapsed échelle diagram from all stars. Vertical dashed lines have been added to guide the eye to the contribution from different angular degrees. For an optimization of the l = 3 signal, modes from $l\ne 3$ will not form a straight ridge in the échelle diagram and will consequently give a wide signal in the collapsed spectrum. An offset has been added to the frequencies for each of the collapsed spectra to align all l = 3 signals at 0.5.

Standard image High-resolution image

For 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 ${\nu }_{\max }$ 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).

Figure 24. Refer to the following caption and surrounding text.

Figure 24. Examples of collapsed échelle diagrams optimized for the detection of l = 4. Shown are the spectra for KIC 6603624 (left; Saxo), KIC 7510397 (middle), and KIC 12069449 (right; 16 Cyg B). These represent some of the cases with the strongest apparent l = 4 signals. An offset has been added to place the position of the expected l = 4 signal at 0.5.

Standard image High-resolution image

5. 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 epsilon to the expected value as a function of ${T}_{\mathrm{eff}}$ (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 ${r}_{\mathrm{01,10,02}}$ and second differences ${{\rm{\Delta }}}_{2}\nu (n,l)$ (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 $\mathrm{Frequency}$ $\mathrm{Amplitude}$ Line width ln K
    $(\mu \mathrm{Hz})$ $(\mathrm{ppm})$ $(\mu \mathrm{Hz})$  
11 1 1351.15${}_{-0.70}^{+0.59}$ 2.2
12 0 1407.23${}_{-1.18}^{+0.95}$ 0.79${}_{-0.13}^{+0.09}$ 2.50${}_{-1.34}^{+2.31}$ >6
12 1 1454.25${}_{-0.70}^{+0.53}$ >6
13 0 1510.10${}_{-0.48}^{+0.70}$ 0.99${}_{-0.15}^{+0.11}$ 2.50${}_{-0.95}^{+2.66}$ >6
13 1 1558.45${}_{-0.42}^{+0.54}$ >6
13 2 1605.68${}_{-0.81}^{+0.74}$ 1.45
14 0 1615.12${}_{-0.29}^{+0.24}$ 1.16${}_{-0.07}^{+0.07}$ 2.60${}_{-0.51}^{+0.67}$ >6
14 1 1664.09${}_{-0.23}^{+0.21}$ >6
14 2 1711.40${}_{-0.60}^{+0.50}$ 3.14
15 0 1720.35${}_{-0.17}^{+0.18}$ 1.46${}_{-0.07}^{+0.06}$ 2.30${}_{-0.37}^{+0.35}$ >6
15 1 1769.65${}_{-0.14}^{+0.15}$ >6
15 2 1816.19${}_{-0.36}^{+0.34}$ >6
16 0 1825.41${}_{-0.13}^{+0.12}$ 1.89${}_{-0.07}^{+0.06}$ 2.19${}_{-0.22}^{+0.34}$ >6
16 1 1873.88${}_{-0.14}^{+0.13}$ >6
16 2 1919.97${}_{-0.26}^{+0.26}$ >6
17 0 1929.05${}_{-0.14}^{+0.12}$ 2.29${}_{-0.06}^{+0.06}$ 2.84${}_{-0.20}^{+0.30}$ >6
17 1 1977.35${}_{-0.12}^{+0.11}$ >6
17 2 2023.80${}_{-0.21}^{+0.22}$ >6
18 0 2032.68${}_{-0.11}^{+0.11}$ 2.77${}_{-0.05}^{+0.07}$ 2.67${}_{-0.18}^{+0.23}$ >6
18 1 2081.57${}_{-0.09}^{+0.09}$ >6
18 2 2128.62${}_{-0.16}^{+0.15}$ >6
19 0 2137.59${}_{-0.09}^{+0.10}$ 3.14${}_{-0.06}^{+0.06}$ 2.50${}_{-0.15}^{+0.21}$ >6
19 1 2186.89${}_{-0.09}^{+0.08}$ >6
19 2 2234.70${}_{-0.16}^{+0.16}$ >6
20 0 2243.42${}_{-0.08}^{+0.08}$ 3.47${}_{-0.06}^{+0.07}$ 2.22${}_{-0.12}^{+0.16}$ >6
19 3 2281.61${}_{-3.36}^{+1.97}$ 3.01
20 1 2293.05${}_{-0.09}^{+0.09}$ >6
20 2 2340.63${}_{-0.16}^{+0.17}$ >6
21 0 2349.64${}_{-0.09}^{+0.08}$ 3.46${}_{-0.06}^{+0.07}$ 2.61${}_{-0.17}^{+0.16}$ >6
20 3 2385.57${}_{-1.16}^{+0.93}$ 3.94
21 1 2399.39${}_{-0.10}^{+0.08}$ >6
21 2 2446.71${}_{-0.15}^{+0.16}$ >6
22 0 2455.69${}_{-0.10}^{+0.11}$ 3.46${}_{-0.06}^{+0.07}$ 3.03${}_{-0.23}^{+0.19}$ >6
21 3 2493.08${}_{-1.64}^{+1.24}$ 3.66
22 1 2505.34${}_{-0.11}^{+0.10}$ >6
22 2 2552.85${}_{-0.24}^{+0.22}$ >6
23 0 2561.29${}_{-0.14}^{+0.15}$ 3.18${}_{-0.06}^{+0.06}$ 4.01${}_{-0.18}^{+0.25}$ >6
22 3 2598.55${}_{-1.66}^{+1.55}$ 3.12
23 1 2611.20${}_{-0.14}^{+0.13}$ >6
23 2 2658.63${}_{-0.32}^{+0.33}$ >6
24 0 2666.49${}_{-0.22}^{+0.24}$ 2.72${}_{-0.04}^{+0.06}$ 5.23${}_{-0.33}^{+0.26}$ >6
24 1 2717.47${}_{-0.17}^{+0.18}$ >6
24 2 2765.05${}_{-0.40}^{+0.40}$ >6
25 0 2773.06${}_{-0.31}^{+0.30}$ 2.41${}_{-0.06}^{+0.05}$ 6.75${}_{-0.51}^{+0.38}$ >6
25 1 2824.15${}_{-0.26}^{+0.27}$ >6
25 2 2872.28${}_{-0.54}^{+0.55}$ >6
26 0 2879.34${}_{-0.50}^{+0.56}$ 1.96${}_{-0.06}^{+0.06}$ 7.60${}_{-0.81}^{+0.63}$ >6
26 1 2931.24${}_{-0.34}^{+0.35}$ >6
26 2 2978.49${}_{-0.84}^{+0.80}$ 3.78
27 0 2987.15${}_{-0.49}^{+0.49}$ 1.70${}_{-0.06}^{+0.04}$ 7.53${}_{-0.74}^{+0.87}$ >6
27 1 3038.67${}_{-0.53}^{+0.51}$ >6
27 2 3084.55${}_{-1.59}^{+1.37}$ 1.45
28 0 3092.80${}_{-0.92}^{+0.88}$ 1.46${}_{-0.06}^{+0.07}$ 8.85${}_{-1.11}^{+1.44}$ >6
28 1 3145.65${}_{-0.62}^{+0.61}$ 4.92
28 2 3194.64${}_{-1.27}^{+1.68}$ 2.18
29 0 3204.41${}_{-0.93}^{+0.88}$ 0.98${}_{-0.10}^{+0.08}$ 5.87${}_{-1.41}^{+2.75}$ 3.81
29 1 3251.96${}_{-1.47}^{+1.71}$ 4.2
29 2 3302.59${}_{-2.57}^{+2.16}$ 1.22
30 0 3314.17${}_{-2.02}^{+2.07}$ 1.12${}_{-0.08}^{+0.07}$ 11.64${}_{-1.65}^{+1.55}$ 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 ${r}_{\mathrm{01,10,02}}(n)$ (Equation (16)) for KIC 6225718 (Saxo2)

$\mathrm{Ratio}\,\mathrm{type}$ n $\mathrm{Ratio}$
r02 14 ${0.0888}_{-0.0155}^{+0.0153}$
r02 15 ${0.0856}_{-0.0115}^{+0.0092}$
r02 16 ${0.0883}_{-0.0091}^{+0.0074}$
r02 17 ${0.0877}_{-0.0066}^{+0.0081}$
r02 18 ${0.0853}_{-0.0075}^{+0.0064}$
r02 19 ${0.0854}_{-0.0049}^{+0.0041}$
r02 20 ${0.0821}_{-0.0056}^{+0.0039}$
r02 21 ${0.0845}_{-0.0039}^{+0.0062}$
r02 22 ${0.0849}_{-0.0069}^{+0.0049}$
r02 23 ${0.0797}_{-0.0062}^{+0.0063}$
r02 24 ${0.0742}_{-0.0087}^{+0.0064}$
r02 25 ${0.0748}_{-0.0101}^{+0.0081}$
r02 26 ${0.0656}_{-0.0121}^{+0.0125}$
r02 27 ${0.0781}_{-0.0135}^{+0.0161}$
r02 28 ${0.0754}_{-0.0204}^{+0.0188}$
r02 29 ${0.0883}_{-0.0257}^{+0.0260}$
r10 12 ${0.0413}_{-0.0202}^{+0.0197}$
r01 13 ${0.0387}_{-0.0208}^{+0.0117}$
r10 13 ${0.0374}_{-0.0152}^{+0.0105}$
r01 14 ${0.0366}_{-0.0083}^{+0.0070}$
r10 14 ${0.0349}_{-0.0060}^{+0.0067}$
r01 15 ${0.0328}_{-0.0042}^{+0.0040}$
r10 15 ${0.0323}_{-0.0040}^{+0.0046}$
r01 16 ${0.0332}_{-0.0043}^{+0.0034}$
r10 16 ${0.0331}_{-0.0034}^{+0.0041}$
r01 17 ${0.0330}_{-0.0053}^{+0.0044}$
r10 17 ${0.0329}_{-0.0053}^{+0.0045}$
r01 18 ${0.0323}_{-0.0043}^{+0.0049}$
r10 18 ${0.0327}_{-0.0027}^{+0.0029}$
r01 19 ${0.0330}_{-0.0031}^{+0.0034}$
r10 19 ${0.0331}_{-0.0053}^{+0.0033}$
r01 20 ${0.0329}_{-0.0041}^{+0.0022}$
r10 20 ${0.0325}_{-0.0051}^{+0.0030}$
r01 21 ${0.0319}_{-0.0042}^{+0.0022}$
r10 21 ${0.0313}_{-0.0046}^{+0.0023}$
r01 22 ${0.0309}_{-0.0025}^{+0.0029}$
r10 22 ${0.0300}_{-0.0038}^{+0.0034}$
r01 23 ${0.0281}_{-0.0033}^{+0.0037}$
r10 23 ${0.0252}_{-0.0039}^{+0.0036}$
r01 24 ${0.0221}_{-0.0052}^{+0.0050}$
r10 24 ${0.0211}_{-0.0053}^{+0.0054}$
r01 25 ${0.0203}_{-0.0077}^{+0.0058}$
r10 25 ${0.0182}_{-0.0072}^{+0.0080}$
r01 26 ${0.0166}_{-0.0076}^{+0.0085}$
r10 26 ${0.0177}_{-0.0087}^{+0.0085}$
r01 27 ${0.0177}_{-0.0118}^{+0.0088}$
r10 27 ${0.0136}_{-0.0117}^{+0.0104}$
r01 28 ${0.0148}_{-0.0137}^{+0.0104}$
r10 28 ${0.0285}_{-0.0138}^{+0.0107}$
r01 29 ${0.0516}_{-0.0198}^{+0.0185}$

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 ${{\rm{\Delta }}}_{2}\nu (n,l)$ (Equation (20)) for KIC 6225718 (Saxo2)

n l ${{\rm{\Delta }}}_{2}\nu $
    $(\mu \mathrm{Hz})$
12 1 1.49${}_{-1.87}^{+1.43}$
13 0 1.64${}_{-1.65}^{+2.97}$
13 1 1.20${}_{-1.23}^{+1.49}$
14 0 0.37${}_{-1.09}^{+0.83}$
14 1 0.13${}_{-0.91}^{+0.54}$
14 2 −1.08${}_{-1.28}^{+1.55}$
15 0 −0.23${}_{-0.41}^{+0.52}$
15 1 −1.27${}_{-0.45}^{+0.34}$
15 2 −0.96${}_{-0.97}^{+0.92}$
16 1 −0.76${}_{-0.32}^{+0.35}$
16 2 0.14${}_{-0.78}^{+0.61}$
16 0 −1.35${}_{-0.40}^{+0.25}$
17 0 0.02${}_{-0.32}^{+0.27}$
17 2 0.97${}_{-0.51}^{+0.57}$
17 1 0.69${}_{-0.20}^{+0.38}$
18 2 1.25${}_{-0.41}^{+0.40}$
18 1 1.13${}_{-0.28}^{+0.21}$
18 0 1.30${}_{-0.29}^{+0.25}$
19 2 −0.15${}_{-0.42}^{+0.39}$
19 1 0.82${}_{-0.19}^{+0.26}$
19 0 0.89${}_{-0.21}^{+0.25}$
20 2 0.13${}_{-0.37}^{+0.41}$
20 3 −1.08${}_{-1.28}^{+1.55}$
20 1 0.15${}_{-0.18}^{+0.24}$
20 0 0.42${}_{-0.23}^{+0.19}$
21 3 −0.96${}_{-0.97}^{+0.92}$
21 1 −0.37${}_{-0.24}^{+0.20}$
21 2 0.16${}_{-0.48}^{+0.33}$
21 0 −0.16${}_{-0.23}^{+0.21}$
22 1 −0.09${}_{-0.27}^{+0.27}$
22 2 −0.42${}_{-0.50}^{+0.63}$
22 0 −0.45${}_{-0.26}^{+0.27}$
23 2 0.43${}_{-0.53}^{+1.06}$
23 0 −0.44${}_{-0.35}^{+0.42}$
23 1 0.38${}_{-0.30}^{+0.37}$
24 0 1.38${}_{-0.61}^{+0.54}$
24 2 0.72${}_{-0.96}^{+1.08}$
24 1 0.40${}_{-0.41}^{+0.47}$
25 0 −0.20${}_{-0.96}^{+0.73}$
25 2 −1.27${}_{-1.17}^{+1.56}$
25 1 0.47${}_{-0.69}^{+0.59}$
26 2 −0.15${}_{-1.96}^{+2.13}$
26 0 1.51${}_{-1.26}^{+1.17}$
26 1 0.36${}_{-0.90}^{+0.85}$
27 0 −1.98${}_{-1.70}^{+1.11}$
27 1 −0.35${}_{-1.34}^{+1.18}$
27 2 4.80${}_{-4.22}^{+2.77}$
28 2 −4.34${}_{-3.40}^{+6.80}$
28 1 −0.51${}_{-2.18}^{+1.68}$
28 0 6.76${}_{-2.89}^{+1.55}$
29 0 −1.46${}_{-2.99}^{+2.12}$

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 ${r}_{\mathrm{01,10,02}}$ as shown in the right panel of Figure 26; (4) a plot of the derived second differences ${{\rm{\Delta }}}_{2}\nu (n,l)$ 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.

Figure 25. Refer to the following caption and surrounding text.

Figure 25.

Example of the peak-bagging fit for KIC 6225718 (Saxo2). The power density spectrum is shown in gray with a 1 μHz Epanechnikov smoothed version overlain in black, and with the fitted power spectrum model given by the red curve. The markers indicate the frequency and angular degree of the fitted modes. (The complete figure set (66 images) is available.)

Standard image High-resolution image
Figure 26. Refer to the following caption and surrounding text.

Figure 26.

Left: example of échelle diagram for KIC 6225718 (Saxo2). The power spectrum used has been background corrected and smoothed by a 1 μHz Epanechnikov filter. The color scale indicates the power to background level, going from white at low level to black at high levels. The number on the right hand of the plot gives the radial order n of the l = 0 modes. Right: frequency difference ratios for KIC 6225718 as a function of the central frequency of the respective ratios (see Equation (16)). (The complete figure set (66 images) is available.)

Standard image High-resolution image
Figure 27. Refer to the following caption and surrounding text.

Figure 27.

Example of the extracted mode amplitudes and line widths for KIC 6225718 (Saxo2). Plotted in black are the amplitudes (left ordinate) and in red the line widths (right ordinate). The vertical dashed line gives the value for ${\nu }_{\max }.$ (The complete figure set (66 images) is available.)

Standard image High-resolution image

6. 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 ${\rm{\Delta }}\nu $ and ${\nu }_{\max }$ 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 ${\rm{\Delta }}\nu $ and ${\nu }_{\max }$ 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 ${\rm{\Delta }}\nu $ 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 ${\nu }_{\max }$ 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 ${\nu }_{\max }.$ We fitted this relation for all stars and were able to derive simple relations between the parameters of the fit and ${T}_{\mathrm{eff}}$ and ${\nu }_{\max }.$ 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 ${T}_{\mathrm{eff}}$ dependence of the line width at ${\nu }_{\max },$ 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 ${T}_{\mathrm{eff}}$.
  • 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 ${T}_{\mathrm{eff}}$, 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 ${T}_{\mathrm{eff}}$ values lower by about 1000–2000 K. Some structure was observed in the visibilities against ${T}_{\mathrm{eff}}$ 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 $x=\tfrac{\pi \nu }{2{\nu }_{\mathrm{nq}}}$, where ${\nu }_{\mathrm{nq}}$ is the Nyquist frequency—this is, however, an imprecise definition, because ${\nu }_{\mathrm{nq}}$ relates to the sampling time whereas the apodization relates to the integration time.

  • 19 

    Expressed in the units ${D}_{n0}^{(1)}=\delta {\nu }_{01}$ and ${D}_{n0}=\delta {\nu }_{02}/3$.

  • 20 

    Note that this ${A}_{\mathrm{dip}}$, where the Lorentzian in Equation (30) is added rather than subtracted, is different from that given by Appourchaux et al. (2014), who instead have ${A}_{\mathrm{dip}}={\rm{\Delta }}{{\rm{\Gamma }}}_{\mathrm{dip}}$.

  • 21 
Please wait… references are loading.
10.3847/1538-4357/835/2/172