How realistic are solar model atmospheres?

¹ Research School of Astronomy and Astrophysics, Australian National University, Cotter Rd., Weston, ACT 2611, Australia
e-mail: tiago.pereira@astro.uio.no
² Institute for Theoretical Astrophysics, University of Oslo, PO Box 1029, 0315 Blindern, Norway
³ Max-Planck-Institut für Astrophysik, Postfach 1317, 85741 Garching b. München, Germany
⁴ JILA, University of Colorado and National Institute of Standards and Technology, 440 UCB, Boulder, CO 80309, USA

Received: 2 February 2013
Accepted: 17 April 2013

Abstract

Context. Recently, new solar model atmospheres have been developed to replace classical 1D local thermodynamical equilibrium (LTE) hydrostatic models and used to for example derive the solar chemical composition.

Aims. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D magnetohydrodynamical (MHD) model (with an imposed 10 mT vertical magnetic field), 1D NLTE and LTE models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & Müller.

Methods. We confronted the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLVs), absolute continuum fluxes, and the wings of hydrogen lines. We also tested the 3D models for the intensity distribution of the granulation and spectral line shapes.

Results. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & Müller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well.

Conclusions. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D hydrostatic models (theoretical and semi-empirical) are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.