A new method to compute limb-darkening coefficients for stellar atmosphere models with spherical symmetry: the space missions TESS, Kepler, CoRoT, and MOST^⋆,⋆⋆

¹ Instituto de Astrofísica de Andalucía CSIC, Apartado 3004, 18080 Granada, Spain
e-mail: claret@iaa.es
² Dept. Física Teórica y del Cosmos, Universidad de Granada, Campus de Fuentenueva s/n, 10871 Granada, Spain

Received: 20 March 2018
Accepted: 15 April 2018

Abstract

Aims. One of the biggest problems we can encounter while dealing with the limb-darkening coefficients for stellar atmospheric models with spherical symmetry is the difficulty of adjusting both the limb and the central parts simultaneously. In particular, the regions near the drop-offs are not well reproduced for most models, depending on Teff, log g, or wavelength. Even if the law with four terms is used, these disagreements still persist. Here we introduce a new method that considerably improves the description of both the limb and the central parts and that will allow users to test models of stellar atmospheres with spherical symmetry more accurately in environments such as exoplanetary transits, eclipsing binaries, etc.

Methods. The method introduced here is simple. Instead of considering all the μ points in the adjustment, as is traditional, we consider only the points until the drop-off (μ_cri) of each model. From this point, we impose a condition I(μ)/I(1) = 0. All calculations were performed by adopting the least-squares method.

Results. The resulting coefficients using this new method reproduce the intensity distribution of the PHOENIX spherical models (COND and DRIFT) quite well for the photometric systems of the space missions TESS, Kepler, CoRoT, and MOST. The calculations cover the following ranges of local gravity and effective temperatures: 2.5 ≤ log g ≤ 6.0 and 1500 K ≤ T_eff ≤ 12 000 K. The new spherical coefficients can easily be adapted to the most commonly used light curve synthesis codes.