Probing 3D and NLTE models using APOGEE observations of globular cluster stars

¹ Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
e-mail: thomas.masseron@iac.es
² Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
³ ELTE Eötvös Loránd University, Gothard Astrophysical Observatory, 9700 Szombathely, Szent Imre H. st. 112, Hungary

Received: 21 September 2020
Accepted: 30 December 2020

Abstract

Context. Hydrodynamical (or 3D) and non-local thermodynamic equilibrium (NLTE) effects are known to affect abundance analyses. However, there are very few observational abundance tests of 3D and NLTE models.

Aims. We developed a new way of testing the abundance predictions of 3D and NLTE models, taking advantage of large spectroscopic survey data.

Methods. We use a line-by-line analysis of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) spectra (H band) with the Brussels Automatic Code for Characterizing High accUracy Spectra (BACCHUS). We compute line-by-line abundances of Mg, Si, Ca, and Fe for a large number of globular cluster K giants in the APOGEE survey. We compare this line-by-line analysis against NLTE and 3D predictions.

Results. While the 1D–NLTE models provide corrections in the right direction, there are quantitative discrepancies between different models. We observe a better agreement with the data for the models including reliable collisional cross-sections. The agreement between data and models is not always satisfactory when the 3D spectra are computed in LTE. However, we note that for a fair comparison, 3D corrections should be computed with self-consistently derived stellar parameters, and not on 1D models with identical stellar parameters. Finally, we focus on 3D and NLTE effects on Fe lines in the H band, where we observe a systematic difference in abundance relative to the value from the optical. Our results suggest that the metallicities obtained from the H band are more accurate in metal-poor giants.

Conclusions. Current 1D–NLTE models provide reliable abundance corrections, but only when the atom data and collisional cross-sections are accurate and complete. Therefore, we call for more atomic data for NLTE calculations. In contrast, we show that 3D corrections in LTE conditions are often not accurate enough, thus confirming that 3D abundance corrections are only valid when NLTE is taken into account. Consequently, more extended self-consistent 3D–NLTE computations need to be made. The method we have developed for testing 3D and NLTE models could be extended to other lines and elements, and is particularly suited for large spectroscopic surveys.