Free Access
Issue
A&A
Volume 597, January 2017
Article Number L3
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201629938
Published online 22 December 2016

© ESO, 2016

1. Introduction

In order to decipher stellar physics and Galactic populations, current stellar surveys such as Gaia-ESO (Gilmore et al. 2012), GALAH (De Silva et al. 2015), LAMOST (Luo et al. 2015) or SDSS/APOGEE (SDSS Collaboration 2016) record several hundred thousand stellar spectra. Due to the large number of observations, this naturally leads to the development of automatic procedures (e.g. Smiljanic et al. 2014; Ness et al. 2015; García Pérez et al. 2016; Masseron et al. 2016) to conduct their analysis. It was quickly realised that one of the most crucial aspect of those surveys resides in evaluating the performances of their internal procedures. The universally adopted method to evaluate the pipeline performances consists of comparing the results to the values of well-established reference stars. There are mainly three kinds of reference stars: the so-called benchmark stars (Heiter et al. 2015), stars in globular and open cluster and stars with asteroseismic data. A recent set of studies (e.g. Mészáros et al. 2013; Holtzman et al. 2015, and reproduced in Fig. 1) have shown a puzzling feature when comparing the surface gravity (log g) spectroscopic and the asteroseismic log g in the APOKASC dataset. Indeed, this figure shows a bimodal discrepancy in surface gravity (log g) obtained by spectroscopic and seismic methods depending on whether the star is at the core helium burning phase or at the H shell burning phase. However, it is not expected that the photosphere is influenced by ongoing internal nucleosynthesis. In other words, no distinct spectroscopic log g is expected between red clump (RC) stars and red giant branch (RGB) stars. Therefore, assuming this log g phenomenon is indeed a photospheric effect, understanding it represents an extraordinary opportunity to unambiguously distinguish the stellar evolutionary status of red giants from spectroscopy alone.

In this Letter, we attempt to narrow down the spectroscopic origin of this effect by constraining and testing unexplored spectroscopic variables such as abundances and isotopic ratios.

thumbnail Fig. 1

Difference in log g derived in the uncalibrated APOGEE DR13 log g and the log g obtained by asteroseismology. The sample is divided in three categories with three different symbols: RC stars as red points, lower RGB as green squares, and blue triangles as upper RGB.

2. Data

To carry out our analysis, we used the APOKASC sample (Pinsonneault et al. 2014). To build up the sample, we cross-matched the APOGEE DR13 uncalibrated stellar parameters from SDSS Collaboration (2016) with the asteroseismic masses and stellar parameters (Teff,log g) from Pinsonneault et al. (2014), as well as the evolutionary status (Elsworth, in prep.). The difference in the log g derived from spectroscopy and seismology as a function of asteroseismic log g for this sample is shown in Fig. 1. We divide the sample into three categories according to asteroseismology: the RC stars, the low RGB stars, defined as RGB stars with log g> 2.1, and upper RGB stars, RGB stars with log g< 2.1. Note that, while we always display these groups, we do not discuss extensively this latter category of stars because those are significantly cooler stars and are possibly affected by other spectroscopic biases.

thumbnail Fig. 2

Comparison of radius obtained by asteroseismology and by parallax measurement for a subsample of our stars.

3. Results and discussion

First, we found evidence that the log g discrepancy is rather inherent to the spectral analysis. While an error in effective temperature could impact the seismic log g determination, the required change in Teff to achieve a 0.2 offset in log g would be unrealistic (~4000 K by extrapolating from the study of Morel et al. 2014). Therefore, another source of systematic error in the determination of the seismic log g has to be explored. Thanks to the recent Gaia release (Gaia Collaboration 2016; Lindegren et al. 2016), we have compared the radius values obtained using asteroseismology and those obtained using parallax. To compute the latter, we use the effective temperature as presented in the previous section as well as the bolometric corrections from Jordi et al. (2010). In Fig. 2, we compare the radius obtained via the two methods and note that both RGB and RC agree well. Therefore, we conclude that the asteroseismic analysis provides at least relatively consistent radii between RGB and RC stars. Thus, although the seismic mass may also be questioned (Constantino et al. 2015; Gaulme et al. 2016), we consider that the log g provided by asteroseismology is not driving the log g discrepancy between RGB and RC stars.

thumbnail Fig. 3

As in Fig. 1, but as a function of stellar mass instead of log g. The lines drawn are least square fit of the data points.

Holtzman et al. (in prep.) illustrated a metallicity dependency in the log g discrepancy. However, metallicity dependent parameter gradients in red giant stars may in fact underline a mass dependent gradient (Masseron et al. 2017). Therefore we argue that it may be more relevant to study the mass dependency rather than metallicity. In Fig. 3, we present the same discrepancy in log g as in Fig. 1 but as a function of stellar mass. While the largest discrepancy between the spectroscopic and seismic (Δlog g) are with the RCs, the most interesting new feature of this diagram is the correlation of Δlog g with stellar mass for low RGB stars. The missing ingredient in the spectral analysis that leads to the discrepancy in log g between RGB and RC stars seems also tied to the only event that low RGB stars have undergone: the first dredge up. According to low-mass stellar evolution theory, the first dredge-up brings the product of H burning to the surface during the main sequence, and would thus affect H, He, Li, Be, B, C, N, 12C /13C, 14N /15N, 16O /17O, and 16O /18O (Iben 1967; Kippenhahn & Weigert 1990; Charbonnel 1994).

thumbnail Fig. 4

Normalised H-band spectrum for a star with Teff = 4500 K, log g = 2.0, [C/H] =−0.3, [N/H] = + 0.4 and solar composition model ([M/H] = 0, [He/H] = 0, [Li/H] = 0, [Be/H] = 0, [B/H] = 0, 12C /13C = 90, 14N /15N = 272, 16O /17O = 2625, 16O /18O = 498; top panel). The difference in normalised flux with the same parameters except, respectively, from upper to lower panel, log g = 2.5, [He/H] = 0.1, [Li/H] =−1.0, [Be/H] =−1.4, [B/H] =−2.7, 12C /13C = 5, 14N /15N = 100 000, 16O /17O = 60, and 16O /18O = 25 000.

While C and N abundances are already taken into account in the APOGEE analysis (with the ASPCAP pipeline, García Pérez et al. 2016), we show in Fig. 4 the impact on the spectrum of the other elements. Not surprisingly, Li and Be do not contribute to the spectrum and B shows only two very weak lines. The isotopic ratios of N and O contribute via the multiple CN, CO and OH molecular features over the wavelength range concerned, albeit relatively insignificantly given that their equilibrium values are very large (14N /15N = 100 000, 16O /17O = 60, and 16O /18O = 25 000, Boothroyd & Sackmann 1999). Therefore, the only remaining parameters that standard analysis did not take into account are He abundance and carbon isotopic ratio.

Helium appears to be a promising candidate as it impacts continuum opacities and mean molecular weight. Consequently, a change in He abundance alters the line profiles of H lines, line strength (and in particular molecular lines), ionisation equilibrium and molecular equilibrium (Boehm-Vitense 1979), especially around 1.5 μm where H opacity reaches a minimum. However, the expected surface abundances along the red giant evolution are relatively small (typically [He/H] ≈ + 0.05− + 0.1, Lagarde et al. 2012, corresponding to Y ≈ 0.27–0.30), yielding a relatively weak contribution compared to a change of log g of 0.5 dex. We note, however, that in the particular case of He, one should also evaluate the change in the atmospheric structure. We do not, however, have the means to achieve that in this work.

In contrast, a similar amplitude in the variation of the spectrum is obtained when changing the 12C /13C ratio from nearly solar values to typical RC values. In fact, because the total number of carbon atoms (N(12C) + N(13C)) needs to be conserved in the spectral synthesis, the different spectrum is an addition of many weaker 12C features and slightly enhanced 13C features. However, those variations in flux do not systematically appear at the same wavelengths between a change in 12C /13C and a change in log g, and to precisely estimate its impact on stellar parameters using only the observation of the spectrum is not straightforward.

thumbnail Fig. 5

Theoretical surface 12C /13C ratio values for the APOKASC sample as a function of their asteroseismic mass derived from the stellar evolution models of Lagarde et al. (2012). Those values are used as input for computation of the spectral grid.

To more quantitatively evaluate the impact of the change of 12C /13C ratio on the stellar parameters, we built a grid of synthetic spectra matching the APOKASC sample using Turbospectrum (Plez 2012). We adopted the same stellar parameters and atomic and molecular linelists as described in Hawkins et al. (2016). We computed the grid for two 12C /13C isotopic ratios: a generic value (15 as assumed for all giants in the ASPCAP grid) and a theoretical value for each star. This theoretical 12C /13C value is derived by interpolating its observed mass and Teff within the theoretical grid of Lagarde et al. (2012). The derived theoretical values for the satrs in our sample are shown in Fig. 5. Note that, because the metallicity of the models are too coarse, we adopted the closest models to the star’s metallicity. Knowing that the main sample covers a metallicity range of −0.8 < [Fe/H] < 0.5, we chose the Z = 0.004 models for stars with [Fe/H] < −0.2 and the Z = 0.014 models for the stars with [Fe/H] > −0.2. This results in a bimodal trend for low RGB and RC stars in Fig. 5 that are nonetheless in the same direction. We then ran the Cannon (a pipeline that employs a data-driven approach to determine stellar parameters and chemical abundances, Ness et al. 2015) on this synthetic grid to derive their stellar parameters using the uncalibrated DR13 stellar parameters and C, N and α element abundances for a sample of 5000 stars as our training set. This training set was employed using the same criteria as in Casey et al. (2016) with an additional criterion that the uncalibrated parameters and C and N abundances must also be known. Once the cannon was trained, the stellar parameters (Teff,log g, [M/H], [α/ Fe], [C/H], and [N/H]) of the spectra could be inferred.

thumbnail Fig. 6

Resulting difference in stellar parameter determination (Teff,log g, [C/H], and [N/H]) obtained by the Cannon on synthetic spectra when assuming a 12C /13C ratio of 15 and the stellar evolution model values described in Lagarde et al. (2012), and illustrated in Fig. 5.

Figure 6 shows the difference in the stellar parameters obtained by the Cannon when changing the 12C /13C ratio. First, it is interesting in this figure to point out the similarities with Fig. 3. In particular, the RC and RGB stars are easily distinguishable. This is expected as they have distinct 12C /13C ratios (Fig. 5). This distinct 12C /13C ratio between the low RGBs and the RCs is explained by the occurrence of non-canonical extra mixing between the two stellar evolutionary stages. During this extra mixing along the upper RGB branch, CN-cycled material is pulled up to the surface decreasing the 12C /13C ratio. It is worth noticing that the extra mixing is more effective at lower masses. However, the overall change between RGB and RC in log g due to the change in 12C /13C ratio is approximately 0.05 dex, which is approximately four time larger in the observed APOGEE spectra (Fig. 3). Therefore, we conclude that the missing spectroscopic factor that could explain the log g problem is probably not 13C. We also demonstrate earlier that it is even less likely to be He, Li, Be, B, nitrogen or oxygen isotopic ratios.

In any case, we observe that a change in 12C /13C ratio impacts not only log g but all stellar parameters. This may indicate that while 13C does contribute to the discrepancy in Δlog g, it is probably not the only parameter. We stress here that our procedure is not the exact pipeline employed by the APOGEE team and thus we would suggest the collaboration considers the isotopic ratio. Nonetheless, our results demonstrate that including 12C /13C ratio in the spectral analysis of red giants is recommended for an improvement in the determination of stellar parameters. Furthermore, when considering, notably, the N abundance, this may lead to systematic effects that can be as large as the quoted precision.

4. Conclusions

We show that the discrepancy between the log g derived spectroscopically in the APOGEE data and the log g derived by

asteroseismology seems related to CN cycling in red giants. If this problem is due to a missing parameter in the spectral analysis, it is very likely to be 13C, but when simulating the 13C effect on stellar parameters, the obtained shift in Δlog g is not as large as the observed one. Naturally, this may be a result of the fact that we use a different procedure than the ASPCAP pipeline. One first, obvious, required future test is to run the APOGEE pipeline with a grid that includes the 12C /13C ratio as an extra variable. Another direction is to measure 12C /13C ratio in a sample of those giants. However, it is difficult to find strong enough lines to measure in the APOGEE spectra. Counterpart optical spectra are thus required.

We believe that, in this letter, we have explored several spectroscopic effects that may explain the presence of an offset in log g. Nevertheless, if the present results are confirmed, this may question assumptions made concerning the yields from stellar evolution and the log g from asteroseismology.

Acknowledgments

This work was partly supported by the European Union FP7 programme through ERC grant number 320360. K.H. thanks the Marshall Scholarship program and the Simons Society of Fellows

References

  1. Boehm-Vitense, E. 1979, ApJ, 234, 521 [NASA ADS] [CrossRef] [Google Scholar]
  2. Boothroyd, A. I., & Sackmann, I.-J. 1999, ApJ, 510, 232 [Google Scholar]
  3. Casey, A. R., Hogg, D. W., Ness, M., et al. 2016, ApJ, submitted [arXiv:1603.03040] [Google Scholar]
  4. Charbonnel, C. 1994, A&A, 282, 811 [NASA ADS] [Google Scholar]
  5. Constantino, T., Campbell, S. W., Christensen-Dalsgaard, J., Lattanzio, J. C., & Stello, D. 2015, MNRAS, 452, 123 [NASA ADS] [CrossRef] [Google Scholar]
  6. De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  7. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. García Pérez, A. E., Allen de Prieto, C., Holtzman, J. A., et al. 2016, AJ, 151, 144 [NASA ADS] [CrossRef] [Google Scholar]
  9. Gaulme, P., McKeever, J., Jackiewicz, J., et al. 2016, ApJ, 832, 121 [NASA ADS] [CrossRef] [Google Scholar]
  10. Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
  11. Hawkins, K., Masseron, T., Jofré, P., et al. 2016, A&A, 594, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Heiter, U., Jofré, P., Gustafsson, B., et al. 2015, A&A, 582, A49 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
  13. Holtzman, J. A., Shetrone, M., Johnson, J. A., et al. 2015, AJ, 150, 148 [NASA ADS] [CrossRef] [Google Scholar]
  14. Iben, Jr., I. 1967, ApJ, 147, 624 [NASA ADS] [CrossRef] [Google Scholar]
  15. Jordi, C., Gebran, M., Carrasco, J. M., et al. 2010, A&A, 523, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Kippenhahn, R., & Weigert, A. 1990, in Stellar Structure and Evolution (Berlin, Heidelberg, New York: Springer Verlag), 192 [Google Scholar]
  17. Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Lindegren, L., Lammers, U., Bastian, U., et al. 2016, A&A, 595, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Luo, A.-L., Zhao, Y.-H., Zhao, G., et al. 2015, RA&A, 15, 1095 [Google Scholar]
  20. Masseron, T., Merle, T., & Hawkins, K. 2016, BACCHUS: Brussels Automatic Code for Characterizing High accUracy Spectra (Astrophysics Source Code Library) [Google Scholar]
  21. Masseron, T., Lagarde, N., Miglio, A., Elsworth, Y., & Gilmore, G. 2017, MNRAS, 464, 3021 [NASA ADS] [CrossRef] [Google Scholar]
  22. Mészáros, S., Holtzman, J., García Pérez, A. E., et al. 2013, AJ, 146, 133 [NASA ADS] [CrossRef] [Google Scholar]
  23. Morel, T., Miglio, A., Lagarde, N., et al. 2014, A&A, 564, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Ness, M., Hogg, D. W., Rix, H.-W., Ho, A. Y. Q., & Zasowski, G. 2015, ApJ, 808, 16 [NASA ADS] [CrossRef] [Google Scholar]
  25. Pinsonneault, M. H., Elsworth, Y., Epstein, C., et al. 2014, ApJS, 215, 19 [NASA ADS] [CrossRef] [Google Scholar]
  26. Plez, B. 2012, Turbospectrum: Code for spectral synthesis (Astrophysics Source Code Library) [Google Scholar]
  27. SDSS Collaboration 2016, ApJ, submitted [arXiv:1608.02013] [Google Scholar]
  28. Smiljanic, R., Korn, A. J., Bergemann, M., et al. 2014, A&A, 570, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Figures

thumbnail Fig. 1

Difference in log g derived in the uncalibrated APOGEE DR13 log g and the log g obtained by asteroseismology. The sample is divided in three categories with three different symbols: RC stars as red points, lower RGB as green squares, and blue triangles as upper RGB.

In the text
thumbnail Fig. 2

Comparison of radius obtained by asteroseismology and by parallax measurement for a subsample of our stars.

In the text
thumbnail Fig. 3

As in Fig. 1, but as a function of stellar mass instead of log g. The lines drawn are least square fit of the data points.

In the text
thumbnail Fig. 4

Normalised H-band spectrum for a star with Teff = 4500 K, log g = 2.0, [C/H] =−0.3, [N/H] = + 0.4 and solar composition model ([M/H] = 0, [He/H] = 0, [Li/H] = 0, [Be/H] = 0, [B/H] = 0, 12C /13C = 90, 14N /15N = 272, 16O /17O = 2625, 16O /18O = 498; top panel). The difference in normalised flux with the same parameters except, respectively, from upper to lower panel, log g = 2.5, [He/H] = 0.1, [Li/H] =−1.0, [Be/H] =−1.4, [B/H] =−2.7, 12C /13C = 5, 14N /15N = 100 000, 16O /17O = 60, and 16O /18O = 25 000.

In the text
thumbnail Fig. 5

Theoretical surface 12C /13C ratio values for the APOKASC sample as a function of their asteroseismic mass derived from the stellar evolution models of Lagarde et al. (2012). Those values are used as input for computation of the spectral grid.

In the text
thumbnail Fig. 6

Resulting difference in stellar parameter determination (Teff,log g, [C/H], and [N/H]) obtained by the Cannon on synthetic spectra when assuming a 12C /13C ratio of 15 and the stellar evolution model values described in Lagarde et al. (2012), and illustrated in Fig. 5.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.