Analytic solutions to the maximum and average exoplanet transit depth for common stellar limb darkening laws

TL;DR

This paper derives analytical formulas to accurately determine the true exoplanet transit depth and radius, accounting for stellar limb darkening effects across various limb darkening laws, improving upon the simple squared radius ratio approximation.

Contribution

It provides new analytical solutions for transit depth overshoot due to limb darkening for arbitrary impact parameters and limb darkening laws, enabling more precise radius measurements without light curve fitting.

Findings

Transit depth overshoot ranges from 15% to 30% depending on star type.

Analytical solutions significantly reduce errors compared to small planet approximation.

Formulas allow direct calculation of planet radius from transit data with high accuracy.

Abstract

The depth of an exoplanetary transit in the light curve of a distant star is commonly approximated as the squared planet-to-star radius ratio, (R_p/R_s)^2. Stellar limb darkening, however, results in significantly deeper transits. Here we derive analytical solutions to the overshoot of the mid-transit depth caused by stellar limb darkening compared to the (R_p/R_s)^2 estimate for arbitrary transit impact parameters. In turn, this allows us to compute the true planet-to-star radius ratio from the transit depth for a given parameterization of a limb darkening law and for a known transit impact parameter. We calculate the maximum emerging specific stellar intensity covered by the planet in transit and derive analytic solutions for the transit depth overshoot. Solutions are presented for the linear, quadratic, square-root, logarithmic, and non-linear stellar limb darkening with arbitrary transit impact parameters. We also derive formulae to calculate the average intensity along the transit chord, which allows us to estimate the actual transit depth (and therefore R_p/R_s) from the mean in-transit flux. The transit depth overshoot of exoplanets compared to the (R_p/R_s)^2 estimate increases from about 15% for A main-sequence stars to roughly 20% for sun-like stars and some 30% for K and M stars. The error in our analytical solutions for R_p/R_s from the small planet approximation is orders of magnitude smaller than the uncertainties arising from typical noise in real light curves and from the uncertain limb darkening. Our equations can be used to predict with high accuracy the expected transit depth of extrasolar planets. The actual planet radius can be calculated from the measured transit depth or from the mean in-transit flux if the stellar limb darkening can be properly parameterized and if the transit impact parameter is known. Light curve fitting is not required.

Tables

Table 1. Table 1: Parameterization of the quadratic limb darkening in the Kepler band (as per Claret & Bloemen, 2011 ) for the main-sequence stars shown in Fig.  3 . All stars have log ⁡ ( g ) = 4.5 𝑔 4.5 \log(g)=4.5 . The corresponding radius overshoot following from Eq. ( 3 ) is shown in the last column.

Spectral Type	$T_{\mathrm{eff}}$	$a$	$b$	$o_{\mathrm{LD}}$
A5	8500	0.2474	0.3026	0.153
F5	6500	0.2987	0.3098	0.178
G2	5750	0.4089	0.2556	0.218
K5	4500	0.6809	0.0645	0.312
M0	4000	0.5444	0.1898	0.271

Equations

$$o_{\mathrm{LC}}=\frac{\delta -(R_{\mathrm{p}}/R_{\mathrm{s}}{)}^2}{(R_{\mathrm{p}}/R_{\mathrm{s}}{)}^2}=\frac{\delta }{(R_{\mathrm{p}}/R_{\mathrm{s}}{)}^2}-1.$$

$${\Bigg{(}}\frac{R_{\rm p}}{R_{\rm s}}{\Bigg{)}}=\sqrt{\frac{\delta}{o_{\rm LC}+1}}\ .$$

$$o_{\mathrm{LD}}=\frac{I_{\mathrm{p}}-⟨I{⟩}_A}{⟨I{⟩}_A}=\frac{I_{\mathrm{p}}}{⟨I{⟩}_A}-1$$

$${\Bigg{(}}\frac{R_{\rm p}}{R_{\rm s}}{\Bigg{)}}=\sqrt{\delta\ \frac{\langle I\rangle_{A}}{I_{\rm p}}}\ .$$

$$\delta =\frac{1-⟨f_{\mathrm{in}}⟩}{⟨I{⟩}_x}.$$

$$\left(\frac{R_{\rm p}}{R_{\rm s}}\right)=\sqrt{{\Big{(}}1-{\langle}f_{\rm in}{\rangle}{\Big{)}}\frac{{\langle}I{\rangle}_{A}}{{\langle}I{\rangle}_{x}\ I_{\rm p}}}\ ,$$

$$I_{\rm lin}(\gamma)=1-u\ {\Big{(}}1-\cos(\gamma){\Big{)}}\ ,$$

$$I_{\mathrm{lin}}(r)=1-u(1-\sqrt{1-r^2}).$$

$$I_{\mathrm{tot},\mathrm{lin}}=$$

$$=$$

$$⟨I_{\mathrm{lin}}{⟩}_A=\frac{I_{\mathrm{tot},\mathrm{lin}}(u)}{\pi }=1-\frac{u}{3},$$

$$⟨I_{\mathrm{lin}}{⟩}_x$$

$$\displaystyle=\frac{1}{\sqrt{1-p^{2}}}\int\displaylimits_{0}^{\sqrt{1-p^{2}}}{\rm d}x\,{\Bigg{(}}1-u\ {\Big{(}}1-\sqrt{1-(p^{2}+x^{2})}\,{\Big{)}}{\Bigg{)}}$$

$$=1-u\left(1-\frac{\pi }{4}\sqrt{1-p^2}\right),$$

$$I_{\mathrm{quad}}(\gamma )$$

$$I_{\mathrm{quad}}(r)$$

$$⟨I_{\mathrm{quad}}{⟩}_A=1-\frac{a}{3}-\frac{b}{6}.$$

$$⟨I_{\mathrm{quad}}{⟩}_x$$

$$\displaystyle\hskip 8.5359pt-b\left(\frac{5}{3}-\frac{\pi}{2\sqrt{1-p^{2}}}-p^{2}{\Big{(}}\frac{2}{3}-\frac{\pi}{2\sqrt{1-p^{2}}}{\Big{)}}\right)\ .$$

$$I_{\mathrm{sqr}}(\gamma )$$

$$I_{\mathrm{sqr}}(r)$$

$$⟨I_{\mathrm{sqr}}{⟩}_A=1-\frac{c}{3}-\frac{d}{5}.$$

$$⟨I_{\mathrm{sqr}}{⟩}_x$$

$$-d\left(1-\frac{\sqrt{\pi }(1-p^2{)}^{1/4}\Gamma \left(\frac{5}{4}\right)}{2\Gamma \left(\frac{7}{4}\right)}\right),$$

$$I_{\log }(\gamma )$$

$$I_{\log }(r)$$

$$⟨I_{\log }{⟩}_A=1+\frac{2}{9}f-\frac{e}{3}.$$

$$⟨I_{\log }{⟩}_x=$$

$$\displaystyle\hskip 8.5359pt-f\frac{\pi}{8}\sqrt{1-p^{2}}\left(1-\ln{\Big{(}}\frac{4}{1-p^{2}}{\Big{)}}\right)\ .$$

$$I_{\mathrm{nl}}(\gamma )$$

$$\displaystyle\hskip 17.64056pt-a_{3}{\Big{(}}1-(\cos(\gamma))^{3/2}{\Big{)}}-a_{4}{\Big{(}}1-(\cos(\gamma))^{2}{\Big{)}}$$

$$I_{\mathrm{nl}}(r)$$

$$\displaystyle\hskip 17.64056pt-a_{3}{\Big{(}}1-(1-r^{2})^{3/4}{\Big{)}}-a_{4}r^{2}\ ,$$

$$⟨I_{\mathrm{nl}}{⟩}_A=1-\frac{a_1}{5}-\frac{a_2}{3}-\frac{3a_3}{7}-\frac{a_4}{2}.$$

$$⟨I_{\mathrm{nl}}{⟩}_x$$

$$-a_2\left(1-\frac{\pi }{4}\sqrt{1-p^2}\right)$$

Full text

*Aims. *We calculate the overshoot of the mid-transit depth caused by stellar limb darkening compared to the $(R_{\rm p}/R_{\rm s})^{2}$ estimate for arbitrary transit impact parameters. In turn, this allows us to compute the true planet-to-star radius ratio from the transit depth for a given parameterization of a limb darkening law and for a known transit impact parameter.

*Methods. *We calculated the maximum emerging specific stellar intensity covered by the planet in transit and derive analytic solutions for the transit depth overshoot. Solutions are presented for the linear, quadratic, square-root, logarithmic, and nonlinear stellar limb darkening with arbitrary transit impact parameters. We also derived formulae to calculate the average intensity along the transit chord, which allows us to estimate the actual transit depth (and therefore $R_{\rm p}/R_{\rm s}$) from the mean in-transit flux.

*Results. *The transit depth overshoot of exoplanets compared to the $(R_{\rm p}/R_{\rm s})^{2}$ estimate increases from about $15$ % for main-sequence stars of spectral type A to roughly $20$ % for sun-like stars and some $30$ % for K and M stars. The error in our analytical solutions for $R_{\rm p}/R_{\rm s}$ from the small planet approximation is orders of magnitude smaller than the uncertainties arising from typical noise in real light curves and from the uncertain limb darkening.

*Conclusions. *Our equations can be used to predict with high accuracy the expected transit depth of extrasolar planets. The actual planet radius can be calculated from the measured transit depth or from the mean in-transit flux if the stellar limb darkening can be properly parameterized and if the transit impact parameter is known. Light curve fitting is not required.

1 Introduction

The planetary radius ($R_{\rm p}$) is one of the key properties that are currently being derived by several exoplanet hunting surveys. The most successful method to determine the radius of an exoplanet is the transit method, which measures the slight decrease of the brightness of a star as the planet traverses the stellar disk as seen from Earth (Struve, 1952).

In a simple picture, the star’s appearance can be modeled as a circle with uniform brightness. Then the ratio of the planetary radius and stellar radius ($R_{\rm s}$) can be estimated from the constant111Even with the neglect of stellar limb darkening, however, the transit would still have an ingress and egress, both of which lead to a gradual transition between the in- and out-of-transit apparent stellar brightness (a mathematical description is given in Appendix A of Heller, 2014). transit depth ($\delta$) via $(R_{\rm p}/R_{\rm s})=\sqrt{\delta}$. In reality, however, stars show center-to-limb brightness variations that affect the estimated planet-to-star radius ratio (Csizmadia et al., 2013).

The common variable used throughout the literature to describe limb darkening is $\mu=\cos(\gamma)$, where $\gamma$ is the angle between the emerging specific intensity of the star and the line of sight of the observer. Early attempts to model stellar limb darkening used equations that are linear in $\mu$ (Milne, 1921) but these became insufficient as computer models could be used to model stellar atmospheres. The quadratic (Manduca et al., 1977), square root (Diaz-Cordoves & Gimenez, 1992), and logarithmic limb darkening laws (Klinglesmith & Sobieski, 1970) provided better agreement with the observations. The advent of exoplanet transit observations (Charbonneau et al., 2000) and new high-accuracy space-based transit photometry (Brown et al., 2001), however, required even better precision.

Claret (2000) presented a four parameter nonlinear limb darkening law that has now widely been adapted in the exoplanet community in combination with proper limb darkening parameters computed with stellar atmosphere models (Claret & Bloemen, 2011). Nevertheless, even though the nonlinear law might be somewhat more precise in certain regimes of the parameter space and for certain stars (Espinoza & Jordán, 2016), the quadratic law is commonly used because Mandel & Agol (2002) provided an analytic solution to the resulting transit light curve for the quadratic law.

Although it is well known that the transit depth is not equal to $(R_{\rm p}/R_{\rm s})^{2}$, no framework exists that conclusively derives the actual deviations from that estimate for stars with limb darkening. Here we derive the correction factors to the $(R_{\rm p}/R_{\rm s})\approx\sqrt{\delta}$ approximation for each of the above-mentioned stellar limb darkening laws.

2 Methods

We present two methods to derive the actual planet-to-star radius ratio for arbitrary parameterizations of the common limb darkening laws. The first method is based on a measurement of the transit depth (Sect. 2.1) and the second method requires a measurement of the arithmetic mean of the in-transit flux (Sect. 2.2).

2.1 Transit depth overshoot

In Fig. 1 we show an example transit light curve using the python implementation222Available at http://www.astro.ucla.edu/~ianc/files as transit.py. of the Mandel & Agol (2002) analytic model for quadratic stellar limb darkening by Ian Crossfield. We arbitrarily chose $a=0.4$ and $b=0.4$, irrespective of the stellar spectral type that this choice could imply, a transit impact parameter of zero, and a circular orbit. The solid line shows the numerical computation and the horizontal dashed line shows the $(R_{\rm p}/R_{\rm s})^{2}$ estimate. The vertical arrow at the bottom of the transit denotes what we refer to as the overshoot ($o$) between the actual transit depth and the $(R_{\rm p}/R_{\rm s})^{2}$ approximation. There are two ways to calculate $o$ and we can show that they are equivalent to very high precision.

2.1.1 Transit depth overshoot in the light curve

The first way to calculate the overshoot uses the light curve. We consider a light curve with an out-of-transit flux normalized to 1, then the transit depth is simply the difference between 1 and the minimum in-transit flux, $\delta=1-f_{\rm min}$. We then define the transit depth overshoot as determined from the light curve as the difference between the transit depth and the $(R_{\rm p}/R_{\rm s})^{2}$ estimate in units of $(R_{\rm p}/R_{\rm s})^{2}$,

[TABLE]

This is the definition of the overshoot shown in Fig. 1 and it illustrates our problem. From an observational perspective, $\delta$ can be readily measured, but $(R_{\rm p}/R_{\rm s})^{2}$ is what we actually want to determine. Hence, if we knew the transit depth overshoot for a given host star of a transiting planet, then we could correct the simple $(R_{\rm p}/R_{\rm s})\approx\sqrt{\delta}$ approximation and derive the true planet-to-star radius ratio as

[TABLE]

Fortunately, there is a second way to calculate the overshoot.

2.1.2 Overshoot of the stellar limb darkening profile

This second approach relates to the star’s emerging specific intensity, $I$, which we simply refer to as intensity from now on. The only simplification that we require is to assume that the intensity profile covered by the planet in transit is constant, which is equivalent to the assumption that the planet is infinitesmal. As shown in Sect. 3.3, the resulting error from this small planet approximation is extremely small.

We consider a star that is transited by a planet (Fig. 2). The minimum in-transit distance between the planet and the center of the stellar disk is referred to as the transit impact parameter $0\,{\leq}\,p\,{\leq}\,1$. The radial coordinate, measured from the disk center in units of $R_{\rm s}$, is $r$ and we orient the abscissa, or $x$ coordinate, of our reference system to be parallel to the transit chord of the planet and to have its origin in the center of the stellar disk, so that $r^{2}=p^{2}+x^{2}$. Introducing $x_{1}$ and $x_{2}$ as the distances traversed by the planet during the first and second halves of its transit, respectively, we find $1=p^{2}+x_{1}^{2}=p^{2}+x_{2}^{2}$.

The key idea of this study is that the overshoot occurs at the minimum radial distance of the planet to the stellar disk center, that is, when $x=0$ and $r=p$. Without loss of generality, we set the intensity at the disk center $I_{0}=1$ to simplify our notation. We only consider radially symmetric intensity profiles $I(r)$ in this paper and so we refer to the intensity covered by the planet at mid-transit as $I(r=p)=I_{\rm p}$.

Irrespective of the actual limb darkening law, we can now define the expected transit depth overshoot as the difference between the intensity at mid-transit and the average intensity across the entire stellar disk area, $\langle I\rangle_{A}$, as

[TABLE]

Replacing in Eq. (2) the overshoot in the light curve $o_{\rm LC}$, which is unknown a priori, with the overshoot of the intensity profile from Eq. (3) gives

[TABLE]

2.2 Mean in-transit flux

Alternatively, we can compute the planet-to-star radius ratio from the average in-transit flux ${\langle}f_{\rm in}{\rangle}$ and from the average intensity sampled by the transit chord ${\langle}I{\rangle}_{x}$. From the mere definition of these two quantities it is clear that

[TABLE]

Plugging this expression for $\delta$ into Eq. (4) we obtain

[TABLE]

which allows us to estimate the actual planet-to-star radius ratio based on the measured mean in-transit flux of a limb-darkened transit light curve and based on the disk average, radial average, and mid-transit intensity of the limb darkening profile.

2.3 Limb darkening laws

Now we set out to derive analytical expressions for the average intensity across the stellar disk $\langle I\rangle_{A}$ and the average intensity along the planetary in-transit path parallel to the $x$ coordinate $\langle I\rangle_{x}$.

2.3.1 Linear limb darkening

We start by rewriting the linear limb darkening law as in Claret (2000), originally derived by (Milne, 1921, Eq. 45 therein), but now as a function of the angle ($\gamma$) between the line of sight and the emerging intensity,

[TABLE]

where $u$ is the limb darkening parameter. With $r=\sin(\gamma)$ as the distance from the star disk center normalized by the stellar radius (see Fig. 2) and using $\cos\left(\,\arcsin(r)\,\right)\leavevmode\nobreak\ =\leavevmode\nobreak\ \sqrt{1-r^{2}}$ we have

[TABLE]

We then calculate the total intensity across the area elements $({\rm d}\varphi\,{\rm d}r\ r)$ of the stellar disk, where $0\,{\leq}\,\varphi\,{\leq}\,2\pi$ is the azimuth angle, as

[TABLE]

The disk-averaged intensity then is

[TABLE]

where the factor $\pi$ in Eq. (11) is the area of the apparent stellar disk with a normalized radius of $r=1$. Using Eq. (11) as $\langle I\rangle_{A}$ and $I_{\rm lin}(r=p)$ from Eq. (8) as $I_{\rm p}$, we can estimate the planet-to-star radius ratio from the observed transit depth via Eq. (4) and based on the linear limb darkening law. Alternatively, we can calculate the transit depth overshoot from Eq. (3).

Moving on to the mean intensity across the transit chord, we use the relationship $x_{1}=\sqrt{1-p^{2}}=-x_{2}$ to parameterize the lower and upper boundaries of integration. Instead of integrating the intensity from ingress (at $x_{2}=-\sqrt{1-p^{2}}$) to egress (at $x_{1}=+\sqrt{1-p^{2}}$), we take advantage of the fact that the average intensity transited by the planet during the first half is the same as during the second half of the transit, which again is equal to the average intensity transited during the entire transit. Hence,

[TABLE]

which can be used as $\langle I\rangle_{x}$, for example, in Eq. (6).

2.3.2 Quadratic limb darkening

Rewriting the quadratic limb darkening law as a function of the angle and the normalized radius to the disk center gives

[TABLE]

where $a$ and $b$ are the two limb darkening parameters. The total disk-integrated intensity $I_{\rm tot,quad}$ can then be calculated as per the right-hand side of Eq. (9) and the disk-averaged intensity then turns out as $I_{\rm tot,quad}/\pi$, giving

[TABLE]

An equivalent expression to Eq. (16) has been presented by Csizmadia et al. (2013, Appendix B therein) to calculate the transit depth as a function of $\mu$. Hippke & Heller (2019) have used this formula for their Transit Least Squares transit detection algorithm (TLS).

The transit chord average of the intensity in the quadratic limb darkening law can be obtained in analogy to the previous section using the right-hand side of Eq. (12). This results in

[TABLE]

2.3.3 Square root limb darkening

The square root limb darkening law as a function $\gamma$ or $r$ appears as

[TABLE]

respectively, where $c$ and $d$ are the corresponding limb darkening parameters. We derive the total disk-integrated intensity $I_{\rm tot,sqr}$ as per Eq. (9) and then divide $I_{\rm tot,sqr}$ by the normalized disk area $\pi$ to obtain the disk-averaged intensity as

[TABLE]

The transit chord average of the intensity in the square root limb darkening then follows from the right-hand side of Eq. (12),

[TABLE]

with $\Gamma\left(\frac{5}{4}\right)=0.90640...$ and $\Gamma\left(\frac{7}{4}\right)=0.91906...$ (see Appendix A).

2.3.4 Logarithmic limb darkening

The logarithmic limb darkening law can be written as a function of $\gamma$ or $r$ in the following ways,

[TABLE]

where $e$ and $f$ are the limb darkening parameters. The total disk-integrated intensity $I_{\rm tot,log}$ then follows via Eq. (9) and the disk-averaged intensity is equal to $I_{\rm tot,log}/\pi$, resulting in

[TABLE]

We note that in contrast to all other limb darkening laws treated in this paper, Eq. (26) can be larger than 1, which means that the star can be brighter on average than in the disk center. In the transit light curve this would manifest itself by a mid-transit depth that is not the deepest point in the transit trough.

The average intensity along the planetary transit chord in the case of logarithmic limb darkening law is computed from the right-hand side of Eq. (12) as

[TABLE]

2.3.5 Nonlinear limb darkening

The nonlinear limb darkening law as a function $\gamma$ and $r$ can be written as

[TABLE]

with $a_{1}$ to $a_{4}$ as the limb darkening parameters. The total disk-integrated intensity for the nonlinear limb darkening law $I_{\rm tot,nl}$ is then obtained from the right-hand side of Eq. (9), and dividing $I_{\rm tot,nl}$ by the normalized apparent disk area $\pi$ gives

[TABLE]

Finally, radial averaging of the specific intensity in the case of the nonlinear limb darkening law is obtained according to Eq. (12) as

[TABLE]

where $\Gamma\left(\frac{9}{4}\right)=1.13300...$. This derivation requires solving an integral that does not have a closed form solution, similar to the case of the square root limb darkening law. We nevertheless find the elementary antiderivative shown above by applying a substitution, see the analogous case for the square root limb darkening law in Appendix A.

2.3.6 Comparing transit depths between limb darkening laws

For the linear limb darkening law, we can derive the single limb darkening parameter $u_{o^{\prime}}$, which produces the same transit depth overshoot as any of the other limb darkening laws. The linear limb darkening law and the quadratic, square, and logarithmic limb darkening laws produce equal overshoots ($o_{\rm LD}=o_{\rm LD}^{\prime}$) if

[TABLE]

Of course, the relations between the other limb darkening parameters can be derived analogously by equating the expressions for $\langle I_{\rm quad}\rangle_{A}$ and $\langle I_{\rm sqr}\rangle_{A}$, $\langle I_{\rm quad}\rangle_{A}$ and $\langle I_{\rm log}\rangle_{A}$, and $\langle I_{\rm quad}\rangle_{A}$ and $\langle I_{\rm log}\rangle_{A}$, respectively, and solving for whichever limb darkening parameters are of interest.

3 Results

3.1 Transit depth overshoot

We apply our framework to the numerically computed light curve in Fig. 1 as an example. First, we measure $\delta\leavevmode\nobreak\ =\leavevmode\nobreak\ 100\,\%\leavevmode\nobreak\ -\leavevmode\nobreak\ 99.9875\,\%=0.0125\,\%\leavevmode\nobreak\ =\leavevmode\nobreak\ 125$ parts per million. As this light curve has been computed using the quadratic limb darkening law and using $a=0.4$ and $b=0.4$, we use the same parameterization of Eq. (16) and calculate $\langle I_{\rm quad}\rangle_{A}\leavevmode\nobreak\ =\leavevmode\nobreak\ 0.8$. Moreover, with $p=0$ in Fig. 1, we have $I_{\rm p}=1$ and therefore Eq. (4) gives the correct value of $(R_{\rm p}/R_{\rm s})=0.01$ used for this simulation. The corresponding overshoot is calculated from Eq. (3) as $o_{\rm LD}=0.25$.

For an analogous application of this procedure to an observed light curve the limb darkening parameters and the transit impact parameter need to be known. The limb darkening parameters can be read from pre-computed tables if the stellar spectral type is known and if the transit impact parameter can be estimated. The latter can be done analytically based on measurements of the transit depth, orbital period, the time between first and fourth contact, and the time between second and third contact (Seager & Mallén-Ornelas, 2003).

Figure 3 illustrates the transit depth overshoot of the quadratic limb darkening law for any combination of $a$ and $b$ and using $p=0$. We use the tables of Claret & Bloemen (2011) to obtain the limb darkening coefficients in the Kepler band for several main-sequence stars with metallicity [Fe/H] = 0, surface gravity $\log(g)=4.5$, and zero microturbulent velocity. The spectral types of these stars are indicated as A5, F5, sun, K5, and M0 in Fig. 3, respectively. As a result, we find that the quadratic limb darkening law predicts an overshoot of $o_{\rm LD}=0.218$ times $(R_{\rm p}/R_{\rm s})^{2}$ for a small planet around a sun-like star. This exact solution can be derived by plugging $a$ and $b$ in Eq. (16) and the resulting $\langle I_{\rm quad}\rangle_{A}$ into Eq. (3), where $I_{\rm p}=1$ for a transit impact parameter of zero. Table 1 lists the overshoots for the other template stars as well, together with the effective temperatures and limb darkening coefficients.

3.2 Mean in-transit flux

For the numerically computed example light curve in Fig. 1, where the simulated planet-to-star ratio is 0.01 and the limb darkening parameters are $a=0.4$ and $b=0.4$ for $\langle I_{\rm quad}\rangle_{A}$ and $\langle I_{\rm quad}\rangle_{x}$ in Eqs. (16) and (17), respectively, we measure ${\langle}f_{\rm in}{\rangle}\leavevmode\nobreak\ =\leavevmode\nobreak\ 0.999891$. Equation (6) then predicts $(R_{\rm p}/R_{\rm s})\leavevmode\nobreak\ =\leavevmode\nobreak\ 0.009996$, the source of the small error in the planetary radius of $4\times 10^{-6}\,R_{\rm s}\approx 2.8$ km being in the ingress and egress, as explained in the next section.

3.3 Small planet approximation

The error in the planetary radius ($\varepsilon_{\rm p}$) from the small planet approximation can be calculated by comparing the planet-to-star radius ratio $(R_{\rm p}/R_{\rm s})_{\rm obs}$ derived from the measured transit depth in a simulated light curve as per Eq. (4) with the actual planet-to-star radius ratio that went into the simulation, $(R_{\rm p}/R_{\rm s})_{\rm sim}$. The simulation of a Jupiter transit with $p=0$ (hence $I_{\rm p}=1$) across a sun-like star, where $(R_{\rm p}/R_{\rm s})_{\rm sim}=0.099386$ and the quadratic limb darkening parameters are {$a=0.4089$, $b=0.2556$}, results in $\delta=1-f_{\rm min}=1.20176\times 10^{-2}$. With $\langle I\rangle_{A}=0.8211$, Eq. (4) then predicts $(R_{\rm p}/R_{\rm s})_{\rm obs}=0.099336$ and the resulting error in the planetary radius of a Jupiter-like planet around a sun-like star is $\varepsilon_{\rm p}=5\leavevmode\nobreak\ {\times}\leavevmode\nobreak\ 10^{-5}\,R_{\odot}\approx 35$ km. For or an Earth-like planet, this same exercise yields an error of $\varepsilon_{\rm p}\leavevmode\nobreak\ =\leavevmode\nobreak\ 4\leavevmode\nobreak\ {\times}\leavevmode\nobreak\ 10^{-8}\,R_{\odot}\approx 28$ m.

The main source of error in $(R_{\rm p}/R_{\rm s})_{\rm obs}$ as predicted from transit depth is in the variation of the intensity over the small area of the stellar disk that is temporarily occulted by the planet. The larger the planet, the larger the variation, and the larger the error.

For both of these test cases, we also calculate the mean in-transit flux ${\langle}f_{\rm in}{\rangle}$ to predict $(R_{\rm p}/R_{\rm s})_{\rm obs}$ via Eq. (6). For the Jupiter-sized test planet, we measure ${\langle}f_{\rm in}{\rangle}=0.9894685$ and with $\langle I\rangle_{A}=\leavevmode\nobreak\ 0.8211$ and $\langle I\rangle_{x}=0.887745$ I obtain $(R_{\rm p}/R_{\rm s})_{\rm obs}\leavevmode\nobreak\ =\leavevmode\nobreak\ 0.0986959$. The resulting deviation from the injected planet-to-star radius ratio is $\varepsilon_{\rm p}=6.9\leavevmode\nobreak\ {\times}\leavevmode\nobreak\ 10^{-4}\,R_{\odot}\approx 480$ km. The same exercise for the Earth-like planet yields an error of $\varepsilon_{\delta}=4.2\times 10^{-6}\,R_{\odot}\approx 2.9$ km.

The main source of error in $R_{\rm p}/R_{\rm s}$ as predicted from the mean flux is in the ingress and egress of the planet, which is not taken into account in our calculations of the radially averaged emerging specific intensity of the star. The larger the planet, the longer the times of ingress and egress, and the larger the error in the resulting $(R_{\rm p}/R_{\rm s})_{\rm obs}$ estimate from the chord-averaged intensity.

These tests show that the effect from the small planet approximation is orders of magnitude smaller than any uncertainties arising from the total noise budget in real light curves. Even high-accuracy space-based observations of the brightest and most photometrically quiet stars (Gilliland et al., 2011) have a noise floor of several times 10 parts per million, which translates into an error in the derived radius ratio of $>10^{-4}$ for $R_{\rm p}/R_{\rm s}=0.1$ and an error $>10^{-3}$ for $R_{\rm p}/R_{\rm s}=0.01$ (the relative error increases for smaller planets). Our small planet approximation produces errors in $R_{\rm p}/R_{\rm s}$ that are orders of magnitude smaller than the error coming from the uncertain limb darkening coefficients (Müller et al., 2013); this is a crucial issue for any fitting of observed light curves with numerical methods.

4 Conclusions

In this paper we present analytical expressions for the overshoot of the transit depth for small planets ($R_{\rm p}\leavevmode\nobreak\ {\ll}\leavevmode\nobreak\ R_{\rm s}$) with arbitrary transit impact parameter for the linear, quadratic, square-root, logarithmic, and nonlinear stellar limb darkening laws. Equation (3) can be used to calculate the overshoot of the stellar emerging intensity with respect to the disk-averaged intensity, which translates into an overshoot of the observed transit depth. Equation (4) gives the actual planet-to-star radius ratio based on the measured transit transit depth ($\delta$) and based on the limb darkening law and parameterization. For sun-like stars, exoplanet transits can be $\sim 20$ % deeper than the $(R_{\rm p}/R_{\rm s})^{2}$ estimate.

We also derive analytical expressions to calculate the maximum transit depth $\delta$ from the mean in-transit flux ${\langle}f_{\rm in}{\rangle}$ for any of the limb darkening laws mentioned above. This approach is based on calculating the average stellar intensity along the transit chord of the planet across the stellar disk ${\langle}I{\rangle}_{x}$, which relates to ${\langle}f_{\rm in}{\rangle}$ and $\delta$ as per Eq. (5). Combined with our results for the overshoot of the transit depth, we derive Eq. (6) which gives an estimate for the actual planet-to-star radius ratio based only on the measured in-transit mean flux and on the parameterization of the respective limb darkening law, which requires knowledge of the transit impact parameter.

These expressions can be used to calculate the expected maximum transit depth for a given planet-star system with any parameterization of the common stellar limb darkening laws. This approach is not meant to replace the numerical fitting of light curves, the latter of which can constrain limb darkening and transit impact parameter. Our approach is certainly less accurate than the use of elliptical integrals (Abubekerov & Gostev, 2013) or than the full modeling of the stellar intensity as a sum of spherical harmonics (Luger et al., 2019), let alone the numerical simulation of transits with model stellar atmospheres (Neilson et al., 2018). But our approach requires no numerical modeling and it could be more intuitive to understand.

Acknowledgements.

The author is thankful to Michael Hippke for posing the problem solved in this study and to an anonymous referee for her/his useful comments. This work was supported by the German space agency (Deutsches Zentrum für Luft- und Raumfahrt) under PLATO Data Center grant 50OO1501. This work made use of NASA’s ADS Bibliographic Services

Appendix A Average intensity of the transit chord for the square root limb darkening law

In Sect. 2.3.3, we present our result for the average intensity sampled by the planet along its transit chord in the case of the square root limb darkening law. Its derivation in analogy to Eq. (12) is a bit more complicated than for the other limb darkening laws because the corresponding integral $\int{\rm d}x\,I_{\rm sqr}(x)$ does not have a closed form integral, hence we need to work on that. First, we separate the integral into its summands,

[TABLE]

and we identify the last summand as the integral that does not have an elementary antiderivative. Hence, we substitute $1-p^{2}-x^{2}\,{\equiv}\,u$, which means

[TABLE]

so that

[TABLE]

Equation (41) then simplifies to

[TABLE]

which is Eq. (22).