Gravity and limb-darkening coefficients for compact stars: DA, DB, and DBA eclipsing white dwarfs star
Context. The distribution of the specific intensity over the stellar disk is an essential tool for modeling the light curves in eclipsing binaries, planetary transits, and stellar diameters through interferometric techniques, line profiles in rotating stars, gravitational microlensing, etc. However, the available theoretical calculations are mostly restricted to stars on the main sequence or the giant branch, and very few calculations are available for compact stars. Aims. The main objective of the present work is to extend these investigations by computing the gravity and limb-darkening coefficients for white dwarf atmosphere models with hydrogen, helium, or mixed compositions (types DA, DB, and DBA). Methods. We computed gravity and limb-darkening coefficients for DA, DB, and DBA white dwarfs atmosphere models, covering the transmission curves of the Sloan, UBVRI, Kepler, TESS, and Gaia photometric systems. Specific calculations for the HiPERCAM instrument were also carried out. For all calculations of the limb-darkening coefficients we used the least-squares method. Concerning the effects of tidal and rotational distortions, we also computed for the first time the gravity-darkening coefficients y(lambda) for white dwarfs using the same models of stellar atmospheres as in the case of limb-darkening. A more general differential equation was introduced to derive these quantities, including the partial derivative ( partial differential ln I-o(lambda)/ partial differential ln g)(Teff). Results. Six laws were adopted to describe the specific intensity distribution: linear, quadratic, square root, logarithmic, power-2, and a more general one with four coefficients. The computations are presented for the chemical compositions log[H/He] = -10.0 (DB), -2.0 (DBA) and He/H = 0 (DA), with log g varying between 5.0 and 9.5 and effective temperatures between 3750 and 100 000 K. For effective temperatures higher than 40 000 K, the models were also computed adopting nonlocal thermal equilibirum (DA). The adopted mixing-length parameters are ML2/alpha = 0.8 (DA case) and 1.25 (DB and DBA). The results are presented in the form of 112 tables. Additional calculations, such as for other photometric systems and/or different values of log[H/He], log g, and T-eff can be performed upon request.© ESO 2020 ; We thank the anonymous referee for helpful comments that have improved the manuscript. The Spanish MEC (AYA2015-71718-R and ESP2017-87676-C5-2-R) is gratefully acknowledged for its support during the development of this work. A.C. also acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). The research leading to these results has received funding from the European Research Council under the European Union's Horizon 2020 research and innovation programme n. 677706 (WD3D). This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France, and of NASA's Astrophysics Data System Abstract Service. ; Peer reviewed