MARCS model atmospheres

TL;DR

The paper reviews the development of the MARCS stellar atmosphere models from 1975 to 2008, discussing their physical data, accuracy, and applications in spectral analysis of stars, especially red giants.

Contribution

It provides a comprehensive overview of the MARCS model evolution, improvements in physical data, and initial spectral comparisons, highlighting the model's accuracy and areas for future work.

Findings

Errors in spectral energy distribution are within a few percent for most models.

High-resolution spectra are planned for future publication.

Scattering in TiO lines affects atmospheric cooling but conflicts with some observations.

Abstract

In this review presented at the Symposium A stellar journey in Uppsala, June 2008, I give my account of the historical development of the MARCS code from the first version published in 1975 and its premises to the 2008 grid. It is shown that the primary driver for the development team is the science that can be done with the models, and that they constantly strive to include the best possible physical data. A few preliminary comparisons of M star model spectra to spectrophotometric observations are presented. Particular results related to opacity effects are discussed. The size of errors in the spectral energy distribution (SED) and model thermal stratification are estimated for different densities of the wavelength sampling. The number of points used in the MARCS 2008 grid (108000) is large enough to ensure errors of only a few K in all models of the grid, except the optically very thin layers of metal-poor stars. Errors in SEDs may reach about 10% locally in the UV. The published sampled SEDs are thus appropriate to compute synthetic broad-band photometry, but higher resolution spectra will be computed in the near future and published as well on the MARCS site (marcs.astro.uu.se). Test model calculations with TiO line opacity accounted for in scattering show an important cooling of the upper atmospheric layers of red giants. Rough estimates of radiative and collisional time scales for electronic transitions of TiO indicate that scattering may well be the dominant mechanism in these lines. However models constructed with this hypothesis are incompatible with optical observations of TiO (Arcturus) or IR observations of OH (Betelgeuse), although they may succeed in explaining H2O line observations. More work is needed in that direction.