Accurate Computation of Light Curves and the Rossiter-McLaughlin Effect in Multi-Body Eclipsing Systems

Abstract

We present here an efficient method for computing the visible flux for each body during a multi-body eclipsing event for all commonly used limb darkening laws. Our approach follows the idea put forth by Pal (2012) to apply Green's Theorem on the limb darkening integral, thus transforming the two-dimensional flux integral over the visible disk into a one-dimensional integral over the visible boundary. We implement this idea through an iterative process which combines a fast method for describing the visible boundary of each body with a fast numerical integration scheme to compute the integrals. For the two-body case, our method compares well in speed with both that of Mandel & Agol (2002) and that of Gimenez (2006a). The strength of the method is that it works for any number of spherical bodies, with a computational accuracy that is adjustable through the use of a tolerance parameter. Most significantly, the method offers two main advantages over previously used techniques: (i) it can employ a multitude of limb darkening laws, including all of the commonly used ones; (ii) it can compute the Rossiter-McLaughlin effect for rigid body rotation with an arbitrary orientation of the rotation axis, using any of these limb darkening laws. In addition, we can compute the Rossiter-McLaughlin effect for stars exhibiting differential rotation, using the quadratic limb darkening law. We provide the mathematical background for the method and explain in detail how to implement the technique with the help of several examples and codes which we make available.

74 pages, 33 figures, published in the Astronomical Journal Changes from original version: correction to Eq. E56, minor corrections to Figures 21, 29, and various minor stylistic corrections

References in corpus (4)

• KOI-126: A Triply-Eclipsing Hierarchical Triple with Two Low-Mass Stars
• Planet-Planet Occultations in TRAPPIST-1 and Other Exoplanet Systems
• A spectral approach to transit timing variations
• KIC 2856960: the impossible triple star

Cited by in corpus (15)

• Analytic Planetary Transit Light Curves and Derivatives for Stars with Polynomial Limb Darkening
• TOI-1338: TESS' First Transiting Circumbinary Planet
• The TESS light curve of AI Phoenicis
• Wobble: A Data-driven Analysis Technique for Time-series Stellar Spectra
• Kepler-1661 b: A Neptune-sized Kepler Transiting Circumbinary Planet around a Grazing Eclipsing Binary
• NEID Rossiter-McLaughlin Measurement of TOI-1268b: A Young Warm Saturn Aligned with Its Cool Host Star
• The Young Planet DS Tuc Ab has a Low Obliquity
• Transmission strings: a technique for spatially mapping exoplanet atmospheres around their terminators
• The K2 M67 Study: Precise Mass for a Turnoff Star in the Old Open Cluster M67
• Then and now: A new look at the eclipse timing variations of hierarchical triple star candidates in the primordial $Kepler$-field, revisited by TESS
• Analytic Light Curve for Mutual Transits of Two Bodies Across a Limb-darkened Star
• Photodynamical Modeling of the Fascinating Eclipses in the Triple-Star System KOI-126
• A differentiable N-body code for transit timing and dynamical modelling -- II. Photodynamics
• Fast Transit Computation Using Tabulated Stellar Intensities
• TIC 172900988: A Transiting Circumbinary Planet Detected in One Sector of TESS Data