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

# The Correlation between Mixing Length and Metallicity on the Giant Branch: Implications for Ages in the Gaia Era

Jamie Tayar¹, Garrett Somers¹,², Marc H. Pinsonneault¹, Dennis Stello³,⁴,⁵, Alexey Mints⁴,⁶, Jennifer A. Johnson¹, O. Zamora⁷,⁸, D. A. García-Hernández⁷,⁸, Claudia Maraston⁹, Aldo Serenelli¹⁰, Carlos Allende Prieto⁷,⁸, Fabienne A. Bastien¹¹,<fnref n="27" />, Sarbani Basu¹², J. C. Bird², R. E. Cohen¹³, Katia Cunha¹⁴, Yvonne Elsworth¹⁵, Rafael A. García¹⁶, Leo Girardi¹⁷, Saskia Hekker⁴,⁶, Jon Holtzman¹⁸, Daniel Huber³,⁴,¹⁹, Savita Mathur²⁰, Szabolcs Mészáros²¹,<fnref n="28" />, B. Mosser²², Matthew Shetrone²³, Victor Silva Aguirre⁴, Keivan Stassun², Guy S. Stringfellow²⁴, Gail Zasowski²⁵, and A. Roman-Lopes²⁶

¹Department of Astronomy, Ohio State University, 140 West 18th Avenue, OH 43210, USA

²Department of Physics and Astronomy, Vanderbilt University, 6301 Stevenson Circle, Nashville, TN 37235, USA

³Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia

⁴Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark

⁵School of Physics, University of New South Wales, NSW 2052, Australia

⁶Max-Planck-Institut für Sonnensystemforschung, Justus-von-Liebig-Weg 3, D-37077 Göttingen, Germany

⁷Instituto de Astrofísica de Canarias (IAC), Vía Lactea s/n, E-38205 La Laguna, Tenerife, Spain

⁸Departamento de Astrofísica, Universidad de La Laguna (ULL), E-38206 La Laguna, Tenerife, Spain

⁹ICG—University of Portsmouth, Burnaby Road, PO1 3FX, Portsmouth, UK

¹⁰Institute of Space Sciences (CSIC-IEEC), Carrer de Can Magrans, Barcelona, E-08193, Spain

¹¹Department of Astronomy and Astrophysics, 525 Davey Lab, The Pennsylvania State University, University Park, PA 16803, USA

¹²Department of Astronomy, Yale University, New Haven, CT 06511, USA

¹³Departamento de Astronomía, Universidad de Concepción, Casilla 160-C, Concepción, Chile

¹⁴Observatório Nacional-MCTI, Brazil

¹⁵School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK

¹⁶Laboratoire AIM, CEA/DRF-CNRS, Université Paris 7 Diderot, IRFU/SAp, Centre de Saclay, F-91191, Gif-sur-Yvette, France

¹⁷Osservatorio Astronomico di Padova-INAF, Vicolo dell'Osservatorio 5, I-35122 Padova, Italy

¹⁸New Mexico State University, Las Cruces, NM 88003, USA

¹⁹SETI Institute, 189 Bernardo Avenue, Mountain View, CA 94043, USA

²⁰Space Science Institute, 4750 Walnut Street Suite 205, Boulder, CO 80301, USA

²¹ELTE Gothard Astrophysical Observatory, H-9704 Szombathely, Szent Imre herceg st. 112, Hungary

²²LESIA, Observatoire de Paris, PSL Research University, CNRS, Université Pierre et Marie Curie, Université Paris Diderot, F-92195 Meudon, France

²³University of Texas at Austin, McDonald Observatory, 32 Fowlkes Road, McDonald Observatory, TX 79734, USA

²⁴Center for Astrophysics and Space Astronomy, University of Colorado, 389 UCB, Boulder, CO 80309, USA

²⁵Department of Physics and Astronomy, JHU, 3400 North Charles Street Baltimore, MD 21218, USA

²⁶Departamento de Física, Facultad de Ciencias, Universidad de La Serena, Cisternas 1200, La Serena, Chile

Received 2016 September 23; revised 2017 March 2; accepted 2017 March 27; published 2017 April 28

## Abstract

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested, and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity-dependent convective mixing length is used, with $\Delta_{\alpha ML}$, $Y_{REC}$ ~ 0.2 per dex in metallicity, a trend inconsistent with the predictions of three-dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of two even at modest deviations from solar metallicity ($[\text{Fe/H}] = -0.5$).

**Key words:** stars: evolution – stars: fundamental parameters

## 1. Introduction

The theory of stellar structure and evolution makes a rich web of predictions about the life histories of stars. Many predictions of this theory have proven accurate; in particular, the agreement between the predicted and actual positions of the core-hydrogen-burning main sequence in the Hertzsprung-Russell diagram was a major triumph for 20th century astrophysics. The situation for more evolved stars, however, has been more challenging to evaluate. The temperature locus of evolved red giants is sensitive to the input physics in

general, and to the efficiency of stellar convection in particular. This efficiency of stellar convection, an inherently three-dimensional process, has typically been parameterized in one-dimensional stellar models as an effective mixing length (Böhm-Vitense 1958). Testing the validity of models in this more evolved regime has historically been difficult, due to the lack of fundamental masses and reliable absolute spectroscopic measurements of temperature, and detailed abundances for large samples of stars, although small samples have generally indicated reasonable agreement or only small discrepancies (e.g., Huber et al. 2012; Takeda et al. 2016).

Correct modeling of the position of the red giant branch on the Hertzsprung-Russell diagram is now of increased importance, in light of the recently released Gaia data (Lindegren et al. 2016).

<footnote n="27">²⁷ Hubble Fellow.</footnote>

<footnote n="28">²⁸Premium Postdoctoral Fellow of the Hungarian Academy of Sciences.</footnote>
Gaia's ability to provide accurate luminosities for many red giants, covering a much wider metallicity and age range than Hipparcos could, means that giants could be used to map galactic formation and evolution, if we were able to correctly interpret their ages. Because of the steepness of the red giant branch in the HR diagram, and the relatively small effect of age, stars with a given luminosity and temperature can have widely different derived ages, depending on the adopted mixing length (Freytag & Salaris 1999).

In this paper, we test the predictions of commonly used stellar models against spectroscopic and asteroseismic data from the APOKASC survey (M. H. Pinsonneault et al. 2017, in preparation). Our test of red giant branch stellar evolution begins with the input physics for the theoretical models. There is a consensus set of physics used in standard stellar models (see Section 2), which leads to robust model predictions for solar-type stars on the main sequence. There are numerous sources of uncertainty in stellar modeling due to the extensive number of physical inputs, including the nuclear reaction rates, convective overshoot, opacities, equation of state, diffusion, and the outer boundary condition. The single largest uncertainty in theoretical predictions about the locus of the red giant branch is the convective efficiency, or mixing length $\alpha$, which is typically calibrated to reproduce the known solar radius and luminosity at the solar age (e.g., Bressan et al. 2012; Choi et al. 2016). There is, however, no strong a priori reason why such a calibration should apply to evolved stars. In fact, three-dimensional convection simulations have indicated that the mixing length should depend on parameters like luminosity, gravity, and metallicity (Trampedach et al. 2014; Magic et al. 2015). Recent observational evidence has also indicated that a solar mixing length is not the best fit to every star (e.g., Bonaca et al. 2012; Metcalfe et al. 2014; Mann et al. 2015; Saio et al. 2015; Wu et al. 2015). For the purposes of this particular study, we treat the mixing length as a free parameter that we use to quantify the metallicity-dependent offset in temperature between theory and data.

## 2. Methods

The first ingredient is a comprehensive spectroscopic data set, here provided by the APOGEE survey (Holtzman et al. 2015; Majewski et al. 2015). This survey employed an automated pipeline analysis (Nidever et al. 2015; García Pérez et al. 2016) of high-resolution H-band spectra obtained for more than 100,000 stars (Zasowski et al. 2013), collected using the 2.5 m Sloan Digital Sky Survey telescope (Gunn et al. 2006). Crucially, this survey included over 6000 red giants with asteroseismic data from the Kepler mission. For these stars, we took the measured frequency of maximum power ($\nu_{max}$) and large frequency spacing ($\Delta\nu$), and used the seismic scaling relations (e.g., Kjeldsen & Bedding 1995) to compute a mass and surface gravity (see Section 3.3 for a discussion of proposed alterations to the scaling relations).

$$ \frac{M}{M_{\odot}} \propto \left( \frac{\nu_{\max}}{\nu_{\max, \odot}} \right)^3 \left( \frac{\Delta\nu}{\Delta\nu_{\odot}} \right)^{-4} \left( \frac{T_{\text{eff}}}{T_{\text{eff}, \odot}} \right)^{3/2} \\ \frac{g}{g_{\odot}} \propto \left( \frac{\nu_{\max}}{\nu_{\max, \odot}} \right) \left( \frac{T_{\text{eff}}}{T_{\text{eff}, \odot}} \right)^{1/2} $$

This allowed us to compare the data to theoretical tracks in the theoretical variables of mass, temperature, and gravity, rather

than observational spaces like color or magnitude (see Section 3).

For this analysis, we used the second version of the combined APOGEE-Kepler data set (APOKASC, M. H. Pinsonneault et al. 2017, in preparation). This is an update of the Pinsonneault et al. (2014) APOKASC catalog containing thousands more stars, improved seismic analysis using the full Kepler data set from a larger number of analysis pipelines (Y. Elsworth et al. 2017, in preparation), and improved spectroscopic parameters from Data Release 13 (DR13, J. Holtzman et al. 2017, in preparation) of the Sloan Digital Sky Survey IV (Blanton et al. 2017). Of the 11876 APOGEE stars with Kepler data, 8887 were not flagged as dwarfs and were checked for solar-like oscillations. Of those stars, 6076 were returned as giants whose seismic frequency pattern reliably indicated whether they were core helium burning (clump, secondary clump) or not (red giant branch, asymptotic giant branch). For our analysis, we use only stars that the seismic analysis indicated are first-ascent red giants in order to avoid uncertainties in mass loss, and therefore initial stellar mass, that arise when stars pass through the tip of the giant branch to the red clump. We used stars with asteroseismic scaling relation masses and surface gravities computed using the $\nu_{max}$ and $\Delta\nu$ returned by the OCT pipeline (Hekker et al. 2010), based on the Kepler Asteroseismic Science Operations Center (KASOC, Handberg & Lund 2014) time series with the appropriate solar values for that pipeline ($\Delta\nu_{\odot} = 135.045, \nu_{max, \odot} = 3139, T_{\text{eff}, \odot} = 5771.8$). However, in Section 3, we confirm that using a different seismic analysis pipeline and various proposed corrections to the scaling relations do not significantly affect our results.

For the following discussion, we restricted the sample to the 3217 stars with OCT scaling relation masses between 0.6 and $2.6 M_{\odot}$, and values available for each seismic and spectroscopic parameter. Our mass limits were chosen to exclude approximately the highest and lowest 0.1% of stars in the APOKASC sample, where different analysis pipelines are most likely to disagree about the stellar properties or fail to report values. We note that the stars in our sample have $log(g)$ values from 1.1 to 3.3, although most of the stars in our sample are on the lower giant branch. The data table and model grid used in this analysis are available online.<fnref n="29" />

Of particular relevance to this paper, post-release analysis of DR13 spectroscopic data found that APOGEE temperatures required a metallicity-dependent correction to match the photometric temperature scale (details available online,<fnref n="30" /> documentation to be published in J. Holtzman et al. 2017, in preparation). A correction was determined using the González Hernández & Bonifacio (2009) J–K color temperature relation in low-extinction fields (see Figure 1); we note that this correction is consistent with or larger than the correction that would have been assumed using clusters, the Ramírez & Meléndez (2005) V–K temperature relation, or angular diameter temperatures collected in Ramírez & Meléndez (2005). We refer the interested reader to J. Holtzman et al. (2017, in preparation) for a more detailed discussion, but assert that the temperature offsets reported in this paper represent a lower limit on the metallicity dependence of the difference between the data and theoretical predictions.

<footnote n="29">²⁹ www.astronomy.ohio-state.edu/~tayar/MixingLength.htm</footnote>

<footnote n="30">³⁰ www.sdss.org/dr13/irspec/parameters</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Difference between the uncorrected Data Release 13 APOGEE temperatures and photometric temperatures computed using the J–K color temperature relations of González Hernández & Bonifacio (2009) for clusters (red stars) and giants in low-extinction fields binned by metallicity (black triangles). The temperature correction fit by J. Holtzman et al. (2017, in preparation) is shown as a black solid line. Note that the [Fe/H] used here represents the bulk metallicity, rather than the specific iron abundance.</figcaption></figure>

There have also been suggestions that the APOGEE metallicity scale needs a correction to match the literature metallicities quoted in Holtzman et al. (2015). Specifically, the APOGEE metallicities of open clusters are, on average, about 0.07 dex more metal-poor than literature values, whereas the globular clusters are, on average, about 0.13 dex more metal-rich. However, given the relatively small number of clusters, it was somewhat difficult to determine the most appropriate form of the correction. If we assume that it is linear, the best fit is $\Delta[\text{Fe/H}] = -0.135[\text{Fe/H}] - 0.050$. Because of our uncertainty on the form of this correction, we chose not to apply it in this paper. However, on a few key plots, we indicate what effect the metallicity correction would have, and show that in all cases it makes the effect we report here stronger and more significant.

These data were compared to two different grids of stellar evolution models (see Table 1 for a summary). The first was a grid run using the Yale Rotating Evolution Code (YREC, Pinsonneault et al. 1989, with updates as discussed in van Saders & Pinsonneault 2012). We ran a grid of masses ($0.6 M_{\odot} - 2.6 M_{\odot}$ in $0.1 M_{\odot}$ increments), metallicities ($[\text{Fe/H}] = -2.0$ to $+0.6$ in steps of 0.2), and $\alpha$-element enhancements ($+0.0, +0.2, \text{and } +0.4$) to generate tracks. These models use a Grevesse & Sauval (1998) elemental mixture, with the updated OPAL equation of state (Rogers et al. 1996; Rogers & Nayfonov 2002), OPAL opacity tables (Iglesias & Rogers 1996) and gray atmospheres. Alpha-enhanced models use alpha-enhanced starting models and opacity tables, but not equations of state. Our models include semiconvection (see Kippenhahn & Weigert 1994), but do not include overshoot or the diffusion of helium or heavy elements. We assumed a Big Bang helium value of $Y = 0.2485$ at $Z = 0$ (Cyburt 2004), and did a linear fit to our solar value of $Y = 0.272683$ at $Z = 0.018011$ (so $Y = 1.3426Z + 0.2485$). We convert Z to metallicity assuming $\log_{10}\left(\frac{Z/X}{Z_{\odot}/X_{\odot}}\right) = [\text{Fe/H}]$ for a solar mixture, with a correction factor for alpha-enhanced models. We compare to other theoretical models in Section 3. In the initial part of this analysis, we use grids with a solar

<table><caption>**Table 1**
Summary of the Input Physics Used in Our Main YREC and PARSEC Model Grids</caption><thead><tr><th>Parameter</th><th>YREC</th><th>PARSEC</th></tr></thead><tbody><tr><td>Atmosphere</td><td>Gray</td><td>Gray</td></tr><tr><td>α-enhancement</td><td>Yes</td><td>No</td></tr><tr><td>Convective Overshoot</td><td>No</td><td>Yes</td></tr><tr><td>Diffusion</td><td>No</td><td>Yes</td></tr><tr><td>Equation of State</td><td>OPAL+SCV</td><td>FREEEOS</td></tr><tr><td>High-Temperature Opacities</td><td>OPAL</td><td>OPAL</td></tr><tr><td>Low-Temperature Opacities</td><td>Ferguson et al. (2005)</td><td>ÆSOPUS (Marigo & Aringer 2009)</td></tr><tr><td>Mixing Length</td><td>1.22, 1.72, 2.22</td><td>1.74</td></tr><tr><td>Mixture</td><td>Grevesse & Sauval (1998)</td><td>Grevesse & Sauval (1998) +Caffau et al. (2011)</td></tr><tr><td>Nuclear Reaction Rates</td><td>Adelberger et al. (2011)</td><td>Bressan et al. (2012)</td></tr><tr><td>Weak Screening</td><td>Salpeter (1954)</td><td>Dewitt et al. (1973) + Graboske et al. (1973)</td></tr><tr><td>BBN He</td><td>0.2485</td><td>0.2485</td></tr><tr><td>Solar X</td><td>0.709306</td><td>0.7091328</td></tr><tr><td>Solar Y</td><td>0.272683</td><td>0.2756272</td></tr><tr><td>Solar Z</td><td>0.018011</td><td>0.01524</td></tr><tr><td>[Fe/H] Range</td><td>-2.0 to +0.6</td><td>-2.2 to +0.5</td></tr><tr><td>[α/Fe] Values</td><td>0.0, +0.2, +0.4</td><td>0.0</td></tr><tr><td>Mass Range</td><td>0.6 M<sub>\odot</sub> - 2.6 M<sub>\odot</sub></td><td>0.1 M<sub>\odot</sub> - 20 M<sub>\odot</sub></td></tr></tbody></table>

mixing length (Böhm-Vitense 1958) of 1.72, chosen to reproduce the solar luminosity ($3.827 \times 10^{33}$ erg s$^{-1}$) and the solar radius ($6.957 \times 10^{10}$ cm) at the solar age (4.57 Gyrs) (Mamajek 2012). For later parts of this study, we interpolate in a grid of models with mixing lengths of 1.22, 1.72, and 2.22. For the analysis in Section 3.2, we varied the initial helium abundance and surface boundary conditions in our model grid to explore their impact on the theoretical uncertainties. We ran models with fixed helium fractions of 0.239, 0.290, and 0.330, in addition to the solar calibrated value of 0.272683. We also created models with different atmospheric surface boundary conditions (defined as the boundary pressure at $\tau = 2/3$): a look-up table from Kurucz (1997) models, and one from Castelli & Kurucz (2004) models, in addition to the gray atmosphere used in our main model grid.

We also compare to PARSEC models (Bressan et al. 2012). These models were run at metallicities from -2.2 to +0.5, and masses between $0.1 M_{\odot}$ and $20 M_{\odot}$, over that full metallicity range, with a mixture that is based on Grevesse & Sauval (1998) with some updates from Caffau et al. (2011). Our PARSEC grid does not include $\alpha$-element enhanced tracks. The PARSEC models use Irwin's FREEEOS$^{31}$ tool to generate equations of state, OPAL high temperature opacities (Iglesias & Rogers 1996), ÆSOPUS low temperature opacities (Marigo & Aringer 2009) and gray atmospheres. These models were run at a solar mixing length of $\alpha_{MLT} = 1.74$ (Böhm-Vitense 1958). They include envelope overshoot as well as core overshoot in massive stars; microscopic diffusion is included in stars with a substantial surface convection zone unless they have a persistent convective core. Solar helium abundance is taken to be 0.276 and these models use $Y = 0.2485 + 1.78Z$.

<footnote>$^{31}$ freeeos.sourceforge.net</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Temperature and gravities of the stars in our sample between 1.1 and 1.3 $M_\odot$ are compared to models of 1.2 $M_\odot$ stars, close to the mean of this subsample. These tracks come from the YREC grid used in this work, and color coding indicates metallicity. The temperature offset is defined as the difference between the corrected APOGEE temperature and the temperature of a model with the measured mass, metallicity, [α/Fe], and surface gravity, and is represented by the horizontal offset between a point and the corresponding line on this plot.</figcaption></figure>

# 3. Analysis

## 3.1. Metallicity-dependent Temperature Offset

After interpolating to the appropriate helium for each metallicity, we perform a cubic spline interpolation in our model grid of mass, [Fe/H], [α/Fe] and log g to infer a predicted $T_{\text{eff}}$ for each star in our sample. This predicted $T_{\text{eff}}$ is then subtracted from the APOKASC $T_{\text{eff}}$ to compute a temperature offset. In Figure 2, we show a subset of our data, chosen to have similar masses and differing compositions. The overlaid tracks represent theoretical expectations; it is visually apparent that the temperatures seen in the data depend less on metallicity than predicted. To quantify this, we rank-ordered the data in metallicity, and averaged the differences between the predicted and actual temperatures in 64 star bins. In Figure 3, we present these binned means as a function of metallicity, along with contours enclosing 68 and 95% of the sample. The differences show clear trends with metallicity that are highly statistically significant. These trends are present both for YREC (top) and PARSEC (bottom). A linear fit gives $\Delta T_{\text{eff,YREC}} = 93.1$ [Fe/H] + 107.5 K and $\Delta T_{\text{eff,PARSEC}} = 127.9$ [Fe/H] + 4.1 K, with theoretical temperatures hotter than APOGEE observations at low metallicity and colder at high metallicity. We suspect that the steeper metallicity dependence of the PARSEC fit results from the use of a solar mixture without α-element enhanced options. Correcting the solar mixture models for the effect of α-element enhancement changes the temperatures by about 20 K on average (see Salaris et al. 1993), and α enhancement correlates with metallicity. We also note that there are statistically significant deviations from a linear fit. For example, there is minimal evidence for a slope using only stars above a metallicity of −0.2. However, given

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Plots of the difference between the temperature measured by APOGEE and the YREC or PARSEC model predictions as a function of metallicity. Contours indicate the extent of 68 and 95% of the stars. The Z values on the top plot assume a solar mixture. Error bars indicate the mean and spread of binned data. The dot-dashed line in the top panel indicates the slope we would have fit using the metallicity correction discussed in Section 2. We note that the temperature offset is best correlated with metallicity, with $\Delta T_{\text{eff,YREC}} = 93.12$ [Fe/H] + 107.50 K, shown as a dashed line in this and following figures. The PARSEC models have a different normalization, but show a similar slope ($\Delta T_{\text{eff,PARSEC}} = 127.88$ [Fe/H] + 4.12 K). It is clear that a single linear fit does not completely capture the variation, and we discuss in the text regions consistent with a flat line. Note that we have terminated the lines at a metallicity of −1.19, as there are fewer than 10 stars in our sample below that metallicity.</figcaption></figure>

the dependence of the slope on the details of the chosen metallicity scale, we caution against over-interpretation and provide only the linear fit. We note that the scatter around the linear relation for the YREC models is approximately Gaussian, with a standard deviation of 41 K, significantly less than the quoted APOGEE temperature uncertainties of 69 K (Alam et al. 2015).

We also show, in Figure 4, the correlations between the temperature offsets and other parameters such as log(g), mass,
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Plots of the difference between the temperature measured by APOGEE and the YREC model predictions as a function of various parameters, with log(g), mass, [α/Fe] and APOGEE temperature clockwise from top left. Points are color coded by metallicity, and contours indicate the extent of 68 and 95% of the sample. We note that the temperature offset is best correlated with metallicity (see Table 2), although correlations appear on some of these plots as a result of stellar evolution and galactic chemical evolution.</figcaption></figure>

temperature, and [α/Fe]. We find that the temperature offset is best correlated with metallicity (see Table 2 for the linear Pearson correlation coefficients). We also see many of the correlations between the stellar properties that we would expect from stellar evolution (e.g., between temperature and gravity) and galactic chemical evolution (e.g., alpha enhanced stars are likely to be old, and therefore have low mass and low metallicity). The presence of the correlation between the temperature offset and the metallicity in two independent sets of models suggests that it may be a generic problem with evolutionary models. However, before we conclude that, we must explore more mundane explanations.

## 3.2. Model Uncertainties

First, we explore whether known uncertainties in the stellar model physics, such as helium content and atmosphere boundary condition, can cause shifts in giant branch temperature on the level of what is observed in the APOKASC data. We show, in Figure 5, the expected scale shifts, as well as the differential offsets with metallicity, caused by choosing a different atmosphere boundary condition (Kurucz or Castelli), or by assuming a different helium abundance. The mean trend in the APOKASC data is shown as a dashed line. It is clear that, although there are many uncertainties in the theoretical models that can cause changes in the temperature of the giant branch locus, these shifts rarely show a strong metallicity dependence and are less than 30 K dex⁻¹, not large enough to explain our offset, which is about three times larger.

## 3.3. Mass Biases

Next, we explore whether this trend could be due to systematic errors in the APOKASC data. Although the seismic masses we use are extremely precise, there have been suggestions that they need to be corrected directly for the effects of metallicity and temperature (White et al. 2011; Epstein et al. 2014; Guggenberger et al. 2016; Sharma et al. 2016). We tested whether the observed correlation between

<caption>**Table 2**</caption>

Linear Pearson Correlation Coefficients between the Stellar Properties (Mass, Gravity, Metallicity, α Element Enhancement, and Temperature) and the Temperature Offset between the Data and the YREC Models (ΔT<sub>eff</sub>)

<table><thead><tr><th></th><th>Mass</th><th>log g</th><th>[Fe/H]</th><th>[α/Fe]</th><th>T<sub>eff</sub></th><th>ΔT<sub>eff</sub></th></tr></thead><tbody><tr><td>Mass</td><td>1.00</td><td>-0.03</td><td>0.31</td><td>-0.51</td><td>0.12</td><td>0.12</td></tr><tr><td>log g</td><td>-0.03</td><td>1.00</td><td>0.21</td><td>-0.27</td><td>0.84</td><td>0.03</td></tr><tr><td>[Fe/H]</td><td>0.31</td><td>0.21</td><td>1.00</td><td>-0.63</td><td>-0.19</td><td>0.20</td></tr><tr><td>[α/Fe]</td><td>-0.51</td><td>-0.27</td><td>-0.63</td><td>1.00</td><td>-0.21</td><td>-0.15</td></tr><tr><td>T<sub>eff</sub></td><td>0.12</td><td>0.84</td><td>-0.19</td><td>-0.21</td><td>1.00</td><td>0.00</td></tr><tr><td>ΔT<sub>eff</sub></td><td>0.12</td><td>0.03</td><td>0.20</td><td>-0.15</td><td>0.00</td><td>1.00</td></tr></tbody></table>

**Note.** A value of 1 indicates a perfect positive correlation, a value of -1 would indicate a perfect negative correlation, and a value of 0 indicates no correlation. We show here that the temperature offset is best correlated with metallicity.

temperature offset and metallicity is a result of using uncorrected scaling relations to determine the mass. For this, we estimated masses from (1) grid-based modeling using Δν computed from radial mode eigenfrequencies and the scaling relation for ν<sub>max</sub> (A. Serenelli et al. 2017, in preparation), (2) masses computed using the scaling relations, but with the grid-interpolated correction to the Δν relation based on radial modes by Sharma et al. (2016), and (3) scaling relation masses using the Guggenberger et al. (2016) analytical correction to the Δν relation. We find that corrections generally have an effect less than 0.1 M<sub>☉</sub> (less than 30 K), and that these corrections do not substantially affect the metallicity dependence of the temperature offsets (see Figure 6). We also show that using a different procedure to measure the global seismic parameters (SYD, Huber et al. 2009) also does not substantially affect the results.

We note that asteroseismic mass estimates are sometimes adjusted to better agree with the properties expected from a grid of stellar models; this is referred to as grid-based modeling. However, we caution that grid-based modeling does not necessarily improve asteroseismic results if the underlying tracks are not correctly located. In the red giant case, naive grid-based modeling will change the inferred posterior temperatures to reduce the discrepancy between the observed properties and the model predictions, dragging the metal-poor stars to hotter posterior temperatures and metal-rich stars to cooler posterior temperatures. These changes average 50–100 K at [Fe/H] = ±0.5. If the quoted temperature uncertainties are reduced to prevent this from occurring, the grid-based modeling will instead alter the masses by about ±0.2 M<sub>☉</sub> at [Fe/H] = ±0.5 to reduce the discrepancy between the data and the models. These changes are completely consistent with our expectations, given the metallicity dependent temperature offset we discuss above, and we therefore chose not to use grid modeling based masses in this work, as it complicates the signal we are examining.

## 3.4. Spectroscopic Uncertainties

In this analysis, we use the updated APOGEE temperatures, which have been corrected to the photometric scale. There should not, therefore, be correlated errors between the temperatures and the metallicities measured by the APOGEE ASPCAP pipeline. We did, however, test whether substantial changes would occur in the measured temperatures and metallicities if we used the APOGEE grid of MARCS atmospheres, rather than the Kurucz atmospheres used in
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Plots of the difference in model temperature over a range of masses and gravities as a function of metallicity for different helium values (top) and atmosphere boundary conditions (bottom). The mean trend line of temperature offset vs. metallicity from Figure 3 is again shown here, as a dashed line.</figcaption></figure>

DR13, for the 27 stars in our sample with $T_{\text{eff}} < 4400$ K, where the MARCS grid is available (see Zamora et al. 2015). We find that this, on average, moves the temperatures about 12 K cooler and makes a star about 0.007 dex more metal-poor. However, these changes are small and not strongly metallicity-dependent; choice of spectroscopic model atmosphere is thus unlikely to contribute to our observed offset.

We have also checked that the temperature offset is still present in the Hawkins et al. (2016) reanalysis of the previous APOKASC data set. This gives us more confidence that it is not some peculiarity of the ASPCAP analysis that is causing the metallicity dependent temperature offset.

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Plots of the offset between the APOGEE temperature and the model temperature for four different mass scales. The top left plot uses scaling relation masses and gravities from the SYD pipeline (Huber et al. 2009), rather than the OCT pipeline. The top right panel uses masses computed from the OCT pipeline values, but with corrections as given in Guggenberger et al. (2016). The bottom left panel uses masses computed using grid-based modeling as in A. Serenelli et al. (2017, in preparation) using the actual APOGEE temperature errors. The bottom right panel uses masses computed using the Sharma et al. (2016) corrections. The temperature offset persists no matter the choice of seismic mass scale. The dashed lines indicate the best fit from Figure 3; solid lines indicate the best fit for each panel. Contours indicate the extent of 68 and 95% of the sample.</figcaption></figure>

One other major source of uncertainty in the spectroscopic parameters is the use of atmosphere models constructed assuming local thermodynamic equilibrium. It is well-known that accounting for the effects of non-local thermodynamic equilibrium (NLTE) can substantially change the measured stellar parameters (see, e.g., Asplund 2005). Investigations by Ruchti et al. (2013), for example, indicated that using an NLTE analysis of optical iron lines could change the measured temperature, metallicity, and gravity in a metallicity dependent way. Although there have been fewer studies of the effects of NLTE corrections to lines in the infrared (Zhang et al. 2016), and it is unclear how the full spectral fitting done for the APOGEE spectra would differ from the analysis of individual lines, we use the Ruchti et al. (2013) study to estimate the kind of parameter changes we might expect if NLTE effects were included. In our analysis, we are calibrating our temperatures to the photometric scale, which should be unaffected by spectroscopic systematics like NLTE effects. Similarly, because we are using seismic gravities, we are not concerned about the effect of NLTE corrections on the measured surface gravity. In terms of metallicity, which has not be calibrated to a fundamental scale, the Ruchti et al. (2013) analysis indicates that the measured metallicities change by less than 0.1 dex for high-metallicity stars like our APOKASC sample. We therefore choose to proceed with the assumption that our measured trend is not caused by spectroscopic uncertainties, but recommend further investigation into the effect of NLTE corrections on the APOGEE analysis.

## 3.5. Comparison with Other Data Sets

### 3.5.1. Star Clusters

Star clusters are a natural laboratory for checking the concordance between theoretical models and stellar data. In
principle, because clusters are coeval and the red giant branch lifetime of a star is relatively short, the red giants in a cluster should represent a homogeneous sample of a single known mass and composition. The offset between the cluster isochrone and the measured temperature at the gravity of each star can therefore be used to measure a temperature offset similar to our procedure for asteroseismic data. This approach is powerful but caution is needed. Some methods for determining color-temperature relationships, for example, or fitting for cluster parameters, involve minimizing deviations between isochrone and data (for example, by applying zero point offsets to align cluster CMDs with theory). Globular clusters are also not simple stellar populations, and deviations from a universal mix may not be captured in theoretical isochrones assuming a solar mixture with alpha enhancement (Beom et al. 2016). With these caveats in mind, there is a large set of APOGEE data for star clusters, and we have used it to test both our temperature scale (Figure 1) and the offsets between cluster data and theory (Figure 7).

For globular clusters, we took cluster members with APOGEE data from M2, M3, M5, M13, M71, and M107 as described in J. Holtzman et al. (2015; 2017, in preparation). We note that this list excludes clusters below a metallicity of −2.0 because of known uncertainties with the spectroscopic gravity corrections in such clusters. We also included only stars with first generation abundance patterns, as the ASPCAP fitting process (García Pérez et al. 2016) yields spurious correlated metallicity and temperature offsets for stars whose base model atmospheres differ significantly from the assumed base model. For the open cluster NGC 2420, we took the members from Souto et al. (2016); M67 members were taken from stars with asteroseismic and APOGEE data as tabulated in Stello et al. (2016, submitted). Finally, for the open clusters NGC 6819 and 6791, we took targets with asteroseismic and spectroscopic data from M. H. Pinsonneault et al. (2017, in preparation). In open clusters where evolutionary states where available, we used the relevant analysis to exclude red clump stars. In the globular clusters, we used the red clump cut defined by Bovy et al. (2014) to limit our sample to shell burning giants. We excluded stars whose measured APOGEE metallicity was more than three sigma away from the cluster mean, as well as stars outside the gravity range of our APOKASC sample (log(g) between 3.3 and 1.1). This limits contamination by AGB stars, and also eliminates any concern that the temperature offset may become gravity-dependent at very low gravities.

Metallicities and α-element enhancements reflect the mean abundances in the chosen samples (see Table 3). In the cases of M67, NGC 6819, and NGC 6791, red giant masses can be accurately predicted by eclipsing binary stars with measured masses on the upper main sequence. We assume a red giant branch mass of 1.63 $M_\odot$ for NGC 2420, consistent with Souto et al. (2016). A red giant branch mass of 0.85 $M_\odot$ was assumed for the old globular clusters. We note that the temperature offsets on the RGB are relatively insensitive to a change in mass (<10 K per 0.05 $M_\odot$), but they are more sensitive to a change in metallicity (~25 K per 0.05 dex). Given our sample sizes, these systematic errors are in most cases more significant than random errors, and their quadrature sum is approximately represented by the size of the diamonds in Figure 7. We assumed that all cluster members had the mean metallicity and [α/Fe] of the cluster and computed the offset between an evolutionary track with the cluster RGB mass and metallicity

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

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

**Figure 7.** Plot of the offset between the measured APOGEE and model temperatures for stars in clusters, using APOGEE metallicities in the top plot and literature metallicities from Holtzman et al. (2015) in the bottom plot. Large diamonds indicate cluster means, and show the approximate size of the systematic error bars given a mass uncertainty of 0.05 $M_\odot$ and a metallicity uncertainty of 0.05 dex. Error bars indicate the standard error on the mean, if it is larger than the systematic uncertainties. Individual stars are represented by smaller points. We note that many of these clusters were also used to check the APOGEE temperature calibration (red stars in Figure 1), so they are consistent with the photometric temperature scale. The solid line indicates the best fit in each plot, treating each cluster as a single point. The dashed line represents the mean trend with metallicity from Figure 3. Note that it is consistent with clusters in the metallicity range spanned by the APOKASC data, but inconsistent with the low-metallicity globular clusters (see text). We therefore caution against extrapolating these results to the low-metallicity regime. Because of the known uncertainties with the APOGEE metallicity scale (see Section 2) we show this figure with both the APOGEE (top) and literature metallicity scales (bottom).

and the individual stellar spectroscopic temperatures at the corrected spectroscopic gravity; the spectroscopic temperatures were corrected for the trends described in Section 2. The results are illustrated in Figure 7. Using the cluster data with APOGEE metallicities, we see no strong evidence for a trend across the full metallicity range, although there is some evidence that the
<table><caption>**Table 3**
Cluster Properties Used in this Paper</caption><thead><tr><th>Cluster</th><th>Mass</th><th>[Fe/H]</th><th>[α/Fe]</th><th>Lit. [Fe/H]</th></tr></thead><tbody><tr><td>M13</td><td>0.85</td><td>-1.49</td><td>0.22</td><td>-1.58</td></tr><tr><td>M2</td><td>0.85</td><td>-1.46</td><td>0.23</td><td>-1.66</td></tr><tr><td>M3</td><td>0.85</td><td>-1.41</td><td>0.18</td><td>-1.50</td></tr><tr><td>M5</td><td>0.85</td><td>-1.24</td><td>0.21</td><td>-1.33</td></tr><tr><td>M107</td><td>0.85</td><td>-0.90</td><td>0.29</td><td>-1.03</td></tr><tr><td>M71</td><td>0.85</td><td>-0.64</td><td>0.18</td><td>-0.82</td></tr><tr><td>N2420</td><td>1.63</td><td>-0.21</td><td>0.01</td><td>-0.13</td></tr><tr><td>M67</td><td>1.36</td><td>-0.01</td><td>0.00</td><td>0.06</td></tr><tr><td>N6819</td><td>1.63</td><td>0.04</td><td>0.00</td><td>0.16</td></tr><tr><td>N6791</td><td>1.15</td><td>0.31</td><td>0.09</td><td>0.37</td></tr></tbody></table>

Note. Literature metallicities come from Holtzman et al. (2015).

intermediate metallicity globular clusters are less offset from the tracks than the lower metallicity clusters. However, if we use literature metallicities from Holtzman et al. (2015), we see a trend that averages about 45 K per dex across the full metallicity range, with clusters in the APOKASC metallicity range (above [Fe/H] ~ -1), consistent with the trend seen in the APOKASC data. Given how dependent this analysis is on the fundamental metallicity scale, something this paper is not equipped to address, we tentatively conclude that the cluster sample is not inconsistent with our trend in the metal-rich domain, but that our data should not be extrapolated into the metal-poor domain without further work on the absolute metallicity calibration of APOGEE in that domain.

### 3.5.2. Hipparcos Giants

Another data set with all the parameters necessary for our analysis is the Massarotti et al. (2008) sample of red giants. These stars have distances from Hipparcos data as well as luminosities computed from their V-band absolute magnitudes combined with bolometric corrections from VandenBerg & Clem (2003). They do have heterogeneously determined temperatures (photometric), metallicities (mostly spectroscopic), and gravities (partially spectroscopic, partially determined through comparison with isochrones). The combination of these parameters allows the calculation of the mass of each star independently of stellar models, and these stars can then be analyzed in exactly the same manner as our APOKASC stars. Their temperature offsets versus metallicity are shown in Figure 8. We use only stars in the same mass range as the APOKASC sample, 0.6–2.6 $M_⊙$, and use the temperature and metallicity dependent gravity cuts from Bovy et al. (2014) to exclude red clump stars in the Massarotti sample, leaving 173 stars. Because these stars only have measured bulk metallicities rather than α-element enhancements, we assume that all stars in the sample have [α/Fe] = 0. Knowing that low-metallicity stars are more likely to be alpha enhanced, we expect the real temperature offsets at low metallicity to be slightly (~15 K) smaller than shown in Figure 8. We note that this slope is significantly larger than the one measured using the APOKASC data, this suggests that there could be correlated metallicity and temperature errors in this sample. We suspect that this could result from the heterogeneous determination of the stellar parameters. Given that the offset seen is much larger than expected, we find this test inconclusive as well and suggest that it should be repeated as radii and spectroscopy for stars in the recent Gaia release (Lindegren et al. 2016) become

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Plot of the offset between the measured and model temperatures for the Massarotti et al. (2008) spectroscopic sample of Hipparcos giants. Error bars represent binned data. The solid line indicates the best fit to the sample, with $\Delta T_{\text{eff}} = 388.3$ [Fe/H] + 3.29 K. The dashed line is our fit to the APOKASC data. We suspect that the difference in slope could be due to the heterogeneous derivation of properties in the Massarotti et al. (2008) work; the offset in the mean temperature difference is likely due to calibration choices.</figcaption></figure>

available. For the purposes of this paper, we will assume that the metallicity-dependent temperature offset measured in the APOKASC sample is real, rather than the result of additional unrecognized systematic errors.

## 4. Results

### 4.1. A Metallicity-dependent Mixing Length

Although there are a variety of model ingredients, we suspect that the assumption of a universal solar calibrated convective efficiency may be causing the observed correlation between temperature offset and metallicity. We explore what metallicity dependent changes to the mixing length α would need to be made to bring the models into agreement with the observations. Specifically, we calculate an “effective mixing length” for each star by comparing the offset between the measured $T_{\text{eff}}$ and the model predicted $T_{\text{eff}}$ for three different mixing lengths. We then quadratically interpolate between them to find the mixing length that produces a temperature offset of 0 K. The mixing length for red giants at solar metallicity is not solar, suggesting a temperature or evolutionary state dependence in the mixing length.

The results are shown in Figure 9, and indicate a metallicity-dependent mixing length. We note that the scatter in effective mixing length (0.065) is about half the uncertainty we would expect given the quoted temperature uncertainties alone (69 K, $\Delta\alpha \sim 0.1$), indicating that there is very little actual scatter in our metallicity-effective mixing length relation. We suggested in Section 3.1 that there were possible nonlinearities in the fit of the temperature offset with metallicity that also show up in the mixing length fit. Except for these possible nonlinearities, there do not seem to be strong residual correlations between the mixing length and the α-element enhancement, mass, gravity, or temperature although we can not conclusively rule out smaller correlations at this point. We therefore report that a mixing length of $\alpha_{\text{ML}} = 0.162$ [Fe/H] + 1.90 would best bring our YREC grid into agreement with the APOKASC data. Preliminary investigations indicate that the change in
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Top: mixing length predicted for the temperatures and gravities of the stars in our sample at the metallicities available from the Magic et al. (2015, α<sub>ML,⊙</sub> = 1.98) (cyan) and Trampedach et al. (2014, α<sub>ML,⊙</sub> = 1.764) (red) simulations. Middle: effective mixing length needed to match the observed data as a function of metallicity for our YREC models (α<sub>ML,⊙</sub> = 1.72). We have divided the data into fifty equally sized bins in metallicity, black error bars represent the mean and standard deviation of each bin. The best fit, α<sub>ML</sub> = 0.1612[Fe/H] + 1.9037, is shown as a solid line; the dot-dashed line indicates the best fit using the metallicity corrections discussed in Section 2. Gray bands indicate the range of predicted mixing lengths for our data set from the simulations. Bottom: mixing length predicted by the simulations (same color coding as top) for stars in our sample with solar metallicity (−0.1 < [Fe/H] < 0.1) compared to the effective mixing length required for these stars (black error bars, 30 equally sized bins).</figcaption></figure>

giant branch temperature with mixing length is also linear and similarly sized in the GARSTEC grid of models (Weiss & Schlattl 2008). We therefore suggest that, although the normalization will be different for different model grids, the slope is unlikely to change appreciably. Although the model atmospheres used to measure the temperatures also assume a mixing length, our tests indicate that the model spectra, and therefore the derived parameters, are not very sensitive to the mixing length assumed. A change in the mixing length of the atmosphere model of the magnitude discussed here will therefore not impact our results.

We note that detailed seismic modeling of dwarfs and subgiants suggests a range of mixing lengths that depend on temperature, composition, and gravity (see e.g., Bonaca et al. 2012; Metcalfe et al. 2014; Creevey et al. 2017; Silva Aguirre et al. 2017). However, none of these studies have more than a single star in the more evolved gravity ranges we are considering. We therefore assume that this difference from the strongly metallicity-dependent temperature offset that we suggest here tentatively supports the idea that our changes in the effective mixing length are actually correcting for metallicity dependent physics errors unique to giants. However, even if this is the case, it is still possible to mitigate the offset between the data and the models in order to compute accurate isochrone masses and ages from Gaia data by calibrating the effective mixing length.

### 4.1.1. Comparison with Convection Simulations

We also show in Figure 9 the range of mixing lengths predicted for the stars in our sample by the Stagger (Magic et al. 2015, α<sub>ML,⊙</sub> = 1.98) and the Trampedach et al. (2014, α<sub>ML,⊙</sub> = 1.764) grids of three-dimensional convection simulations. For each grid, we use the stars in our sample within 0.1 dex of the metalliticies at which the mixing lengths were provided, and we use the temperatures and gravities of those stars to compute a predicted mixing length (top panel). For comparison with the empirical results (middle panel), we mark in gray the full range of results for the Trampedach case. For the Magic sample, we show in gray a mean fit and a two sigma spread. We remind the reader that the Magic et al. (2015) grid uses a substantially different solar mixing length than the YREC grid does, and thus the agreement at solar metallicity might be a coincidence.

We note that the theoretical predictions from the simulations are inconsistent with the range of effective mixing lengths we fit using one-dimensional stellar models. Specifically, to explain the observed temperatures of the stars in our sample would require a much stronger metallicity dependence than the simulations predict (middle panel), at least in some regimes. We also compare the predictions of the simulations as a function of gravity for stars around solar metallicity (−0.1 < [Fe/H] < 0.1) to our effective mixing lengths (bottom panel). It is clear that the simulations predict a trend with gravity that is not observed in our sample, and that the spread of mixing lengths at fixed gravity is larger in the simulations than our sample would require to explain the observed effective temperatures. Disagreements with the predicted trend direction and spread have also been observed for dwarf stars whose mixing lengths were obtained from detailed seismic modeling (Silva Aguirre et al. 2017). Together,
these discrepancies could indicate either that there is physics not captured or correctly converted to the 1D-mixing length proxy from the 3D-simulations, or that the temperature offset between the data and the models is not entirely due to a change in convection properties on the giant branch, but rather due to other metallicity dependent errors in our stellar evolution model physics.

## 4.2. Implications for Age Measurements

This correlation between metallicity and temperature offset complicates the interpretation of the recent Gaia data. Combining a photometric temperature (and possibly metallicity) with the luminosity inferred from the parallax will not allow for accurate isochrone fitting of masses, and therefore ages, of red giants using standard stellar models. Given our observed temperature offset, we would expect isochrone masses computed from uncorrected models to be off by as much as $0.2 M_\odot$ (~50 K) at [Fe/H] of $\pm 0.5$, even in the optimistic case where the solar metallicity tracks have been properly calibrated. This corresponds to ages being incorrect by about a factor of two, which is several gigayears in the solar mass case. Moreover, low-metallicity stars will be inferred as less massive (and therefore older), whereas high-metallicity stars will be inferred as more massive (therefore younger), giving a shallower age-metallicity relation in the local galaxy. We show in Figure 10 the expected difference between a star's real age and the age that would be inferred using a solar mixing length isochrone, for stars whose actual and inferred masses fall within our grid of models, and whose actual age is less than 15 Gyr. We note that the age errors are particularly large for low-mass stars, and can become larger than five gigayears for low-mass, low-metallicity stars. We therefore emphasize the need for isochrones to be made at a variety of mixing lengths, which can be calibrated to correct for these effects.

## 4.3. Implications for Nucleosynthesis

If the temperature offset is caused by a real change in convection, this would have important implications for the stellar nucleosynthesis models of low- and intermediate-mass stars ($1 < M < 10 M_\odot$). This is because theoretical stellar nucleosynthesis models usually assume the same solar calibrated mixing length parameter for different metallicities and evolutionary stages, including the red giant branch and asymptotic giant branch phases (see e.g., Lugaro et al. 2012; Karakas & Lattanzio 2014). For example, a higher mixing length parameter produces a higher horizontal branch temperature and the predicted yields for massive AGB and super-AGB stars may vary significantly as a result (for some elements, up to a factor of 3; see e.g., Doherty et al. 2014).

## 4.4. Implications for Stars in the Instability Strip

We have noted here that changes in the mixing length are required to match giant branch observations. Previous work on subgiant stars has indicated a definite, although possibly smaller, metallicity dependence for the mixing length (Bonaca et al. 2012; Metcalfe et al. 2014). We suspect, therefore, that there is likely a metallicity dependence for the mixing length of all stars between the main sequence and the main giant branch locus, including those in the RR Lyrae and Cepheid instability strips. Changing the stellar mixing length affects the minimum masses of stars that develop the Cepheid loop, with the reduced

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

**Figure 10.** Top: the difference between the real age, assuming our mixing length formula, and the age that would have been inferred from a fixed solar mixing length isochrone of the correct metallicity for a 100 $L_\odot$ star (colors). We have removed points that required extrapolation outside our grid, and excluded stars whose actual age would be greater than 15 Gyrs, because they should not exist in the real Gaia sample. Bottom: a similar plot in fractional age, with contours and colors indicating percent and fractional errors respectively. We have truncated the color scheme at age errors of 100% for clarity. Above about two solar masses, low-metallicity stars can still be crossing the Hertzsprung gap as subgiants evolving toward the red giant branch locus at a 100 $L_\odot$, which reduces the impact of the giant branch temperature offset. Because of these large errors, ages should not be computed for Gaia stars using solar mixing length isochrones.

mixing length at low metallicity making it harder to reach the instability strip, and hence increasing the minimum Cepheid mass. Similarly, the reduced mixing length at low metallicity will make it hard for low-mass stars to reach the RR Lyrae instability strip, and hence will require even more mass loss than commonly invoked to explain the CMD of low-metallicity globular clusters.

## 4.5. Implications for Galaxies

The physical properties of extra-galactic unresolved stellar populations like star clusters and galaxies are obtained by comparing observational spectrophotometry with so-called evolutionary population synthesis models (e.g., Maraston 1998, 2005; Bruzual & Charlot 2003; Conroy et al. 2009, etc.). These models, which provide integrated spectral energy distributions and
mass-to-light ratios as a function of time, metallicity, initial mass function, etc., are based on stellar evolutionary tracks. Hence, they are affected by the input physics and the unknowns of stellar evolution. Relevant to this work, Maraston (2005) describes how the mixing length adopted in different evolutionary calculations affect the temperature-luminosity relation of the red giant branch (their Figure 9). Because the red giant branch contributes a large part of the bolometric energy (~30%–40%, Maraston 1998, 2005) of old stellar populations ($t >$ a few Gyr), the adopted mixing length impacts population synthesis models. For example, the theoretical $V-K$ color of the models can be affected by ~0.25 mag (their Figure 27). When population synthesis models are fit to data, the derived physical properties will depend on the model input physics. In Goddard et al. (2017), we find that MaNGA galaxies that are fit with a warmer red giant branch are older by ~2.5 Gyr around absolute ages of 10 Gyr than galaxies fit to a cooler giant branch. Given the possible offset from solar mixing length, the possible nonlinearity of the mixing length we propose here, and the correlation between metallicity and galaxy mass, we suggest that more work needs to be done to understand the impact of our temperature offset trend on our understanding of galaxies. Galaxy ages and metallicities are key to understanding galaxy formation and evolution in the universe, hence it is important to realize that the calibration of the mixing length carries cosmological implications.

**5. Why Was This Not Noticed Before?**

Given the strength of this effect, it seems surprising that it was not previously recognized. In fact, offsets from isochrones have been noticed both in metal-poor globular clusters (e.g., Brasseur et al. 2010; Cohen et al. 2015) and metal-rich bulge stars (Ness et al. 2013). Additionally, work on open clusters has indicated that, depending on the choice of models, isochrones can be too blue in some bands and too red in others (see e.g., Hayes et al. 2015). However, because each of these regimes is usually analyzed individually, using each author's choice of models and colors, this likely made it more difficult to identify a systematic offset between the models and the data. We suspect that the calibration of the color-temperature relations and the reddening assumptions could also be masking the offset. In low-metallicity clusters, for example, where the model giant branch was bluer than the data, it is likely that increased reddening or higher metallicity were assumed during the isochrone fitting, in order to make the models better match the data. When combined with the fact that most authors studying star clusters do qualitative, rather than quantitative, assessments of the accuracy of their giant branch locations after fitting to the location of the main sequence or main sequence turnoff, it is somewhat less surprising that these offsets of only about 100 K were not noticed previously.

We also suggest that the lack of attention paid to these offsets could be due partially to the fact that the effect is much more obvious in the theoretical gravity-temperature plane than in the observational plane of color and magnitude. One can see this in, for example, the recent work by Smiljanic et al. (2016) on the metal-rich cluster Trumpler 20. Their Figure 2 shows an isochrone that is reasonably close to the data in a color-magnitude diagram, but it is clear in their Figure 3 that the same isochrone is over a hundred Kelvin too cool at fixed gravity.

**6. Summary**

When comparing the APOKASC data to grids of stellar models, we see a clear, metallicity-dependent offset in temperature. Although the exact offset depends on the calibration of the data and the stellar models used, the trend with metallicity is consistent with $\Delta T_{\text{eff}} \gtrsim 100$ K dex$^{-1}$ for the models and calibrations we tested, with $\Delta T_{\text{eff,YREC}} = 93.1[\text{Fe/H}] + 107.5$ K for the YREC grid of models. We find no residual dependence of the offset with mass, surface gravity, temperature, or [$\alpha$/Fe] ratio. The offset is difficult to study in optical spectra and cluster data and we suggest that more work should be done in both regimes to determine whether this effect is in fact real and present in all data sets. The models can be brought into agreement with the data if a metallicity-dependent mixing length is used, with the required correction for the YREC models being $\alpha_{\text{ML,YREC}} = 0.161[\text{Fe/H}] + 1.90$. We note that this predicts a non-solar mixing length for solar metallicity giants and does not agree with the temperature, gravity, and metallicity trends predicted by three-dimensional convection simulations. We assert that quantification of this metallicity-dependent temperature offset is particularly important for the calculation of ages from the recent Gaia data, as uncorrected isochrones give ages that are wrong by about a factor of two, even at modest (0.5 dex) deviations from solar metallicity.

We would like to thank J. Choi, C. Conroy, and A. Dotter for helpful discussions that led to discoveries of systematic errors in the APOGEE DR13 temperatures. We would also like to thank J. Zinn, M. Ness, and C. Hayes for pointing out useful references.

Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS web site is www.sdss.org.

SDSS-IV is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, the Chilean Participation Group, the French Participation Group, Harvard-Smithsonian Center for Astrophysics, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatory of China, New Mexico State University, New York University, University of Notre Dame, Observatório Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.

J.T., M.H.P., and J.A.J. acknowledge support from NASA grant NNX15AF13G. J.A.J. also acknowledges support from NSF grant AST-1211673. The research leading to the presented results has received funding from the European Research
Council under the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no 338251 (StellarAges). G.S. and J.C.B. acknowledge the support of the Vanderbilt Office of the Provost through the Vanderbilt Initiative in Data-intensive Astrophysics (VIDA). Szabolcs Mészáros has been supported by the Premium Postdoctoral Research Program of the Hungarian Academy of Sciences, and by the Hungarian NKFI Grants K-119517 of the Hungarian National Research, Development and Innovation Office. D.A.G.H. was funded by the Ramón y Cajal fellowship number RYC-2013-14182. D.A.G.H. and O.Z. acknowledge support provided by the Spanish Ministry of Economy and Competitiveness (MINECO) under grant AYA-2014-58082-P. F.B. is supported by NASA through Hubble Fellowship grant \#HST-HF2-51335 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. R.A.G. and B.M. acknowledge the support of CNES. S.B. acknowledges support from NSF grant AST-1514676 and NASA grant NNX16A109G. Funding for *Kepler* Discovery Mission is provided by NASA's Science Mission Directorate. Funding for the Stellar Astrophysics Foundation (Grant DNRF106). V.S.A. acknowledges support from VILLUM FONDEN (research grant 10118). R.E.C. acknowledges funding through Gemini-CONICYT Project 32140007. A.S. acknowledges support from grants 2014-SGR-1458 (Generalitat de Catalunya), ESP2014-56003-R and ESP2015-66134-R (MINECO). 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. S.M. would like to acknowledge support from NASA grants NNX12AE17G and NNX15AF13G and NSF grant AST-1411685. C.A.P. is thankful to the Spanish Ministry of Economy and Competitiveness (MINECO) for funding (grant AYA2014-56359-P).

## References

Adelberger, E. G., García, A., Robertson, R. G. H., et al. 2011, *RvMP*, 83, 195

Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, *ApJS*, 219, 12

Asplund, M. 2005, *ARA&A*, 43, 481

Beom, M., Na, C., Ferguson, J. W., & Kim, Y.-C. 2016, *ApJ*, 826, 155

Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, *ApJ*, submitted (arXiv:1703.00052)

Böhm-Vitense, E. 1958, *ZAp*, 46, 108

Bonaca, A., Tanner, J. D., Basu, S., et al. 2012, *ApJL*, 755, L12

Bovy, J., Nidever, D. L., Rix, H.-W., et al. 2014, *ApJ*, 790, 127

Brasseur, C. M., Stetson, P. B., VandenBerg, D. A., et al. 2010, *AJ*, 140, 1672

Bressan, A., Marigo, P., Girardi, L., et al. 2012, *MNRAS*, 427, 127

Bruzual, G., & Charlot, S. 2003, *MNRAS*, 344, 1000

Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, *SoPh*, 268, 255

Castelli, F., & Kurucz, R. L. 2004, *arXiv:astro-ph/0405087*

Choi, J., Dotter, A., Conroy, C., et al. 2016, *ApJ*, 823, 102

Cohen, R. E., Hempel, M., Mauro, F., et al. 2015, *AJ*, 150, 176

Conroy, C., Gunn, J. E., & White, M. 2009, *ApJ*, 699, 486

Creevey, O., Metcalfe, T. S., Schultheis, M., et al. 2017, A&A, in press (arXiv:1612.08990)

Cyburt, R. H. 2004, *PhRvD*, 70, 023505

Dewitt, H. E., Graboske, H. C., & Cooper, M. S. 1973, *ApJ*, 181, 439

Doherty, C. L., Gil-Pons, P., Lau, H. H. B., Lattanzio, J. C., & Siess, L. 2014, *MNRAS*, 437, 195

Epstein, C. R., Elsworth, Y. P., Johnson, J. A., et al. 2014, *ApJL*, 785, L28

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

Freytag, B., & Salaris, M. 1999, *ApJL*, 513, L49

García Pérez, A. E., Allende Prieto, C., Holtzman, J. A., et al. 2016, *AJ*, 151, 144

Goddard, D., Thomas, D., Maraston, C., et al. 2017, *MNRAS*, 466, 4731

González Hernández, J. I., & Bonifacio, P. 2009, *A&A*, 497, 497

Graboske, H. C., Dewitt, H. E., Grossman, A. S., & Cooper, M. S. 1973, *ApJ*, 181, 457

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

Guggenberger, E., Hekker, S., Basu, S., & Bellinger, E. 2016, *MNRAS*, 460, 4277

Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, *AJ*, 131, 2332

Handberg, R., & Lund, M. N. 2014, *MNRAS*, 445, 2698

Hawkins, K., Masseron, T., Jofré, P., et al. 2016, *A&A*, 594, A43

Hayes, C. R., Friel, E. D., Slack, T. J., & Boberg, O. M. 2015, *AJ*, 150, 200

Hekker, S., Broomhall, A.-M., Chaplin, W. J., et al. 2010, *MNRAS*, 402, 2049

Holtzman, J. A., Shetrone, M., Johnson, J. A., et al. 2015, *AJ*, 150, 148

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

Huber, D., Stello, D., Bedding, T. R., et al. 2009, *CoAst*, 160, 74

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

Karakas, A. I., & Lattanzio, J. C. 2014, *PASA*, 31, e030

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

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

Kurucz, R. L. 1997, in IAU Symp. 189, ed. T. R. Bedding, A. J. Booth, & J. Davis (Cambridge: Cambridge Univ. Press), 217

Lindegren, L., Lammers, U., Bastian, U., et al. 2016, *A&A*, 595, A4

Lugaro, M., Doherty, C. L., Karakas, A. I., et al. 2012, *M&PS*, 47, 1998

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

Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2015, *arXiv:1509.05420*

Mamajek, E. E. 2012, *ApJL*, 754, L20

Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, *ApJ*, 804, 64

Maraston, C. 1998, *MNRAS*, 300, 872

Maraston, C. 2005, *MNRAS*, 362, 799

Marigo, P., & Aringer, B. 2009, *A&A*, 508, 1539

Massarotti, A., Latham, D. W., Stefanik, R. P., & Fogel, J. 2008, *AJ*, 135, 209

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

Ness, M., Freeman, K., Athanassoula, E., et al. 2013, *MNRAS*, 430, 836

Nidever, D. L., Holtzman, J. A., Allende Prieto, C., et al. 2015, *AJ*, 150, 173

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

Pinsonneault, M. H., Kawaler, S. D., Sofia, S., & Demarque, P. 1989, *ApJ*, 338, 424

Ramírez, I., & Meléndez, J. 2005, *ApJ*, 626, 465

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

Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, *ApJ*, 456, 902

Ruchti, G. R., Bergemann, M., Serenelli, A., Casagrande, L., & Lind, K. 2013, *MNRAS*, 429, 126

Saio, H., Wood, P. R., Takayama, M., & Ita, Y. 2015, *MNRAS*, 452, 3863

Salaris, M., Chieffi, A., & Straniero, O. 1993, *ApJ*, 414, 580

Salpeter, E. E. 1954, *AuJPh*, 7, 373

Sharma, S., Stello, D., Bland-Hawthorn, J., Huber, D., & Bedding, T. R. 2016, *ApJ*, 822, 15

Silva Aguirre, V., Lund, M. N., Antia, H. M., et al. 2017, *ApJ*, 835, 173

Smiljanic, R., Franciosini, E., Randich, S., et al. 2016, *A&A*, 591, A62

Souto, D., Cunha, K., Smith, V., et al. 2016, *ApJ*, 830, 35

Takeda, Y., Tajitsu, A., Sato, B., et al. 2016, *MNRAS*, 457, 4454

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

van Saders, J. L., & Pinsonneault, M. H. 2012, *ApJ*, 746, 16

VandenBerg, D. A., & Clem, J. L. 2003, *AJ*, 126, 778

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

White, T. R., Bedding, T. R., Stello, D., et al. 2011, *ApJL*, 742, L3

Wu, X. S., Alexeeva, S., Mashonkina, L., et al. 2015, *A&A*, 577, A134

Zamora, O., García-Hernández, D. A., Allende Prieto, C., et al. 2015, *AJ*, 149, 181

Zasowski, G., Johnson, J. A., Frinchaboy, P. M., et al. 2013, *AJ*, 146, 81

Zhang, J., Shi, J., Pan, K., Allende Prieto, C., & Liu, C. 2016, *ApJ*, 833, 137