Analytical transit light curves for arbitrary power-law limb darkening: a unified framework

Abstract

We present a unified analytical framework for computing exoplanet transit light curves under arbitrary real power-law limb darkening $I(\mu )=I_0{\mu }^{\alpha }$, $\alpha >-2$, eliminating the two-decade restriction to integer-polynomial forms. Our central result is an exact closed-form expression for the normalised stellar flux in terms of Appell's bivariate hypergeometric function $F_1$, valid for all real $\alpha$ and all geometric transit configurations. Three mutually equivalent formulations -- a geometric kernel representation, a Riemann-Liouville fractional-calculus framework, and a hypergeometric-integrand representation -- provide complementary physical insights and independent computational routes, with inter-method agreement at $|\Delta \mathcal{F}|\lesssim 10^{-17}$. The framework encompasses the square-root law ($\alpha =1/2$) essential for M-dwarf characterisation, half-integer powers required by Claret's four-parameter law, and arbitrary real values enabling empirical fitting, while recovering integer-polynomial elliptic-integral expressions as special cases. Validation at 40-digit precision across five geometric regions and eight $\alpha$ values confirms a $\sim 10^4\times$ computational speedup over Monte Carlo integration and thirteen orders of magnitude accuracy advantage. Linear superposition enables exact treatment of arbitrary multi-parameter limb-darkening laws without special-case logic.

18 pages, Accepted version at Monthly Notices of the Royal Astronomical Society journal

References in corpus (14)

• The Transiting Exoplanet Survey Satellite
• Analytic Lightcurves for Planetary Transit Searches
• The James Webb Space Telescope
• A temperate rocky super-Earth transiting a nearby cool star
• STARRY: Analytic Occultation Light Curves
• Early Release Science of the exoplanet WASP-39b with JWST NIRSpec PRISM
• Stellar Limb-Darkening Coefficients for CoRot and Kepler
• Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H
• Limb darkening and exoplanets: testing stellar model atmospheres and identifying biases in transit parameters
• Analytic Planetary Transit Light Curves and Derivatives for Stars with Polynomial Limb Darkening
• The Stagger-grid: A grid of 3D stellar atmosphere models - IV. Limb darkening coefficients
• Properties of analytic transit light curve models
• Light curve modelling for mutual transits
• CHEOPS observations of the HD 108236 planetary system: A fifth planet, improved ephemerides, and planetary radii