# Seismic inference of 57 stars using full-length Kepler data sets

Orlagh Louise Creevey, Orlagh Creevey, Travis Metcalfe, David Salabert, Savita Mathur, Mathias Schultheis, Rafael García, Frédéric Thévenin, Michaël Bazot, Haiying Xu

► To cite this version:

Orlagh Louise Creevey, Orlagh Creevey, Travis Metcalfe, David Salabert, Savita Mathur, et al.. Seismic inference of 57 stars using full-length Kepler data sets. EPJ Web of Conferences, 2017, 160, pp.03007. 10.1051/epjconf/201716003007. hal-03428801

HAL Id: hal-03428801

https://univ-cotedazur.hal.science/hal-03428801v1

Submitted on 24 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
# Seismic inference of 57 stars using full-length *Kepler* data sets

Orlagh Creevey<fnref n="1" />,*, Travis S. Metcalfe<fnref n="2" />,<fnref n="3" />, David Salabert<fnref n="4" />,**, Savita Mathur<fnref n="2" />, Mathias Schultheis<fnref n="1" />, Rafael A. García<fnref n="4" />, Frédéric Thévenin<fnref n="1" />, Michaël Bazot<fnref n="5" />, and Haiying Xu<fnref n="6" />

**Abstract.** We present stellar properties of 57 stars from a seismic inference using full-length data sets from *Kepler* (mass, age, radius, distances). These stars comprise active stars, planet-hosts, solar-analogs, and binary systems. We validate the distances derived from the astrometric Gaia-Tycho solution. Ensemble analysis of the stellar properties reveals a trend of mixing-length parameter with the surface gravity and effective temperature. We derive a linear relationship with the seismic quantity $\langle r_{02} \rangle$ to estimate the stellar age. Finally, we define the stellar regimes where the Kjeldsen et al (2008) empirical surface correction for 1D model frequencies is valid.

## 1 Introduction

Time-series analysis of the full-length *Kepler* data sets of solar-like main sequence and subgiants stars is presented in [1]. They identify modes of oscillations for 66 stars. Using the *Asteroseismic Modeling Portal*, AMP<fnref n="1" />, we analyse the individual frequency data from [1] for 57 of these stars and we supplement them with spectroscopic data from [2–7]. In this version of AMP, AMP 1.3, we fit the frequency separation ratios $r_{01}$ and $r_{02}$ [8] to determine the optimal models for each star in our sample.

The optimization method in AMP is genetic algorithm (GA) which efficiently samples the full parameter ranges without imposing constraints on unknown parameters, such as the mixing-length parameter or the initial chemical composition. The results of the GA is a dense clustering of models (1000s) around the optimal parameters. We analyse the distribution of these models to determine the stellar properties, such as mass, radius and age, along with their uncertainties. Details of the methods can be found in [9–11], with the most recent reference containing the tables of stellar parameters for the sample of 57 stars presented here. The results are validated using solar data and independent determinations of luminosity, radii, and ages.

In these proceedings we analyse the derived stellar properties of our sample. We compare their predicted distances with the recent solution provided by the Gaia-Tycho analysis (TGAS, [12], Sect. 2). In Section 3 we derive expressions for estimating the mixing-length parameter and

the stellar age based on observed properties. Finally, in Section 4 we explore the range of parameters where the correction proposed by [13] to mimic the so-called surface effect related to seismic data is useful.

<figure><img src="image_2.png" /><figcaption>**Figure 1:** HR diagram showing the position of the 57 *Kepler* stars used in this work. Evolutionary tracks for solar-metallicity models with 1.0, 1.2, and 1.4 M⊙ stellar masses are shown.</figcaption></figure>

## 2 Asteroseismic distances

We use the stellar luminosity, *L*, constrained from the asteroseismic analysis to compute the stellar distance, as a parallax. The model surface gravity and the observed $T_{eff}$ and [Fe/H] are used to derive the amount of interstellar absorption between the top of the Earth's atmosphere and the

<footnote n="1">¹Université Côte d'Azur, Laboratoire Lagrange, CNRS, Blvd de l'Observatoire, CS 34229, 06304 Nice cedex 4, France</footnote>

<footnote n="2">²Center for Extrasolar Planetary Systems, Space Science Institute, 4750 Walnut St. Suite 205, Boulder CO 80301, USA</footnote>

<footnote n="3">³Visiting Scientist, National Solar Observatory, 3665 Discovery Dr., Boulder CO 80303, USA</footnote>

<footnote n="4">⁴Laboratoire AIM, CEA/DRF-CNRS, Université Paris 7 Diderot, IRFU/SAp, Centre de Saclay, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="5">⁵Center for Space Science, NYUAD Institute, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, UAE</footnote>

<footnote n="6">⁶Computational & Information Systems Laboratory, NCAR, P.O. Box 3000, Boulder CO 80307, USA</footnote>

<footnote>*e-mail: ocreevey@oca.eu</footnote>

<footnote>**Visiting Scientist at 1</footnote>

<footnote n="1">¹https://amp.phys.au.dk/</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2:** Comparison of the parallaxes derived in this work with those from the TGAS solution. Error bars are TGAS errors.</figcaption></figure>

star, $A_{Ks}$, by applying the isochrone method described in [14]. Here, this subscript $Ks$ refers to the 2MASS $Ks$ filter [15]. We compute the bolometric correction $BC_{Ks}$, using $BC_{Ks} = 4.514650 - 0.000524 T_{eff}$ [16] where the solar bolometric magnitude is 4.72 mag. With $L$, $Ks$, $A_{Ks}$, and $BC_{Ks}$ we derive the distance to each of the stars in our sample, and then its parallax.

The TGAS catalogue of stellar parallaxes was recently made available through the first Gaia Data Release [17]. A comparison of the parallaxes we derive and these new values is shown in Figure 2. Apart from a few outliers there is an overall good agreement between the independent methods. However, the parallaxes that we derive are systematically larger, with a mean difference of +0.7 mas, just over $2\sigma$ the current uncertainties on the TGAS parallaxes. The cause of this systematic difference is currently under investigation. Nevertheless it is worthy to note that such comparisons of independent methods may allow us to reveal errors in the extinction measurements, the bolometric corrections, the $T_{eff}$ or the underlying assumptions in the models which compute the stellar luminosities.

# 3 Trends in stellar properties

Performing a homogenous analysis on a large sample allows us to check for trends in some stellar parameters, and compare them to trends derived or established by other methods. We performed this check for two parameters: the mixing-length parameter and the stellar age.

## 3.1 The mixing-length parameter versus $T_{eff}$ and $\log g$

The mixing-length parameter, $\alpha$, is usually calibrated for a solar model and then applied to all models for a set range of masses and metallicities. However, several authors have shown that this approach is not correct [18, 19]. The values of $\alpha$ resulting from a GA offer an optimal approach to effectively test and subsequently constrain this parameter, since the only assumption is that $\alpha$ lies between 1 and 3.

The distribution of $\alpha$ with $\log g$ and $T_{eff}$ (colour-coded, see figure caption) is shown in Fig. 3. We see that for a

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

<figure><img src="image_5.png" /><figcaption>**Figure 3:** Top: Distribution of the $\log g$ and $\alpha$ for our sample. The colour coding is as follows: red – $T_{eff} < 5600$ K, yellow – $5600 K < T_{eff} < 6000$ K, green – $6000 K < T_{eff} < 6300$ K, blue – $T_{eff} > 6300$ K. Lower: The distribution of age and $\langle r_{02} \rangle$.</figcaption></figure>

given value of $\log g$, the value of $\alpha$ has an upper limit, which we can approximate by $\alpha < 1.65 \log g - 4.75$, and this is denoted by the dashed line in the figure. A regression analysis considering the model values of $\log g$, $\log T_{eff}$ and [M/H] yields:

$$ \alpha = 5.972778 + 0.636997 \log g - 1.799968 \log T_{eff} + 0.040094 [\text{M/H}], \quad (1) $$

with a mean and rms of the residual to the fit of -0.01 ± 0.15. This equation yields a value of $\alpha = 2.03$ for the known solar properties.

These results agree in part with those derived by [20], using full 3D radiative hydrodynamic calculations for convective envelopes. These authors also found that $\alpha$ increases with $\log g$ and decreases with $T_{eff}$. Our fit indicates a very small dependence on metallicity while their results finds an opposite and more significant trend with this parameter. Our sample, however, does not span a very large range in [M/H], and the low coefficient is consistent with zero within the error bars.

## 3.2 Age and $\langle r_{02} \rangle$

The $r_{02}$ frequency separation ratios are effective at probing the gradients near the core of the star [8]. As the core is most sensitive to nuclear processing, the $r_{02}$ are a diagnostic of the evolutionary state of the star. Using theoretical models [21] showed a relationship between the mean
value of $r_{02}$ and the stellar age. That relationship was recently used by [22] to estimate the age of KIC 7510397 (HIP 93511). Figure 3 shows the distribution of $\langle r_{02} \rangle$ versus the derived age for the sample of stars studied here. A linear fit to these data leads to the following estimate of the stellar age, $\tau$, based on $\langle r_{02} \rangle$

$$ \tau = 17.910 - 193.918 \langle r_{02} \rangle, \quad (2) $$

This is, of course, only valid for the range covered by our sample. Note that when inserting the solar value of $\langle r_{02} \rangle = 0.068$ µHz, Eq. 2 yields an age of 4.7 Gyr, in agreement with the Sun's age as determined by other means.

## 4 Surface Effects

The comparison of the observed oscillation frequencies from the Sun and other stars with model frequencies calculated from 1D stellar models reveals a systematic error in the models which increases with frequency. These are known as surface effects, which arise from incomplete modelling of the near-surface layers and the use of an adiabatic treatment on the stellar oscillation modes. Efforts to improve the stellar modelling is underway, but the application of improved models on a large scale is still out of reach. To alleviate this problem, several authors have advocated for the use of combination frequencies which are insensitive to this systematic offset [8], hence the exclusive use of $r_{01}$ and $r_{02}$ in the AMP 1.3 methodology. However, since individual frequencies contain more information than $r_{01}$ and $r_{02}$, some authors have derived simple prescriptions to mitigate the surface effect. One such parametrization is that of [13] who suggest a simple correction to the individual frequencies $\delta v_{n,l}$ in the form of a power law, namely

$$ \delta v_{n,l} = a_0 \left( \frac{v_{n,l}^{\text{obs}}}{v_{\text{max}}} \right)^b \quad (3) $$

where $b = 4.82$ is a fixed value, calibrated by a solar model, $a_0$ is computed from the differences between the observed and model frequencies [10, 23], $v_{n,l}^{\text{obs}}$ is the observed frequency $n, l$ mode and $v_{\text{max}}$ is the frequency corresponding to the highest amplitude mode, see [1].

The AMP 1.3 methodology uses exclusively $r_{01}$ and $r_{02}$ as the seismic constraints, hence our results are insensitive to the surface effects. Using the model frequencies of the best-matched models, we can test how useful Eq. 3 is for different stellar regimes along the main sequence and early sub-giant phase. The interest in testing this becomes apparent when we consider that for many stars we do not have a high enough precision on the frequencies, or a sufficiently large range of radial orders, to use $r_{01}$ and $r_{02}$ to effectively constrain the stellar modelling. This is the case for some ground-based observations and for some stars that will be observed by TESS [24], and PLATO [25], where a limited time series of only one to two months may only be available.

In order to test where Eq. 3 is useful, we calculated the residual between the observed oscillation frequencies

and the model frequencies corrected for the surface term (Eq. 3): $q_{n,l} = v_{n,l}^{\text{obs}} - v_{n,l}^{\text{mod}} + \delta v_{n,l}$. Then we defined the metric $Q$ as the median of the square root of the squared residuals:

$$ Q = \text{median} \left\{ \sqrt{q_{n,l}^2} \right\}, \quad (4) $$

for all $n$ and all $l$ defined in the region of $0.7 \le v_{n,l}^{\text{obs}}/v_{\text{max}} \le 1.3$. We find values of $Q$ that vary between 0.3 and 10 and this variation is anti-correlated with $a_0$: as $a_0$ becomes more negative the match with the model becomes worse.

We define a subsample of 44 stars with the best fits to $r_{01}$ and $r_{02}$. This subsample is shown in Fig. 4 as open circles for the observed values of $\langle \Delta v \rangle$ and $T_{\text{eff}}$, and the derived values of mass and radius. In these figures the stars represented by the dark blue filled circles have $Q \le 1.0$ and those represented by the light blue filled circles have $1.0 < Q < 1.2$. It is evident from the figures that there exists regions of the parameter space where the surface correction proposed by [13] adequately corrects for the surface term. These regions are highlighted by the dashed lines, defined more conservatively here as $\log g > 4.2$, $T_{\text{eff}} < 6250$ K, $\langle \Delta v \rangle > 80\mu\text{Hz}$ and $v_{\text{max}} > 1700\mu\text{Hz}$. In physical properties this corresponds to a star with $R < 1.6 R_{\odot}$, $M < 1.3 M_{\odot}$, and $L < 2.5 L_{\odot}$, with no apparent evidence of the age or the metallicity playing any role.

## 5 Conclusions

In these proceedings we used the derived properties of 57 *Kepler* stars presented in [11] to predict the star's parallaxes, to derive an expression relating the mixing-length parameter and the age to observed properties, and to explore the regions of parameters where the proposed surface correction to individual frequencies by [13] (Eq. 3) is useful. For this latter point, the aim of defining valid regions is to use the correction for data sets where the frequency precision or the range of radial orders is not sufficient for $r_{01}$ and $r_{02}$ to constrain the stellar models. This will be the case for some stars that will be observed in low ecliptic latitudes with TESS and for the *step-and-stare* pointing phases of the PLATO mission.

The parallaxes that we derived in this work were compared to the parallaxes derived from the Tycho-Gaia solution. We found in general a good agreement between our results. The mean differences between our results (ours – theirs) is +0.7 mas and the cause of this error is currently under study. Nevertheless we validate the new parallaxes from [12]. This comparison of parallaxes also highlights the ability of these external measurements to allow us to investigate the sources of systematic errors. We believe that in providing a prior on the luminosity, we may even be able to constrain the initial helium abundance in the star. We look forward to the wealth of forthcoming asteroseismic data on bright nearby stars with TESS and PLATO, and the forthcoming Gaia parallaxes with $\mu$s precision.

## References

[1] Lund, M. N., Silva Aguirre, V., Davies, G. R., et al., ApJ, **835**, 172 (2017)
<figure><img src="image_6.png" /><figcaption>**Figure 4:** Distribution of observed (left) and derived parameters (right) for a subsample of our stars (open circles). The dark/light blue filled dots represent the stars with $Q \le 1.0/1.2$. The prescription for correcting model frequencies is valid within the region delimited by dashed lines.</figcaption></figure>

[2] I. Ramírez, J. Meléndez, M. Asplund, A&A, **508**, L17 (2009)

[3] L.A. Buchhave, D.W. Latham, A. Johansen, M. Bizarro, G. Torres, J.F. Rowe, N.M. Batalha, W.J. Borucki, E. Brugamyer, C. Caldwell et al., Nature, **486**, 375 (2012)

[4] M.H. Pinsonneault, D. An, J. Molenda-Žakowicz, W.J. Chaplin, T.S. Metcalfe, H. Bruntt, ApJS, **199**, 30 (2012), 1110.4456

[5] D. Huber, W.J. Chaplin, J. Christensen-Dalsgaard, R.L. Gilliland, H. Kjeldsen, L.A. Buchhave, D.A. Fischer, J.J. Lissauer, J.F. Rowe, R. Sanchis-Ojeda et al., ApJ, **767**, 127 (2013), 1302.2624

[6] W.J. Chaplin, S. Basu, D. Huber, A. Serenelli, L. Casagrande, V. Silva Aguirre, W.H. Ball, O.L. Creevey, L. Gizon, R. Handberg et al., ApJS, **210**, 1 (2014), 1310.4001

[7] M.H. Pinsonneault, Y. Elsworth, C. Epstein, S. Hekker, S. Mészáros, W.J. Chaplin, J.A. Johnson, R.A. García, J. Holtzman, S. Mathur et al., ApJS, **215**, 19 (2014), 1410.2503

[8] I.W. Roxburgh, S.V. Vorontsov, A&A, **411**, 215 (2003)

[9] T.S. Metcalfe, P. Charbonneau, Journal of Computational Physics **185**, 176 (2003), astro-ph/0208315

[10] T.S. Metcalfe, O.L. Creevey, J. Christensen-Dalsgaard, ApJ, **699**, 373 (2009), 0903.0616

[11] O. Creevey, T.S. Metcalfe, D. Salabert, S. Mathur, M. Schultheis, R.A. García, F. Thévenin, M. Bazot, Z. Haiying, A&A, subm. (2017)

[12] D. Michalik, L. Lindegren, D. Hobbs, A&A, **574**, A115 (2015), 1412.8770

[13] H. Kjeldsen, T.R. Bedding, J. Christensen-Dalsgaard, ApJL, **683**, L175 (2008), 0807.1769

[14] M. Schultheis, G. Zasowski, C. Allende Prieto, F. Anders, R.L. Beaton, T.C. Beers, D. Bizyaev, C. Chiappini, P.M. Frinchaboy, A.E. García Pérez et al., AJ, **148**, 24 (2014), 1405.2180

[15] M.F. Skrutskie, R.M. Cutri, R. Stiening, M.D. Weinberg, S. Schneider, J.M. Carpenter, C. Beichman, R. Capps, T. Chester, J. Elias et al., AJ, **131**, 1163 (2006)

[16] P. Marigo, L. Girardi, A. Bressan, M.A.T. Groenewegen, L. Silva, G.L. Granato, A&A, **482**, 883 (2008), 0711.4922

[17] Gaia Collaboration, ArXiv e-prints (2016), 1609.04153

[18] A. Bonaca, J.D. Tanner, S. Basu, W.J. Chaplin, T.S. Metcalfe, M.J.P.F.G. Monteiro, J. Ballot, T.R. Bedding, A. Bonanno, A.M. Broomhall et al., ApJL, **755**, L12 (2012), 1207.2765

[19] O.L. Creevey, F. Thévenin, T.S. Boyajian, P. Kervella, A. Chiavassa, L. Bigot, A. Mérand, U. Heiter, P. Morel, B. Pichon et al., A&A, **545**, A17 (2012), 1207.5954

[20] Z. Magic, A. Weiss, M. Asplund, A&A, **573**, A89 (2015), 1403.1062

[21] Y. Lebreton, J. Montalbán, *Stellar ages from asteroseismology*, in *The Ages of Stars*, edited by E.E. Majek, D.R. Soderblom, R.F.G. Wyse (2009), Vol. 258 of *IAU Symposium*, pp. 419–430, 0811.2908

[22] T. Appourchaux, H.M. Antia, W. Ball, O. Creevey, Y. Lebreton, K. Verma, S. Vorontsov, T.L. Campante, G.R. Davies, P. Gaulme et al., A&A, **582**, A25 (2015)

[23] T.S. Metcalfe, O.L. Creevey, G. Doğan, S. Mathur, H. Xu, T.R. Bedding, W.J. Chaplin, J. Christensen-Dalsgaard, C. Karoff, R. Trampedach et al., ApJS, **214**, 27 (2014), 1402.3614

[24] G.R. Ricker, J.N. Winn, R. Vanderspek, D.W. Latham, G. Á. Bakos, J.L. Bean, Z.K. Berta-Thompson, T.M. Brown, L. Buchhave, N.R. Butler et al., Journal of Astronomical Telescopes, Instruments, and Systems **1**, 014003 (2015)

[25] H. Rauer, C. Catala, C. Aerts, T. Appourchaux, W. Benz, A. Brandeker, J. Christensen-Dalsgaard, M. Deleuil, L. Gizon, M.J. Goupil et al., Experimental Astronomy **38**, 249 (2014), 1310.0696