<figure><img src="image_1.png" /></figure>

# Standing on the Shoulders of Dwarfs: the *Kepler* Asteroseismic LEGACY Sample.
## II. Radii, Masses, and Ages

Víctor Silva Aguirre<fnref n="1" />, Mikkel N. Lund<fnref n="1" />,<fnref n="2" />, H. M. Antia<fnref n="3" />, Warrick H. Ball<fnref n="4" />,<fnref n="5" />, Sarbani Basu<fnref n="6" />, Jørgen Christensen-Dalsgaard<fnref n="1" />,
Yveline Lebreton<fnref n="7" />,<fnref n="8" />, Daniel R. Reese<fnref n="1" />,<fnref n="2" />,<fnref n="9" />, Kuldeep Verma<fnref n="1" />,<fnref n="3" />, Luca Casagrande<fnref n="10" />, Anders B. Justesen<fnref n="1" />, Jakob R. Mosumgaard<fnref n="1" />,
William J. Chaplin<fnref n="1" />,<fnref n="2" />, Timothy R. Bedding<fnref n="1" />,<fnref n="11" />, Guy R. Davies<fnref n="1" />,<fnref n="2" />, Rasmus Handberg<fnref n="1" />, Günter Houdek<fnref n="1" />, Daniel Huber<fnref n="1" />,<fnref n="11" />,<fnref n="12" />,
Hans Kjeldsen<fnref n="1" />, David W. Latham<fnref n="13" />, Timothy R. White<fnref n="1" />, Hugo R. Coelho<fnref n="2" />, Andrea Miglio<fnref n="2" />, and Ben Rendle<fnref n="2" />

Received 2016 September 12; revised 2016 November 12; accepted 2016 November 19; published 2017 January 30

## Abstract

We use asteroseismic data from the *Kepler* satellite to determine fundamental stellar properties of the 66 main-sequence targets observed for at least one full year by the mission. We distributed tens of individual oscillation frequencies extracted from the time series of each star among seven modeling teams who applied different methods to determine radii, masses, and ages for all stars in the sample. Comparisons among the different results reveal a good level of agreement in all stellar properties, which is remarkable considering the variety of codes, input physics, and analysis methods employed by the different teams. Average uncertainties are of the order of ~2% in radius, ~4% in mass, and ~10% in age, making this the best-characterized sample of main-sequence stars available to date. Our predicted initial abundances and mixing-length parameters are checked against inferences from chemical enrichment laws $\Delta Y/\Delta Z$ and predictions from 3D atmospheric simulations. We test the accuracy of the determined stellar properties by comparing them to the Sun, angular diameter measurements, *Gaia* parallaxes, and binary evolution, finding excellent agreement in all cases and further confirming the robustness of asteroseismically determined physical parameters of stars when individual frequencies of oscillation are available. Baptised as the *Kepler* dwarfs LEGACY sample, these stars are the solar-like oscillators with the best asteroseismic properties available for at least another decade. All data used in this analysis and the resulting stellar parameters are made publicly available for the community.

**Key words:** asteroseismology – stars: fundamental parameters – stars: oscillations

**Supporting material:** machine-readable tables

## 1. Introduction

Asteroseismology has changed the way we determine properties of stars. The detection of oscillations at the surface is a window to their interiors because they allow us to determine their physical properties with a precision that is otherwise extremely hard to achieve in field stars. These properties, ages in particular, are crucial to our understanding of different fields of astrophysics (such as Galactic Archaeology, see, e.g., Miglio et al. 2013; Casagrande et al. 2016). Stellar age determinations, by traditional means, can be very uncertain (see, e.g., Soderblom 2010) and most age estimates, even when they claim to be precise, are not necessarily accurate.

The most widely used method to obtain ages of stars is to compare their observed atmospheric properties with those predicted by theoretical evolutionary sequences (e.g., Edvardsson et al. 1993; Ng & Bertelli 1998; Pont & Eyer 2004) or theoretical isochrones in the case of star clusters. This approach works reasonably well for clusters as long as the color–magnitude

diagram is well defined and contains stars in all stages of evolution, in particular, the main sequence, the main-sequence turnoff, and the red-giant branch. However, the technique of evolutionary-track or isochrone fitting does not work very well for single stars. More sophisticated statistical treatments, particularly Bayesian estimates (e.g., Jørgensen & Lindegren 2005; Takeda et al. 2007), are needed to obtain unbiased ages and, even then, uncertainties can be large. Estimating stellar ages also becomes increasingly problematic for stars that are more distant and hence fainter: reddening is uncertain, and the distances to the stars are not yet well known (though *Gaia* is drastically changing this picture).

The CoRoT and Kepler missions have completely changed the field of asteroseismology. Prior to these missions, solar-like oscillations for stars other than the Sun were observed in only a small number of them (see, e.g., Chaplin & Miglio 2013 for a review). Space-based photometry has now made asteroseismic analyses almost routine. The Kepler mission detected solar-like

<footnote n="1">¹Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark</footnote>

<footnote n="2">²School of Physics and Astronomy, University of Birmingham, Edgbaston Park Road, West Midlands, Birmingham, B15 2TT UK</footnote>

<footnote n="3">³Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India</footnote>

<footnote n="4">⁴Institut für Astrophysik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, D-37077, Göttingen, Germany</footnote>

<footnote n="5">⁵Max-Planck-Institut für Sonnensystemforschung, Justus-von-Liebig-Weg 3, D-37077 Göttingen, Germany</footnote>

<footnote n="6">⁶Department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101, USA</footnote>

<footnote n="7">⁷Observatoire de Paris, GEPI, CNRS UMR 8111, F-92195 Meudon, France</footnote>

<footnote n="8">⁸Institut de Physique de Rennes, Université de Rennes 1, CNRS UMR 6251, F-35042 Rennes, France</footnote>

<footnote n="9">⁹LESIA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06,
Univ. Paris Diderot, Sorbonne Paris Cité, F-92195 Meudon, France</footnote>

<footnote n="10">¹⁰Research School of Astronomy and Astrophysics, Mount Stromlo Observatory, The Australian National University, ACT 2611, Australia</footnote>

<footnote n="11">¹¹Sydney Institute for Astronomy (SIFA), School of Physics, University of Sydney, NSW 2006, Australia</footnote>

<footnote n="12">¹²SETI Institute, 189 Bernardo Avenue, Mountain View, CA 94043, USA</footnote>

<footnote n="13">¹³Harvard-Smithsonian Center for Astrophysics, 60 Garden Street Cambridge, MA 02138, USA</footnote>
oscillations in hundreds of main-sequence and subgiant stars, and the global asteroseismic parameters of these targets have been used to determine properties such as mass and radius (Chaplin et al. 2014). While the mass and radius estimates of these stars are reliable, their ages are less so due to the lack of sensitivity of these global properties to the deep stellar layers, where most of the evolutionary change takes place. To overcome this and to determine robust asteroseismic ages for these stars, it is necessary to extract and model the full spectrum of individual oscillation modes.

Silva Aguirre et al. (2015) and Davies et al. (2016) analyzed and modeled a sample of 33 exoplanet hosts for which the *Kepler* mission had obtained asteroseismic data, including individual oscillation frequencies. Those results have informed other studies of exoplanet systems (e.g., Van Eylen & Albrecht 2015; Bonfanti et al. 2016; Jontof-Hutter et al. 2016) and of the relation between stellar rotation and age (van Saders et al. 2016). In this paper, we perform a similar analysis to that of Silva Aguirre et al. (2015) for a set of 66 main-sequence stars observed by *Kepler* that are not known to be exoplanet hosts. Instead, the target stars are selected as those with the highest signal-to-noise ratios and the most precisely determined oscillation frequencies. In recent years, there have been similar studies using smaller samples and observations of shorter durations (e.g., Appourchaux et al. 2012; Mathur et al. 2012; Metcalfe et al. 2014). The data set used for this work has been extracted from the full available set of observations made by the *Kepler* mission, which for most stars means frequencies obtained from a four-year time series. Thus, the targets comprising this study, which we baptise as the *Kepler* dwarfs LEGACY sample, have the best asteroseismic data available among solar-like stars for at least another decade (i.e., until the PLATO 2.0 mission, Rauer et al. 2014).

The companion paper by Lund et al. (2017) describes the details of the data analysis and the frequencies that have been used in the modeling presented in this study. The rest of this paper is organized as follows. In Section 2, we describe the selection criteria for choosing the sample of stars that have been analyzed. We also discuss the ancillary data such as effective temperatures and metallicity needed for the analysis. The different techniques used in modeling the stars are described in Section 3. Section 4 is devoted to the results of the analysis, and tests of the derived properties using independent constraints are given in Section 5. We summarize the results and mention some of their applications in Section 6.

## 2. The LEGACY Sample

Our goal was to select main-sequence stars of the highest asteroseismic data quality present in the *Kepler* sample. A thorough description of the selection process and the main observational properties of our targets are given in the accompanying paper by Lund et al. (2017). We give a brief account of the most important points in this section.

The 66 stars comprising the LEGACY sample were chosen from more than 500 main-sequence and subgiant targets in which *Kepler* detected oscillations (Chaplin et al. 2014). We selected all targets that had more than one year of short-cadence observations, and where inspection of the power spectrum did not reveal any clear signature of bumped $\ell = 1$ modes. These are modes of mixed character that behave like pressure modes in the envelope and buoyancy modes in the core, and their appearance is related to the end of the core-hydrogen burning phase (Aizenman et al. 1977).

Figure 1 shows an asteroseismic Hertzsprung-Russell diagram (HRD), where the luminosity on the y-axis has been replaced by an asteroseismic property known as the average large frequency separation. This quantity decreases as stars evolve, since it is approximately proportional to the square root of the mean stellar density (Ulrich 1986; Gough 1987):

$$ \langle \Delta\nu \rangle \propto \left( \frac{M}{M_{\odot}} \right)^{1/2} \left( \frac{R}{R_{\odot}} \right)^{-3/2} . \qquad (1) $$

The large frequency separation in Figure 1 has been determined for the data using a Gaussian-weighted linear fit to the $\ell = 0$ frequencies as a function of radial order as described in Lund et al. (2017). Tracks at solar metallicity are shown in the figure as a reference, depicting the evolution of solar-type stars in this parameter space. Our sample predominantly comprises stars hotter and more evolved than the Sun, an expected bias since the amplitudes of oscillations scale approximately with stellar luminosity (Kjeldsen & Bedding 1995), making their detection easier in that region of the asteroseismic HRD. The atmospheric composition of our targets spans values from super-solar metallicity at $[\text{Fe/H}] = 0.35$ to metal-poor stars of $[\text{Fe/H}] = -1.00$, making this sample the most extended one in terms of heavy-element abundances, where a homogeneous asteroseismic analysis is feasible using individual frequencies of oscillation. A summary of the global asteroseismic quantities and atmospheric properties of the sample as compiled for this modeling effort is given in Appendix B.

## 3. Asteroseismic Modeling Pipelines

Seven independent pipelines were used to determine the stellar properties of the LEGACY sample. These pipelines use a variety of evolutionary and pulsation codes and consider different sets of input physics when computing theoretical evolutionary sequences and oscillation frequencies. Moreover, each pipeline uses its own preferred set of asteroseismic observables and statistical approach to analyze the goodness of fit of the models and extract the stellar properties of each target. In the following, we give a brief general description of the common characteristics in the process of estimating stellar properties from asteroseismic modeling, and give the specific details of each pipeline in an itemized list below.

First, one must generate a set of stellar models for different masses and chemical compositions using a stellar evolution code and set of input physics. Next, for many of these models, theoretical oscillation frequencies are computed using a standard (adiabatic) pulsation code. The actual number of evolutionary models and frequencies calculated is either defined in advance (when pipelines use precomputed grids of models) or determined during the fitting process (when pipelines chose masses and composition as part of the optimization). In both cases, the net result is stellar structures with atmospheric properties and theoretically computed frequencies that can be compared to those obtained from observations.

All pipelines use the atmospheric parameters $T_{\text{eff}}$ and $[\text{Fe/H}]$ as constraints but there are variations in the way asteroseismic data are used when determining stellar properties. One method is to fit the individual oscillation frequencies, $\nu_i(n)$, using a suitable prescription to correct for the surface effect (see, e.g., Kjeldsen et al. 2008; Ball & Gizon 2014; Sonoi et al. 2015, for a physical motivation and description of these corrections). Another method is to fit combinations of frequencies, normally
<figure><img src="image_2.png" /><figcaption>Figure 1. Position of the LEGACY sample (filled circles) and the Sun in an asteroseismic HRD, color coded according to their surface chemical composition. The typical error bar in effective temperature and large frequency separation is shown in black. Also depicted are representative stellar evolution tracks at solar metallicity spanning a mass range between 0.8 and 1.4 M$_{\odot}$ in steps of 0.1 M$_{\odot}$ (from Silva Aguirre et al. 2015). As a side note, we mention that the two lowest mass tracks do not evolve further into the red-giant phase because the computations were stopped at an age of 16 Gyr.</figcaption></figure>

the frequency separation ratios, defined as (Roxburgh & Vorontsov 2003)

$$r_{02}(n) = \frac{d_{02}(n)}{\Delta\nu_1(n)} \qquad (2)$$

$$r_{01}(n) = \frac{d_{01}(n)}{\Delta\nu_1(n)}, \quad r_{10}(n) = \frac{d_{10}(n)}{\Delta\nu_0(n+1)}. \qquad (3)$$

Here, $\Delta\nu_t(n) = \nu_{n,t} - \nu_{n-1,t}$ is the large separation between modes of the same angular degree and consecutive overtone, $d_{02}(n) = \nu_{n,0} - \nu_{n-1,0}$ is the small frequency separation, and $d_{01}(n)$ and $d_{10}(n)$ are the five-point small frequency separations:

$$d_{01}(n) = \frac{1}{8}(\nu_{n-1,0} - 4\nu_{n-1,1} + 6\nu_{n,0} - 4\nu_{n,1} + \nu_{n+1,0}) \qquad (4)$$

$$d_{10}(n) = -\frac{1}{8}(\nu_{n-1,1} - 4\nu_{n,0} + 6\nu_{n,1} - 4\nu_{n+1,0} + \nu_{n+1,1}). \qquad (5)$$

The frequency separation ratios have been shown to be mostly sensitive to the deep layers of the star (see, e.g., Otí Floranes et al. 2005; Roxburgh 2005; Silva Aguirre et al. 2011a), effectively diminishing the impact of poor modeling of the stellar outer layers in 1D evolutionary codes as well as the requirement of line-of-sight velocity corrections (see Davies et al. 2014). Due to the strong correlations between combinations including five individual frequencies, the ratios are customarily written and reproduced as one unique set of observables, called $r_{010}$ (see, e.g., Silva Aguirre et al. 2013):

$$r_{010} = \{r_{01}(n), r_{10}(n), r_{01}(n+1), r_{10}(n+1), ...\}. \qquad (6)$$

The pipelines participating in this study use individual frequencies, or combinations of frequencies, or both, following their own preference.

The process of determining stellar properties can be computationally expensive since some methods produce individual grids for each of the fitted targets that can require calculating thousands of models and frequencies for a given set of input physics. The inclusion of additional physical processes such as microscopic

diffusion or overshooting effectively increases the grid-dimensions to explore and can lead to a rapid increase in the number of models required to properly sample the parameter space. One alternative used by some pipelines to reduce the computational load is to use the average large frequency separation $\langle\Delta\nu\rangle$ and the frequency of maximum power in the oscillation spectrum $\nu_{max}$ as physically motivated proxies to constrain the parameter space that needs to be covered by the models. Since these quantities are sensitive to the mean stellar density and surface gravity, respectively (Ulrich 1986; Brown et al. 1991), they can restrict the mass and radius combinations that need to be explored.

All modeling teams were provided with the global asteroseismic parameters ($\langle\Delta\nu\rangle$ and $\nu_{max}$), the individual frequencies, frequency ratios, and correlations as determined by Lund et al. (2017), as well as the atmospheric properties $T_{eff}$ and [Fe/H] compiled from the literature. For most of our targets, the effective temperature and composition were taken from the spectroscopic analysis made by Buchhave & Latham (2015), and for the remainder we opted for other sources compiled from the literature. Table 3 reports the values given to the modeling teams, including the references for the atmospheric properties used.

In the following, we describe the main characteristics of the seven pipelines employed in the analysis and enclose a summary of them in Table 1.

**AIMS:** The “Asteroseismic Inference on a Massive Scale” pipeline is described in detail for the first time in this paper, and we give a description of it in Appendix A. Briefly, it uses the grid of MESA models from Coelho et al. (2015) and Delaunay tessellation for linear interpolation across the grid. It then applies an MCMC algorithm to find a representative set of models that satisfy the seismic and classic constraints. Stellar properties and uncertainties are obtained by averaging and calculating the standard deviations of the properties from the set of models.

**ASTFIT:** The “ASTEC FITting” method uses the ASTEC evolutionary code (Christensen-Dalsgaard 2008b) coupled with ADIPLS (Christensen-Dalsgaard 2008a) for computation of theoretical pulsation frequencies in a grid of stellar models. The initial chemical composition is determined by exploring a range of values for the galactic enrichment law, and it includes the effects of microscopic diffusion for low-mass stars and three values of the mixing-length efficiency $\alpha_{MLT} = 1.5, 1.8, 2.1$. The statistical methods applied to determine the stellar properties are described in Appendix A1 of Silva Aguirre et al. (2015).

**BASTA:** The “BAyesian STellar Algorithm” uses precomputed grids of evolutionary models and performs a global search for the optimal solution using the Bayesian approach described by Silva Aguirre et al. (2015). The current implementation considers microscopic diffusion for masses below $1.2 M_{\odot}$, the NACRE rates with the updated $^{14}$N(p, $\gamma$)$^{15}$O reaction from Formicola et al. (2004), and overshooting in the exponential decay formulation of Freytag et al. (1996) using the calibrated efficiency determined by Magic et al. (2010). The original grids have been extended in this work to include lower metallicities as required to match the spectroscopically determined surface composition of some targets.

**C2kSMO:** The “Cesam2k Stellar Model Optimization” pipeline (C2kSMO) uses the Cesam2k evolutionary code coupled to the Liège Oscillations Code and the procedures described by Lebreton & Goupil (2014). The input physics considered includes
<table><caption>**Table 1**
Summary of Codes, Input Physics, and Optimization Methods Applied by Each Pipeline</caption><thead><tr><th></th><th>AIMS</th><th>ASTFIT</th><th>BASTA</th><th>C2kSMO</th><th>GOE</th><th>V&A</th><th>YMCM</th></tr></thead><tbody><tr><td>Models</td><td>MESA</td><td>ASTEC</td><td>GARSTEC</td><td>Cesam2k</td><td>MESA</td><td>MESA</td><td>YREC</td></tr><tr><td>Frequencies</td><td>InversionKit</td><td>ADIPLS</td><td>ADIPLS</td><td>LOSC</td><td>ADIPLS</td><td>ADIPLS</td><td>AB94</td></tr><tr><td>Solar mixture</td><td>GN93</td><td>GN93</td><td>GS98</td><td>GN93</td><td>GS98</td><td>GS98</td><td>GS98</td></tr><tr><td>Opacities</td><td>OPAL96+JF05</td><td>OPAL96+JF05</td><td>OPAL96+JF05</td><td>OPAL96+JF05</td><td>OPAL96<br>+JF05</td><td>OP05+JF05</td><td>OPAL96<br>+JF05</td></tr><tr><td>EOS</td><td>OPAL05</td><td>OPAL05</td><td>OPAL05</td><td>OPAL05</td><td>OPAL05</td><td>OPAL05</td><td>OPAL05</td></tr><tr><td>Nuclear<br>reactions</td><td>NACRE</td><td>NACRE</td><td>NACRE</td><td>NACRE</td><td>NACRE</td><td>NACRE</td><td>Solar fusion</td></tr><tr><td>Atmosphere</td><td>Eddington gray</td><td>Eddington gray</td><td>Eddington gray</td><td>Eddington gray</td><td>Eddington gray</td><td>Eddington gray</td><td>Eddington<br>gray</td></tr><tr><td>Diffusion</td><td>No</td><td>MP93, &le;1.1 M<sub>&odot;</sub></td><td>T94, &le;1.2 M<sub>&odot;</sub></td><td>MP93</td><td>T94</td><td>T94, &le;1.35 M<sub>&odot;</sub></td><td>No</td></tr><tr><td>Overshoot</td><td>Yes</td><td>No</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td></tr><tr><td>Convection</td><td>CG68</td><td>MTL58</td><td>MLT12</td><td>CGM96</td><td>CG68</td><td>MLT58</td><td>MTL58</td></tr><tr><td>&alpha;<sub>conv</sub></td><td>1.8</td><td>1.5, 1.8, 2.1</td><td>1.791</td><td>Variable</td><td>Variable</td><td>Variable</td><td>Variable</td></tr><tr><td>&Delta;Y/&Delta;Z</td><td>2.0</td><td>1.0-2.0</td><td>1.4</td><td>Variable</td><td>Variable</td><td>Variable</td><td>Variable</td></tr><tr><td>Fitted data</td><td>&nu;<sub>i</sub>(n)</td><td>&nu;<sub>i</sub>(n)</td><td>r<sub>010</sub>, r<sub>02</sub>(n)</td><td>r<sub>010</sub>, r<sub>02</sub>(n),<br>&nu;<sub>0</sub>(n<sub>min</sub>)</td><td>&nu;<sub>i</sub>(n)</td><td>&lang;&Delta;&nu;&rang;, r<sub>02</sub>,<br>r<sub>01</sub>(n), r<sub>10</sub>(n + 3)</td><td>&nu;<sub>i</sub>(n),<br>r<sub>010</sub>, r<sub>02</sub>(n)</td></tr><tr><td>Surface<br>correction</td><td>BG14</td><td>SC</td><td>None</td><td>TS15</td><td>BG14</td><td>HK08</td><td>BG14</td></tr><tr><td>Optimization</td><td>MCMC</td><td>&chi;<sup>2</sup> minimization</td><td>Bayesian</td><td>Levenberg-<br>Marquardt</td><td>Downhill<br>simplex</td><td>&chi;<sup>2</sup> minimization</td><td>Monte Carlo</td></tr><tr><td>Reference</td><td>Appendix A</td><td>VSA15</td><td>VSA15</td><td>LG14</td><td>TA15</td><td>TA15</td><td>VSA15</td></tr></tbody></table>

**References.** Evolutionary models: ASTEC (Christensen-Dalsgaard 2008b), Cesam2k (Morel & Lebreton 2008), GARSTEC (Weiss & Schlattl 2008), MESA (Paxton et al. 2011), YREC (Demarque et al. 2008). Pulsation frequencies: AB94 (Antia & Basu 1994), ADIPLS Christensen-Dalsgaard (2008a), LOSC Scuflaire et al. (2008). Solar mixture: GN93 (Grevesse & Noels 1993), GS98 (Grevesse & Sauval 1998). Opacities: OPAL96 (Iglesias & Rogers 1996), OP05 (Badnell et al. 2005), JF05 (Ferguson et al. 2005). Equation of state (EOS): OPAL05 (Rogers & Nayfonov 2002). Nuclear reactions: Solar Fusion (Adelberger et al. 1998), NACRE (Angulo et al. 1999). Convection theory: MLT58 (Böhm-Vitense 1958), CG68 (Cox & Giuli 1968), MLT12 (Kippenhahn et al. 2012), CGM96 (Canuto et al. 1996). Diffusion: T94 (Thoul et al. 1994), MP93 (Michaud & Proffitt 1993). Surface corrections: HK08 (Kjeldsen et al. 2008), SC (solar scaled, see Appendix A1 in Silva Aguirre et al. 2015), BG14 (Ball & Gizon 2014), TS15 (Sonoi et al. 2015). Pipelines: VSA15 (Silva Aguirre et al. 2015), LG14 (Lebreton & Goupil 2014), TA15 (Appourchaux et al. 2015). See the text for details.

atomic diffusion, convective core overshooting with an efficiency of $d_{\text{ov}} = 0.15 \times \min(H_p, R_{\text{convcore}})$, convection being treated under the Canuto et al. (1996) formalism, and the NACRE rates with the updated $^{14}$N($p, \gamma$)$^{15}$O reaction from Imbriani et al. (2005). An initial approach to the solution is performed by finding the best fit to the observed values of $T_{\text{eff}}$, [Fe/H], and $\log(g)$ (the latter obtained from the global asteroseismic parameter $\nu_{\text{max}}$), using a chemical enrichment law $\Delta Y/\Delta Z$ from a solar calibration computed with the same input physics. The resulting age, mass, and initial chemical composition are used as starting values in a Levenberg-Marquardt minimization that adjusts free parameters in the modeling in order to minimize the merit function (see Miglio & Montalbán 2005). These free model parameters are the age and mass of the star, the initial helium content and initial ($Z/X)_0$ ratio, and the convective efficiency $\alpha_{\text{CGM}}$. The observables fitted by this optimization process are [Fe/H], $T_{\text{eff}}$, $\log(g)$, the lowest observed radial mode, and the frequency separation ratios $r_{010}$ and $r_{02}$. The error bars on the free parameters are obtained as the diagonal coefficients of the inverse of the Hessian matrix, while uncertainties in the other parameters are obtained following Equation (10) in Deheuvels et al. (2016).

**GOE:** The setting of the “GOEttingen” pipeline is very similar to that described by Appourchaux et al. (2015) and Reese et al. (2016). Since more massive stars are included in the LEGACY sample, the Ledoux (1947) criterion for convection was used with an extra free parameter for the convective overshooting described according to the exponential model of

Freytag et al. (1996). The NACRE thermonuclear reactions are used with the updated rates in $^{14}$N($p, \gamma$)$^{15}$O from Imbriani et al. (2005). The free adjustable parameters to find the best-fitting model were age, initial metallicity, initial helium, the mixing length, the overshooting efficiency, and the two terms in the Ball & Gizon (2014) surface correction.

**V&A:** This approach uses a combination of the MESA and ADIPLS codes to compute evolutionary models and pulsation frequencies with the same input physics as listed in Table 3 of Appourchaux et al. (2015). Diffusion of helium and heavy elements is typically not included for stars with masses greater than 1.35 $M_{\odot}$, while overshoot is taken into account with an exponentially decaying efficiency for stars in the relevant mass regime (above 1.10 $M_{\odot}$). For thermonuclear reactions the NACRE compilation is used with the exception of updated rates in the $^{14}$N($p, \gamma$)$^{15}$O reaction from Imbriani et al. (2005). The properties of stars are determined via $\chi^2$ minimization over a suitable grid of models. This is the only pipeline that uses the large frequency separation as an asteroseismic constraint, obtained by a linear least-squares fit to the observed frequencies as a function of radial order.

**YMCM:** The “Yale Monte Carlo Method” is applied as described by Silva Aguirre et al. (2015), with the exceptions of a variable value of the mixing-length parameter and the use of the Ball & Gizon (2014) formulation to correct for surface effects. The Solar fusion (Adelberger et al. 1998) cross-section for $^{14}$N($p, \gamma$)$^{15}$O has been updated according to Formicola et al. (2004).
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Distribution of uncertainties in density, radius, mass, and age for all stars in the sample as determined by each pipeline. When asymmetric error bars are reported, they are added in quadrature for the figure. Values for the C2kSMO density uncertainties are beyond the scale of the plot and therefore not shown. See the text for details.</figcaption></figure>

# 4. Results

## 4.1. Precision of the Stellar Properties

We determined stellar properties for the 66 stars in the LEGACY sample using the pipelines described in the previous section. Considering the variety of methods and input physics involved in the process, we choose to report the results from all of pipelines in complete and homogeneous sets given in Appendix C.

Figure 2 shows the uncertainty distributions obtained by each pipeline for the LEGACY sample. Overall, the results agree at the level of the size of their formal uncertainties, with some exceptions for given parameters. In terms of formal errors in density the BASTA results are of the order of ~2.5%, a factor of two larger than the majority of pipelines, while uncertainties from C2kSMO have median values at the 16% level and are much larger than those of any other pipeline. The reasons for these discrepancies are related to the details of each individual method: BASTA is the only pipeline to fit only frequency ratios, which are less sensitive to the mean density of the star than individual frequencies (see Silva Aguirre et al. 2015 for a discussion). On the other hand, the local minimization Levenberg-Marquardt algorithm used by C2kSMO does not determine uncertainties for density, so these are derived by

(uncorrelated) error propagation from mass and radius, which results in much higher formal errors than all other methods.

In Figure 2, the results for radius and mass also reveal some differences between the pipelines: C2kSMO has much larger radius uncertainties and formal mass errors only slightly smaller than the rest of the methods, resulting in the large density fractional uncertainties. Age distributions show some pipelines with fractional errors below 10% (C2kSMO and GOE), while the rest encompass values between 10% and 30%, which is in much better agreement with the level of age uncertainties determined from hare-and-hounds exercises (Reese et al. 2016) and those normally obtained from individual frequency fitting of Kepler main-sequence targets (e.g., Mathur et al. 2012; Silva Aguirre et al. 2013; Lebreton & Goupil 2014; Lund et al. 2014; Silva Aguirre et al. 2015).

The distributions in age show a tail toward large fractional error values (up to ~60%), much higher than the expected formal uncertainties from asteroseismic determinations. These are dominated by the youngest stars in our sample that have statistical errors of similar magnitude than the oldest stars. A more physically motivated scale to compare the age uncertainty distributions is given by the Terminal Age Main Sequence (TAMS), which we define as the point where the remaining central hydrogen content in a stellar model reaches $1 \times 10^{-5}$.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Age uncertainty distributions normalized to the TAMS age (solid) or the seismically determined age (dotted) of the target. See the text for details.</figcaption></figure>

Five pipelines were able to produce a proxy for the TAMS values of our targets by evolving their best-fitting model beyond the main-sequence phase. Figure 3 shows the uncertainty distributions normalized by current and TAMS age, where the tail at large fractional uncertainties is clearly reduced by the change in scale while the peak of the distribution remains very close to its original value. The latter is a consequence of an asteroseismic detection bias, which favors more evolved stars, as the oscillation amplitudes scale proportionally with stellar luminosity (e.g., Kjeldsen & Bedding 1995).

In summary, stellar properties of our targets are determined to varying degrees of precision by the pipelines, with median uncertainties ranging from 0.5% to 2.6% in density (excluding the C2kSMO results, see above); 1.3%-4.2% in radius; 2.3%-4.5% in mass; and 6.7%-20% in age. This level of precision is comparable to previous asteroseismic studies of smaller samples using individual frequencies of oscillations and makes this set the best-characterized main-sequence sample from the nominal *Kepler* mission.

## 4.2. Comparison Between Pipelines

The stellar properties predicted by each pipeline represent an estimate of the true physical parameters of our sample based on a given set of data. For the large majority of our targets, we lack an independent verification of these stellar properties that would allow us to test the validity of our results (but see Section 5 below for some exceptions). It is customary in these cases to compare the results across methods by selecting a reference pipeline and plot the fractional differences in a given quantity as a function of the result predicted by the reference. Since a priori we cannot consider any set of results to be an accurate representation of the underlying stellar properties of the sample, this approach can introduce spurious correlations misleading the analysis (see, e.g., Ludbrook 1997). An alternative is to use the Altman–Blandt method for analyzing differences (Bland & Altman 1986, 1995), consisting of calculating the actual difference between the results of a pipeline and the reference and plotting them against the corresponding mean values. Figures 4–7 depict these differences for density, radius, mass, and age using the BASTA results as a reference. In each panel, we include the results of computing the Pearson

product-moment correlation coefficient *r* of the mean values, as well as the *p*-value of a one-sample *t*-test of the weighted mean of the differences.

If there is no proportional bias between the results of a pipeline and the reference, the regression of the difference on means should have a slope of zero. For a sample of the size considered here, it translates into an absolute value of $|r| > 0.25$ for the bias to be significant at the 5% level. Except for the AIMS comparison, *r* values in Figure 4 reveal no proportional bias in density as a function of its mean value, a reassuring result considering that the large frequency separation is reproduced by all methods when fitting individual frequencies and it scales to good approximation with the square root of the mean stellar density (see Equation (1)).

In the comparison plots, we include dashed lines showing the mean and twice the standard deviation of the differences (weighted by the uncertainties). If there were no fixed biases between the results, then the mean of the differences should be zero. To test this, we performed a one-sample *t*-test of the weighted mean of the difference in each stellar property, where resulting *p*-values below 0.05 reject the $H_0$ hypothesis of zero mean differences and reveal a fixed bias between the methods. Thus, we find no such biases in the comparisons with ASTFIT, and YMCM, while the other four pipelines show a statistically significant (and positive) fixed bias.

To gain further understanding of these results, we followed up with a similar analysis of the results for radius and mass, shown in Figures 5 and 6. There is a clear proportional bias in radius between BASTA and ASTFIT (increasing), and between BASTA and V&A (decreasing). Considering that this correlation was not observed in any pipeline for the density results, it must be compensated in the mass determinations as confirmed by Figure 6, where the same two pipelines show a proportional bias according to their *r* values. On the other hand, the proportional bias in density between BASTA and AIMS is a result of a similar bias in radius, as the mass comparison between these pipelines shows no traces of such bias.

The results of the *t*-test for radius revealed a fixed bias for C2kSMO and YMCM, both cases predicting larger radii than BASTA. In the case of C2kSMO, a similar behavior was already present in the density: since this pipeline extracts this quantity from propagation of the mass and radius determinations, we expected one of these quantities to be the underlying effect causing the bias. The mass comparison does not show a fixed bias between C2kSMO and BASTA.

Masses determined with YMCM are the only ones showing a statistically significant fixed bias (higher values), while the ASTFIT and V&A results hint toward proportional biases in mass (see Figure 6). One would naively expect a similar behavior in the age determination, since mass is the main property controlling the main-sequence lifetime of a star. Figure 7 reveals that there are no correlations as a function of age for any of the pipeline results and also that the fixed bias between BASTA and YMCM is not significant. However, only C2kSMO and YMCM do not show a significant fixed bias, while all other pipelines do predict such a difference. The lack of correspondence between masses and ages is a result of the differences in evolutionary and pulsation codes employed by the methods, as well as the included input physics (see Table 1, also the discussion in Section 4 of Silva Aguirre et al. 2015).

In summary, the results from different pipelines show some fixed and proportional biases for given stellar properties, but
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Densities determined from each pipeline compared to those obtained with BASTA. The abscissa shows the average density value calculated between the pipeline and BASTA, while the ordinate gives the density differences in the sense (pipeline–BASTA). Solid lines mark the zero level to guide the eye, while dashed lines correspond to the mean and twice the standard deviation of the differences weighted by their uncertainties (added in quadrature). Pearson correlation coefficients *r*- and *p*-values are also given. See the text for details.</figcaption></figure>

overall there is an excellent level of agreement considering the variety of evolutionary and pulsation codes, input physics employed, as well as the differences in the seismic quantities fitted by each pipeline. This heterogeneity across methods as seen from Table 1 makes it difficult to isolate the unique physical culprit of a particular outlier in our sample. Thus, we have focused our comparison on general trends, with descriptions of individual cases when possible, and note that detailed discussions of the impact of changing the input physics on asteroseismically derived properties are given by Lebreton & Goupil (2014) and Silva Aguirre et al. (2015).

### 4.3. Initial Helium Abundance and Mixing-length Parameter

Four of the pipelines used in our analysis include initial abundances and efficiency of convection as free parameters that are determined during the optimization process. We can therefore compare the resulting values with expectations based on chemical enrichment of the Galaxy and predictions of 3D hydrodynamical simulation of atmospheres.

One of the most poorly constrained ingredients in modeling solar-type stars is the helium content because its abundance cannot be determined from spectroscopic observations. The transformation from the observed metallicity [Fe/H] to

fractional abundances X, Y, Z is achieved by assuming a relation between heavy elements and helium content known as the chemical enrichment law $\Delta Y/\Delta Z$. The existence and value of such a relation have long been debated, and current understanding favors the range between $1 \le \Delta Y/\Delta Z \le 3$ as obtained from different sets of indicators (see, e.g., Lebreton et al. 1999; Jimenez 2003; Balser 2006; Casagrande et al. 2007; Serenelli & Basu 2010, and references therein). The relation is anchored to a known point in the helium-heavy elements plane, normally the values obtained from Standard Big Bang Nucleosynthesis (SBBN) calculations $Z_0 = 0$ and $Y_0 = 0.248$ (Steigman 2010).

Figure 8 shows the initial helium and heavy-element abundance for all targets as predicted by the four pipelines freely varying the composition as part of the optimization. All panels show cases where $Y_{ini}$ is below the measurement from SBBN, a problem commonly appearing in asteroseismic results fitting individual oscillation frequencies (see, e.g., Mathur et al. 2012; Metcalfe et al. 2014; Silva Aguirre et al. 2015).

A linear fit to the results from each pipeline (weighted by the uncertainties when available) shows slopes compatible with the broad range of $1 \le \Delta Y/\Delta Z \le 3$ for the C2kSMO and GOE pipelines, while the fit to the V&A predictions returns a slightly lower slope. The YMCM abundances cluster systematically
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Same as Figure 4 but for radius determinations.</figcaption></figure>

below the $\Delta Y/\Delta Z = 1$ line and in combination with sub-SBBN helium for its most metal-poor stars predict the highest slope of all targets.

Despite the aforementioned differences, the overall asteroseismic predictions are broadly consistent with the slopes commonly adopted as galactic chemical enrichment laws. The recent measurement of surface helium abundance using asteroseismology in the binary 16 Cyg A&B (Verma et al. 2014b) and the detection of glitches associated with the He II ionization zone in all stars from our sample (K. Verma et al. 2016, in preparation) open the possibility for detailed studies that can potentially put firm constraints on $Y_{ini}$, and it will be pursued in a future study.

Another poorly constrained parameter that has a huge impact in the modeling results of solar-type stars is the efficiency of convection, parameterized in both the mixing length (Böhm-Vitense 1958) and the full spectrum of turbulent eddies (Canuto et al. 1996) theories. Normally, this efficiency is calibrated to the Sun and used across the HRD, but asteroseismic fitting allows us to obtain a determination of this parameter by including it in the optimization process. These predictions have not been tested to date due to the lack of benchmarks for comparison, but a new approach is now possible by extracting information about convective efficiency from large grids of 3D hydrodynamical simulations of stellar atmospheres. These simulations follow the properties of convection in the outer stellar layers without the need of free parameters and can be used to calibrate the expected

mixing-length parameter by matching a 1D envelope model to a suitable average of the 3D simulations (see Trampedach et al. 2014; Magic et al. 2015, for a description of the procedure).

In particular, the Stagger collaboration has produced a grid of more than 200 simulations covering a wide range in metallicity, effective temperature, and surface gravity. An equivalent mixing-length parameter based on these 3D simulations is available as a function of these properties (Magic et al. 2015), and we make a first comparison of their results to our determinations in Figure 9. For a given pipeline, we show, as a function of $T_{eff}$, the ratio between the mixing-length parameter of each star and the value for the Sun obtained from a solar calibration. In the 3D hydrodynamical case, the solar value is determined from an average solar simulation. The rationale behind presenting the comparison in terms of ratios in the convective efficiency is to minimize as much as possible the effects of the different input physics used in each pipeline and the atmospheric simulations. We caution that the comparison with C2kSMO is shown only for guidance because this pipelines uses the Canuto et al. (1996) theory for convection for which the efficiency value is not directly comparable to the mixing-length one.

None of the pipelines in our study follow the decreasing trend in convective efficiency as a function of $T_{eff}$ predicted by the simulations. In fact, two of the pipelines are consistent with a relatively flat relation (C2kSMO and V&A), while the GOE and YMCM pipelines seem to suggest the opposite behavior as the
<figure><img src="image_7.png" /><figcaption>Figure 6. Same as Figure 4 but for mass determinations.</figcaption></figure>

observed in the 3D case. It is worth noticing that V&A predicts mixing-length parameters systematically lower than the solar one. The results from 3D simulations in terms of an HRD-position-dependent convective efficiency are starting to be adopted by the stellar community in evolutionary calculations (see, e.g., Salaris & Cassisi 2015; Mosumgaard et al. 2016, for the first examples), and our sample of stars would be perfectly suited to test the impact of including results from 3D simulations on asteroseismically derived properties.

**5. Accuracy of Stellar Properties**

In the previous section, we compared the results from the different pipelines assuming that the real values of the predicted stellar properties were not known, and therefore each determination of a given quantity was treated as an independent method of measurement. While this is true for the large majority of our targets, in some cases additional constraints on the stellar properties are available from interferometry, parallaxes, or binarity. The following sections put our derived quantities under scrutiny using these additional data.

**5.1. The Sun**

Our own Sun is the star we have best characterized in the Universe and displays a rich spectrum of stochastically excited modes of oscillations. To test the accuracy of our pipelines, we

included in the target list circulated among all modelers a set of frequencies determined from solar data degraded in quality to match that of the *Kepler* mission (see Lund et al. 2016, submitted, for details). The true identity of this star was concealed from the modelers, and the results presented here include the same relevant input physics as given in Table 1 and fitting methods described in Section 3. Besides the frequencies, we provided a temperature of T<sub>eff</sub> = 5752 ± 77 K and a surface composition of [Fe/H] = −0.02 ± 0.15 dex determined from a random realization of Gaussian noise perturbations of the solar values and having uncertainties representative of the rest of the *Kepler* LEGACY sample. The results are given in Table 2 and are shown in Figure 10, showing that, in general, most methods return values compatible with the solar properties. However, some discrepant results appear that deserve further discussion.

The mass reported by C2kSMO is more than 3σ away from the Sun despite determining values comparable to other pipelines that agree with solar properties within one standard deviation. The reason for this is the smaller uncertainties in this parameter returned by this method, a feature that was already visible in relation to Figure 2. A similar issue is found in the derived mass and radius from AIMS, where the results lie around 2σ away from the true values. It is interesting to see the impact of diffusion on the solar age: while two methods (C2kSMO and YMCM) predict central masses higher than 1 M<sub>☉</sub>,
<figure><img src="image_8.png" /><figcaption>Figure 7. Same as Figure 4 but for age determinations.</figcaption></figure>

the age of the latter is not lower than the solar one due to the lack of microscopic diffusion in the modeling (see Table 1). A test including this effect in the YMCM modeling results in an estimated age of $\tau = 4.83$ Gyr, in better agreement with the solar value of 4.57 Gyr.

The solar properties determined by the V&A method favor a sub-solar mass and a surface helium abundance higher than that measured by helioseismology (see Table 2). This anti-correlation between mass and helium is expected from homology relations (e.g., Weiss et al. 2005) and has been reported before in asteroseismic modeling methods using local optimization (e.g., Mathur et al. 2012; Metcalfe et al. 2014, see also Section 4.3 in Lebreton & Goupil 2014 and Section 4.4 in Silva Aguirre et al. 2015).

According to homology relations, there is a strong dependence of the stellar luminosity on the mean molecular weight. Since the luminosity is related to the effective temperature and radius, once these are determined, there is a narrow range of combinations of mass and helium abundance that can comply with this constraint. The V&A results for the Sun reveal a bimodal distribution in mass, and the most probable of these peaks is reported in Table 2. However, the secondary solution is in better agreement with the solar properties, yielding a mass of $1.036 M_\odot$ and surface helium of $Y_{sup} = 0.237$. We note that when comparing the two solutions, the secondary one cannot be ruled out at the $2\sigma$ level, and analysis of the amplitude of

frequency glitches in the spectrum, as done by Verma et al. (2014a, 2014b), favors the high-mass results. This type of study can help breaking the degeneracy between mass and helium abundance and will be considered it a subsequent paper.

Finally, we turn to parameters that have been successfully determined from helioseismic inferences, namely the depth of the convective envelope and the helium surface abundance (see Christensen-Dalsgaard 2002 for a review). In both cases, all but the V&A results give agreement within $1\sigma$ to the helioseismic values, giving us confidence about the stellar parameters determined from these techniques for stars like the Sun.

## 5.2. Interferometry and Photometry

Independent radius determinations provide another check on the stellar properties derived from asteroseismology. Results from interferometric campaigns with CHARA combined with Hipparcos parallaxes have started to appear for targets where oscillations have been detected by Kepler, but currently amount to a handful in the main-sequence and red-giant phases (Huber et al. 2012; White et al. 2013; Johnson et al. 2014). In fact, only five of our targets have published interferometric radii that we can compare to the determinations from asteroseismology.

Figure 11 shows the results for the binary 16 Cyg A&B (White et al. 2013) and KIC 8006161 (Huber et al. 2012). For
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Initial helium abundance as a function of initial heavy-element abundance for the pipelines that do not constrain the chemical composition with a galactic enrichment law. Two dashed gray lines depict different slopes in the ΔY/ΔZ relation: 3.0 (upper), 1.0 (lower). The dotted gray line shows the primordial helium abundance predicted by SBBN (Steigman 2010), while the solid blue lines present a linear least-square fit to the results from each pipeline, characterized by the slope and intercept given in each panel. Median uncertainties are plotted to avoid clutter (black circles, except for V&A that does not provide them for these parameters). See the text for details.</figcaption></figure>

all three there is very good agreement between interferometric and asteroseismic radii, except for the *AIMS* results, which are more than $1\sigma$ away from the observations in two cases. For the other two stars in our sample that have interferometric radii, the comparison is shown in Figure 12. There is a clear difference between the radius from CHARA and the asteroseismically determined ones, the latter being systematically smaller. However, these stars have the smallest angular diameters measured in the Huber et al. (2012) sample (0.289 and 0.231 mas), making them more prone to systematic errors in the adopted calibrator diameters than for larger target stars, which are better resolved. In light of this, we restrict the comparison with interferometry to stars with angular diameters larger than >0.3 mas and also consider another independent radius determination coming from the InfraRed Flux Method (IRFM).

In the implementation considered here, the IRFM used several photometric bands to determine the bolometric flux of a star and its effective temperature in a self-consistent manner (see Casagrande et al. 2010, 2014a, and Section 5.3 below for details). A natural by-product of this procedure is the stellar angular diameter, which, combined with parallax measurements, yields stellar radius as depicted in Figure 12.

The agreement with asteroseismic radii is better than the comparison with interferometry, suggesting that the differences with the interferometric results indeed come from calibration problems. Additionally, the interferometric effective temperatures for these two targets are lower (by 1–2$\sigma$) than the IRFM and spectroscopic temperatures presented in this work and from Bruntt et al. (2012), translating into systematic differences of the order of ~3%. Only additional data from CHARA or another interferometer in the future will yield further information about these discrepancies, but we can conclude that the current available observations support the accuracy of asteroseismic radius determinations for these five targets. We note in passing that the IRFM results for KIC 8006161 and 16 Cyg A&B also agree well with the predictions from all pipelines.

## 5.3. Parallaxes

Using parallaxes to test the accuracy of asteroseismic results is a less direct method of comparison than, e.g., radius determination from interferometry. The measured parallaxes can accurately determine the distance to our targets, which in turn can be predicted from asteroseismic radii if an estimate of
<figure><img src="image_10.png" /><figcaption>**Figure 9.** Mixing-length parameter as a function of effective temperature for the LEGACY sample. Plotted on the Y-axis is the ratio between the α parameter of the targets and the solar mixing length determined for each pipeline (blue circles) or from 3D hydrodynamical simulations (gray circles). Dashed lines represent a linear fit weighted by the uncertainties when available. To avoid clutter, median uncertainties from each pipeline are plotted (black circles). The C2kSMO results are computed with the Canuto et al. (1996) formalism of convection. See the text for details.</figcaption></figure>

<table><caption>**Table 2**

Solar Properties Determined from Each Pipeline</caption><thead><tr><th></th><th>Mass (M<sub>☉</sub>)</th><th>Radius (R<sub>☉</sub>)</th><th>Age (Myr)</th><th>Luminosity (L<sub>☉</sub>)</th><th>Density (g cm<sup>-3</sup>)</th><th>Y<sub>sup</sub></th><th>R<sub>BCZ</sub> (R<sub>☉</sub>)</th></tr></thead><tbody><tr><td>AIMS</td><td>0.979 ± 0.013</td><td>0.992 ± 0.005</td><td>4840 ± 367</td><td>1.060 ± 0.048</td><td>1.413 ± 0.004</td><td>...</td><td>...</td></tr><tr><td>ASTFIT</td><td>0.986 ± 0.023</td><td>0.994 ± 0.008</td><td>4686 ± 393</td><td>0.972 ± 0.052</td><td>1.411 ± 0.003</td><td>0.249 ± 0.009</td><td>0.71 ± 0.007</td></tr><tr><td>BASTA</td><td>0.978<sup>+0.039</sup><sub>-0.030</sub></td><td>0.993<sup>+0.012</sup><sub>-0.012</sub></td><td>4852<sup>+1181</sup><sub>-1069</sub></td><td>0.976<sup>+0.054</sup><sub>-0.052</sub></td><td>1.411<sup>+0.021</sup><sub>-0.022</sub></td><td>0.247<sup>+0.012</sup><sub>-0.01</sub></td><td>0.713<sup>+0.009</sup><sub>-0.009</sub></td></tr><tr><td>C2kSMO</td><td>1.021 ± 0.003</td><td>1.006 ± 0.010</td><td>4331 ± 85.</td><td>1.084 ± 0.048</td><td>1.412 ± 0.048</td><td>0.245 ± 0.003</td><td>0.715 ± 0.004</td></tr><tr><td>GOE</td><td>0.997 ± 0.006</td><td>0.995 ± 0.018</td><td>4859 ± 128</td><td>0.947 ± 0.041</td><td>1.412 ± 0.002</td><td>0.234 ± 0.009</td><td>0.720 ± 0.003</td></tr><tr><td>V&A</td><td>0.927 ± 0.030</td><td>0.973 ± 0.015</td><td>4621 ± 200</td><td>0.937</td><td>1.418 ± 0.006</td><td>0.277</td><td>0.725</td></tr><tr><td>YMCM</td><td>1.037<sup>+0.031</sup><sub>-0.047</sub></td><td>1.012<sup>+0.005</sup><sub>-0.005</sub></td><td>5297<sup>+350</sup><sub>-350</sub></td><td>1.008<sup>+0.043</sup><sub>-0.042</sub></td><td>1.406<sup>+0.001</sup><sub>-0.001</sub></td><td>0.248<sup>+0.01</sup><sub>-0.01</sub></td><td>0.716<sup>+0.003</sup><sub>-0.003</sub></td></tr></tbody></table>

**Note.** See the text for details. The surface helium abundance ($Y_{sup}$) and radius of the base of the convective envelope ($R_{BCZ}$) are not determined by AIMS, while V&A reports only values from the best-fit model (and thus no uncertainties) for these parameters and the luminosity. For reference, the solar age is determined to be 4.57 Gyr (Bahcall et al. 1995), while helioseismology predicts $Y_{sup} = 0.248 \pm 0.003$ (Basu 1998) and $R_{BCZ} = 0.713 \pm 0.001 R_{☉}$ (Christensen-Dalsgaard et al. 1991).
(This table is available in its entirety in machine-readable form.)

the stellar angular diameter is available. Now that the Gaia mission has provided parallaxes for all stars in common with the Hipparcos and Tycho-2 catalogs (Lindegren et al. 2016), we can compare our distance determinations for 64 of the targets appearing in the Gaia data release 1 (DR1, Gaia Collaboration et al. 2016).

We base our astroseismic distance determinations on a method that couples radii to IRFM results in a self-consistent manner, as described by Silva Aguirre et al. (2011b, 2012). Briefly, the bolometric flux (and therefore angular diameter) was estimated from Tycho and 2MASS photometry using the spectroscopic metallicity of each star and a reddening value
<figure><img src="image_11.png" /><figcaption>**Figure 10.** Radius, mass, and age of the Sun determined by the pipelines for the concealed solar data set. See the text for details.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11.** Comparison of interferometric radii from Huber et al. (2012; KIC 8006161) and White et al. (2013; 16 Cyg A&B, KIC 12069424 and KIC 12069449) with those determined seismically by each pipeline. Gray shaded areas and horizontal dashed lines depict the standard deviation reported in the interferometric measurements.</figcaption></figure>

given by the distance dependent map of Amôres & Lépine (2005). We then interpolated using the asteroseismic gravity determined by each pipeline to extract a value of angular diameter and estimate its distance using the radius. We have assumed a uniform uncertainty in angular diameter of 3% (see Casagrande et al. 2014b, for further details on the implementation adopted). The resulting asteroseismic distances have median uncertainties of the order of 4.5%, while the median astrometric parallax uncertainty from the Gaia DR1 is close to 1%.

The resulting comparison is presented in Figure 13, where we plot the distribution of normalized differences between the asteroseismically determined distance transformed to parallax and the measurements from Gaia. There is an overall good level of agreement between all seismic pipelines, but an offset at the level of $1\sigma$ appears when comparing to the results of DR1 (with the astrometric parallax being systematically smaller than those predicted by asteroseismology). A similar result has

<figure><img src="image_13.png" /><figcaption>**Figure 12.** Same as Figure 11 showing the interferometric radius from Huber et al. (2012; gray shaded region and solid lines) and the radius determined using the IRFM formulation from Casagrande et al. (2014a; violet shaded region and dashed lines). See the text for details.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 13.** Comparison of parallaxes determined from asteroseismic distances and those measured by the Gaia satellite. Each panel shows the distribution of the difference between both determinations normalized by the individual uncertainties. Gray shaded area and vertical dashed lines mark the $1\sigma$ region. Left: using parallaxes directly from Gaia DR1. Right: applying the Stassun & Torres (2016) correction to the Gaia parallaxes. See the text for details.</figcaption></figure>

recently been found by Stassun & Torres (2016) who compared Gaia distances with independent measurements from eclipsing binaries, and found a constant offset of -0.25 mas in the sense of Gaia parallaxes being too small. The right panel in Figure 13 shows the comparison of our distances with those corrected by the constant offset found by Stassun & Torres (2016), where it is clear that the results agree much better if this systematic shift is applied. We note that 20 of our stars are also present in the Hipparcos catalog and had asteroseismic distances available (from Silva Aguirre et al. 2012), which were recently compared to Gaia results by De Ridder et al. (2016). The authors conclude that there is excellent agreement between the new astrometric and asteroseismic parallaxes, but there is a small offset pointing toward Gaia parallaxes being too small, though within the uncertainties in the determinations (similar to what is found here, see Figure 2 in De Ridder et al. 2016).
Although it is tempting to corroborate the conclusion of Gaia predicting too long distances drawn by Stassun & Torres (2016) from our sample, we note that there are other factors that could lie at the foundation of the systematic offset observed in the distances. For instance, our results and the properties of the eclipsing binaries used by Stassun & Torres (2016) depend on the adopted temperature scale. If our implementation of the IRFM were to return hotter temperatures and therefore smaller angular diameters, that would result in smaller radii and therefore larger distances and could potentially solve this discrepancy. Further analysis goes beyond the scope of this paper and will be presented for larger samples in upcoming studies (D. Huber et al. 2017, in preparation; C. Sahlholdt et al. 2017, in preparation). We also mention that an uncertainty of 0.3 mas is quoted in the Gaia DR1 as systematic, which is of the same order as the one found here (Gaia Collaboration et al. 2016).

Despite the overall good agreement between asteroseismic and astrometric parallaxes, some outliers remain beyond the 3σ level. Reasons that could explain this discrepancy are contaminated photometry compromising the IRFM determinations of angular diameters, inaccurate reddening estimates in the 3D dust maps, or problems with the astrometric solution used to derive parallaxes. The cases of KIC 8379927 and KIC 7510397 are examples of the latter because they are known binary stars with their fluxes affected by the presence of a companion. Similarly, the targets KIC 10454113 and KIC 9025370 have only recently been identified as spectroscopic binaries based on HARPS-N data (P. E. Nissen 2017, private communication). There are two other binary systems in our parallax sample, namely 16 Cyg A&B (KIC 12069424 and KIC 12069449) and the common proper motion pair KIC 9139151 and KIC 9139163 (Halbwachs 1986). However, the components of these systems have wide orbital separations and can be individually resolved in photometry, so there is no reason to doubt the reliability of their parallaxes and IRFM estimates of angular diameters.

There seems to be photometric contamination in KIC 1435467 (a rotationally variable star, see Ceillier et al. 2015) and KIC 7940546 (a magnetically active star presenting flares in its light curve, see Balona 2015). In the same vein, KIC 12317678 is an interesting case as it shows an infrared excess in its photometry attributed to circumstellar matter (McDonald et al. 2012) that could affect the fluxes used by the IRFM to extract the angular diameter. On the other hand, the parallax listed in Gaia DR1 is not consistent within the uncertainties with the value measured by the Hipparcos satellite. We expect that future data releases from Gaia will help clarify this situation. KIC 4914923 and KIC 9965715 also present deviations beyond the 3σ level that are currently attributed to problems in the reddening estimation. The remaining stars, which account for 55 out of the 64 targets with Gaia parallaxes, have no warning flag in their photometry and show an excellent agreement in distances, confirming the reliability of asteroseismically determined radii.

## 5.4. Binaries

One final consistency check on the seismic results can be extracted from information on known binaries in our sample. Two stars in our sample have a binary companion that has either no oscillations detected (KIC 8379927) or only global seismic parameters have been extracted in one of the

<figure><img src="image_15.png" /><figcaption>**Figure 14.** Ages from binary stars with both components included in our sample. Plotted in the x-axis are the age differences for both components normalized by the uncertainties, while the average age of the components is plotted in the y-axis. Gray shaded region and vertical dashed lines mark the 1σ level.</figcaption></figure>

components (KIC 7510397, Appourchaux et al. 2015). We are therefore left with two systems: KIC 12069424–12069449 (16 Cyg A&B) and KIC 9139151–9139163. Unfortunately, there are no dynamical masses derived for either of these pairs, but the fact that each system was presumably formed at the same time allows us to put strong constraints on the age of the systems.

None of the pipelines participating in this work forced the ages of these binaries to match, and Figure 14 shows the results for the two systems included in our sample. All methods return compatible results for 16 Cyg A&B, while GOE is the only outlier in the results for the other two stars (at the 1.5σ level). All considered, age determinations from all pipelines are consistent for the binaries in our sample and provide further support for the accuracy of asteroseismically determined ages from individual frequencies in main-sequence stars.

## 6. Conclusions and Outlook

We have presented an asteroseismic analysis for the *Kepler* dwarfs LEGACY sample, comprising main-sequence stars observed for more than a year for which individual oscillation frequencies have been detected. Using seven independent modeling pipelines, we have determined precise physical properties for the 66 stars in the set, with average uncertainties of the order of ~2% in radius, ~4% in mass, and ~10% in age. Comparison of the results across pipelines shows a good level of agreement in the derived properties, with differences at the level expected from the variety of codes and input physics that were used in the present analysis.

Our results were subject to tests using available independent measurements of stellar properties such as those of the Sun, radii from interferometry and the IRFM, and distances from parallaxes. We found that the asteroseismically inferred properties can reproduce the large majority of independent constraints within the uncertainties, while the deviating results can be understood in terms of stellar evolution relations (such as the mass-initial helium abundance anti-correlation) and different issues with the complementary data needed to estimate the desired property (such as contaminated photometric fluxes).
These comparisons serve as confirmation of the accuracy of
asteroseismic properties determined by fitting individual oscilla-
tions frequencies (or combinations thereof).

The stellar and oscillation mode properties presented here and in the accompanying paper by Lund et al. (2016, submitted) constitute the largest sample of asteroseismically analyzed main-sequence stars using individual oscillation frequencies. Adding the stars studied by Davies et al. (2016) and Silva Aguirre et al. (2015) gives 99 stars that comprise the best-characterized set of main-sequence solar-like oscillators observed by the *Kepler* mission. We have compared our determined initial abundances with expectations from galactic chemical enrichment, and our obtained efficiencies of convection with predictions from 3D hydrodynamical simulations. In both cases, we found interesting features that motivate further study, such as sub-SBBN initial helium abundances and the non-dependance of the convective efficiency on stellar effective temperature. We make these results available for the community and hope they will provide the basis of extensive work going beyond our analysis, such as studies of convective envelopes (Mazumdar et al. 2014), surface helium abundance (Verma et al. 2014b), core convection (e.g., Deheuvels et al. 2016), inversion of the stellar interior (e.g., Buldgen et al. 2015a, 2015b), amplitudes and line widths (e.g., G. Houdek et al. 2016, in preparation), average 3D simulations (e.g., Sonoi et al. 2015; Ball et al. 2016), and many others.

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 Jens
Jessen-Hansen for inspiration to the paper title. Funding for the
Stellar Astrophysics Centre is provided by The Danish National
Research Foundation (Grant DNRF106). The research was
supported by the ASTERISK project (ASTERoseismic Investiga-
tions with SONG and *Kepler*) funded by the European Research
Council (Grant agreement no.: 267864). V.S.A. and T.R.W.
acknowledge support from VILLUM FONDEN (research grant
10118). M.N.L. acknowledges the support of The Danish
Council for Independent Research—Natural Science (Grant
DFF-4181-00415). W.H.B. acknowledges research funding by
Deutsche Forschungsgemeinschaft (DFG) under grant SFB
963/1 “Astrophysical flow instabilities and turbulence”
(Project A18). S.B. is partially supported by NSF grant AST-
1514676 and NASA grant NNX16A109G. D.R.R. was funded
by the European Community’s Seventh Framework Pro-
gramme (FP7/2007-2013) under grant agreement No. 312844
(SPACEINN), thereby allowing him to develop the AIMS
pipeline. W.J.C. and G.R.D. acknowledge the support of the
UK Science and Technology Facilities Council (STFC). 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. This work
has made use of data from the European Space Agency (ESA)
mission *Gaia* (http://www.cosmos.esa.int/gaia), processed by
the *Gaia* Data Processing and Analysis Consortium (DPAC,
http://www.cosmos.esa.int/web/gaia/dpac/consortium).
Funding for the DPAC has been provided by national
institutions, in particular, the institutions participating in the
*Gaia* Multilateral Agreement.

# Appendix A
## The AIMS Pipeline

The “Asteroseimic Inference on a Massive Scale” (AIMS)
pipeline was recently developed at the University of Birming-
ham in the context of the SpaceInn network. This pipeline
relies on various basic components:

1. a precalculated grid of models. In the present paper, the grid of main-sequence models from Coelho et al. (2015) was used. This grid of models was calculated via the MESA (Paxton et al. 2013, 2015) stellar evolution code and ranges from 0.8 $M_\odot$ to 1.5 $M_\odot$ in steps of 0.01 $M_\odot$, and from -0.8 to 0.8 dex in steps of 0.1 dex for [Fe/H]. The Grevesse & Noels (1993) value of $(Z/X)_\odot = 0.0245$ was used to convert [Fe/H] to $Z/X$. The enrichment law used $\Delta Y/\Delta Z = 2$ (e.g., Chiosi & Matteucci 1982) and a primordial helium abundance $Y_p = 0.2484$ (Cyburt et al. 2003). Convection was based on the standard mixing-length theory, using a solar-calibrated parameter. No diffusion or rotational mixing was included, while overshoot was taken into account as an adiabatic extension of the mixed region with an efficiency of $d_{\text{ov}} = 0.20 \times H_p$. The NACRE thermonuclear reactions are used with the updated rate in $^{14}$N($p, \gamma$)$^{15}$O from Imbriani et al. (2005).

2. the MCMC algorithm from Foreman-Mackey et al. (2013). This provides a way of approximating the probability distribution function in stellar parameter space that results from the provided seismic and classic constraints.

3. linear interpolation within the grid of models using a Delaunay tessellation in stellar parameter space (except for the time dimension, which is dealt with separately). Hence, AIMS can explore points in parameter space which lie between the stellar models from the grid. The Delaunay tessellation is calculated thanks to the qhull<sup>14</sup> library (via python's numpy library). The advantage of using a Delaunay tessellation is that the grid of models does not necessarily have to be structured.

The AIMS pipeline is freely available<sup>15</sup> along with a more
detailed description than what is provided here.

Appendix B
Modeling Input Parameters

In order to determine the stellar properties of our targets,
each modeling team used the individual frequencies and
combinations determined by Lund et al. (2017) in the manner
described in Section 3. The necessary complementary atmo-
spheric parameters $T_{\text{eff}}$ and [Fe/H] have been compiled from
different sources in the literature and listed in Table 3, together
with the global asteroseismic quantities.

Simple estimates of $\langle \Delta\nu \rangle$ were obtained using the method of White et al. (2011), that is, from a weighted fit of the asymptotic function ($\nu_0(n) \approx (n + 1/2 + \epsilon)\Delta\nu$) to the radial mode frequencies from the peak-bagging as a function of radial order $n$. The weights were given by a Gaussian with an FWHM of $0.25\nu_{\text{max}}$. Estimates from a more elaborate method including deviations from the simple asymptotic relation used here can be found in Lund et al. (2016, submitted). The estimate of $\nu_{\text{max}}$

<footnote n="14"><sup>14</sup> See http://www.qhull.org.</footnote>

<footnote n="15">¹⁵ See http://bison.ph.bham.ac.uk/spaceinn/aims.</footnote>
<table><caption>**Table 3**
Global Asteroseismic and Atmospheric Properties of Our Sample</caption><thead><tr><th>KIC</th><th>ν<sub>max</sub> (μHz)</th><th>⟨Δν⟩ (μHz)</th><th>T<sub>eff</sub> (K)</th><th>[Fe/H] (dex)</th><th>Reference</th></tr></thead><tbody><tr><td>1435467</td><td>1406.7<sup>+6.3</sup><sub>-8.4</sub></td><td>70.369<sup>+0.034</sup><sub>-0.033</sub></td><td>6326 ± 77</td><td>0.01 ± 0.1</td><td>1</td></tr><tr><td>2837475</td><td>1557.6<sup>+8.2</sup><sub>-9.2</sub></td><td>75.729<sup>+0.041</sup><sub>-0.042</sub></td><td>6614 ± 77</td><td>0.01 ± 0.1</td><td>1</td></tr><tr><td>3427720</td><td>2737.0<sup>+10.7</sup><sub>-17.7</sub></td><td>120.068<sup>+0.031</sup><sub>-0.032</sub></td><td>6045 ± 77</td><td>-0.06 ± 0.1</td><td>1</td></tr><tr><td>3456181</td><td>970.0<sup>+8.3</sup><sub>-5.9</sub></td><td>52.264<sup>+0.041</sup><sub>-0.039</sub></td><td>6384 ± 77</td><td>-0.15 ± 0.1</td><td>1</td></tr><tr><td>3632418</td><td>1166.8<sup>+3.0</sup><sub>-3.8</sub></td><td>60.704<sup>+0.019</sup><sub>-0.018</sub></td><td>6193 ± 77</td><td>-0.12 ± 0.1</td><td>1</td></tr><tr><td>3656476</td><td>1925.0<sup>+7.0</sup><sub>-6.3</sub></td><td>93.194<sup>+0.018</sup><sub>-0.020</sub></td><td>5668 ± 77</td><td>0.25 ± 0.1</td><td>1</td></tr><tr><td>3735871</td><td>2862.6<sup>+16.6</sup><sub>-26.5</sub></td><td>123.049<sup>+0.047</sup><sub>-0.046</sub></td><td>6107 ± 77</td><td>-0.04 ± 0.1</td><td>1</td></tr><tr><td>4914923</td><td>1817.0<sup>+6.3</sup><sub>-5.2</sub></td><td>88.531<sup>+0.019</sup><sub>-0.019</sub></td><td>5805 ± 77</td><td>0.08 ± 0.1</td><td>1</td></tr><tr><td>5184732</td><td>2089.3<sup>+4.4</sup><sub>-4.1</sub></td><td>95.545<sup>+0.024</sup><sub>-0.023</sub></td><td>5846 ± 77</td><td>0.36 ± 0.1</td><td>1</td></tr><tr><td>5773345</td><td>1101.2<sup>+5.7</sup><sub>-6.6</sub></td><td>57.303<sup>+0.030</sup><sub>-0.027</sub></td><td>6130 ± 84</td><td>0.21 ± 0.09</td><td>6</td></tr><tr><td>5950854</td><td>1926.7<sup>+21.9</sup><sub>-20.4</sub></td><td>96.629<sup>+0.102</sup><sub>-0.107</sub></td><td>5853 ± 77</td><td>-0.23 ± 0.1</td><td>1</td></tr><tr><td>6106415</td><td>2248.6<sup>+4.6</sup><sub>-3.9</sub></td><td>104.074<sup>+0.023</sup><sub>-0.026</sub></td><td>6037 ± 77</td><td>-0.04 ± 0.1</td><td>1</td></tr><tr><td>6116048</td><td>2126.9<sup>+5.5</sup><sub>-5.0</sub></td><td>100.754<sup>+0.017</sup><sub>-0.017</sub></td><td>6033 ± 77</td><td>-0.23 ± 0.1</td><td>1</td></tr><tr><td>6225718</td><td>2364.2<sup>+4.9</sup><sub>-4.6</sub></td><td>105.695<sup>+0.018</sup><sub>-0.017</sub></td><td>6313 ± 76</td><td>-0.07 ± 0.1</td><td>1</td></tr><tr><td>6508366</td><td>958.3<sup>+4.6</sup><sub>-3.6</sub></td><td>51.553<sup>+0.046</sup><sub>-0.047</sub></td><td>6331 ± 77</td><td>-0.05 ± 0.1</td><td>1</td></tr><tr><td>6603624</td><td>2384.0<sup>+5.4</sup><sub>-5.6</sub></td><td>110.128<sup>+0.012</sup><sub>-0.012</sub></td><td>5674 ± 77</td><td>0.28 ± 0.1</td><td>1</td></tr><tr><td>6679371</td><td>941.8<sup>+5.1</sup><sub>-5.0</sub></td><td>50.601<sup>+0.029</sup><sub>-0.029</sub></td><td>6479 ± 77</td><td>0.01 ± 0.1</td><td>1</td></tr><tr><td>6933899</td><td>1389.9<sup>+3.9</sup><sub>-3.6</sub></td><td>72.135<sup>+0.018</sup><sub>-0.018</sub></td><td>5832 ± 77</td><td>-0.01 ± 0.1</td><td>1</td></tr><tr><td>7103006</td><td>1167.9<sup>+7.2</sup><sub>-6.9</sub></td><td>59.658<sup>+0.029</sup><sub>-0.030</sub></td><td>6344 ± 77</td><td>0.02 ± 0.1</td><td>1</td></tr><tr><td>7106245</td><td>2397.9<sup>+24.0</sup><sub>-28.7</sub></td><td>111.376<sup>+0.063</sup><sub>-0.061</sub></td><td>6068 ± 102</td><td>-0.99 ± 0.19</td><td>4</td></tr><tr><td>7206837</td><td>1652.5<sup>+10.6</sup><sub>-11.7</sub></td><td>79.131<sup>+0.037</sup><sub>-0.039</sub></td><td>6305 ± 77</td><td>0.10 ± 0.1</td><td>1</td></tr><tr><td>7296438</td><td>1847.8<sup>+8.5</sup><sub>-12.6</sub></td><td>88.698<sup>+0.040</sup><sub>-0.036</sub></td><td>5775 ± 77</td><td>0.19 ± 0.1</td><td>1</td></tr><tr><td>7510397</td><td>1189.1<sup>+3.4</sup><sub>-4.4</sub></td><td>62.249<sup>+0.020</sup><sub>-0.020</sub></td><td>6171 ± 77</td><td>-0.21 ± 0.1</td><td>1</td></tr><tr><td>7680114</td><td>1709.1<sup>+7.1</sup><sub>-6.5</sub></td><td>85.145<sup>+0.039</sup><sub>-0.043</sub></td><td>5811 ± 77</td><td>0.05 ± 0.1</td><td>1</td></tr><tr><td>7771282</td><td>1465.1<sup>+27.0</sup><sub>-18.7</sub></td><td>72.463<sup>+0.069</sup><sub>-0.079</sub></td><td>6248 ± 77</td><td>-0.02 ± 0.1</td><td>1</td></tr><tr><td>7871531</td><td>3455.9<sup>+19.3</sup><sub>-26.5</sub></td><td>151.329<sup>+0.025</sup><sub>-0.023</sub></td><td>5501 ± 77</td><td>-0.26 ± 0.1</td><td>1</td></tr><tr><td>7940546</td><td>1116.6<sup>+3.3</sup><sub>-3.6</sub></td><td>58.762<sup>+0.029</sup><sub>-0.029</sub></td><td>6235 ± 77</td><td>-0.20 ± 0.1</td><td>1</td></tr><tr><td>7970740</td><td>4197.4<sup>+21.2</sup><sub>-18.4</sub></td><td>173.541<sup>+0.060</sup><sub>-0.068</sub></td><td>5309 ± 77</td><td>-0.54 ± 0.1</td><td>1</td></tr><tr><td>8006161</td><td>3574.7<sup>+11.4</sup><sub>-10.5</sub></td><td>149.427<sup>+0.015</sup><sub>-0.014</sub></td><td>5488 ± 77</td><td>0.34 ± 0.1</td><td>1</td></tr><tr><td>8150065</td><td>1876.9<sup>+38.1</sup><sub>-32.4</sub></td><td>89.264<sup>+0.134</sup><sub>-0.121</sub></td><td>6173 ± 101</td><td>-0.13 ± 0.15</td><td>4</td></tr><tr><td>8179536</td><td>2074.9<sup>+13.8</sup><sub>-12.0</sub></td><td>95.090<sup>+0.058</sup><sub>-0.054</sub></td><td>6343 ± 77</td><td>-0.03 ± 0.1</td><td>1</td></tr><tr><td>8228742</td><td>1190.5<sup>+3.4</sup><sub>-3.7</sub></td><td>62.071<sup>+0.022</sup><sub>-0.021</sub></td><td>6122 ± 77</td><td>-0.08 ± 0.1</td><td>1</td></tr><tr><td>8379927</td><td>2795.3<sup>+6.0</sup><sub>-5.7</sub></td><td>120.288<sup>+0.017</sup><sub>-0.018</sub></td><td>6067 ± 120</td><td>-0.10 ± 0.15</td><td>2</td></tr><tr><td>8394589</td><td>2396.7<sup>+10.5</sup><sub>-9.4</sub></td><td>109.488<sup>+0.034</sup><sub>-0.035</sub></td><td>6143 ± 77</td><td>-0.29 ± 0.1</td><td>1</td></tr><tr><td>8424992</td><td>2533.7<sup>+27.0</sup><sub>-28.1</sub></td><td>120.584<sup>+0.062</sup><sub>-0.064</sub></td><td>5719 ± 77</td><td>-0.12 ± 0.1</td><td>1</td></tr><tr><td>8694723</td><td>1470.5<sup>+3.7</sup><sub>-4.1</sub></td><td>75.112<sup>+0.019</sup><sub>-0.021</sub></td><td>6246 ± 77</td><td>-0.42 ± 0.1</td><td>1</td></tr><tr><td>8760414</td><td>2455.3<sup>+9.1</sup><sub>-8.3</sub></td><td>117.230<sup>+0.022</sup><sub>-0.018</sub></td><td>5873 ± 77</td><td>-0.92 ± 0.1</td><td>1</td></tr><tr><td>8938364</td><td>1675.1<sup>+5.2</sup><sub>-5.8</sub></td><td>85.684<sup>+0.018</sup><sub>-0.020</sub></td><td>5677 ± 77</td><td>-0.13 ± 0.1</td><td>1</td></tr><tr><td>9025370</td><td>2988.6<sup>+20.0</sup><sub>-16.9</sub></td><td>132.628<sup>+0.030</sup><sub>-0.024</sub></td><td>5270 ± 180</td><td>-0.12 ± 0.18</td><td>3</td></tr><tr><td>9098294</td><td>2314.7<sup>+9.2</sup><sub>-10.4</sub></td><td>108.894<sup>+0.023</sup><sub>-0.022</sub></td><td>5852 ± 77</td><td>-0.18 ± 0.1</td><td>1</td></tr><tr><td>9139151</td><td>2690.4<sup>+14.5</sup><sub>-9.0</sub></td><td>117.294<sup>+0.031</sup><sub>-0.032</sub></td><td>6302 ± 77</td><td>0.10 ± 0.1</td><td>1</td></tr><tr><td>9139163</td><td>1729.8<sup>+6.2</sup><sub>-5.9</sub></td><td>81.170<sup>+0.042</sup><sub>-0.036</sub></td><td>6400 ± 84</td><td>0.15 ± 0.09</td><td>6</td></tr><tr><td>9206432</td><td>1866.4<sup>+10.3</sup><sub>-14.9</sub></td><td>84.926<sup>+0.046</sup><sub>-0.051</sub></td><td>6538 ± 77</td><td>0.16 ± 0.1</td><td>1</td></tr><tr><td>9353712</td><td>934.3<sup>+11.1</sup><sub>-8.3</sub></td><td>51.467<sup>+0.091</sup><sub>-0.104</sub></td><td>6278 ± 77</td><td>-0.05 ± 0.1</td><td>1</td></tr><tr><td>9410862</td><td>2278.8<sup>+31.2</sup><sub>-16.6</sub></td><td>107.390<sup>+0.050</sup><sub>-0.053</sub></td><td>6047 ± 77</td><td>-0.31 ± 0.1</td><td>1</td></tr><tr><td>9414417</td><td>1155.3<sup>+6.1</sup><sub>-4.6</sub></td><td>60.115<sup>+0.024</sup><sub>-0.024</sub></td><td>6253 ± 75</td><td>-0.13 ± 0.1</td><td>7</td></tr><tr><td>9812850</td><td>1255.2<sup>+9.1</sup><sub>-7.0</sub></td><td>64.746<sup>+0.067</sup><sub>-0.068</sub></td><td>6321 ± 77</td><td>-0.07 ± 0.1</td><td>1</td></tr><tr><td>9955598</td><td>3616.8<sup>+21.2</sup><sub>-29.6</sub></td><td>153.283<sup>+0.029</sup><sub>-0.032</sub></td><td>5457 ± 77</td><td>0.05 ± 0.1</td><td>1</td></tr><tr><td>9965715</td><td>2079.3<sup>+9.2</sup><sub>-10.4</sub></td><td>97.236<sup>+0.041</sup><sub>-0.042</sub></td><td>5860 ± 180</td><td>-0.44 ± 0.18</td><td>3</td></tr><tr><td>10068307</td><td>995.1<sup>+2.8</sup><sub>-2.7</sub></td><td>53.945<sup>+0.019</sup><sub>-0.020</sub></td><td>6132 ± 77</td><td>-0.23 ± 0.1</td><td>1</td></tr><tr><td>10079226</td><td>2653.0<sup>+47.7</sup><sub>-44.3</sub></td><td>116.345<sup>+0.059</sup><sub>-0.052</sub></td><td>5949 ± 77</td><td>0.11 ± 0.1</td><td>1</td></tr><tr><td>10162436</td><td>1052.0<sup>+4.0</sup><sub>-4.2</sub></td><td>55.725<sup>+0.035</sup><sub>-0.039</sub></td><td>6146 ± 77</td><td>-0.16 ± 0.1</td><td>1</td></tr><tr><td>10454113</td><td>2357.2<sup>+8.2</sup><sub>-9.1</sub></td><td>105.063<sup>+0.031</sup><sub>-0.033</sub></td><td>6177 ± 77</td><td>-0.07 ± 0.1</td><td>1</td></tr><tr><td>10516096</td><td>1689.8<sup>+4.6</sup><sub>-5.8</sub></td><td>84.424<sup>+0.022</sup><sub>-0.025</sub></td><td>5964 ± 77</td><td>-0.11 ± 0.1</td><td>1</td></tr><tr><td>10644253</td><td>2899.7<sup>+21.3</sup><sub>-22.8</sub></td><td>123.080<sup>+0.056</sup><sub>-0.055</sub></td><td>6045 ± 77</td><td>0.06 ± 0.1</td><td>1</td></tr><tr><td>10730618</td><td>1282.1<sup>+14.6</sup><sub>-12.7</sub></td><td>66.333<sup>+0.061</sup><sub>-0.064</sub></td><td>6150 ± 180</td><td>-0.11 ± 0.18</td><td>3</td></tr><tr><td>10963065</td><td>2203.7<sup>+6.7</sup><sub>-6.3</sub></td><td>103.179<sup>+0.027</sup><sub>-0.027</sub></td><td>6140 ± 77</td><td>-0.19 ± 0.1</td><td>1</td></tr><tr><td>11081729</td><td>1968.3<sup>+11.0</sup><sub>-12.6</sub></td><td>90.116<sup>+0.048</sup><sub>-0.047</sub></td><td>6548 ± 82</td><td>0.11 ± 0.1</td><td>1</td></tr></tbody></table>
<table><caption>**Table 3**
(Continued)</caption><thead><tr><th>KIC</th><th>ν<sub>max</sub> (μHz)</th><th>⟨Δν⟩ (μHz)</th><th>T<sub>eff</sub> (K)</th><th>[Fe/H] (dex)</th><th>Reference</th></tr></thead><tbody><tr><td>11253226</td><td>1590.6<sup>+10.6</sup><sub>-6.8</sub></td><td>76.858<sup>+0.026</sup><sub>-0.030</sub></td><td>6642 ± 77</td><td>-0.08 ± 0.1</td><td>1</td></tr><tr><td>11772920</td><td>3674.7<sup>+55.1</sup><sub>-36.1</sub></td><td>157.746<sup>+0.032</sup><sub>-0.033</sub></td><td>5180 ± 180</td><td>-0.09 ± 0.18</td><td>3</td></tr><tr><td>12009504</td><td>1865.6<sup>+7.7</sup><sub>-6.2</sub></td><td>88.217<sup>+0.026</sup><sub>-0.025</sub></td><td>6179 ± 77</td><td>-0.08 ± 0.1</td><td>1</td></tr><tr><td>12069127</td><td>884.7<sup>+10.1</sup><sub>-8.0</sub></td><td>48.400<sup>+0.048</sup><sub>-0.048</sub></td><td>6276 ± 77</td><td>0.08 ± 0.1</td><td>1</td></tr><tr><td>12069424</td><td>2188.5<sup>+4.6</sup><sub>-3.0</sub></td><td>103.277<sup>+0.021</sup><sub>-0.020</sub></td><td>5825 ± 50</td><td>0.10 ± 0.026</td><td>5</td></tr><tr><td>12069449</td><td>2561.3<sup>+5.0</sup><sub>-5.6</sub></td><td>116.929<sup>+0.012</sup><sub>-0.013</sub></td><td>5750 ± 50</td><td>0.05 ± 0.021</td><td>5</td></tr><tr><td>12258514</td><td>1512.7<sup>+3.3</sup><sub>-2.9</sub></td><td>74.799<sup>+0.016</sup><sub>-0.015</sub></td><td>5964 ± 77</td><td>-0.00 ± 0.1</td><td>1</td></tr><tr><td>12317678</td><td>1212.4<sup>+5.5</sup><sub>-4.9</sub></td><td>63.464<sup>+0.025</sup><sub>-0.024</sub></td><td>6580 ± 77</td><td>-0.28 ± 0.1</td><td>1</td></tr></tbody></table>

Note. All values of the average large frequency separation ⟨Δν⟩ and the frequency of maximum power νmax have been determined by Lund et al. (2017). Reference column indicates the source of the atmospheric parameters: (1) Buchhave & Latham (2015), (2) Pinsonneault et al. (2012), (3) Pinsonneault et al. (2014), (4) Casagrande et al. (2014b), (5) Ramírez et al. (2009), (6) Chaplin et al. (2014), (7) Huber et al. (2013).

(This table is available in its entirety in machine-readable form.)

was obtained from a Gaussian fit to the radial mode amplitudes from the peak-bagging.

## Appendix C
### Catalogs of Stellar Properties

The resulting stellar properties determined by each pipeline are published in the online version of this paper. A description of all fields available is given in Table 4. For the pipelines that do not return asymmetric uncertainties in the stellar properties, we report the standard deviation in the positive and negative uncertainties column.

<table><caption>**Table 4**
Stellar Properties for the LEGACY Dwarfs Sample
Determined by each Pipeline</caption><thead><tr><th>Field</th><th>Description</th></tr></thead><tbody><tr><td>KIC</td><td>Kepler Input Catalog Identifier</td></tr><tr><td>Mass</td><td>Mass in solar units</td></tr><tr><td>sMassP</td><td>Positive mass uncertainty in solar units</td></tr><tr><td>sMassM</td><td>Negative mass uncertainty in solar units</td></tr><tr><td>Rad</td><td>Radius in solar units</td></tr><tr><td>sRadP</td><td>Positive radius uncertainty in solar units</td></tr><tr><td>sRadM</td><td>Negative radius uncertainty in solar units</td></tr><tr><td>Grav</td><td>Surface gravity in dex</td></tr><tr><td>sGravP</td><td>Positive surface gravity uncertainty in dex</td></tr><tr><td>sGravM</td><td>Negative surface gravity uncertainty in dex</td></tr><tr><td>Age</td><td>Age in units of Gyr</td></tr><tr><td>sAgeP</td><td>Positive age uncertainty in units of Gyr</td></tr><tr><td>sAgeM</td><td>Negative age uncertainty in units of Gyr</td></tr><tr><td>Lum</td><td>Luminosity in solar units</td></tr><tr><td>sLumP</td><td>Positive luminosity uncertainty in solar units</td></tr><tr><td>sLumM</td><td>Negative luminosity uncertainty in solar units</td></tr><tr><td>Rho</td><td>Density in g cm<sup>-3</sup></td></tr><tr><td>sRhoP</td><td>Positive density uncertainty in g cm<sup>-3</sup></td></tr><tr><td>sRhoM</td><td>Negative density uncertainty in g cm<sup>-3</sup></td></tr><tr><td>Dist</td><td>Distance in pc</td></tr><tr><td>sDistP</td><td>Positive distance uncertainty in pc</td></tr><tr><td>sDistM</td><td>Negative distance uncertainty in pc</td></tr><tr><td>Xini</td><td>Fractional initial hydrogen abundance</td></tr><tr><td>sXiniP</td><td>Positive fractional initial hydrogen abundance uncertainty</td></tr><tr><td>sXiniM</td><td>Negative fractional initial hydrogen abundance uncertainty</td></tr><tr><td>Yini</td><td>Fractional initial helium abundance</td></tr><tr><td>sYiniP</td><td>Positive fractional initial helium abundance uncertainty</td></tr><tr><td>sYiniM</td><td>Negative fractional initial helium abundance uncertainty</td></tr></tbody></table>

<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>Field</th><th>Description</th></tr></thead><tbody><tr><td>Xsup</td><td>Fractional surface hydrogen abundance</td></tr><tr><td>sXsupP</td><td>Positive fractional surface hydrogen abundance uncertainty</td></tr><tr><td>sXsupM</td><td>Negative fractional surface hydrogen abundance uncertainty</td></tr><tr><td>Ysup</td><td>Fractional surface helium abundance</td></tr><tr><td>sYsupP</td><td>Positive fractional surface helium abundance uncertainty</td></tr><tr><td>sYsupM</td><td>Negative fractional surface helium abundance uncertainty</td></tr><tr><td>Xcen</td><td>Fractional central hydrogen abundance</td></tr><tr><td>sXcenP</td><td>Positive fractional central hydrogen abundance uncertainty</td></tr><tr><td>sXcenM</td><td>Negative fractional central hydrogen abundance uncertainty</td></tr><tr><td>Ycen</td><td>Fractional central helium abundance</td></tr><tr><td>sYcenP</td><td>Positive fractional central helium abundance uncertainty</td></tr><tr><td>sYcenM</td><td>Negative fractional central helium abundance uncertainty</td></tr><tr><td>MCcore</td><td>Mass coordinate of the convective core edge</td></tr><tr><td>sMCcoreP</td><td>Positive fractional uncertainty in mass coordinate of convective core edge</td></tr><tr><td>sMCcoreM</td><td>Negative fractional uncertainty in mass coordinate of convective core edge</td></tr><tr><td>Rbce</td><td>Radius coordinate of the base of the convective envelope</td></tr><tr><td>sRbceP</td><td>Positive fractional uncertainty in radius coordinate of the base of the convective envelope</td></tr><tr><td>sRbceM</td><td>Negative fractional uncertainty in radius coordinate of the base of the convective envelope</td></tr><tr><td>α<sub>conv</sub></td><td>Convective efficiency</td></tr><tr><td>sα<sub>convP</sub></td><td>Positive fractional uncertainty in convective efficiency</td></tr><tr><td>sα<sub>convM</sub></td><td>Negative fractional uncertainty in convective efficiency</td></tr><tr><td>TAMS</td><td>Terminal age main sequence</td></tr></tbody></table>

Note. Asymmetric uncertainties are given when available, otherwise the standard deviation of the property is given in the positive and negative uncertainty fields. Density uncertainties for the C2kSMO results are not reported, while convective efficiencies and TAMS values are not available for all pipelines (see Section 4.1 for details).

(This table is available in its entirety in machine-readable form.)

## References

Adelberger, E. G., Austin, S. M., Bahcall, J. N., et al. 1998, RvMP, 70, 1265
Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41
Amôres, E. B., & Lépine, J. R. D. 2005, AJ, 130, 659
Angulo, C., Arnould, M., Rayet, M., et al. 1999, NuPhA, 65, 3
Antia, H. M., & Basu, S. 1994, A&AS, 107, 421
Appourchaux, T., Antia, H. M., Ball, W., et al. 2015, A&A, 582, A25
Appourchaux, T., Chaplin, W. J., García, R. A., et al. 2012, A&A, 543, A54

Badnell, N. R., Bautista, M. A., Butler, K., et al. 2005, MNRAS, 360, 458

Bahcall, J. N., Pinsonneault, M. H., & Wasserburg, G. J. 1995, *RvMP*, 67, 781

Ball, W. H., Beeck, B., Cameron, R. H., & Gizon, L. 2016, A&A, 592, A159

Ball, W. H., & Gizon, L. 2014, A&A, 568, A123

Balona, L. A. 2015, MNRAS, 447, 2714

Balser, D. S. 2006, AJ, 132, 2326

Basu, S. 1998, MNRAS, 298, 719

Bland, J. M., & Altman, D. G. 1986, *Lancet* (London, England), 1, 307

Bland, J. M., & Altman, D. G. 1995, *Lancet* (London, England), 346, 1085

Böhm-Vitense, E. 1958, ZA, 46, 108

Bonfanti, A., Ortolani, S., & Nascimbeni, V. 2016, A&A, 585, A5

Brown, T. M., Gilliland, R. L., Noyes, R. W., & Ramsey, L. W. 1991, ApJ, 368, 599

Bruntt, H., Basu, S., Smalley, B., et al. 2012, MNRAS, 423, 122

Buchhave, L. A., & Latham, D. W. 2015, ApJ, 808, 187

Buldgen, G., Reese, D. R., & Dupret, M. A. 2015a, A&A, 583, A62

Buldgen, G., Reese, D. R., Dupret, M. A., & Samadi, R. 2015b, A&A, 574, A42

Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550

Casagrande, L., Flynn, C., Portinari, L., Girardi, L., & Jimenez, R. 2007, MNRAS, 382, 1516

Casagrande, L., Portinari, L., Glass, I. S., et al. 2014a, MNRAS, 439, 2060

Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., & Asplund, M. 2010, A&A, 512, 54

Casagrande, L., Silva Aguirre, V., Schlesinger, K. J., et al. 2016, MNRAS, 455, 987

Casagrande, L., Silva Aguirre, V., Stello, D., et al. 2014b, ApJ, 787, 110

Ceillier, T., van Saders, J., García, R. A., et al. 2015, MNRAS, 456, 119

Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1

Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353

Chiosi, C., & Matteucci, F. M. 1982, A&A, 105, 140

Christensen-Dalsgaard, J. 2002, *RvMP*, 74, 1073

Christensen-Dalsgaard, J. 2008a, Ap&SS, 316, 113

Christensen-Dalsgaard, J. 2008b, Ap&SS, 316, 13

Christensen-Dalsgaard, J., Gough, D. O., & Thompson, M. J. 1991, ApJ, 378, 413

Coelho, H. R., Chaplin, W. J., Basu, S., et al. 2015, MNRAS, 451, 3011

Cox, J.-P., & Giuli, R. T. 1968, Principles of Stellar Structure (New York: Gordon and Breach)

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2003, PhLB, 567, 227

Davies, G. R., Aguirre, V. S., Bedding, T. R., et al. 2016, MNRAS, 456, 2183

Davies, G. R., Handberg, R., Miglio, A., et al. 2014, MNRAS Letters, 445, L94

Deheuvels, S., Brandão, I., Silva Aguirre, V., et al. 2016, A&A, 589, A93

Demarque, P., Guenther, D. B., Li, L. H., Mazumdar, A., & Straka, C. W. 2008, Ap&SS, 316, 31

De Ridder, J., Molenberghs, G., Eyer, L., & Aerts, C. 2016, A&A, 595, L3

Edvardsson, B., Andersen, J., Gustafsson, B., et al. 1993, A&A, 275, 101

Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Formicola, A., Imbriani, G., Costantini, H., et al. 2004, PhLB, 591, 61

Freytag, B., Ludwig, H.-G., & Steffen, M. 1996, A&A, 313, 497

Gaia Collaboration G., Brown, A. G. A., Vallenari, A., et al. 2016, arXiv:1609.04172

Gough, D. 1987, Natur, 326, 257

Grevesse, N., & Noels, A. 1993, in 35ème cours de perfectionnement de l'Association Vandoise des Chercheurs en Physique, ed. B. Hauck, S. Paltani, & D. Raboud (Lausanne: AVCP), 205

Grevesse, N., & Sauval, A. J. 1998, SSRv, 85, 161

Halbwachs, J. L. 1986, A&AS, 66, 131

Huber, D., Chaplin, W. J., Christensen-Dalsgaard, J., et al. 2013, ApJ, 767, 127

Huber, D., Ireland, M. J., Bedding, T. R., et al. 2012, ApJ, 760, 32

Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943

Imbriani, G., Costantini, H., Formicola, A., et al. 2005, EPJA, 25, 455

Jimenez, R. 2003, Sci, 299, 1552

Johnson, J. A., Huber, D., Boyajian, T., et al. 2014, ApJ, 794, 15

Jontof-Hutter, D., Ford, E. B., Rowe, J. F., et al. 2016, ApJ, 820, 39

Jørgensen, B. R., & Lindegren, L. 2005, A&A, 436, 127

Kippenhahn, R., Weigert, A., & Weiss, A. 2012, Stellar Structure and Evolution (Berlin: Springer)

Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87

Kjeldsen, H., Bedding, T. R., & Christensen-Dalsgaard, J. 2008, ApJL, 683, L175

Lebreton, Y., & Goupil, M. J. 2014, A&A, 569, A21

Lebreton, Y., Perrin, M. N., Cayrel, R., Baglin, A., & Fernandes, J. 1999, A&A, 350, 587

Ledoux, W. P. 1947, ApJ, 105, 305

Lindegren, L., Lammers, U., Bastian, U., et al. 2016, arXiv:1609.04303

Ludbrook, J. 1997, *Clinical and experimental pharmacology & physiology*, 24, 193

Lund, M. N., Lundkvist, M., Aguirre, V. S., et al. 2014, A&A, 570, A54

Lund, M. N., Silva Aguirre, V., Davies, G. R., et al. 2017, ApJ, 835, 172

Magic, Z., Serenelli, A., Weiss, A., & Chaboyer, B. 2010, ApJ, 718, 1378

Magic, Z., Weiss, A., & Asplund, M. 2015, A&A, 573, A89

Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152

Mazumdar, A., Monteiro, M. J. P. F. G., Ballot, J., et al. 2014, ApJ, 782, 18

McDonald, I., Zijlstra, A. A., & Boyer, M. L. 2012, MNRAS, 427, 343

Metcalfe, T. S., Creevey, O. L., Doğan, G., et al. 2014, ApJS, 214, 27

Michaud, G., & Profitt, C. R. 1993, in IAU Coll. Proc. 137, Inside the Stars, ed. A. Baglin & W. W. Weiss (Cambridge: Cambridge Univ. Press), 246

Miglio, A., Chiappini, C., Morel, T., et al. 2013, MNRAS, 429, 423

Miglio, A., & Montalbán, J. 2005, A&A, 441, 615

Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61

Mosumgaard, J. R., Aguirre, V. S., Weiss, A., Christensen-Dalsgaard, J., & Trampedach, R. 2016, arXiv:1610.07323v1

Ng, Y. K., & Bertelli, G. 1998, A&A, 329, 943

Otú Floranes, H., Christensen-Dalsgaard, J., & Thompson, M. J. 2005, MNRAS, 356, 671

Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3

Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4

Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15

Pinsonneault, M. H., An, D., Molenda-Zakowicz, J., et al. 2012, ApJS, 199, 30

Pinsonneault, M. H., Elsworth, Y., Epstein, C., et al. 2014, ApJS, 215, 19

Pont, F., & Eyer, L. 2004, MNRAS, 351, 487

Ramírez, I., Meléndez, J., & Asplund, M. 2009, A&A, 508, L17

Rauer, H., Catala, C., Aerts, C., et al. 2014, ApJ, 38, 249

Reese, D. R., Chaplin, W. J., Davies, G. R., et al. 2016, A&A, 592, A14

Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064

Roxburgh, I. W. 2005, A&A, 434, 665

Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215

Salaris, M., & Cassisi, S. 2015, A&A, 577, A60

Scuflaire, R., Montalbán, J., Théado, S., et al. 2008, Ap&SS, 316, 149

Serenelli, A. M., & Basu, S. 2010, ApJ, 719, 865

Silva Aguirre, V., Ballot, J., Serenelli, A. M., & Weiss, A. 2011a, A&A, 529, 63

Silva Aguirre, V., Basu, S., Brandão, I. M., et al. 2013, ApJ, 769, 141

Silva Aguirre, V., Casagrande, L., Basu, S., et al. 2012, ApJ, 757, 99

Silva Aguirre, V., Chaplin, W. J., Ballot, J., et al. 2011b, ApJL, 740, L2

Silva Aguirre, V., Davies, G. R., Basu, S., et al. 2015, MNRAS, 452, 2127

Soderblom, D. R. 2010, ARA&A, 48, 581

Sonoi, T., Samadi, R., Belkacem, K., et al. 2015, A&A, 583, A112

Stassun, K. G., & Torres, G. 2016, ApJL, 831, L6

Steigman, G. 2010, JCAP, 2010, 29

Takeda, G., Ford, E. B., Sills, A., et al. 2007, ApJS, 168, 297

Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828

Trampedach, R., Stein, R. F., Christensen-Dalsgaard, J., Nordlund, Å., & Asplund, M. 2014, MNRAS, 445, 4366

Ulrich, R. K. 1986, ApJL, 306, L37

Van Eylen, V., & Albrecht, S. 2015, ApJ, 808, 126

van Saders, J. L., Ceillier, T., Metcalfe, T. S., et al. 2016, Natur, 529, 181

Verma, K., Antia, H. M., Basu, S., & Mazumdar, A. 2014a, ApJ, 794, 114

Verma, K., Faria, J. P., Antia, H. M., et al. 2014b, ApJ, 790, 138

Weiss, A., Hillebrandt, W., Thomas, H.-C., & Ritter, H. 2005, Cox & Giulì's Principles of Stellar Structure (2nd ed.; New York: Gordon and Breach)

Weiss, A., & Schlattl, H. 2008, Ap&SS, 316, 99

White, T. R., Bedding, T. R., Stello, D., et al. 2011, ApJ, 743, 161

White, T. R., Huber, D., Maestro, V., et al. 2013, MNRAS, 433, 1262