Issue 
A&A
Volume 573, January 2015



Article Number  A90  
Number of page(s)  10  
Section  Stellar atmospheres  
DOI  https://doi.org/10.1051/00046361/201423804  
Published online  23 December 2014 
The Staggergrid: A grid of 3D stellar atmosphere models
IV. Limb darkening coefficients^{⋆}
^{1}
MaxPlanckInstitut für Astrophysik,
KarlSchwarzschildStr. 1,
85741
Garching,
Germany
email:
magic@mpagarching.mpg.de
^{2}
Research School of Astronomy & Astrophysics, Cotter Road,
Weston
ACT
2611,
Australia
^{3}
Laboratoire Lagrange, UMR 7293, CNRS, Observatoire de la Côte
d’Azur, Université de Nice SophiaAntipolis, Nice, France
Received: 14 March 2014
Accepted: 9 June 2014
Aims. We compute the emergent stellar spectra from the UV to far infrared for different viewing angles using realistic 3D model atmospheres for a large range in stellar parameters to predict the stellar limb darkening.
Methods. We have computed full 3D LTE synthetic spectra based on 3D radiative hydrodynamic atmosphere models from the Staggergrid in the ranges: T_{eff} from 4000 to 7000 K, log g from 1.5 to 5.0, and [Fe/H], from −4.0 to +0.5. From the resulting intensities, we derived coefficients for the standard limb darkening laws considering a number of oftenused photometric filters. Furthermore, we calculated theoretical transit light curves, in order to quantify the differences between predictions by the widely used 1D model atmosphere and our 3D models.
Results. The 3D models are often found to predict steeper darkening towards the limb compared to the 1D models, mainly due to the temperature stratifications and temperature gradients being different in the 3D models compared to those predicted with 1D models based on the mixing length theory description of convective energy transport. The resulting differences in the transit light curves are rather small; however, these can be significant for highprecision observations of extrasolar transits, and are able to lower the residuals from the fits with 1D limb darkening profiles.
Conclusions. We advocate the use of the new limb darkening coefficients provided for the standard fourparameter nonlinear power law, which can fit the limb darkening more accurately than other choices.
Key words: convection / hydrodynamics / radiative transfer / stars: atmospheres / binaries: eclipsing / planetary systems
Full Table A.1 and the grid of spectra are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/573/A90, as well as at http://www.staggerstars.net
© ESO, 2014
1. Introduction
The emergent intensity across the surface of latetype stars diminishes gradually from the center of the stellar disk towards the edge (limb), since the optical depth depends on the angle of view. Rays crossing the stellar photosphere near the limb reach optical depth unity in layers at higher altitude and at typically lower densities and temperatures than rays crossing the stellar disk near the center. The intensity is very sensitive to the temperature, therefore, one observes darker brightness from the limb, which emerges from higher and cooler regions of the stellar atmosphere. This effect is known as limb darkening (Gray 2005). An accurate knowledge of the surface brightness distribution is essential for the analysis of light curves from stars with transiting objects in the line of sight, such as exoplanets and eclipsing stellar companions in binary systems. Furthermore, the precise determination of stellar angular diameters with stellar interferometry relies also on the theoretical limb darkening predictions (Davis et al. 2000). The variation in surface intensity with angular distance from the stellar disk center is usually expressed in the form of limb darkening laws (Claret 2000). Multiple functional basis have been used in the past, from simple linear to higher order nonlinear laws, in order to fit the surface brightness variations predicted by theoretical model atmospheres leading to socalled limb darkening coefficients (LDC). For example, the individual shape of a light curve for transiting exoplanets is important, because it contains information about the structure of the external layers of the occulted stellar object (e.g., Southworth 2008). The observed light curves are interpreted by comparisons with theoretical transit light curves that are based on limb darkening predictions arising from model atmospheres. More accurate theoretical atmosphere models will reduce the uncertainties in the comparison, and thereby improving the quality of the analysis in favor of other transitparameters like planettostar ratio or the inclination of the orbit. Also, the goodness of the transmissionspectroscopy of exoplanet atmospheres relies on the underlying theoretical atmospheres of the host stars (e.g., Seager & Sasselov 2000).
The first estimates of the intensity variation over the disk were performed with a simple linear law (Milne 1921). However, with theoretical 1D model atmospheres it was shown that a linear law is insufficient to describe the limb darkening of a real star adequately (e.g., van Hamme 1993). Then, various alternatives with a twoparameter law was introduced starting from a quadratic, over square root to a logarithmic, and finally an exponential law (e.g. see DiazCordoves et al. 1995; Claret et al. 1995). These restricted functional bases are only marginally accurate for a certain range in effective temperatures, therefore, Claret (2000) introduced a new nonlinear power law with four coefficients, which is powerful enough to fit the LDC for a broad range in stellar parameters, while conserving the flux to a high accuracy. Later on, limb darkening variations were fitted and provided for the community derived from extensive grids with the latest model atmospheres (e.g, MARCS, ATLAS and PHOENIX) for several broad band filters, e.g. the SDSS (Claret 2004), Kepler and CoRoT (Sing 2010). An extensive comparison of the various limb darkening laws has been performed by Southworth (2008). All of these developments revealed that a wellconsidered choice of an appropriate functional basis is mandatory for a precise description of the intensity variations.
The next step in improving the systematic errors prevailing in the predicted limb darkening laws was yielded in the underlying model atmospheres, since the limb darkening is mainly determined by the temperature gradient (see Knutson et al. 2007; Hayek et al. 2012). Therefore, flaws in the theoretical atmospheric temperature stratification will directly propagate into the predicted limb darkening. The hydrostatic 1D models make use of several simplifications, the most prominent one being the use of the mixing length theory to account for convective energy transport (BöhmVitense 1958). Cool latetype stars feature a convective envelope, thereby convective motions are present in the thin photospheric transition region due to overshooting of convective flows. These stars exhibit a typical granulation pattern in its emergent intensity due to inhomogeneities arising from the asymmetric up and downflowing stellar plasma. Therefore, only 3D atmosphere models are able to predict these properties accurately. With the advent of 3D atmosphere modeling (Nordlund 1982), which solves from firstprinciple the hydrodynamic equations coupled with a realistic radiative transfer, the deficiencies of the 1D models were revealed and quantified (e.g., Nordlund et al. 2009, and references therein). Comparisons of the 3D models with the Sun showed that these models can predict accurately the intensity distribution (Pereira et al. 2013), while 1D models overestimate the limb darkening of our resolved host star. Bigot et al. (2006) studied the limb darkening of α Centauri B by comparing its interferometrically observed visibility curves with theoretical predictions. The latter is sensitive to the limb darkening, and they found an significant improvement with the predictions from 3D models. Furthermore, Hayek et al. (2012) showed on the basis of the extremely accurately measured light curves of the transiting exoplanet HD 209458 that the intrinsic residuals of the 1D models can be resolved with the more realistic 3D model atmospheres. The largest differences were found close to the limb, hence during the ingress and egress of the transition. With 1D model predictions, the wellstudied closeorbit Jupiterlike transit planet HD 209458 exhibited priorly systematic residuals due to the simplified treatment of convection leading to insufficient temperature stratifications (see Knutson et al. 2007). For another wellstudied star, Procyon, and four K giants Chiavassa et al. (2010, 2012) could also find the limb darkening and stellar diameter predictions to be coherent with independent asteroseismic observations (see also Allende Prieto et al. 2002; Aufdenberg et al. 2005).
After the first detection of a Jupiterlike extra solar planet through radial velocity detection (Mayor & Queloz 1995), five years later, eventually a transiting exoplanet around a solarlike star was also found (Charbonneau et al. 2000). These spectacular landmark discoveries triggered literally a goldrush in the hunt for new exoplanets. With advanced satellite missions, like Kepler and CoRoT, nowadays up 1491 transiting extra solar planets have been detected (Wright et al. 2011). Both of the mentioned satellite missions operate in the visible spectral range, therefore, the effects of limb darkening are strong. These sophisticated observations evoke rightfully a demand in more accurate theoretical limb darkening predictions. In order to fulfill this call, we present in this work LDCs derived from realistic full 3D synthetic spectra based on a comprehensive grid of 3D RHD atmosphere models.
In Sects. 2 and 3, we explain the methods we utilized to obtain the LDC. Subsequently, the resulting theoretical limb darkening variations (Sect. 4) and transit light curves (Sect. 5) are presented and discussed. We compare our results with previous predictions from 1D ATLAS models in Sect. 6. Finally, we conclude our findings in Sect. 7.
2. 3D atmosphere models and 3D synthetic spectra
We have computed a large grid of realistic 3D atmosphere models (see Magic et al. 2013a, hereafter Paper I). We employed the Staggercode, a stateoftheart (magneto)hydrodynamic code that solves the timedependent equations for conservation of mass, momentum and energy. In the optically thin regime, the code solves the radiative transfer for the vertical direction and eight inclined rays along long characteristics (Nordlund 1982; Stein & Nordlund 1998). The original large set of wavelength points for the opacity sampling data are grouped together into 12 opacity bins in the socalled opacity binning method (Nordlund 1982; Skartlien 2000) while solving the radiative transfer. The Staggercode utilizes a realistic EOS (Mihalas 1970) and continuum and line opacities (Kurucz 1979 with subsequent updates; Gustafsson et al. 2008). We perform socalled “boxinthestar” simulations, where only a small representative volume is considered that accommodates the top of the convection zone and the extended photosphere. The vertical directions feature open boundaries, while the horizontal ones are periodic. The numerical resolution of the geometrical mesh is 240^{3}. The horizontal mesh is equidistant, while the vertical depth scale is optimized to resolve the photospheric transition region with an enhanced resolution, thereby exploiting the given resolution best possibly. The Staggergrid covers a wide range in stellar parameters with effective temperatures from 4000 K to 7000 K in steps of 500 K, surface gravities from 1.5 to 5.0 in steps of 0.5 and metallicities from −4.0 to +0.5 in steps of 1.0 below −1.0, and steps of 0.5 above. The range in surface gravity and metallicity for the Staggergrid covers the range of the planet hosting stars, while the effective temperature is slightly smaller. For further details on the model atmospheres we refer to Paper I. Furthermore, in Magic et al. (2013b) we explain in detail the horizontal averaging methods for the ⟨3D⟩ models. Furthermore, additional results on the Staggergrid are given in Magic et al. (2015, 2014) and Magic & Asplund (2014).
Based on the 3D RHD models from the Staggergrid, we computed a comprehensive library of full 3D synthetic spectra. Therefore, we used the Optim3Dcode (see Chiavassa et al. 2009, 2010, for further details), which is a postprocessing 3D radiative transfer code that assumes local thermodynamic equilibrium (). The code considers the realistic velocity field due to convective motions present in the 3D RHD simulations, thereby taking Doppler broadening and shifts into account. Optim3D employs pretabulated extinction coefficients that are the same as in the MARCS and Stagger codes, thereby accounting for continuous and sampled line opacities (Gustafsson et al. 2008). We assumed the latest solar composition by Asplund et al. (2009) consistently in Optim3D as in the 3D RHD simulations performed with the Staggercode. In contrast to the RHDcode, the spectral synthesis code computes the large number of wavelength points with N_{λ} = 105 767 explicitly, thereby raising the computational costs enormously. We achieve a wavelength resolution with a constant sampling rate of λ/ Δλ = 20 000, however, we cover a broad range with λ = 1010.0 to . We apply N_{t} = 10 snapshots for each simulation, while we keep the horizontal mesh resolution fixed with N_{xy} = 240. Four (equidistant) azimutal φangles are considered (N_{φ} = 4), while the centertolimb resolution is resolved higher with nine μangles besides the diskcenter (μ = 0.01, 0.05, 0.10, 0.20, 0.30, 0.50, 0.70, 0.80, 0.90, 1.00). After carrying out tests with more μangles, we found the ten chosen angles being sufficient to resolve the limb darkening accurately at lower computational costs. The strongest decline in the limb darkening is usually found towards the limb, therefore, we decided to resolve the limb more instead of equidistant μangles. As shown by Hayek et al. (2012) the numerical resolution of our 3D RHD models and the resolution for the spectral flux computations are sufficient to predict realistic observed limb darkening laws accurately. The synthetic spectra will be discussed in a separate work (Chiavassa et al., in prep.).
We show an overview of spatially resolved intensity maps with different inclined μ angles for a selection of distinct stellar parameters (see Fig. 1). These exhibit the typical granulation pattern of cool stars due to convection. The bright, bulk regions are the hotter upflowing granules, which are interspersed with the dark intergranular downdrafts. From the diskcenter towards the limb, the brightness is diminishing significantly, and the intensity contrast is also slightly dropping. Moreover, one can also obtain that bright features are often highly angle dependent.
In Fig. 2, we show a subset of the resulting averaged synthetic fluxes in the range for a number of dwarfs with solar metallicity. Furthermore, we show also the continuum fluxes as well, and one can discern spectral absorption features, the prominent one being the Balmer lines (indicated in the figure). For higher T_{eff} the continuum flux is increasing, while individual spectral absorption features are changing as well.
Fig. 1 Emergent monochromatic intensity at 500 nm shown for eight μ angles: 0.2, 0.3, 0.5, 0.7, 0.8, 0.9, 1.0 (from top to down, respectively) and for a selection of stars: mainsequence, turnoff, Kgiant, Kdwarf (from left to right, respectively) with solar metallicity. Furthermore, we indicated the normalized mean intensity and the intensity contrast (both in percent). 

Open with DEXTER 
Fig. 2 Spectral energy distribution vs. wavelength for models with log g = 4.5, [Fe/H] = 0.0 and various T_{eff} (red/orange colors). The continuum fluxes are also shown (solid black lines). 

Open with DEXTER 
3. Deriving the limb darkening
The variation of the inclination of the line of sight from the diskcenter is parameterized with the projected polar angle, μ = cosθ, where θ is the angle between the line of sight and the direction of the emergent radiation. Therefore, the diskcenter is depicted with μ = 1, while μ = 0 is the limb. The limb darkening law is expressed as the variation in intensity with μangle that is normalized to the diskcenter, i.e. . The resulting monochromatic intensity depends on the horizontal position x and y, the viewing angles φ and μ and the time, t, thus I_{λ}(x,y,μ,φ,t). In order to yield the mean monochromatic intensity ⟨I_{λ}⟩(μ) from the latter, we average the intensity first spatially, then over the azimutal angles, and finally over all time steps, i.e. As next, we compute the inclinationdependent total emergent intensity I(μ) by integrating the mean monochromatic intensity ⟨I_{λ}⟩ over all wavelength points with Then, the total surface brightness variation can be easily derived by normalizing the angular intensities with the diskcenter value, I_{μ}/I_{1}, and we can fit the various (biparametric) functional bases, which are the linear (Eq. (1)), quadratic (Eq. (2)), square root (Eq. (3)), logarithmic (Eq. (4)) and exponential (Eq. (5)) limb darkening law. Also the threeparameter nonlinear limb darkening law, (6)introduced by Sing (2010) is also considered. However, we recommend the use of the standard fourparameter nonlinear functional basis (Eq. (7)) introduced by Claret (2000), which is the default limb darkening law in the present study. The fourparameter power law is the fourth order Taylorseries expansion in μ^{1/2} given by (7)This functional basis conserves the flux to better than 0.05% (see Claret 2000). In order to fit the LDC, we applied the LevenbergMarquardt leastsquare minimization, since Claret (2000) showed that this fitting method performs best.
Precision of the functional fits for the different limb darkening laws (see text for details).
To illustrate the performance of the individual functional basis we show in Fig. 3 (top panel) the limb darkening laws of the solar simulation seen in the Kepler filter. The twocoefficient laws are obviously rather inadequate and show the largest deviations at the limb (μ ~ 0.0), in particular the linear, quadratic, logarithmic and exponential law (Eqs. (1), (2), (4) and (5)), while the square root law (Eq. (3)) exhibits a rather good match (it is already known that the square root law performs better for hotter stars, while for cooler stars the quadratic law is better see Claret 2000). The threeparameter functional basis (Eq. (6)) is performing well, however, it mismatches the limb rather significantly. The standard fourparameter nonlinear power law (Eq. (7)) is an excellent functional basis, and due to its versatility the fits result in extreme small residuals. In order to depict the precision of the individual laws quantitatively, we list the average χ^{2}, average maximal relative deviation and average location of the latter from all stellar models in Table 1. The fourparameter law performs for all stellar parameters significantly much better than any other limb darkening law, therefore, we will discuss subsequently the latter only.
Fig. 3 Solar intensity distribution vs. μangle for different limb darkening laws in the Kepler filter (top) and different bands from the Johnson filters (bottom panel). The 3D results are indicated with filled circles. In the top panel the blue and red lines are shifted by 0.1 and 0.2 respectively, while the black lines are unshifted. 

Open with DEXTER 
We consider a number of broad band filters b by convolving the response function S_{λ}, which considers the transmission of the filter b, with the integration of the intensity, We applied multiple standard broad band filters taken from the SYNPHOT package ^{1}, which comprises Bessel (JHK), Johnson (UBVRI) and Strömgren (uvby). Additionally we considered individual important instruments with CoRoT, Kepler ^{2}, Mauna Kea (JHKLM), SDSS (ugriz) and HST (ACS, STIS). Furthermore, the complete grid of spectra will be made available online, such that limb darkening from different filters or wavelength bands can be derived. In Fig. 3 (bottom panel) we show the five different Johnson filters for the solar model. The brightness distribution and its curvature are becoming more enhanced towards higher wavelength from the ultraviolet to the infrared, which is a general feature for all stellar parameters. The optical depth and the temperature gradient are dependent on the considered wavelength, since radiation at higher wavelength is emerging from higher geometrical depth.
4. Limb darkening
Fig. 4 Diskcenter normalized intensity, , against the inclination, μ, in the Kepler filter showing the limb darkening for different stellar parameters. Furthermore, we illustrate also the 1D ATLAS predictions for comparison (dashed lines), and the deviations Δ = 1D−3D (dashed dotted lines), which are enhanced by a factor of 2 and shifted by +0.2 for clarity. 

Open with DEXTER 
The emergent radiative intensity of stars decreases almost linear monotonically from the center to the edge, until it dropsoff sharply close to the limb (μ ~ 0.1), which is known as limb darkening and can be observed in the Sun. The radiation at diskcenter emerges from lower depths, while towards the limb one observes light from higher layers, where the temperature stratification has dropped very quickly, hence lower (darker) brightness (I ∝ T^{4}). We mention that the limb darkening has obvious boundary constraints being that the intensity is maximal with I(μ = 1) = 1 and minimal with I(μ = 0) ≃ 0. Furthermore, we note briefly that the variation of the limb darkening with stellar parameter depends also on the considered filter, and the results can differ significantly (Fig. 3).
In Fig. 4, we show an overview of different limb darkening variations for various stellar parameters. One finds in general that the variations are rather smooth and systematic, and the largest differences are given close to the limb (μ ~ 0.2), however, the variations with T_{eff} and log g are distinctive for different metallicities. For higher T_{eff} the decline in brightness exhibits a more pronounced convex curvature between μ = 0.2−0.7, while at cooler ones it is close to a linear drop. Therefore, hotter models end up with brighter intensities towards the limb around μ ~ 0.2 and a steeper drop beyond. For instance for dwarfs with solarmetallicity we find I_{0.2}/I_{1} = 0.4−0.6 for T_{eff} = 4500−7000 K. In order to illustrate the relative brightening at the limb for hotter dwarfs with solar metallicity, we show synthetic stellar disks in Fig. 5 with higher T_{eff}. Towards giants (lower log g) the brightness is higher than for dwarfs, however the changes are more subtle compared to the effective temperature. At lower metallicity, we find the differences with T_{eff} between metalpoor dwarfs (log g = 4.5,[Fe/H] = −3) being distinctively smaller, so that a pronounced curvature is given even for the coolest effective temperature (compare bottom panels in Fig. 4). In fact, we find basically no increase with T_{eff} around μ ~ 0.2 (I_{0.2}/I_{1} = 0.59−0.61). As we show further down, the Tinsensitivity of the limb darkening at [Fe/H] = −3 arises due to the temperature gradient. Claret (2000) had also found an enhanced curvature at lower metallicity with 1D models. The hottest and most metalpoor dwarfs are the brightest at the edge, and the sharp drop is the steepest. Furthermore, we note that the centertolimb variation curves for metalpoor models cross the corresponding ones for solarmetallicity models with otherwise the same stellar parameters (T_{eff}), with the exception of the models at high T_{eff}, however, this is not the case towards higher T_{eff}. On the other hand, the metalpoor giants are more similar to the solarmetallicity case, and exhibit also a clear T_{eff} sensitivity.
Fig. 5 Synthetic stellar disks for stars with log g = 4.5, [Fe/H] = 0, and increasing T_{eff} from 4500 K to 6500 K (from left to right) seen in the Kepler filter, and including a transiting exoplanet with p = 0.05 with progressing transit phase. Note the brighter limb towards higher T_{eff}. 

Open with DEXTER 
As next, we want to discuss the individual coefficients for the fourparameter limb darkening law (Eq. (7)). Therefore, we show in Fig. 6 the four coefficients for different T_{eff} and [Fe/H] and log g = 4.5,2.0. Towards higher T_{eff} both coefficients, a_{2} and a_{4}, are increasing until 6000 K, then above they decrease, while the coefficients a_{1} and a_{3} vary the opposite at solar metallicity. For lower [Fe/H] the Tdependence is inverted for a_{k} and the even coefficients, a_{2} and a_{4}, are decreasing, and the odd ones, a_{1} and a_{3}, are increasing. Another aspect worthy of attention is the correlation between the coefficients with halfinteger exponents in μ (a_{1} and a_{3}), and integer one (a_{2} and a_{4}) with T_{eff} and [Fe/H] (compare left with right panels in Fig. 6). In the Kepler filter we find the correlations to amount with C[a_{1},a_{3}] = 0.82 and C[a_{2},a_{4}] = 0.91 for all stellar parameters. Furthermore, the halfinteger exponents anticorrelate with the integer ones (compare top with bottom panels in Fig. 6). While the distant coefficients are less anticorrelated with C[a_{1},a_{4}] = −0.76, we find the anticorrelation between the successive coefficients being much tighter with C[a_{1},a_{2}] = −0.95 and C[a_{3},a_{4}] = −0.99. Claret (2000) noted also correlations between the coefficients. The coefficients a_{k} of the fourparameter law can be decomposed and considered individually. Then the integer exponents with the coefficients a_{2} and a_{4} are leading to a linear and quadratic polynomial respectively, which are describing the general slope of the limb darkening. On the other hand, the halfinteger exponents with the coefficients a_{1} and a_{3} are square root like functions, and are responsible for the curvature towards the limb.
The coefficient u for the linear limb darkening law (Eq. (1)) is rather crude, however, it has the major advantage of simplicity, since the linear law reduces the complex shape of the limb darkening into a single value. We display in Fig. 7u against T_{eff} for different stellar parameters (top panel). A larger value in u relates to a steeper drop in intensity and indicate lower brightness at the limb, and a lower value for u results in brighter limbs (see Eq. (1)). The coefficient u is mainly sensitive to the effective temperature and is decreasing with higher T_{eff}, and for different log g and [Fe/H] the differences are rather small. As a remark we note that the linear coefficients a and c from Eqs. (2) and (3) are similar, while the quadratic and square root coefficients, b and d, behave oppositely and increase for higher T_{eff} (not shown). The values for u ranges from 0.56 to 0.77 in the Kepler filter, which would yield I = 0.44 and 0.23 at the limb (μ = 0) respectively (the global range in u is from 0.03 to 1.05).
With a linear approximation of the Planck function one can derive the slope of the linear limb darkening law, u (e.g. Gray 2005; Hayek et al. 2012), which depends primarily on the temperature gradient at the optical surface and is given by (8)The approximation implies that a steeper temperature gradient at the optical surface will lead to stronger (steeper) limb darkening (larger u). In Fig. 7 (middle panel) we show also u′_{500 nm} considered at the optical surface (τ = 1) vs. T_{eff}, and we find u′_{500 nm} to correlate well with u (top panel). Also, similar as given in u, we find u′_{500 nm} being Tinsensitive for dwarf models with very low metallicity ([Fe/H] = −3), which arises from the temperature gradient term dT/dlog τ_{λ} in Eq. (8). The intensity is given by the source function and the lost radiation, which is stated by the radiative transfer eq., dI_{λ}/dτ_{λ} = S_{λ}−I_{λ}. Under the assumption of LTE, the source function can be approximated with the Planck function at the local temperature, . Therefore, the variation of the intensity with μ is sensitive to the temperature structure and in particular the temperature gradient. In Fig. 7 we show also the temperature gradient, ∇ = dlnT/dlnp_{tot}, considered at the optical surface (bottom panel). For lower metallicity the range in temperature gradient is enhanced, which is very similar to the intensity contrast. We had already mentioned the enhancement of the intensity contrast and temperature gradient at lower metallicity in Paper I. We find the reason for the enhancement being the lack of metals that are usually the most important electron donors for the formation of H^{−}, which is the dominating opacity source (see Nordlund & Dravins 1990). Therefore, in metalpoor models, the main contribution of electrons arises from the ionization of hydrogen. This is the reason for the strong enhancement of the of temperature gradient towards T_{eff} = 6000 K with [Fe/H] = −3, which is the reason for the Tinsensitivity of u and the limb darkening that we found for metalpoor dwarfs (see Fig. 4).
Fig. 6 Four coefficients for the fourparameter nonlinear limb darkening law (Eq. (7)) predicted in the Kepler filter from 3D RHD models for different stellar parameters. 

Open with DEXTER 
Fig. 7 Coefficient u of the linear limb darkening law (Eq. (1)) in the Kepler filter (top) and the approximation from Eq. (8) (middle) and temperature gradient at the optical surface (bottom panel) vs. T_{eff} for log g = 2.0, 4.5 and different [Fe/H]. 

Open with DEXTER 
5. Transit light curves
Fig. 8 Left figure: transit light curve vs. transit phase with p = 0.1 in the Kepler filter for different stellar parameters. The predictions from the 1D ATLAS models are also included (dashed lines). Right figure: relative deviations in the transit light curve with p = 0.1 between 3D atmosphere models and 1D ATLAS models given in %. The difference is δ = 1D/3D−1. 

Open with DEXTER 
During exosolar planet transits the planet eclipses its host star in the line of sight to earth, thereby diminishing the emergent intensity and leaving a characteristic imprint in the observed light curve. Theoretical light curve predictions are described by two main parameters, namely the ratio of the radii of planet and parent star, p = r_{planet}/r_{star}, and the normalized separation of the centers, z = d/r_{star}, with d being the centertocenter distance of the two occulting bodies. For computing the theoretical transit light curves, we used the publicly available code by Mandel & Agol (2002). In Fig. 8a, we show light curves with p = 0.1, i.e. for the case that the radius of the planet is 1/10 of the star, in the Kepler filter for the same LDC as discussed in Sect. 4. Then one can extract that the dwarfs with solarmetallicity exhibit for higher T_{eff} light curves with a more rectangular shaped feature, where the transit center is shallower and the width of the transit is broader than lower T_{eff}. The giant models show similar features with T_{eff} for both [Fe/H]. On the other hand, the metalpoor dwarfs are similarly boxshaped and do not alter much for lower effective temperatures. The limb of the intensity distribution will shape the ingress and egress of the light curve, e.g. when the limb darkening would be a step function (a_{2} ≃ 0. and a_{k} = 0.), then the light curve would be entirely rectangular. A straight linear dropping intensity distribution (a_{2} = 1 and a_{k} = 0.) results in a more elliptical shape with a narrow width, while a curved square root drop (a_{1} = 1 and a_{k} = 0.) would lead to an evenly circular shaped light curve. Therefore, we caution for the use of limb darkening laws that are incapable of rendering the dropoff at the limb, such as the bi and threeparametric laws, since these will introduce inevitably systematic errors in the theoretical transit light curves. Mandel & Agol (2002) found in a comparison between the quadratic and the fourparameter nonlinear power law (with p = 0.1) differences by 3%! Furthermore, the depth of the transit light curve depends primarily on the ratio of planet to host p. For larger p the light curves are increasingly deeper at the center, since more stellar light is effectively blocked during the transition due to a larger surface ratio of the planet. The limb darkening is also sensitive to the considered wavelengthregime (see Fig. 3). Therefore, the transit light curve will be different depending on the actual considered broad band filter, which samples its specific range in λ. We find in general that the light curve is towards lower wavelength (ultraviolet) systematically more convex shaped with a deeper center and more slender width, while towards higher wavelength (infrared) the light curve is more boxshaped with shallower centers and broader width. We note that a multiband photometry approach states a solution to this issue (see Knutson et al. 2007).
Fig. 9 Overview of the maximal relative difference in the limb darkening compared between 1D and 3D against T_{eff} with the Kepler filter. The relative deviations are retrieved by Δ = 1D−3D. 

Open with DEXTER 
6. Comparison with results from 1D models
The ATLAS models are the widest applied 1D atmosphere models for retrieving LDC, since its grid covers the broadest range in stellar parameters leading to a seamless coverage. The differences between the 1D and 3D models arise mainly from the differences in the temperature stratification, in particular, the temperature gradients.
For the comparison, we show in Fig. 4 also the limb darkening derived from 1D MLT models from the ATLAS grid (see Sing 2010). To ensure consistency the fourparameter nonlinear laws for 1D models are also shown in the Kepler filter. The 1D dwarf models with solar metallicity exhibit similar increasing curvatures for higher T_{eff} as given by the 3D models, only the 1D models are slightly brighter than the 3D, except for the coolest one (4500 K). In the case of giants the 1D limb darkening is distributed at much lower and also more linear intensities. The limb of the metalpoor 1D models lacks of a similar smooth sharp dropoff that is given in the 3D models, instead they depict a rather discontinuous behavior at the limb, which is doubtful to be correct. For metalpoor models it is known that the enforcement of radiative equilibrium is leading to an overestimation of the temperature stratification in the upper layers due to lack of spectral line absorption (see Asplund et al. 1999; Collet et al. 2007). The radiation at the limb is emerging from higher layers, therefore, the deviations at the limb are consistent (see also Paper I).
The largest deviations are usually found at the limb, since some of the 1D models depict an almost linear run at the edge, while the 3D models exhibit a comparably smooth dropoff at the limb. In Fig. 9, we show the deviations at the limb, ΔI_{0}/I_{1}, between the 1D and 3D limb darkening predictions. The differences increase significantly towards lower metallicity. Therefore, we advise the use of 3D limb darkening coefficients in particular for metalpoor stars.
We display the relative differences in the transit light curve between the 3D and 1D predictions in Fig. 8b. Similar to the findings of Hayek et al. (2012), we find a characteristic shape in the residuals with two extrema that are usually of opposite sign. One being at the ingress and egress phase, which arises from differences at the limb, and the other one taking place at the diskcenter. We found in Sect. 4 the 1D models being brighter than the 3D limb darkening, except for the giants at solarmetallicity. The transit light curves for the 1D models are similarly brighter at the diskcenter and dimmer at the transit phase edges for all models, only for the solarmetallicity giants it is the opposite. However, the maximal residuals are relatively small with a global maximum of ~0.04% in the transit light curve.
Since the differences in the transit light curves between 3D and 1D are rather small, the expected differences in transit parameters are also small. Hayek et al. (2012) determined a larger (smaller) radius of the exoplanet HD 209458b (HD 189733b) by 0.2% in the wavelength region between 2900 and 5700 . Hence, a reanalysis of the observed transit light curves with our LDC might change the transit parameters less than a percent.
Limb darkening coefficients derived from 3D RHD models in the Kepler filter for different stellar parameters (Cols. 1–3).
7. Conclusions
We derived on the basis of the Staggergrid, a large grid of 3D RHD atmosphere models, the limb darkening coefficients for the various biparametric and nonlinear limb darkening laws. The fourparameter nonlinear power law introduced by Claret (2000) is the only limb darkening law that is sufficiently versatile to express the intensity distribution with an excellent accuracy, while all other limb darkening laws are insufficient, in particular at the limb. Therefore, we recommend the use of the fourparameter functional basis only, in particular for the comparison with highprecision measurements in the hunt of extrasolar planets. We discussed the limb darkening in the Kepler filter for various stellar parameters, and outlined systematical variations that exposed the complex changes of the brightness distribution, in particular with the effective temperature. We compared also our new LDC with predictions from widely used 1D ATLAS models, and the largest differences are given towards the limb. The 1D models are often brighter than 3D predictions, only for giant models with solarmetallicity we find opposite differences. Furthermore, we displayed the systematic (anti)correlations between the coefficients a_{k} between halfinteger and integer exponents of the fourparameter law. We found the coefficient of linear limb darkening law, u, to scale with the temperature gradient and the Planck function. Theoretical transit light curves indicate similar systematical differences between 1D and 3D as the limb darkening variations implied, which are relatively small. However, as observations indicate (Knutson et al. 2007; Hayek et al. 2012), these can be measured with highprecision observations. Therefore, we advise to use of the new LDC.
Acknowledgments
We acknowledge access to computing facilities at the Rechenzentrum Garching (RZG) of the Max Planck Society and at the Australian National Computational Infrastructure (NCI) where the simulations were carried out. Also, we acknowledge the Action de Recherche Concertée (ARC) grant provided by the Direction générale de l’Enseignement non obligatoire et de la Recherche scientifique – Direction de la Recherche scientifique – Communauté française de Belgique, and the F.R.S. FNRS. Remo Collet is the recipient of an Australian Research Council Discovery Early Career Researcher Award (project number DE120102940).
References
 Allende Prieto, C., Asplund, M., García López, R. J., & Lambert, D. L. 2002, ApJ, 567, 544 [NASA ADS] [CrossRef] [Google Scholar]
 Asplund, M., Nordlund, Å., Trampedach, R., & Stein, R. F. 1999, A&A, 346, L17 [NASA ADS] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Aufdenberg, J. P., Ludwig, H.G., & Kervella, P. 2005, ApJ, 633, 424 [NASA ADS] [CrossRef] [Google Scholar]
 Bigot, L., Kervella, P., Thévenin, F., & Ségransan, D. 2006, A&A, 446, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 BöhmVitense, E. 1958, Z. Astrophys., 46, 108 [Google Scholar]
 Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, L45 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Chiavassa, A., Plez, B., Josselin, E., & Freytag, B. 2009, A&A, 506, 1351 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Chiavassa, A., Collet, R., Casagrande, L., & Asplund, M. 2010, A&A, 524, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chiavassa, A., Bigot, L., Kervella, P., et al. 2012, A&A, 540, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A. 2000, A&A, 363, 1081 [NASA ADS] [Google Scholar]
 Claret, A. 2004, A&A, 428, 1001 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A., DiazCordoves, J., & Gimenez, A. 1995, A&AS, 114, 247 [NASA ADS] [Google Scholar]
 Collet, R., Asplund, M., & Trampedach, R. 2007, A&A, 469, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Davis, J., Tango, W. J., & Booth, A. J. 2000, MNRAS, 318, 387 [NASA ADS] [CrossRef] [Google Scholar]
 DiazCordoves, J., Claret, A., & Gimenez, A. 1995, A&AS, 110, 329 [NASA ADS] [Google Scholar]
 Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres (Cambridge: Cambridge University Press) [Google Scholar]
 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hayek, W., Sing, D., Pont, F., & Asplund, M. 2012, A&A, 539, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, T. M., & Gilliland, R. L. 2007, ApJ, 655, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Kurucz, R. L. 1979, ApJS, 40, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Magic, Z. & Asplund, M. 2014, A&A, submitted [arXiv:1405.7628] [Google Scholar]
 Magic, Z., Collet, R., Asplund, M., et al. 2013a, A&A, 557, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Magic, Z., Collet, R., Hayek, W., & Asplund, M. 2013b, A&A, 560, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Magic, Z., Collet, R., & Asplund, M. 2014, A&A, submitted [arXiv:1403.6245] [Google Scholar]
 Magic, Z., Weiss, A., & Asplund, M. 2015, A&A, 573, A89 (Paper III) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [NASA ADS] [CrossRef] [Google Scholar]
 Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Mihalas, D. 1970, Stellar atmospheres, Series of Books in Astronomy and Astrophysics (San Francisco: Freeman) [Google Scholar]
 Milne, E. A. 1921, MNRAS, 81, 361 [NASA ADS] [CrossRef] [Google Scholar]
 Nordlund, A. 1982, A&A, 107, 1 [NASA ADS] [Google Scholar]
 Nordlund, A., & Dravins, D. 1990, A&A, 228, 155 [Google Scholar]
 Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2 [Google Scholar]
 Pereira, T. M. D., Asplund, M., Collet, R., et al. 2013, A&A, 554, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seager, S., & Sasselov, D. D. 2000, ApJ, 537, 916 [NASA ADS] [CrossRef] [Google Scholar]
 Sing, D. K. 2010, A&A, 510, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Skartlien, R. 2000, ApJ, 536, 465 [NASA ADS] [CrossRef] [Google Scholar]
 Southworth, J. 2008, MNRAS, 386, 1644 [NASA ADS] [CrossRef] [Google Scholar]
 Stein, R. F., & Nordlund, A. 1998, ApJ, 499, 914 [NASA ADS] [CrossRef] [Google Scholar]
 van Hamme, W. 1993, AJ, 106, 2096 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2011, PASP, 123, 412 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Table with the limb darkening coefficients
In Table A.1 we listed a subset of the limb darkening coefficients for F, G, and K mainsequence stars. The full table is available at the CDS.
All Tables
Precision of the functional fits for the different limb darkening laws (see text for details).
Limb darkening coefficients derived from 3D RHD models in the Kepler filter for different stellar parameters (Cols. 1–3).
All Figures
Fig. 1 Emergent monochromatic intensity at 500 nm shown for eight μ angles: 0.2, 0.3, 0.5, 0.7, 0.8, 0.9, 1.0 (from top to down, respectively) and for a selection of stars: mainsequence, turnoff, Kgiant, Kdwarf (from left to right, respectively) with solar metallicity. Furthermore, we indicated the normalized mean intensity and the intensity contrast (both in percent). 

Open with DEXTER  
In the text 
Fig. 2 Spectral energy distribution vs. wavelength for models with log g = 4.5, [Fe/H] = 0.0 and various T_{eff} (red/orange colors). The continuum fluxes are also shown (solid black lines). 

Open with DEXTER  
In the text 
Fig. 3 Solar intensity distribution vs. μangle for different limb darkening laws in the Kepler filter (top) and different bands from the Johnson filters (bottom panel). The 3D results are indicated with filled circles. In the top panel the blue and red lines are shifted by 0.1 and 0.2 respectively, while the black lines are unshifted. 

Open with DEXTER  
In the text 
Fig. 4 Diskcenter normalized intensity, , against the inclination, μ, in the Kepler filter showing the limb darkening for different stellar parameters. Furthermore, we illustrate also the 1D ATLAS predictions for comparison (dashed lines), and the deviations Δ = 1D−3D (dashed dotted lines), which are enhanced by a factor of 2 and shifted by +0.2 for clarity. 

Open with DEXTER  
In the text 
Fig. 5 Synthetic stellar disks for stars with log g = 4.5, [Fe/H] = 0, and increasing T_{eff} from 4500 K to 6500 K (from left to right) seen in the Kepler filter, and including a transiting exoplanet with p = 0.05 with progressing transit phase. Note the brighter limb towards higher T_{eff}. 

Open with DEXTER  
In the text 
Fig. 6 Four coefficients for the fourparameter nonlinear limb darkening law (Eq. (7)) predicted in the Kepler filter from 3D RHD models for different stellar parameters. 

Open with DEXTER  
In the text 
Fig. 7 Coefficient u of the linear limb darkening law (Eq. (1)) in the Kepler filter (top) and the approximation from Eq. (8) (middle) and temperature gradient at the optical surface (bottom panel) vs. T_{eff} for log g = 2.0, 4.5 and different [Fe/H]. 

Open with DEXTER  
In the text 
Fig. 8 Left figure: transit light curve vs. transit phase with p = 0.1 in the Kepler filter for different stellar parameters. The predictions from the 1D ATLAS models are also included (dashed lines). Right figure: relative deviations in the transit light curve with p = 0.1 between 3D atmosphere models and 1D ATLAS models given in %. The difference is δ = 1D/3D−1. 

Open with DEXTER  
In the text 
Fig. 9 Overview of the maximal relative difference in the limb darkening compared between 1D and 3D against T_{eff} with the Kepler filter. The relative deviations are retrieved by Δ = 1D−3D. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.