Analytical Transit Light Curves for Arbitrary Power-Law Limb Darkening: A Unified Framework

TL;DR

This work resolves the longstanding restriction to integer-power limb-darkening in exoplanet transit modeling by introducing a unified analytical framework for arbitrary real power α>−2. Through three equivalent formulations—the geometric kernel, the fractional-calculus approach, and the hypergeometric-integrand representation—the authors derive exact solutions that reproduce Mandel & Agol’s elliptic-integral results when α is integer, while providing closed forms for central transits and finite elliptic representations for half-integer exponents. The framework enables linear superposition for multi-parameter limb-darkening laws and delivers nine orders of magnitude improvement in accuracy and ~10^4× speedup over Monte Carlo methods, facilitating real-time Bayesian inference for JWST, TESS, and future facilities. The combination of mathematical rigor, numerical stability, and computational efficiency offers a powerful tool for utilizing physically motivated limb-darkening profiles in high-precision transit analyses, with immediate observational impact and broad theoretical insights through connections to fractional calculus and special functions. These advances push transit modeling toward physical optimization, enabling more accurate planetary radii, atmospheric characterizations, and tests of stellar atmosphere theories across a wide range of α, including non-integer exponents central to M-dwarfs and Claret-type multi-parameter laws.

Abstract

We present the first exact analytical solutions for exoplanet transit light curves under arbitrary real power-law limb darkening, $I(μ) = I_0μ^α$ with $α> -2$, removing the two-decade restriction to integer-polynomial forms. This generalization addresses a fundamental barrier: physically motivated stellar atmosphere models require non-integer exponents -- most critically the square-root law ($α=1/2$) essential for M-dwarf characterization and half-integer powers in Claret's four-parameter law -- yet existing frameworks support only integer powers, forcing reliance on numerical integration or polynomial approximations that introduce systematic errors at the $\sim 10^{-3}$ level. We develop three mathematically equivalent formulations: (1) a geometric kernel representation expressing blocked flux as a one-dimensional radial integral; (2) a fractional-calculus framework yielding exact hypergeometric solutions via Riemann-Liouville operators; and (3) a numerically stable hypergeometric-integrand form with guaranteed convergence. All three reproduce established elliptic-integral expressions for integer $α$, confirming our framework as a true generalization. For central transits, we obtain the exact closed form $\mathcal{F} = (1-p^2)^{1+α/2}$. Linear superposition enables exact treatment of multi-parameter laws without special-case logic. Validation using 45-digit precision arithmetic demonstrates inter-method agreement at $|Δ\mathcal{F}| \lesssim 10^{-13}$ across all geometric regions, including transcendental exponents ($π/2$, $\sqrt{2}$, $e-2$). Compared with Monte Carlo integration, our method achieves nine orders of magnitude improvement in accuracy with $\sim 10^4\times$ computational speedup, enabling efficient Bayesian inference for JWST, TESS, and next-generation facilities.

Figures (6)

• Figure 1: Inter-method agreement for 80+ test cases spanning all geometric regions and eight $\alpha$ values. Maximum discrepancy $<3\times 10^{-13}$ at machine epsilon. Uniform distribution across flux range confirms stability for deep transits and grazing configurations.
• Figure 2: Normalized flux versus $\alpha$ for central transit ($z=0$, $p=0.3$). Closed form $\mathcal{F}=(1-p^2)^{1+\alpha/2}$ (solid) matches both implementations perfectly across $-1.5 \le \alpha \le 3.5$ with discrepancies $<10^{-13}$.
• Figure 3: Geometric regions with validation test overlay. Background colors: R1 (light blue), R2 (yellow), R3a (green), R3b (purple), R3c (orange). Bold curves mark boundaries. Marker sizes proportional to $|\Delta\mathcal{F}|$. Uniform tiny markers confirm machine-epsilon accuracy across all regions including boundaries.
• Figure 4: Complete transit light curves for $p=0.1$ across $z \in [-1.5,1.5]$ and $\alpha \in \{-1,0,0.5,1,2\}$. Curves deepen monotonically with increasing limb darkening. Square-root law ($\alpha=0.5$, green) critical for M-dwarfs lies between uniform and linear. Smooth ingress/egress with no contact-point artifacts.
• Figure 5: Verification of $\alpha$-differentiation identity for $(p,z)=(0.3,0.7)$. Both sides of $\partial_\alpha\ln\mathcal{F} = (1/2)[\langle\ln(1-\rho^2)\rangle_{\text{vis}}+2/(\alpha+2)]$ overlay perfectly across $-1.5 \le \alpha \le 3.5$. Finite difference (blue) matches direct quadrature (red) with residual $<10^{-12}$.

Theorems & Definitions (18)

• Theorem 1: Triple Equivalence
• Definition 2: Riemann-Liouville fractional integral
• Lemma 3: Radial fractional fundamental theorem
• proof
• Theorem 4: Fractional boundary reduction in polar coordinates
• proof
• Corollary 5: Power-law limb darkening
• Remark 6: Origin of hypergeometric parameters
• Theorem 7: Finite termination for half-integers
• proof